AJTMH Transactions of the Royal Society of Tropical Medicine and Hygiene
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am. J. Trop. Med. Hyg., 73(4), 2005, pp. 749-752
Copyright © 2005 by The American Society of Tropical Medicine and Hygiene

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by MARSHALL, J. C.
Right arrow Articles by CACCONE, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by MARSHALL, J. C.
Right arrow Articles by CACCONE, A.
Related Collections
Right arrow Mosquitoes
Right arrow Vector Biology
Right arrow Malaria
Right arrow Medical Entomology

SHORT REPORT


PHYLOGENETIC RELATIONSHIPS OF THE ANTHROPOPHILIC PLASMODIUM FALCIPARUM MALARIA VECTORS IN AFRICA

JONATHON C. MARSHALL, JEFFREY R. POWELL*, AND ADALGISA CACCONE
Department of Ecology and Evolutionary Biology, and Yale Institute for Biospherics Studies–Molecular Systematics and Conservation Biology Laboratory, Yale University, New Haven, Connecticut

 

ABSTRACT

Malaria kills more than one million people a year, and understanding the historical association between its most notorious causative agent, Plasmodium falciparum, and its mosquito vectors is important in fighting the disease. We present a phylogenetic analysis of a number of species within the mosquito subgenus Cellia based on a selection of mitochondrial and nuclear genes. Although some of these relationships have been estimated in other studies, generally few species were included and/or statistical support at many nodes was low. Here we include two additional species of anthropophilic P. falciparum malaria vectors and reanalyze these relationships using a Bayesian method that allows us to simultaneously incorporate different models of evolution. We report data that indicate a paraphyletic relationship between five anthropophilic African mosquito vectors. Such a relationship suggests that these species can serve as independent natural experiments for anopheline immunologic responses to regular, prolonged contact with P. falciparum.


Malaria remains one of the world’s major health problems even after years of research and public health efforts. Malaria kills more than one million people each year and its affects cost sub-Saharan Africa an estimated US$12 billion in gross domestic product.1 Human malaria is caused by four separate species of parasitic protozoan within the genus Plasmodium. A number of anopheline mosquito species serve as disease vectors and infect humans with Plasmodium sporozoites during blood meals. One major problem in controlling malaria has been the enormous variation of transmission patterns.2 Transmission rates may vary widely in entomologic inoculation rates (0.01–1,000 infected bites per person per year) with contact occurring throughout the entire year or only during selected months.3 Additionally, these factors may vary from year to year and between geographically proximate localities.3 The most lethal form of malaria is caused by Plasmodium falciparum, which is particularly dangerous when carried by anthropophilic species of mosquitoes.

Five species of anopheline mosquitoes serve as the major vectors for P. falciparum in continental sub-Saharan Africa. All five species belong to the subgenus Cellia3,4 with the two most widely studied and most notorious being Anopheles gambiae and An. arabiensis. The three remaining species, An. funestus, An. nili, and An. moucheti, have also received significant recent attention as major vectors.5,6 These five species are the only anophelines in sub-Saharan Africa that frequently transmit P. falciparum malaria. However, it is not necessarily the case that other species of anophelines cannot transmit P. falciparum, but no other species are so anthropophilic and therefore encounter P. falciparum on a regular, if not constant, basis.

Prolonged, consistent contact between P. falciparum and these anthropophilic mosquitoes over long periods of time gives rise to distinct selective pressures on the organism’s immune systems. To use this in the fight against malaria, an accurate understanding of this association requires some knowledge of their shared history. For instance, if all five anthropophilic African mosquito species were each other’s closest relatives (i.e., form a monophyletic clade), one might assume speciation occurred from a highly anthropophilic ancestor, and thus the close association between mosquito, P. falciparum, and humans arose only once and was retained through the several speciation events. Conversely, one could infer that this close, prolonged association arose independently several times if phylogenetic studies show these mosquito species to be non-sister taxa separated by less anthropophilic and/or non-P. falciparum vectors. Thus, phylogenies can help us to understand the historic patterns of association between mosquito vector, P. falciparum, and humans.

