A major genetic locus controlling natural Plasmodium falciparum infection is shared by East and West African Anopheles gambiae
© Riehle et al; licensee BioMed Central Ltd. 2007
Received: 02 January 2007
Accepted: 06 July 2007
Published: 06 July 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.
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.
Mapping of the EA_Pfin1 locus.
Position on 2L
# of Alleles
Allele Linked to High Oocyst #
Allele Linked to Low Oocyst #
Average # of Oocysts with Susceptible Allele
Average # of Oocysts with Resistant Allele
Parasite strain transcendence. Comparison of nominal and genome wide p values for linkage at H796 before and after removing mosquitoes fed on a parasite isolate that did not infect any mosquitoes in the pedigree.
Only infective parasites
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].
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.
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.