Impact of natural selection on global patterns of genetic variation and association with clinical phenotypes at genes involved in SARS-CoV-2 infection.

Human genomic diversity has been shaped by both ancient and ongoing challenges from viruses. The current coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a devastating impact on population health. However, genetic diversity and evolutionary forces impacting host genes related to SARS-CoV-2 infection are not well understood. We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (angiotensin converting enzyme 2 [ACE2], transmembrane protease serine 2 [TMPRSS2], dipeptidyl peptidase 4 [DPP4], and lymphocyte antigen 6 complex locus E [LY6E]). We analyzed data from 2,012 ethnically diverse Africans and 15,977 individuals of European and African ancestry with electronic health records and integrated with global data from the 1000 Genomes Project. At ACE2, we identified 41 nonsynonymous variants that were rare in most populations, several of which impact protein function. However, three nonsynonymous variants (rs138390800, rs147311723, and rs145437639) were common among central African hunter-gatherers from Cameroon (minor allele frequency 0.083 to 0.164) and are on haplotypes that exhibit signatures of positive selection. We identify signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage compared with the chimpanzee genome. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19. Our study provides insights into global variation at host genes related to SARS-CoV-2 infection, which have been shaped by natural selection in some populations, possibly due to prior viral infections.

Genetic origins of lactase persistence and the spread of pastoralism in Africa.

In humans, the ability to digest lactose, the sugar in milk, declines after weaning because of decreasing levels of the enzyme lactase-phlorizin hydrolase, encoded by LCT. However, some individuals maintain high enzyme amounts and are able to digest lactose into adulthood (i.e., they have the lactase-persistence [LP] trait). It is thought that selection has played a major role in maintaining this genetically determined phenotypic trait in different human populations that practice pastoralism. To identify variants associated with the LP trait and to study its evolutionary history in Africa, we sequenced MCM6 introns 9 and 13 and ~2 kb of the LCT promoter region in 819 individuals from 63 African populations and in 154 non-Africans from nine populations. We also genotyped four microsatellites in an ~198 kb region in a subset of 252 individuals to reconstruct the origin and spread of LP-associated variants in Africa. Additionally, we examined the association between LP and genetic variability at candidate regulatory regions in 513 individuals from eastern Africa. Our analyses confirmed the association between the LP trait and three common variants in intron 13 (C-14010, G-13907, and G-13915). Furthermore, we identified two additional LP-associated SNPs in intron 13 and the promoter region (G-12962 and T-956, respectively). Using neutrality tests based on the allele frequency spectrum and long-range linkage disequilibrium, we detected strong signatures of recent positive selection in eastern African populations and the Fulani from central Africa. In addition, haplotype analysis supported an eastern African origin of the C-14010 LP-associated mutation in southern Africa.

Identifying Darwinian selection acting on different human APOL1 variants among diverse African populations.

Disease susceptibility can arise as a consequence of adaptation to infectious disease. Recent findings have suggested that higher rates of chronic kidney disease (CKD) in individuals with recent African ancestry might be attributed to two risk alleles (G1 and G2) at the serum-resistance-associated (SRA)-interacting-domain-encoding region of APOL1. These two alleles appear to have arisen adaptively, possibly as a result of their protective effects against human African trypanosomiasis (HAT), or African sleeping sickness. In order to explore the distribution of potential functional variation at APOL1, we studied nucleotide variation in 187 individuals across ten geographically and genetically diverse African ethnic groups with exposure to two Trypanosoma brucei subspecies that cause HAT. We observed unusually high levels of nonsynonymous polymorphism in the regions encoding the functional domains that are required for lysing parasites. Whereas allele frequencies of G2 were similar across all populations (3%-8%), the G1 allele was only common in the Yoruba (39%). Additionally, we identified a haplotype (termed G3) that contains a nonsynonymous change at the membrane-addressing-domain-encoding region of APOL1 and is present in all populations except for the Yoruba. Analyses of long-range patterns of linkage disequilibrium indicate evidence of recent selection acting on the G3 haplotype in Fulani from Cameroon. Our results indicate that the G1 and G2 variants in APOL1 are geographically restricted and that there might be other functional variants that could play a role in HAT resistance and CKD risk in African populations.