The subgenus Cellia contains five taxonomic series with more than 20 species in each (except the Cellia Series).4 Some work has been done in investigating the relationships between and within series; however, much work is needed. For instance, most phylogenetic analyses are characterized by insufficient sampling within or between many series712 or poorly supported phylogenetic trees at the nodes of interest.9,12 We do not criticize these studies because their main objectives were generally not to test the monophyly and relationships of the taxonomic series within Cellia. In one of the more extensive studies, Sallum and others estimated phylogenetic relationships among 32 species of mosquitoes within the subfamily Anophelinae based on portions of the mitochondrial genes COI and COII, the nuclear 18S small subunit ribosomal DNA (rDNA), and the expansion D2 region of the nuclear large subunit 28S rDNA gene.12 In that study, they sampled nine species within the subgenus Cellia, the Asiatic species An. stephensi, An. dirus, An. farauti, An. minimus A, An. subpictus, and An. sundaicus, as well as the African species An. arabiensis, An. gambiae, and An. funestus. Similar to the Africa species, Asiatic species are known to be found infective with P. falciparum.1317 However, unlike the African species, few would be considered highly anthropophilic18 and regular carriers of P. falciparum3,19 and thus not subjected to a constant association with P. falciparum over evolutionary time. As one example, the anthropophilic behavior exhibited by An. minimus A is dependent on the availability of cattle and thus the contact with P. falciparum is sporadic at best.20

In our study, we have compiled a new data set by generating sequence data for two additional anthropophilic African mosquitoes, An. moucheti and An. nili, which we added to a subset of species from the study by Sallum and others.12 Our data set consists of all nine Cellia species from the study by Sallum and others,12 two additional species, An. moucheti and An. nili, and four outgroup species, Uranotaenia lowii, Aedeomyia squamipennis, An. marajoara, and An. atropos, which were selected from various regions of the best estimated phylogenetic tree of Sallum and others.12 GenBank numbers for these taxa are found in Sallum et al.12 Although the addition of only two species from separate taxonomic series (Figure 1Go) presents only a slightly more complete sampling from other studies, our study will be the first phylogenetic analysis of all five African P. falciparum vector species, as well as the first phylogenetic analysis specifically aimed at a rigorous test of monophyly and association patterns of these vectors through alternate topology tests and ancestral state reconstruction.



View larger version (40K):
[in this window]
[in a new window]
 
    FIGURE 1. Phylogenetic relationships for various mosquito species within in the subgenus Cellia estimated from mitochondrial DNA (mt DNA) (COI and COII) and ribosomal DNA (rDNA) (18S and 28S) sequence data. At each node, posterior probabilities from the Bayesian analysis are located above the line and bootstrap values from the parsimony analysis are located below the line. Only bootstrap values and posterior probabilities ≥ 50% are shown. OG = outgroup taxa; An. = Anopheles; Aed. = Aedes; Ur. = Uranotaenia; * = close, continual association between members of the mosquito species and Plasmodium falciparum; {Delta} = likely origin of close mosquito-P. falciparum association. a, Taxonomic series within the subgenus Cellia with corresponding species used in this study. b, Phylogram representing the best phylogenetic estimate and branch lengths from a Bayesian analysis of the combined rDNA (18S and 28S) sequences, c, Parsimony bootstrap consensus tree illustrating the low phylogenetic signal for the mtDNA sequence. d, Bayesian tree estimated from a data set containing all mtDNA and rDNA sequence data.

 
