- Open Access
A major genetic locus controlling natural Plasmodium falciparum infection is shared by East and West African Anopheles gambiae
Malaria Journal volume 6, Article number: 87 (2007)
Genetic linkage mapping identified a region of chromosome 2L in the Anopheles gambiae genome that exerts major control over natural infection by Plasmodium falciparum. This 2L Plasmodium-resistance interval was mapped in mosquitoes from a natural population in Mali, West Africa, and controls the numbers of P. falciparum oocysts that develop on the vector midgut. An important question is whether genetic variation with respect to Plasmodium-resistance exists across Africa, and if so whether the same or multiple geographically distinct resistance mechanisms are responsible for the trait.
To identify P falciparum resistance loci in pedigrees generated and infected in Kenya, East Africa, 28 microsatellite loci were typed across the mosquito genome. Genetic linkage mapping was used to detect significant linkage between genotype and numbers of midgut oocysts surviving to 7–8 days post-infection.
A major malaria-control locus was identified on chromosome 2L in East African mosquitoes, in the same apparent position originally identified from the West African population. Presence of this resistance locus explains 75% of parasite free mosquitoes. The Kenyan resistance locus is named EA_Pfin1 (E ast A frica_ P lasmodium f alciparum Infection In tensity).
Detection of a malaria-control locus at the same chromosomal location in both East and West African mosquitoes indicates that, to the level of genetic resolution of the analysis, the same mechanism of Plasmodium-resistance, or a mechanism controlled by the same genomic region, is found across Africa, and thus probably operates in A. gambiae throughout its entire range.
The mosquito Anopheles gambiae is the major African vector of human malaria caused by Plasmodium falciparum. More than 300 million cases of acute malaria occur globally each year resulting in over one million deaths with greater than 90% of these deaths occurring in children from Sub-Saharan Africa . In malaria endemic countries where the poorest sections of society are most affected, inadequate public health initiatives and the selection of insecticide resistance hinder the efforts to control malaria transmission. Recent proposals to curb the malaria burden involve more technological approaches including the generation and introduction of transgenic mosquitoes . However, all approaches have their own logistical, political and ethical limitations.
Here the genetics of natural vector resistance to P. falciparum is dissected, with the ultimate goal to utilize successful existing natural Plasmodium resistance mechanisms for control of human malaria transmission. Knowledge of the prevalence, strength and genomic location of loci that control mosquito resistance to malaria parasites in vector populations could aid future efforts to control parasite transmission. Genetically resistant and susceptible mosquitoes exist in nature, and phenotypic variation can segregate as a simple Mendelian trait . A recent comprehensive study reported the prevalence, strength and genomic location of natural P. falciparum resistance loci in A. gambiae in Mali, West Africa . The loci with the strongest phenotypic effect on parasite numbers clustered to form a Plasmodium-resistance island (PRI) on chromosome 2L. It would be valuable to know whether the genetic mechanism of the PRI identified in Mali is specific to that population, suggesting that there might be many geographically local A. gambiae resistance mechanisms across Africa, or conversely whether the mechanism identified by the PRI exists in A. gambiae throughout Africa. A major barrier to gene flow in A. gambiae occurs between West and East African populations , and thus comparing these populations should offer the most stringent test of the geographic generality of genetic loci controlling resistance mechanisms.
The source of the pedigrees was previously reported . Briefly, isofemale pedigrees were initiated from a laboratory colony (established in Mbita, Kenya in 2001, approximately three years prior to use in this study) and raised to the F1 generation at which point they were challenged on gametocyte containing blood. The same pedigrees were previously used to investigate the role of genotype by genotype interactions in controlling the parasite load of vector mosquitoes . Thus the parasite challenge was heterogeneous: each F1 pedigree was split into 3–5 groups that were each fed on independent natural gametocytemic blood samples. Unfed and partially fed females were excluded from further analysis. Seven to eight days after the infective blood meal, mosquito midguts were dissected and stained with 2% mercurochrome, and the numbers of oocysts were counted by light microscopy. Mosquito carcasses were saved for later DNA extraction.
DNA was isolated from mosquitoes using Qiagen DNeasy (Qiagen, CA, USA). Microsatellite genotyping was performed multiplexed as described previously . Fragments were separated on an ABI3100 Genetic Analyzer and sized using GeneMapper Version 3.5 (both from Applied Biosystems, CA, USA). All automated genotype calls were visually verified. Mapping was carried out in those pedigrees deemed potentially informative using previously described guidelines; ≥20 females with ≥30% infection prevalence . Six pedigrees, named 57, 92, 141, 142, 146 and 1102 met the typing criteria. Twenty-eight microsatellite loci were screened in these informative pedigrees and as many as 15 additional loci were typed in some, but not all pedigrees (Figure 1). Many of these additional markers were added to increase the signal resolution around marker H325. All primer sequences are available elsewhere [4, 7] except 2L.24D left 5' TGCCAATCAATCAGTGTGCT 3' and right 5' CACTAGCAACGGCACTACCC 3' and 2L.26D left 5' GAAGTGGAAGAACACGCTCA and right 5' CACACGATTGCAGATGAGTT 3'. Physical positions given are based on Ensembl A. gambiae genome sequence version 39, converted to the inverted 2La/a form by use of inversion breakpoint coordinates from .
Linkage analysis methods are described in detail elsewhere . Despite the fact that each pedigree was divided into groups that were fed on different gametocyte carriers, in the current study all individuals within a pedigree were first analysed together. In one pedigree yielding a linkage signal, the infections from different parasite isolates were also analysed individually. In brief, at every marker, the pedigree members were split into two samples, using an allele or genotype observed at that marker as the selection criterion. At each locus, the distribution of oocyst numbers for each allele/genotype was compared against the pooled oocyst distribution for all other alleles/genotypes. Permutations were used to establish empirical genome-wide p-values. Wilcoxon-Mann-Whitney (WMW), a nonparametric statistical test, was used to determine statistical significance.
A novel and important aspect of this mapping methodolgy is that it examines each locus in isolation from all other loci and does not require knowledge of relative marker position. Thus, the result is robust to the presence of segregating chromosomal inversions, which are common in A. gambiae, and can alter the physical relationship between markers in different individuals within a pedigree. For example, another recent mapping study of P. falciparum resistance in other Kenyan pedigrees  used a mapping approach in which valid linkage results are dependent on correct knowledge of marker order, but marker order was not empirically determined. That study assumed a marker order  that predated release of the A. gambiae genome sequence. In fact, marker order in the genome sequence differs markedly from that in [7, 9] (e.g., chromosome 2L ref [7, 9] order: AG2H46-AG2H197-AG2H786-AG2H135-AG2H603-AG2H787-AG2H325; genome order and current study: AG2H46-AG2H197-AG2H135-AG2H786-AG2H325-AG2H787-AG2H603).
The statistical analysis was run using both the program OUTBRED_LINES as previously described  and a new implementation of the same algorithm in the statistical analysis environment R , called OUTBRED_LINES-R. The motivation was to avoid use of approximations employed in OUTBRED_LINES and instead use the built-in exact WMW tests. However, very similar results were obtained from OUTBRED_LINES and OUTBRED_LINES-R. Exact p-values from the R implementation OUTBRED_LINES-R are presented throughout the manuscript.
Six Kenyan pedigrees met the minimum informativeness criteria for genotyping (≥20 mosquitoes surviving to dissection, with infection prevalence ≥ 30%), and these pedigrees were completely genotyped (Figure 2). One of the pedigrees (n = 25, pedigree #57 in ) yielded a significant locus for P. falciparum resistance (nominal p = 0.002, genome-wide p = 0.040). The locus displays significant linkage to 5 nearby markers (genome wide p <0.05; from left to right; 2L.H13808912, H796, 2L.17686896, 2L.25C, 2L.25A; Figure 1, Table 1). The Kenyan resistance locus is here named EA_Pfin1 (E ast A frica_ P lasmodium f alciparum Infection In tensity), to distinguish it from the previously named Pfin1 locus, which it overlaps, identified in the chromosome 2L PRI of West African A. gambiae [3, 4]. This mapping with greater marker density shows solid statistical support for the EA_Pfin1 resistance locus from multiple genetic markers (Figure 3). The detection of maximum linkage over a large physical region (~16 Mb) is due to the limited number of recombinants expected in such a small pedigree (n = 25). Despite small pedigree size, the EA_Pfin1 locus had high explanatory power for infection phenotype: the marker allele linked to resistance explained 75% of completely uninfected mosquitoes in the pedigree, and the susceptible allele explained 89% of oocyst-containing mosquitoes (Figure 4). The other five informative pedigrees gave no significant linkage signal (discussed below).
Careful examination of the EA_Pfin1 pedigree indicates that one of the three parasite isolates used failed to infect any mosquitoes in the pedigree, and infected only five mosquitoes of 76 total assayed across all pedigrees . Due to the lack of infectivity of the parasites in this blood sample for the EA_Pfin1 pedigree, the data from mosquitoes fed on this isolate were removed and the pedigree was reanalysed. Without these six mosquitoes, the remaining 19 mosquitoes still give significant linkage signal at the same microsatellite loci despite the small sample size (Table 2). This result indicates that the linkage signal did not depend solely or disproportionately upon defeat by the mosquitoes of one of the three parasite isolates fed to them. Indeed, even with loss of 24% of the pedigree sample size by removal of the six mosquitoes, the linkage signal was still comparable to that from the entire pedigree. The surviving linkage signal in the 19 mosquitoes was due to mosquitoes in the pedigree whose parasite challenge was divided between two different gametocyte carriers. Analysis of the remaining two groups fed on either gametocyte carrier alone did not yield a linkage signal (not shown).
Thus, the detected linkage signal in the entire pedigree or any subset did not result from a strong protective effect by mosquitoes against any single parasite isolate. Instead, detection of linkage required the larger sample size obtained by combining groups of mosquitoes that fed on at least two different parasite isolates, and consequently the signal must have originated in the pedigree's consistent genetic response to at least two different parasite isolates. Taken together, these results indicate that Plasmodium resistance controlled by the A. gambiae EA_Pfin1 locus can be independent of and transcend parasite genotype.
A P. falciparum resistance locus, EA_Pfin1, is reported in East African A. gambiae that resides within the chromosomal Plasmodium-resistance Island (PRI) detected in West African A. gambiae. The phenotype is also the same: parasite elimination prior to oocyst development. The vast majority (>75%) of mosquitoes in the EA_Pfin1 pedigree carrying the resistance allele completely blocked parasite development, and thus were parasite-free. This is comparable to the phenotypic effect size of resistance loci in the 2L PRI in West Africa . Due to their chromosomal co-localization and similarity of phenotype, EA_Pfin1 and the 2L PRI in West Africa likely control the same underlying Plasmodium resistance mechanism.
Geographic generality of resistance
The geographic generality of resistance mechanisms indicates that the PRI-controlled mechanism of resistance is not a nuance of the pedigrees derived from the West Africa population initially studied. Another recent study in Kenya  detected an effect on P. falciparum infection outcome, but physical positions of the markers and therefore the mapped loci in that study were not determined and thus it is not possible to resolve whether the effect was controlled by the same genomic region as reported here. Most parsimoniously, a common resistance mechanism across the African continent suggests a relatively ancient origin for natural P. falciparum resistance rather than the evolution of distinct resistance mechanisms multiple times. The ancestral PRI-controlled mechanism might have been selected in an A. gambiae ancestor by exposure to a non-Plasmodium pathogen, or might have newly evolved in A. gambiae soon after the advent of the P. falciparum host-pathogen system.
Factors on chromosome 2L control malaria parasite infection
Genetic factors on chromosome 2L have previously been shown to affect malaria infection both in nature  and in a laboratory system [12, 13]. In addition, the form of the polymorphic 2La inversion is also correlated with adaptation to climate , biting and resting behavior , and morphology . The previous mapping pedigrees from West Africa were all fixed inverted homozygotes (2La/a) for the 2La inversion . A molecular diagnostic assay designed to determine the orientation of the breakpoints of the 2La inversion  indicated that the EA_Pfin1 pedigree in the current study is similarly fixed for the inverted chromosome (data not shown).
Therefore, the data show that at least one strong resistance locus segregates within the inverted form of 2La across Africa. Without homozygous 2La standard form (2La+/+) pedigrees or pedigrees with the 2La inversion segregating (which would include heterozygotes and both homozygotes) the data cannot address whether there is linkage between inversion state and parasite infection loads. Given the lack of any pedigrees in the West African studies fixed for or segregating the 2La+/+ standard form arrangement, and the colony origin and small size of the current East African pedigrees, it is also impossible to know whether this PRI-related mechanism of P. falciparum resistance is unique to the inverted form or common across forms. In a confusing turn of nomenclature, the inverted arrangement of 2La is actually ancestral, and the standard form is derived . Thus, it is formally possible that the monophyletic inversion event leading to the "standard" chromosome (2La+/+) could have captured a susceptible allele from the PRI. Another mapping study done in East Africa cannot be compared because it did not report the inversion state of 2La .
Resistance alleles maintained in colonies
Prior linkage mapping [3, 4] of P. falciparum resistance loci used F1 or F2 pedigrees that were only a generation or two removed from the wild, while the current study uses pedigrees generated from a colony established several years earlier. It is important to emphasize the importance of detecting resistance loci in a natural population versus observation in a colony, where selection and genetic drift can have a strong influence on the prevalence of susceptibility/resistance alleles. Repeated observations in pedigrees derived from the natural mosquito population  is a reliable proxy for prevalence of resistance/susceptibility in the field. In contrast, breeding in a colony can amplify or diminish the prevalence of the relevant alleles. Thus, it is only in light of the previous West African results [3, 4], as well as a consistent observation from an independent effort , that it can be concluded that the chromosome 2L PRI is indeed a common resistance mechanism in the Kenyan field population.
Detection of P. falciparum resistance alleles in a colony-derived pedigree demonstrates that despite the large loss of genetic diversity that accompanies colonization , and the absence of any putative parasite-imposed selective pressure for many generations, alleles for resistance and susceptibility to parasite infection can still segregate within a colony. The ability of resistant and susceptible alleles to persist in a colony setting is advantageous for future mechanistic studies of natural genetic resistance factors.
Parasite strain transcendence
The current mapping study allowed testing for effects across parasite isolates within a mosquito pedigree. Detection of linkage in a pedigree in which different mosquitoes were fed on multiple distinct parasite isolates is consistent with the fact that the allele for resistance in the EA_Pfin1 mapping pedigree controls an effect that transcends parasite strain. This is based on the assumption of non-relatedness between parasite isolates, an assumption which is consistent with the previous analysis of the same infected pedigrees . Parasite strain transcendence or independence is evident from the genetic mapping in West Africa , in which multiple independent isofemale pedigrees over multiple years, each fed on blood from a different gametocyte carrier, nevertheless yielded loci in the same 2La genomic region.
Here it is established that the genomic region of A. gambiae chromosome 2L, containing strong Plasmodium-resistance loci originally identified in pedigrees from a West Africa population of vector mosquitoes, is not a geographically local phenomenon. At the level of resolution of the genetic mapping, the same genetic factor also exists in East African A. gambiae. The intriguing prospect of shared P. falciparum resistance loci across Africa makes it a promising approach to consider taking advantage of widespread natural Plasmodium resistance mechanisms for the control of malaria transmission [18, 19].
World Health Organization. [http://www.who.int/tdr/diseases/malaria]
Christophides GK: Transgenic mosquitoes and malaria transmission. Cell Microbiol. 2005, 7 (3): 325-333. 10.1111/j.1462-5822.2005.00495.x.
Niare O, Markianos K, Volz J, Oduol F, Toure A, Bagayoko M, Sangare D, Traore SF, Wang R, Blass C, Dolo G, Bouare M, Kafatos FC, Kruglyak L, Toure YT, Vernick KD: Genetic loci affecting resistance to human malaria parasites in a West African mosquito vector population. Science. 2002, 298 (5591): 213-216. 10.1126/science.1073420.
Riehle MM, Markianos K, Niare O, Xu J, Li J, Toure AM, Podiougou B, Oduol F, Diawara S, Diallo M, Coulibaly B, Ouatara A, Kruglyak L, Traore SF, Vernick KD: Natural malaria infection in Anopheles gambiae is regulated by a single genomic control region. Science. 2006, 312 (5773): 577-579. 10.1126/science.1124153.
Lehmann T, Licht M, Elissa N, Maega BT, Chimumbwa JM, Watsenga FT, Wondji CS, Simard F, Hawley WA: Population Structure of Anopheles gambiae in Africa. J Hered. 2003, 94 (2): 133-147. 10.1093/jhered/esg024.
Lambrechts L, Halbert J, Durand P, Gouagna LC, Koella JC: Host genotype by parasite genotype interactions underlying the resistance of anopheline mosquitoes to Plasmodium falciparum. Malar J. 2005, 4: 3-10.1186/1475-2875-4-3.
Zheng L, Benedict MQ, Cornel AJ, Collins FH, Kafatos FC: An integrated genetic map of the African human malaria vector mosquito, Anopheles gambiae. Genetics. 1996, 143: 941-952.
Sharakhov IV, White BJ, Sharakhova MV, Kayondo J, Lobo NF, Santolamazza F, Della Torre A, Simard F, Collins FH, Besansky NJ: Breakpoint structure reveals the unique origin of an interspecific chromosomal inversion (2La) in the Anopheles gambiae complex. Proc Natl Acad Sci U S A. 2006, 103 (16): 6258-6262. 10.1073/pnas.0509683103.
Menge DM, Zhong D, Guda T, Gouagna L, Githure J, Beier J, Yan G: Quantitative trait loci controlling refractoriness to Plasmodium falciparum in natural Anopheles gambiae mosquitoes from a malaria-endemic region in western Kenya. Genetics. 2006, 173 (1): 235-241. 10.1534/genetics.105.055129.
R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. [http://www.r-project.org]
Petrarca V, Beier JC: Intraspecific chromosomal polymorphism in the Anopheles gambiae complex as a factor affecting malaria transmission in the Kisumu area of Kenya. Am J Trop Med Hyg. 1992, 46 (2): 229-237.
Vernick KD, Collins FH: Association of a Plasmodium-refractory phenotype with an esterase locus in Anopheles gambiae. Am J Trop Med Hyg. 1989, 40 (6): 593-597.
Crews-Oyen AE, Kumar V, Collins FH: Association of two esterase genes, a chromosomal inversion, and susceptibility to Plasmodium cynomolgi in the African malaria vector Anopheles gambiae. Am J Trop Med Hyg. 1993, 49 (3): 341-347.
Toure YT, Petrarca V, Traore SF, Coulibaly A, Maiga HM, Sankare O, Sow M, Di Deco MA, Coluzzi M: The distribution and inversion polymorphism of chromosomally recognized taxa of the Anopheles gambiae complex in Mali, West Africa. Parassitologia. 1998, 40 (4): 477-511.
Rishikesh N, Di Deco MA, Petrarca V, Coluzzi M: Seasonal variations in indoor resting Anopheles gambiae and Anopheles arabiensis in Kaduna, Nigeria. Acta Trop. 1985, 42 (2): 165-170.
Petrarca V, Sabatinelli G, Di Deco MA, Papakay M: The Anopheles gambiae complex in the Federal Islamic Republic of Comoros (Indian Ocean): some cytogenetic and biometric data. Parassitologia. 1990, 32 (3): 371-380.
Norris DE, Shurtleff AC, Toure YT, Lanzaro GC: Microsatellite DNA polymorphism and heterozygosity among field and laboratory populations of Anopheles gambiae ss (Diptera: Culicidae). J Med Entomol. 2001, 38 (2): 336-340.
Hahn MW, Nuzhdin SV: The fixation of malaria refractoriness in mosquitoes. Curr Biol. 2004, 14 (7): R264-265. 10.1016/j.cub.2004.03.018.
Hurd H, Taylor PJ, Adams D, Underhill A, Eggleston P: Evaluating the costs of mosquito resistance to malaria parasites. Evolution Int J Org Evolution. 2005, 59 (12): 2560-2572.
This work was supported by grants from NIAID (KDV), NHGRI (KM), and was aided by resources of the University of Minnesota Supercomputing Institute Computational Genetics Laboratory. The assistance of Dr. Louis Gouagna and the staff of the Malaria Department at ICIPE-Mbita is acknowledged for their assistance during the generation and experimental infection of the pedigrees, and for helpful comments on the manuscript. Informed consent was obtained for studies on humans after the nature and possible consequences of the studies were explained.
This work was a collaborative effort involving all authors. LL initiated and infected the Kenyan pedigrees and JK provided DNA samples from them. MMR performed genotyping assays and statistical analysis. KM performed statistical analyses. IS and AX typed the pedigree yielding a linkage signal for the 2La inversion. MMR, KM, and KV discussed analysis of results and interpreted their meaning. Both MMR and KM wrote initial drafts of the manuscript. The final draft of the ms was read and approved by all authors.
About this article
Cite this article
Riehle, M.M., Markianos, K., Lambrechts, L. et al. A major genetic locus controlling natural Plasmodium falciparum infection is shared by East and West African Anopheles gambiae. Malar J 6, 87 (2007) doi:10.1186/1475-2875-6-87
- Resistance Locus
- Parasite Isolate
- Gametocyte Carrier
- Geographic Generality