Convergent adaptation of human lactase persistence in Africa and Europe.

A SNP in the gene encoding lactase (LCT) (C/T-13910) is associated with the ability to digest milk as adults (lactase persistence) in Europeans, but the genetic basis of lactase persistence in Africans was previously unknown. We conducted a genotype-phenotype association study in 470 Tanzanians, Kenyans and Sudanese and identified three SNPs (G/C-14010, T/G-13915 and C/G-13907) that are associated with lactase persistence and that have derived alleles that significantly enhance transcription from the LCT promoter in vitro. These SNPs originated on different haplotype backgrounds from the European C/T-13910 SNP and from each other. Genotyping across a 3-Mb region demonstrated haplotype homozygosity extending >2.0 Mb on chromosomes carrying C-14010, consistent with a selective sweep over the past approximately 7,000 years. These data provide a marked example of convergent evolution due to strong selective pressure resulting from shared cultural traits-animal domestication and adult milk consumption.

Origin and differential selection of allelic variation at TAS2R16 associated with salicin bitter taste sensitivity in Africa.

Bitter taste perception influences human nutrition and health, and the genetic variation underlying this trait may play a role in disease susceptibility. To better understand the genetic architecture and patterns of phenotypic variability of bitter taste perception, we sequenced a 996 bp region, encompassing the coding exon of TAS2R16, a bitter taste receptor gene, in 595 individuals from 74 African populations and in 94 non-Africans from 11 populations. We also performed genotype-phenotype association analyses of threshold levels of sensitivity to salicin, a bitter anti-inflammatory compound, in 296 individuals from Central and East Africa. In addition, we characterized TAS2R16 mutants in vitro to investigate the effects of polymorphic loci identified at this locus on receptor function. Here, we report striking signatures of positive selection, including significant Fay and Wu's H statistics predominantly in East Africa, indicating strong local adaptation and greater genetic structure among African populations than expected under neutrality. Furthermore, we observed a "star-like" phylogeny for haplotypes with the derived allele at polymorphic site 516 associated with increased bitter taste perception that is consistent with a model of selection for "high-sensitivity" variation. In contrast, haplotypes carrying the "low-sensitivity" ancestral allele at site 516 showed evidence of strong purifying selection. We also demonstrated, for the first time, the functional effect of nonsynonymous variation at site 516 on salicin phenotypic variance in vivo in diverse Africans and showed that most other nonsynonymous substitutions have weak or no effect on cell surface expression in vitro, suggesting that one main polymorphism at TAS2R16 influences salicin recognition. Additionally, we detected geographic differences in levels of bitter taste perception in Africa not previously reported and infer an East African origin for high salicin sensitivity in human populations.

Directional selection at the pfmdr1, pfcrt, pfubp1, and pfap2mu loci of Plasmodium falciparum in Kenyan children treated with ACT.

BACKGROUND: The efficacy of artemisinin-based combination therapy (ACT) for Plasmodium falciparum malaria may be threatened by parasites with reduced responsiveness to artemisinins. Among 298 ACT-treated children from Mbita, Kenya, submicroscopic persistence of P. falciparum on day 3 posttreatment was associated with subsequent microscopically detected parasitemia at days 28 or 42. METHODS: DNA sequences of resistance-associated parasite loci pfcrt, pfmdr1, pfubp1, and pfap2mu were determined in the Mbita cohort before treatment, on days 2 and 3 after initiation of treatment, and on the day of treatment failure. RESULTS: Parasites surviving ACT on day 2 or day 3 posttreatment were significantly more likely than the baseline population to carry the wild-type haplotypes of pfcrt (CVMNK at codons 72-76; P < .001) and pfmdr1 (NFD at codons 86, 184, 1246; P < .001). In contrast, variant alleles of the novel candidate resistance genes pfap2mu (S160N/T; P = .006) and pfubp-1 (E1528D; P < .001) were significantly more prevalent posttreatment. No genetic similarities were found to artemisinin-tolerant parasites recently described in Cambodia. CONCLUSIONS: Among treated children in western Kenya, certain P. falciparum genotypes defined at pfcrt, pfmdr1, pfap2mu, and pfubp1 more often survive ACT at the submicroscopic level, and contribute to onward transmission and subsequent patent recrudescence.

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