Sequence data were generated as in the study by Sallum and others.12 GenBank accession numbers for the sequences are DQ069719–DQ069726. All sequences were aligned in Clustal X 1.8121 and adjusted by visual inspection in Mac-Clade 4.06.22 Mitochondrial genes were translated into amino acids to aid in alignment and then trimmed to newly generated sequence lengths (666 basepairs for COI and 603 base-pairs for COII). We took a conservative approach to the alignment of the nuclear rDNA and used only sequence that was easily and unambiguously aligned. This reduced the total number of characters used for each gene to 475 for 18S and 197 for 28S. Modeltest23 was used to select models of evolution (under the Akaike Information Criterion) for the following data sets (selected models of nucleotide evolution are in parentheses): 18S gene (TIM + I + G), 28S gene (TrN + I), 18S + 28S genes (TrN + I + G), first codon position of the COI gene (TrN + I + G), second codon position of the COI gene (TVM + I), third codon position of the COI gene (TrN + I + G), first codon position of the COII gene (GTR + I), second codon position of the COII gene (TrN + I), and third codon position of the COII gene (TIM + G). Phylogenetic relationships were estimated using unweighted parsimony (MP) and Bayesian analyses. For the parsimony analyses, tree space was searched in PAUP*24 using heuristic searches with 200 random additions and the tree bisection-reconnection branch swapping method. Nodal support was estimated using 1,000 bootstrap replicates. For Bayesian analyses, combined data sets were partitioned by gene and/or codon position and estimated in MRBAYES.25 Bayesian tree estimates were generated from 3,000,000 generations with burn-ins ranging from 10,000 to 50,000 generations (depending on the data set). Phylogenetic estimates were done for each gene separately, all mitochondrial genes combined, all rDNA genes combined, and lastly, for all genes combined.

Figure 1Go shows to which taxonomic series within the subgenus Cellia each species belongs, as well as summarizes the most important phylogenetic results of our study. We believe that the phylogenetic tree based on a combined 18S and 12S data set (Figure 1BGo) represents the best estimate of relationships. Figure 1BGo shows An. nili, An. farauti, and An. dirus forming a monophyletic Neomyzomyia clade (posterior probability [PP] = 72), as well as An. funestus, An. moucheti, and An. minimus forming a monophyletic Myzomyia clade (PP = 100, bootstrap [BP] value = 71). Similar to Sallum and others,12 we found little overall phylogenetic signal within the COI and COII mitochondrial genes (Figure 1CGo), with exception of some weak support (BS value = 54) for the monophyly of the Pyretophorus Series (which was paraphyletic in the rDNA analysis), and strong support (BS value = 100) for sister-taxon status of An. arabiensis with An. gambiae and An. subpictus with An. sundaicus. The mitochondrial DNA (mtDNA) data also supported (BS value = 100) the monophyly of all anopheline taxa. However, we should also note the surprising sister relationship recovered between An. marajoara (an outgroup taxon) and An. dirus. This result and lack of resolution call into question the utility of the COI and COII data sets for estimating the phylogenetic relationships within the subgenus Cellia.

Combining data sets can recover hidden signal because, unlike phylogenetic noise, true signal is not randomly distributed across data sets.26 We tested for hidden signal by performing a combined analysis with all available data. Figure 1DGo shows the results of our combined analysis. Not surprisingly, the addition of 1,269 characters (from a total of 1941 characters) from the COI and COII data sets caused the Pyretophorus Series to form a monophyletic clade. However, the addition of these data broke apart the monophyly of both the Myzomyia and Neomyzomyia Series. The mtDNA sequence makes up the majority (65%) of the combined data set and given the poor overall signal contained in mitochondrial sequences, we tend to give more weight to the relationships estimated with the rDNA data set. The exception may be the monophyly of the Pyretophorus Series because this relationship was strongly supported by the mtDNA and further investigation, possibly with different genes, is needed to clarify this question. To date, several studies, including this one, have supported the monophyly of most taxonomic series,7,8,12 whereas other studies have supported the paraphyly of various series.7,9 A much better sampling of both taxa and genes is needed to truly test the monophyly of these taxonomic series.

We have presented best phylogenetic estimates for the relationships between some species within the subgenus Cellia (Figure 1Go). However, a best phylogenetic estimate does not always equate to a significantly better phylogenetic estimate when compared with an alternative topology.27 We used maximum likelihood methods to test our best estimate (Figure 1BGo) against some possible alternative topologies. We generated three alternative topologies by estimating likelihood trees (using the TrN + I + G model of evolution) from the combined rDNA data sets with the following constraints: 1) all highly anthropophilic African species form a monophyletic clade, 2) all species from the Pyretophorus Series form a monophyletic clade, and 3) all relationships illustrated by Sallum and others12 (their Figure 5) are maintained. Likelihood scores for our best estimate tree and the three constrained trees were as follows: our best estimate tree (which was identical to the tree estimated by Bayesian methods) = –1,839.29, the African monophyly tree = –1,857.05, the Pyretophorus monophyly tree = –1841.90, and the relationship tree of Sallum and others = –1845.1121. We performed the conservative Shimodaira-Hasegawa test28 to compare each of the alternative constrained trees with our best estimate. This test showed that our best estimate tree was preferred to the alternative at the following P values: African monophyly tree, P < 0.05; Pyretophorus monophyly, P = 0.20; and relationships of Sallum and others, P = 0.12. Due to the conservative nature of this test, we are confident that our best estimate tree is significantly better than the African monophyly tree and slightly better than the estimate of Sallum and others. However, our analysis only provides minor evidence against the monophyly of the Pyretophorus Series and as such further studies are needed.

This study shows that the constant association between P. falciparum and its African mosquito vectors probably arose independently several different times in the past. However, it is also possible that this association arose in a common ancestral species and was then lost or altered (a shift toward more zoophilic or exophagic behavior) in some of the later evolved species. To test between these two alternative hypotheses, we reconstructed the most parsimonious character states (highly anthropophilic and encounter P. falciparum on a constant basis versus limited P. falciparum contact) for ancestral taxa based on our best estimate tree using Mesquite 1.05.29 Weighting both character states equally, we found that the regular association between mosquito and P. falciparum has arisen two or three times in the past with equal probability. However, if we assume that gaining the association is even slightly more probable than escaping it, it may have occurred as many as three or four times independently. Also, An. funestus belongs to a group of at least five morphologically similar species (not included in this study) of which An. funestus in the only human Plasmodium vector,3 indicating that P. falciparum acquisition in An. funestus and An. moucheti were probably independent events. Figure 1BGo shows the preferred scenario (four acquisitions) and the likely positions on the tree branches where long-term associations with P. falciparum likely began.

The conclusion that the five African vectors began regular, constant encounters with P. falciparum at least three times independently is not surprising. Most researchers are aware these vectors are member of taxonomic series and species complexes with non-vector species. Our study provides addition evidence, within a more rigorous framework of statistical tests, that the close association of P. falciparum and several of its mosquito vectors has occurred multiple times in the evolutionary past. This is important because as such, each instance will serve as a natural replicate experiment. Since each African vector species has gone through thousands of years of co-evolution with P. falciparum independently, researchers can contrast immune responses in the vectors and determine which immune genes have responded most to this association and possibly which are most closely related to resistance. This knowledge could be very important in the construction of genetically modified mosquitoes and in the battle against malaria.


Received February 5, 2005. Accepted for publication May 24, 2005.

Acknowledgments: We thank Ralph E. Harbach for useful discussions on phylogenetic relationships among anopheline mosquito species; Vincent Medjibe for supplying the An. nili and An. moucheti samples collected in Bayanga, Central African Republic; and Carlo Costantini (Centre National de Recherche et Formation Sur le Paludisme, Ougadouguo, Burkina Faso) for identifying them. We also thank three anonymous reviewers for their very useful and insightful comments on an earlier draft of this report.