Malaria has been a very strong selection pressure in recent human evolution, particularly in Africa. Of the one million deaths per year due to malaria, more than 90% are in sub-Saharan Africa, a region with high levels of genetic variation and population substructure. However, there have been few studies of nucleotide variation at genetic loci that are relevant to malaria susceptibility across geographically and genetically diverse ethnic groups in Africa. Invasion of erythrocytes by Plasmodium falciparum parasites is central to the pathology of malaria. Glycophorin A (GYPA) and B (GYPB), which determine MN and Ss blood types, are two major receptors that are expressed on erythrocyte surfaces and interact with parasite ligands. We analyzed nucleotide diversity of the glycophorin gene family in 15 African populations with different levels of malaria exposure. High levels of nucleotide diversity and gene conversion were found at these genes. We observed divergent patterns of genetic variation between these duplicated genes and between different extracellular domains of GYPA. Specifically, we identified fixed adaptive changes at exons 3-4 of GYPA. By contrast, we observed an allele frequency spectrum skewed toward a significant excess of intermediate-frequency alleles at GYPA exon 2 in many populations; the degree of spectrum distortion is correlated with malaria exposure, possibly because of the joint effects of gene conversion and balancing selection. We also identified a haplotype causing three amino acid changes in the extracellular domain of glycophorin B. This haplotype might have evolved adaptively in five populations with high exposure to malaria.

Limited evidence for adaptive evolution and functional effect of allelic variation at rs702424 in the promoter of the TAS2R16 bitter taste receptor gene in Africa.

Bitter taste perception, mediated by receptors encoded by the TAS2R loci, has important roles in human health and nutrition. Prior studies have demonstrated that nonsynonymous variation at site 516 in the coding exon of TAS2R16, a bitter taste receptor gene on chromosome 7, has been subject to positive selection and is strongly correlated with differences in sensitivity to salicin, a bitter anti-inflammatory compound, in human populations. However, a recent study suggested that the derived G-allele at rs702424 in the TAS2R16 promoter has also been the target of recent selection and may have an additional effect on the levels of salicin bitter taste perception. Here, we examined alleles at rs702424 for signatures of selection using Extended Haplotype Homozygosity (EHH) and FST statistics in diverse populations from West Central, Central and East Africa. We also performed a genotype-phenotype analysis of salicin sensitivity in a subset of 135 individuals from East Africa. Based on our data, we did not find evidence for positive selection at rs702424 in African populations, suggesting that nucleotide position 516 is likely the site under selection at TAS2R16. Moreover, we did not detect a significant association between rs702424 alleles and salicin bitter taste recognition, implying that this site does not contribute to salicin phenotypic variance. Overall, this study of African diversity provides further information regarding the genetic architecture and evolutionary history of a biologically-relevant trait in humans.

Residual Plasmodium falciparum parasitemia in Kenyan children after artemisinin-combination therapy is associated with increased transmission to mosquitoes and parasite recurrence.

BACKGROUND: Parasite clearance time after artemisinin-based combination therapy (ACT) may be increasing in Asian and African settings. The association between parasite clearance following ACT and transmissibility is currently unknown. METHODS: We determined parasite clearance dynamics by duplex quantitative polymerase chain reaction (qPCR) in samples collected in the first 3 days after treatment of uncomplicated malaria with ACT. Gametocyte carriage was determined by Pfs25 quantitative nucleic acid sequence-based amplification assays; infectiousness to mosquitoes by membrane-feeding assays on day 7 after treatment. RESULTS: Residual parasitemia was detected by qPCR in 31.8% (95% confidence interval [CI], 24.6-39.8) of the children on day 3 after initiation of treatment. Residual parasitemia was associated with a 2-fold longer duration of gametocyte carriage (P = .0007), a higher likelihood of infecting mosquitoes (relative risk, 1.95; 95% CI, 1.17-3.24; P = .015), and a higher parasite burden in mosquitoes (incidence rate ratio, 2.92; 95% CI, 1.61-5.31; P < .001). Children with residual parasitemia were also significantly more likely to experience microscopically detectable parasitemia during follow-up (relative risk, 11.25; 95% CI, 4.08-31.01; P < .001). CONCLUSIONS: Residual submicroscopic parasitemia is common after ACT and is associated with a higher transmission potential. Residual parasitemia may also have consequences for individual patients because of its higher risk of recurrent parasitemia.

Malaria transmission after artemether-lumefantrine and dihydroartemisinin-piperaquine: a randomized trial.

BACKGROUND: Artemisinin-based combination therapy (ACT) reduces the potential for malaria transmission, compared with non-ACTs. It is unclear whether this effect differs between ACTs. METHODS: A total of 298 children (age, 6 months to 10 years) with uncomplicated falciparum malaria were randomized to artemether-lumefantrine (AL; n = 153) or dihydroartemisinin-piperaquine (DP; n = 145) in Mbita, a community in western Kenya. Gametocyte carriage was determined by molecular methods on days 0, 1, 2, 3, 7, 14, 28, and 42 after treatment initiation. The gametocyte infectiousness to mosquitoes was determined by mosquito-feeding assays on day 7 after beginning therapy. RESULTS: The cumulative risk of recurrent parasitemia on day 42 after initiation of treatment, unadjusted by polymerase chain reaction findings, was 20.7% (95% confidence interval [CI], 14.4-28.2) for AL, compared with 3.7% (95% CI, 1.2-8.5) for DP (P < .001). The mean duration of gametocyte carriage was 5.5 days (95% CI, 3.6-8.5) for AL and 15.3 days (95% CI, 9.7-24.2) for DP (P = .001). The proportion of mosquitoes that became infected after feeding on blood from AL-treated children was 1.88% (43 of 2293), compared with 3.50% (83 of 2371) for those that fed on blood from DP-treated children (P = .06); the oocyst burden among mosquitoes was lower among those that fed on blood from AL-treated children (P = .005) CONCLUSIONS: While DP was associated with a longer prophylactic time after treatment, gametocyte carriage and malaria transmission to mosquitoes was lower after AL treatment. CLINICAL TRIALS REGISTRATION: NCT00868465.

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.

The genetic and evolutionary basis of gene expression variation in East Africans.

BACKGROUND: Mapping of quantitative trait loci (QTL) associated with molecular phenotypes is a powerful approach for identifying the genes and molecular mechanisms underlying human traits and diseases, though most studies have focused on individuals of European descent. While important progress has been made to study a greater diversity of human populations, many groups remain unstudied, particularly among indigenous populations within Africa. To better understand the genetics of gene regulation in East Africans, we perform expression and splicing QTL mapping in whole blood from a cohort of 162 diverse Africans from Ethiopia and Tanzania. We assess replication of these QTLs in cohorts of predominantly European ancestry and identify candidate genes under selection in human populations. RESULTS: We find the gene regulatory architecture of African and non-African populations is broadly shared, though there is a considerable amount of variation at individual loci across populations. Comparing our analyses to an equivalently sized cohort of European Americans, we find that QTL mapping in Africans improves the detection of expression QTLs and fine-mapping of causal variation. Integrating our QTL scans with signatures of natural selection, we find several genes related to immunity and metabolism that are highly differentiated between Africans and non-Africans, as well as a gene associated with pigmentation. CONCLUSION: Extending QTL mapping studies beyond European ancestry, particularly to diverse indigenous populations, is vital for a complete understanding of the genetic architecture of human traits and can reveal novel functional variation underlying human traits and disease.

Adaptive differentiation of Plasmodium falciparum populations inferred from single-nucleotide polymorphisms (SNPs) conferring drug resistance and from neutral SNPs.

BACKGROUND: Theoretical and experimental data support the geographic differentiation strategy as a valuable tool for detecting loci under selection. In the context of Plasmodium falciparum malaria, few populations have been studied, with limited genomic coverage. METHODS: We examined geographic differentiation in P. falciparum populations on the basis of 12 single-nucleotide polymorphisms (SNPs) in 4 genes encoding drug resistance determinants, 5 SNPs in 2 genes encoding antigens, and a set of 17 putatively neutral SNPs dispersed on 13 chromosomes. We sampled 326 parasite isolates representing 7 P. falciparum populations from regions with varied levels of malaria transmission (Gabon, Kenya, Madagascar, Mali, Mayotte, Haiti, and the Philippines). RESULTS: Frequencies of drug resistance alleles varied considerably among populations (mean F(ST), 0.52). In contrast, allele frequencies varied significantly less for antigenic and neutral SNPs (mean F(ST), 0.16 and 0.24, respectively). This contrasting pattern was more pronounced when only the African populations were considered. Signature of selection was detected for most of the resistant SNPs but not for the antigenic SNPs. CONCLUSION: These data further validate the utility of geographic differentiation for identifying loci under strong positive selection, such as drug resistance loci. This study also provides frequencies of molecular makers of resistance in some overlooked populations.

Impact of natural selection on global patterns of genetic variation, and association with clinical phenotypes, at genes involved in SARS-CoV-2 infection.

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (ACE2, TMPRSS2, DPP4, and LY6E). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2, we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations.

BACKGROUND: Africa is the origin of modern humans within the past 300 thousand years. To infer the complex demographic history of African populations and adaptation to diverse environments, we sequenced the genomes of 92 individuals from 44 indigenous African populations. RESULTS: Genetic structure analyses indicate that among Africans, genetic ancestry is largely partitioned by geography and language, though we observe mixed ancestry in many individuals, consistent with both short- and long-range migration events followed by admixture. Phylogenetic analysis indicates that the San genetic lineage is basal to all modern human lineages. The San and Niger-Congo, Afroasiatic, and Nilo-Saharan lineages were substantially diverged by 160 kya (thousand years ago). In contrast, the San and Central African rainforest hunter-gatherer (CRHG), Hadza hunter-gatherer, and Sandawe hunter-gatherer lineages were diverged by ~ 120-100 kya. Niger-Congo, Nilo-Saharan, and Afroasiatic lineages diverged more recently by ~ 54-16 kya. Eastern and western CRHG lineages diverged by ~ 50-31 kya, and the western CRHG lineages diverged by ~ 18-12 kya. The San and CRHG populations maintained the largest effective population size compared to other populations prior to 60 kya. Further, we observed signatures of positive selection at genes involved in muscle development, bone synthesis, reproduction, immune function, energy metabolism, and cell signaling, which may contribute to local adaptation of African populations. CONCLUSIONS: We observe high levels of genomic variation between ethnically diverse Africans which is largely correlated with geography and language. Our study indicates ancient population substructure and local adaptation of Africans.

Correction to: African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations.

Following publication of the original article [1], a typographical error in the formula for calculating di in the "Scans for local adaptation" subsection in the Method section, was identified. The correct formula should be.

Increased Plasmodium falciparum gametocyte production in mixed infections with P. malariae.

Plasmodium falciparum and P. malariae occur endemically in many parts of Africa. Observations from malariotherapy patients suggest that co-infection with P. malariae may increase P. falciparum gametocyte production. We determined P. falciparum gametocyte prevalence and density by quantitative nucleic acid sequence-based amplification (QT-NASBA) after antimalarial treatment of Kenyan children with either P. falciparum mono-infection or P. falciparum and P. malariae mixed infection. In addition, we analyzed the relationship between mixed species infections and microscopic P. falciparum gametocyte prevalence in three datasets from previously published studies. In Kenyan children, QT-NASBA gametocyte density was increased in mixed species infections (P = 0.03). We also observed higher microscopic prevalences of P. falciparum gametocytes in mixed species infections in studies from Tanzania and Kenya (odds ratio = 2.15, 95% confidence interval = 0.99-4.65 and 2.39, 1.58-3.63) but not in a study from Nigeria. These data suggest that co-infection with P. malariae is correlated with increased P. falciparum gametocytemia.

A randomized trial to monitor the efficacy and effectiveness by QT-NASBA of artemether-lumefantrine versus dihydroartemisinin-piperaquine for treatment and transmission control of uncomplicated Plasmodium falciparum malaria in western Kenya.

BACKGROUND: Many countries have implemented artemisinin-based combination therapy (ACT) for the first-line treatment of malaria. Although many studies have been performed on efficacy and tolerability of the combination arthemeter-lumefantrine (AL) or dihydroartemisinin-piperaquine (DP), less is known of the effect of these drugs on gametocyte development, which is an important issue in malaria control. METHODS AND RESULTS: In this two-arm randomized controlled trial, 146 children were treated with either AL or DP. Both groups received directly observed therapy and were followed for 28 days after treatment. Blood samples were analysed with microscopy and NASBA. In comparison with microscopy NASBA detected much more gametocyte positive individuals. Moreover, NASBA showed a significant difference in gametocyte clearance in favour of AL compared to DP. The decline of parasitaemia was slower and persistence or development of gametocytes was significantly higher and longer at day 3, 7 and 14 in the DP group but after 28 days no difference could be observed between both treatment arms. CONCLUSION: Although practical considerations could favour the use of one drug over another, the effect on gametocytogenesis should also be taken into account and studied further using molecular tools like NASBA. This also applies when a new drug is introduced. TRIAL REGISTRATION: Current controlled trials ISRCTN36463274.

Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection.

Submicroscopic Plasmodium falciparum gametocytemia (<5,000 gametocytes/mL) is common and may result in mosquito infection. We assessed the relation between gametocyte density and mosquito infection under experimental and field conditions using real-time quantitative nucleic acid sequence-based amplification (QT-NASBA) for gametocyte quantification. Serial dilutions of NF54 P. falciparum gametocytes showed a positive association between gametocyte density and the proportion of infected mosquitoes (beta=6.1; 95% confidence interval [CI], 2.7-9.6; P=0.001). Successful infection became unlikely below an estimated density of 250-300 gametocytes/mL. In the field, blood samples of 100 naturally infected children showed a positive association between gametocyte density and oocyst counts in mosquitoes (beta=0.38; 95% CI, 0.14-0.61; P=0.002). The relative contribution to malaria transmission was similar for carriers with submicroscopic and microscopic gametocytemia. Our results show that transmission occurs efficiently at submicroscopic gametocyte densities and that carriers harboring submicroscopic gametocytemia constitute a considerable proportion of the human infectious reservoir.

Comparison of microsatellite and antigen-coding loci for differentiating recrudescing Plasmodium falciparum infections from reinfections in Kenya.

We have compared the ability of five Plasmodium falciparum microsatellites and three antigen-coding loci to differentiate recrudescence from reinfection. We used 133 pairs of P. falciparum-infected blood samples collected during in vivo drug efficacy trials from three sites in Kenya with different malaria endemicities. There were no significant differences between the marker subsets in their ability to discriminate recrudescences from new infections across the three sites. Overall, microsatellite loci revealed significantly higher expected heterozygosity and multiplicity of infection levels than antigen-coding loci. The mean expected heterozygosity across all loci in the three populations was significantly higher with microsatellites (0.70, 0.78 and 0.79) than antigen-coding loci (0.53, 0.60 and 0.62) for Mwea, Tiwi and Bondo areas, respectively. These observations can be explained by three non-exclusive hypotheses: (i) microsatellites are more polymorphic than antigenic loci; (ii) partially immune hosts remove certain parasites from infections on the basis of their antigenic alleles; and/or (iii) recombination occurs in vitro or in vivo with microsatellites.