Financial support: Jonathon C. Marshall is supported by a National Institutes of Health (NIH) multidisciplinary parasitology training grant (PI: Dr. Diane McMahon-Pratt, Yale School of Public Health #2T32AI07404). Jeffrey R. Powell is supported by NIH grant RO1 AI46018.

* Address correspondence to Jeffrey R. Powell, Department of Ecology and Evolutionary Biology, Yale University, PO Box 8105, New Haven, CT 06520-8105. E-mail: jeffrey.powell{at}yale.edu Back

Authors’ addresses: Jonathon C. Marshall, Department of Ecology and Evolutionary Biology, Yale University, PO Box 8105, New Haven, CT 96520-8105, Telephone: 203-432-3886, Fax: 203-432-6066, E-mail: jonathan.marshall{at}yale.edu. Jeffrey R. Powell, ESC 170, Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT 96520-8105, Telephone: 203-432-3887, Fax: 203-432-6066, E-mail: jeffrey.powell{at}yale.edu. Adalgisa Caccone, ESC 140, Yale Institute for Biospherics Studies–Molecular Systematics and Conservation Biology Laboratory, Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT 96520-8105, Telephone: 203-432-5259, Fax: 203-432-7394, E-mail: adalgisa.caccone{at}yale.edu.

 

REFERENCES

  1. Mocumbi P, 2004. Plague of my people. Nature 430: 925.[Medline]
  2. Coluzzi M, 1984. Heterogeneities of the malaria vectorial system in tropical Africa and their significance in malaria epidemiology and control. Bull World Health Organ 62: 107–113.
  3. Fontenille D, Simard F, 2004. Unraveling complexities in human malaria transmission dynamics in Africa through a comprehensive knowledge of vector populations. Comp Immunol Microbial Infect Dis 27: 357–375.
  4. Harbach RE, 2004. The classification of genus Anopheles (Diptera: Culicidae): a working hypothesis of phylogenetic relationships. Bull Entomol Res 94: 537–553.[Medline]
  5. Antonio-Nkondjio C, Awono-Ambene P, Toto JC, Meunier JY, Zabaze-Kemleu S, Nyambam R, Wondji CS, Tchuinkam T, Fontenille D, 2002. High malaria transmission intensity in a village close to Yaounde, the capital city of Cameroon. J Med Entomol 39: 350–355.[ISI][Medline]
  6. Fontenille D, Lochouran L, Diagne N, Sokhna C, Lemasson JJ, Diatta M, Konate L, Faye F, Rogier C, Trape JK, 1997. High annual and seasonal variation in malaria transmission by anophelines and vector species composition in Dielmo, a holoendemic area of Senegal. Am J Trop Med Hyg 56: 247–253.
  7. Foley DH, Bryan JH, Yeates D, Saul A, 1998. Evolution and systematics of Anopheles: insights from a molecular phylogeny of Australasian mosquitoes. Mol Phylogenet Evol 9: 262–275.[ISI][Medline]
  8. Anthony TG, Harbach RE, Kitching IJ, 1999. Phylogeny of the Pyretophorus Series of Anopheles subgenus Cellia (Diptera: Culicidae). Syst Entomol 24: 193–205.
  9. Sallum MAM, Schultz TR, Wilkerson RC, 2000. Phylogeny of Anophelinae (Diptera Culicidae) based on morphological characters. Ann Entomol Soc Am 93: 745–775.
  10. Krzywinski J, Wilkerson RC, Beasansky NJ, 2001. Evolution of mitochondrial and ribosomal gene sequences in Anophelinae (Diptera: Culicidae): implications for phylogeny reconstruction. Mol Phylogenet Evol 18: 479–487.[Medline]
  11. Krzywinski J, Wilkerson RC, Beasansky NJ, 2001. Toward understanding Anophelinae (Diptera, Culicidae) phylogeny: insights from nuclear single-copy genes and the weight of evidence. Syst Biol 50: 540–556.[ISI][Medline]
  12. Sallum MAM, Schultz TR, Foster PG, Aonstein K, Wirtz RA, Wilkerson RC, 2002. Phylogeny of anophelinae (Diptera: Culicidae) based on nuclear ribosomal and mitochondrial DNA sequences. Syst Entomol 27: 361–382.
  13. Prakash A, Bhattacharyya DR, Mohapatra PK, Mahanta J, 2004. Role of the prevalent Anopheles species in the transmission of Plasmodium falciparum and P. vivax in Assam state, northeastern India. Ann Trop Med Parasitol 98: 559–568.[Medline]
  14. Cooper RD, Frances SP, 2002. Malaria vectors on Buka and Bougainville islands, Papua New Guinea. J Am Mosq Control Assoc 18: 100–106.[Medline]
  15. Fryauff DJ, Leksana B, Masbar S, Wiady I, Sismadi P, Susanti AI, Nagesha HS, Syafruddin A, Atmosoedjono S, Bangs MJ, Baird JK, 2002. The drug sensitivity and transmission dynamics of human malaria on Nias Island, North Sumatra, Indonesia. Ann Trop Med Parasitol 96: 447–462.[ISI][Medline]
  16. Sahu SS, 1998. Comparative susceptibility of Anopheles subpictus from fresh and brackish water areas to Plasmodium falciparum infection. Acta Trop 70: 1–7.[Medline]
  17. Sumodan PK, Kumar A, Yadav RS, 2004. Resting behavior and malaria vector incrimination of Anopheles stephensi in Goa, India. J Am Mosq Control Assoc 20: 317–318.[Medline]
  18. Verandi A, Donnelly MJ, Rowland M, Townson H, 2002. Isolation and characterization of microsatellite loci in the mosquito Anopheles stephensi Liston (Diptera: Culicidae). Mol Ecol Notes 2: 488–490.
  19. Fontenille D, Lochouarn L, 1999. The complexity of the malaria vectorial system in Africa. Parassitologia 41: 267–271.[Medline]
  20. Trung HD, Van Bortel W, Sochantha T, Keokenchanh K, Briet OJT, Coosemans M, 2005. Behavioral heterogeneity of Anopheles species in ecologically different localities in Southeast Asia: a challenge for vector control. Trop Med Int Health 10: 251–262.[Medline]
  21. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG, 1997. The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 24: 4876–4882.
  22. Maddison DR, Maddison WP, 2000. MacClade 4. Analysis of Phylogeny and Character Evolution. Sunderland, MA: Sinauer Associates.
  23. Posada D, Crandall KA, 1998. ModelTest: testing the model of DNA substitution. Bioinformatics 14: 817–818.[Abstract/Free Full Text]
  24. Swofford DL, 1999. PAUP: Phylogenetic Analysis Using Parsimony (* and Other Methods). Version 4. Sunderland, MA: Sinauer Associates.
  25. Huelsenbeck JP, Ronquist F, 2001. MrBayes: Bayesian inference of phylogeny. Bioinformatics 17: 754–755.[Abstract/Free Full Text]
  26. Gatesy J, O’Grady P, Baker R, 1999. Corroboration among data sets in simultaneous analysis: hidden support for phylogenetic relationships among higher level artiodactyl taxa. Cladistics 15: 271–313.
  27. Goldman N, Anderson JP, Rodrigo AG, 2000. Likelihood-based tests of topologies in phyloegenetics. Syst Biol 49: 652–670.[ISI][Medline]
  28. Shimodaira H, Hasegawa M, 1999. Multiple comparisons of log-likelihoods with applications to phylogenetic tree inference. Mol Biol Evol 16: 1114–1116.[ISI]
  29. Maddison WP, Maddison DR, 2004. Mesquite: A Modular System for Evolutionary Analysis. Version 1.05. Available from http://mesquiteproject.org




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by MARSHALL, J. C.
Right arrow Articles by CACCONE, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by MARSHALL, J. C.
Right arrow Articles by CACCONE, A.
Related Collections
Right arrow Mosquitoes
Right arrow Vector Biology
Right arrow Malaria
Right arrow Medical Entomology


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS