African origin of the malaria parasite Plasmodium vivax

Abstract

Plasmodium vivax is the leading cause of human malaria in Asia and Latin America but is absent from most of central Africa due to the near fixation of a mutation that inhibits the expression of its receptor, the Duffy antigen, on human erythrocytes. The emergence of this protective allele is not understood because P. vivax is believed to have originated in Asia. Here we show, using a non-invasive approach, that wild chimpanzees and gorillas throughout central Africa are endemically infected with parasites that are closely related to human P. vivax. Sequence analyses reveal that ape parasites lack host specificity and are much more diverse than human parasites, which form a monophyletic lineage within the ape parasite radiation. These findings indicate that human P. vivax is of African origin and likely selected for the Duffy-negative mutation. All extant human P. vivax parasites are derived from a single ancestor that escaped out of Africa.

ARTICLE

Received 11 Sep 2013 | Accepted 29 Jan 2014 | Published 21 Feb 2014

Weimin Liu1, Yingying Li1, Katharina S. Shaw2, Gerald H. Learn1, Lindsey J. Plenderleith3, Jordan A. Malenke1, Sesh A. Sundararaman1,4, Miguel A. Ramirez1, Patricia A. Crystal1, Andrew G. Smith1, Frederic Bibollet-Ruche1, Ahidjo Ayouba5, Sabrina Locatelli5, Amandine Esteban5, Fatima Mouacha5, Emilande Guichet5, Christelle Butel5, Steve Ahuka-Mundeke5,6, Bila-Isia Inogwabini7, Jean-Bosco N. Ndjango8, Sheri Speede9, Crickette M. Sanz10,11, David B. Morgan11,12, Mary K. Gonder13, Philip J. Kranzusch14, Peter D. Walsh15, Alexander V. Georgiev16,w, Martin N. Muller17, Alex K. Piel15,18, Fiona A. Stewart15, Michael L. Wilson19, Anne E. Pusey20, Liwang Cui21, Zenglei Wang21, Anna Farnert22, Colin J. Sutherland23, Debbie Nolder23, John A. Hart24, Terese B. Hart24,Paco Bertolani25, Amethyst Gillis26, Matthew LeBreton26, Babila Tafon27, John Kiyang28, Cyrille F. Djoko26, Bradley S. Schneider26, Nathan D. Wolfe26, Eitel Mpoudi-Ngole29, Eric Delaporte5, Richard Carter30,31, Richard L. Culleton32, George M. Shaw1,4, Julian C. Rayner33, Martine Peeters5, Beatrice H. Hahn1,4 & Paul M. Sharp3,31

Plasmodium vivax is the leading cause of human malaria in Asia and Latin America but is absent from most of central Africa due to the near xation of a mutation that inhibits the expression of its receptor, the Duffy antigen, on human erythrocytes. The emergence of this protective allele is not understood because P. vivax is believed to have originated in Asia. Here we show, using a noninvasive approach, that wild chimpanzees and gorillas throughout central Africa are endemically infected with parasites that are closely related to human P. vivax. Sequence analyses reveal that ape parasites lack host specicity and are much more diverse than human parasites, which form a monophyletic lineage within the ape parasite radiation. These ndings indicate that human P. vivax is of African origin and likely selected for the Duffy-negative mutation. All extant human P. vivax parasites are derived from a single ancestor that escaped out of Africa.

1 Department of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA. 2 Department of Biological Sciences, Columbia University, New York, New York 10027, USA. 3 Institute of Evolutionary Biology, University of Edinburgh, Edinburgh EH9 3JT, UK. 4 Department of Microbiology, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA. 5 Unit Mixte Internationale 233 Institut de Recherche pour le Dveloppement and University of Montpellier 1, 34394 Montpellier, France. 6 Institut National de Recherche Biomedicale, BP 1197 Kinshasa, Democratic Republic of the Congo. 7 Department of Aquatic Sciences and Assessment, Swedish University of Agricultural Sciences, SE-75007 Uppsala, Sweden. 8 Department of Ecology and Management of Plant and Animal Resources, Faculty of Sciences, University of Kisangani, BP 2012 Kisangani, Democratic Republic of the Congo. 9 Sanaga-Yong Chimpanzee Rescue Center, IDA-Africa, Portland, Oregon 97204, USA. 10 Department of Anthropology, Washington University, St Louis, Missouri 63130, USA. 11 Wildlife Conservation Society, Congo Program, BP 14537 Brazzaville, Republic of the Congo. 12 Lester E. Fisher Center for the Study and Conservation of Apes, Lincoln Park Zoo, Chicago, Illinois 60614, USA. 13 Department of Biological Sciences, University at Albany, State University of New York, Albany, New York 12222, USA. 14 Department of Molecular and Cell Biology, University of California, Berkeley, California 94720, USA. 15 Division of Biological Anthropology, University of Cambridge, Cambridge CB2 1QH, UK. 16 Department of Human Evolutionary Biology, Harvard University, Boston, Massachusetts 02115, USA. 17 Department of Anthropology, University of New Mexico, Albuquerque, New Mexico 87131, USA. 18 Department of Anthropology, University of California, San Diego, California 92093, USA. 19 Department of Anthropology, University of Minnesota, Minneapolis, Minnesota 55455, USA. 20 Department of Evolutionary Anthropology, Duke University, Durham, North Carolina 27708, USA. 21 Department of Entomology, Pennsylvania State University, University Park, Pennsylvania 16802, USA. 22 Infectious Disease Unit, Karolinska Institute, Department of Medicine Solna, Karolinska University Hospital, Stockholm SE-17176, Sweden. 23 Public Health England Malaria Reference Laboratory, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London WC1E 7HT, UK. 24 Lukuru Wildlife Research Foundation, Tshuapa-Lomami-Lualaba Project, BP 2012 Kinshasa, Democratic Republic of the Congo. 25 Leverhulme Centre for Human Evolutionary Studies, University of Cambridge, Cambridge, CB2 1QH, UK. 26 Global Viral Forecasting Initiative, San Francisco, California 94104, USA. 27 Ape Action Africa, BP 20072 Yaounde, Cameroon. 28 Limbe Wildlife Centre, PO Box 878 Limbe, Cameroon. 29Institut de Recherches Mdicales et dtudes des Plantes Mdicinales Prvention du Sida au Cameroun, Centre de Recherche Mdicale, BP 906 Yaound, Republic of Cameroon.

30 Institute of Immunology and Infection Research, University of Edinburgh, Edinburgh EH9 3JT, UK. 31 Centre for Immunity, Infection and Evolution, University of Edinburgh, Edinburgh EH9 3JT, UK. 32 Malaria Unit, Institute of Tropical Medicine, Nagasaki University, Nagasaki 852-8523, Japan. 33 Malaria Programme, Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, UK. w Present address: Institute for Mind and Biology, The University of Chicago,

Chicago, Illinois 60637, USA. Correspondence and requests for materials should be addressed to P.M.S. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2014 Macmillan Publishers Limited. All rights reserved.

DOI: 10.1038/ncomms4346

African origin of the malaria parasite Plasmodium vivax

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346

Of the ve Plasmodium species known to cause malaria in humans, P. vivax is the most widespread1. Although highly prevalent in Asia and Latin America, P. vivax is

thought to be absent from west and central Africa due to the near xation of a mutation that causes the Duffy-negative phenotype in indigenous African people1,2. The Duffy antigen receptor for chemokines (DARC) is used by P. vivax merozoites to invade red blood cells3. Since the absence of DARC on the surface of erythrocytes confers protection against P. vivax malaria, this parasite has long been suspected to be the agent that selected for this mutation2,4,5. However, this hypothesis has been difcult to reconcile with the proposed evolutionary origin of P. vivax68. The closest relative of P. vivax is believed to be P. cynomolgi9, which infects macaques. These two parasites form a lineage within a clade comprised of at least seven other Plasmodium species, all of which infect primates found only in Asia. The consensus view has thus been that P. vivax emerged in Southeast Asia following the cross-species transmission of a macaque parasite69. Under this scenario, the Duffy-negative mutation prevalent in west central African people was selected by another unidentied pathogen4, and its high frequency prevented P. vivax from entering central Africa.

Recently, P. vivax-like parasites have been identied in a limited number of African apes1013 and some mosquitoes (Anopheles species) trapped in their vicinity13,14. Molecular characterization of these parasites showed that they are very similar to, but apparently distinct from, human P. vivax13. These ndings raised the possibility that a sylvatic P. vivax reservoir exists in wild-living apes or other African primate species. However, since only very few of these sylvatic parasites have been identied, information concerning their geographic distribution, host species association, prevalence and relationship to humanP. vivax is lacking. In fact, most evidence of this natural P. vivax reservoir has come from pets and apes in wildlife rescue centres. Since captive apes can become infected with Plasmodium species that do not normally infect them in their natural habitat10,15, studies of wild-living populations are essential.

African apes are highly endangered and live in remote forest regions, rendering invasive screening for infectious agents both impractical and unethical. As an alternative, we have developed methods that permit the detection and amplication of pathogen-specic nucleic acids from ape faecal DNA1619. This approach enabled us to trace the origins of human immunodeciency virus type 1 (HIV-1) to chimpanzees (Pan troglodytes) in west central Africa18 and to identify the precursor of human P. falciparum in western gorillas (Gorilla gorilla)10. Here, we used a similar approach to investigate the molecular epidemiology of P. vivax infection in wild-living apes. Screening more than 5,000 faecal samples from 78 remote forest sites, we tested ape communities throughout central Africa. Since conventional PCR analysis is error prone and has the potential to confound phylogenetic analyses10, parasite sequences were generated using single-genome amplication (SGA), which eliminates Taq polymerase-induced recombination and nucleotide substitutions from nished sequences20,21.

In this study, we show that western (G. gorilla) and eastern gorillas (G. beringei) and chimpanzees (P. troglodytes), but not bonobos (P. paniscus), are endemically infected with P. vivax-like parasites, and that infection rates in wild ape communities are similar to those in human populations with stable parasite transmission. Analysing over 2,600 SGA-derived mitochondrial, apicoplast and nuclear sequences, we also show that ape parasites are considerably more diverse than human parasites and do not cluster according to their host species. In contrast, human parasites form a monophyletic lineage that falls within the radiation of the ape parasite sequences. These results indicate that

human P. vivax arose from within a Plasmodium species that infect wild-living chimpanzees and gorillas, and that all extant human P. vivax parasites evolved from a single ancestor that spread out of Africa.

ResultsHost species association and distribution of ape P. vivax. Using PCR primers designed to amplify P. vivax mitochondrial (mt) DNA, we screened 5,469 faecal specimens from ape communities sampled at 78 forest sites throughout sub-Saharan Africa (Fig. 1). Except for 196 samples from habituated apes, all other specimens were derived from non-habituated communities (Supplementary Table 1). Ape species and subspecies were identied by faecal mtDNA analysis10,18,19. A subset of specimens was also subjected to microsatellite analyses to estimate the number of sampled individuals (Supplementary Table 1). Targeting a 297-base-pair (bp) mtDNA fragment (Supplementary Fig. 1a), we foundP. vivax-like sequences in faecal DNA from western gorillas (G. gorilla), eastern gorillas (G. beringei) and chimpanzees (P. troglodytes), but not from bonobos (P. paniscus) (Table 1). Infections were most common in central chimpanzees (P. t. troglodytes) and western lowland gorillas (G. g. gorilla), with infected individuals identied at 76% of eld sites, including six locations where P. vivax was found in both of these species (Fig. 1). Ape P. vivax was also endemic in eastern chimpanzees (P. t. schweinfurthii) and eastern lowland gorillas (G. b. graueri), with infected apes documented at 38% of eld sites. Despite this wide geographic distribution (Fig. 1), the proportion of ape faecal samples that contained P. vivax-like sequences at any given eld site was low: among 2,871 chimpanzee and 1,844 gorilla samples that were analysed, only 45 and 32 were found to be PCR positive, respectively (Table 1). Correcting for specimen degradation and redundant sampling, and taking into account the sensitivity of the non-invasive diagnostic test, we estimated the proportion ofP. vivax-sequence-positive individuals for each eld site (Supplementary Table 1). The resulting values of 48% (Table 1) were lower than prevalence rates previously determined for P. falciparum-like (Laverania) parasites in wild apes10, but they were very similar to P. vivax parasite rates reported for endemically infected human populations1. In humans, point estimates of patent blood infection rarely exceed 7%, even in hyperendemic areas, and a parasite rate of greater than 1% indicates stable transmission1.

Since human P. vivax can induce dormant liver infections, we considered the possibility that ape parasite DNA might be excreted into faeces in the absence of a productive blood-stage infection and thus inate our infection rate estimates. To examine this, we compared the sensitivity of PCR-based parasite detection in blood and faecal samples from captive chimpanzees housed at a wildlife rescue centre (SY). Importantly, these chimpanzees were kept in outside enclosures immediately adjacent to the habitat of wild apes and were thus exposed to the same mosquito populations. Although blood and faecal samples were not matched, 11 of 48 chimpanzees (23%) were found to be P. vivax positive by blood analysis, as compared with 1 of 68 chimpanzees(1.5%) by faecal analysis (Supplementary Table 2). Thus, faecalP. vivax detection is considerably less sensitive than blood detection, most probably because of lower parasite loads, and may underestimate the number of infected apes by an order of magnitude. This likely explains our failure to detect ape P. vivax in wild-living NigeriaCameroonian chimpanzees (P. t. ellioti) and Cross River gorillas (G. g. diehli), for which only very few faecal samples were available (Table 1). Indeed, we subsequently conrmed P. vivax infection by blood analysis in ve Nigeria Cameroonian chimpanzees that were sampled in captivity

2 NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346 ARTICLE

Nigeria

Cameroon

South Sudan

TK MN

BK VM MD

Central African Republic

MK WC

KM

MF

YW LH

MG

NT

YB WE

MP

NL

DD DG

DP BQ

EK LM MM

LI

MT

SY

BD AN

MO

KG

AL

KO

UB

SL

GB

PO

MS

BY

DS

BM BB

NK LB GT

CP

MB

Democratic Republic of the Congo

KK

Congo

Lulonga

Salonga

IS GO

MU

EP BFKS

AZ

AM

PA

WL

DJ

BL

WB

Republic of the Congo

KA

BA

ON

LA

Uganda

KR

BJ

TL

WA

BI

OP

Gabon

LU

IK BN

BR

WK

KE

Rwanda

Tanzania

ML LK

Burundi

GM

UG

Figure 1 | Geographic distribution of P. vivax in wild-living apes. Field sites are shown in relation to the ranges of three subspecies of the common chimpanzee (P. t. ellioti, magenta; P. t. troglodytes, red; and P. t. schweinfurthii, blue), western (G. gorilla, yellow) and eastern (G. beringei, light green) gorillas, as well as bonobos (P. paniscus, green). Circles, squares and hexagons identify eld sites where wild-living chimpanzees, gorillas or both species were sampled, respectively. Ovals indicate bonobo sampling sites. Triangles denote the location of wildlife rescue centres (see Supplementary Table 1 for a list of all eld sites and their two-letter codes). Forested areas are shown in dark green, while arid and semiarid areas are depicted in yellow and brown, respectively. Major lakes and major rivers are shown in blue. Dashed white lines indicate national boundaries. Sites where ape P. vivax was detected are highlighted in yellow, with red lettering indicating that both chimpanzees and gorillas were infected.

Table 1 | Magnitude of the sylvatic P. vivax reservoir.

Species/subspecies Field sites tested*

Field sites positivew

Faecal samples tested

Faecal samples positivez

Infection rate (CI)y

Nigeria-Cameroon chimpanzee (P. t. ellioti) 14 0 126 0 0% (05%) Central chimpanzee (P. t. troglodytes) 25 11 1,130 25 8% (610%) Eastern chimpanzee (P. t. schweinfurthii) 28 10 1,615 20 4% (37%) Cross River gorilla (G. g. diehli) 2 0 80 0 0% (08%) Western lowland gorilla (G. g. gorilla) 22 14 1,575 30 7% (59%) Eastern lowland gorilla (G. b. graueri) 4 1 189 2 4% (19%) Bonobo (P. paniscus) 8 0 754 0 0% (01%)

CI, condence intervals.*Field sites are listed in Supplementary Table 1 and their locations are shown in Fig. 1.

wField sites where sylvatic P. vivax was found are highlighted in Fig. 1.

zFaecal samples were tested for P. vivax mitochondrial DNA by diagnostic PCR; all amplicons were sequence conrmed.yApe P. vivax infection rates were estimated based on the proportion of PCR-positive samples, correcting for specimen degradation, redundant sampling and the sensitivity of the diagnostic test. Brackets indicate 95% CI. Since faecal P. vivax detection is less sensitive than blood detection, the values represent minimum estimates.

(Supplementary Table 3), and although not tested in this study, western chimpanzees (P. t. verus) have previously been shown to carry this parasite in the wild12. Thus, all four subspecies of the common chimpanzee as well as western and eastern gorillas are infected with P. vivax, indicating the existence of a substantial sylvatic reservoir.

Given the widespread infection of both chimpanzee and gorilla populations, the fact that over 700 bonobo faecal samples from eight different collection sites were P. vivax negative came as a surprise (Fig. 1). Since wild-living bonobos also lackP. falciparum-related parasites10, yet are susceptible to infection with human P. falciparum in captivity21, this nding may reect a

paucity of transmitting mosquito vectors. In humans, a mutation (T to C, at position -33) in the GATA-1 transcription-factor-binding site within the promoter region of the DARC gene22 yields resistance to P. vivax infection, but sequence analysis of the same region in 134 ape samples, including 28 from bonobos, indicated that none had this substitution (Supplementary Fig. 2). In addition, all ape DARC genes analysed encoded a blood-group antigen with aspartic acid at amino acid 42, rather than the glycine found in the protective Fya allele in humans (Supplementary Fig. 2)23.

Finally, we asked whether other primates in central Africa might harbour P. vivax-like parasites. Using the same

NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346

P. vivax-specic PCR primers, we screened 998 blood samples from 16 Old World monkey species that had previously been collected for molecular epidemiological studies of simian immunodeciency viruses24,25. Testing samples from 11 different locations in southern Cameroon and the western parts of the Democratic Republic of the Congo (DRC), we failed to detect P. vivax infection in any of the animals tested (Supplementary Fig. 3). Although 501 of the 998 blood samples(50.2%) yielded a PCR amplicon, all of these represented Hepatocystis spp. infections as determined by sequence analysis (Supplementary Table 4). Thus, we found no evidence for aP. vivax reservoir in these African monkey species.

SGA of P. vivax sequences. To examine the evolutionary relationships of ape and human parasites, we amplied the completeP. vivax mitochondrial genome in three partially overlapping fragments (Supplementary Fig. 1a). This was done using SGA followed by direct amplicon sequencing, which eliminates Taq polymerase-induced recombination and nucleotide substitution errors, and provides a proportional representation of the parasite sequences that are present in vivo10,20. Alignment of these sequences revealed two single-nucleotide variants (SNVs) that distinguished all ape from all human parasites (Supplementary Fig. 1b; a third previously proposed SNV13 was polymorphic among the ape samples in our dataset). We thus designed PCR primers to amplify a fragment (fragment D) that included both SNVs on the same SGA amplicon (Supplementary Fig. 1). Although only a subset of P. vivax-positive faecal samples yielded this larger mtDNA fragment, we were able to generate fragment D sequences from 22 chimpanzees and 9 gorillas, 17 of which were sampled in the wild (Supplementary Table 3). Since most database sequences are derived by conventional PCR approaches, we also used SGA to amplify fragment D sequences from the blood of P. vivax-infected humans to produce Taq polymerase error-free sequences20. These samples included 94 international travellers, who had acquired P. vivax while visiting malaria endemic areas, as well as 25 P. vivax-infected individuals from China, Thailand, Myanmar and India, who sought treatment for clinical malaria (Supplementary Table 5), and thus provide a globally representative sample of human P. vivax infections.

Phylogenetic analysis of all SGA-derived P. vivax mtDNA sequences showed that the ape parasites formed two distinct clades (Fig. 2). One clade, represented by sequences from just two chimpanzee samples (termed BQptt392 and DGptt540; see legend of Fig. 2 for an explanation of sample nomenclature), was almost

MOptt51114

KApts1680

SYptt1

EKggg514

BQptt392

BYptt4890

EKggg1179

BQptt194

BQptt392

SYptt43

DGptt540

SYptt45

DGptt540

MOpte51058

KApts1680

DGptt540

LIpte50577

LIpte50577

SYptt31

DPptt143

EKggg1179

SYptt1

SYptt41

BApts1285

SYptt13

DPptt143

EKggg1179

SYpte56

MBptt323

LBggg1222

SYptt47

LIpte50579

SYptt17

LIpte50577

MSggg7183

EKggg1179

North Korean

Mauritania I

n = 4

n = 4 n = 2

n = 2

n = 2

n = 2

n = 50

n = 5

India VII

84 1 Salvador I

n = 2

n = 2

Brazil I

n = 9

n = 17

Figure 2 | Evolutionary relationships of ape and human P. vivax parasites in mitochondrial gene regions. The phylogenetic positions of mitochondrial fragment D (2,524 bp; Supplementary Fig. 1a) sequences from ape and human P. vivax strains are shown in relation to human and macaque parasite reference sequences. All sequences were generated by SGA20, except for human (Salvador I, India VII, Mauritania I, North Korean and Brazil I) and simian reference strains from the database (see Supplementary Tables 68 for GenBank accession numbers). Ape sequences are colour coded, with capital letters indicating the eld site (Fig. 1) and lower case letters denoting species and subspecies origin (ptt: P. t. troglodytes, red; pte: P. t. ellioti, orange; pts: P. t. schweinfurthii, blue; ggg: G. g. gorilla, green). Human sequences are depicted by haplotype (rectangles) and labelled according to their geographic origin in Oceania (light grey), Africa (white), South and Central America (black), South and Southeast Asia (striped) and the Middle East (dark grey). Haplotypes that include more than one sequence are indicated, with the numbers of sequences listed to the right. A second lineage of related parasite sequences from chimpanzee samples DGptt540 and BQptt392 likely represents a new Plasmodium species. The tree was inferred using maximum likelihood methods56. Numbers above and below nodes indicate bootstrap values (Z70%) and Bayesian posterior probabilities (Z0.95), respectively (the scale bar represents ve nucleotide substitutions).

n = 3

n = 3

SYptt32

EKggg1179

EKggg1179

EKggg1179

SYptt32

DGptt540

EKggg1179

EKggg1179

NKggg1167

BQptt392

EKggg514

DSggg86

LBggg1222

DSggg83

CPggg1426

EKggg1179

DSggg83

CPggg1444

5 nt

LBggg1222

86 1 LBggg1222

EKggg1179

LBggg1222

EKggg514

NKggg1167

EKggg1179

96 KGggg1246

1 NKggg1167

100

1

P. simiovale

P. cynomolgi

EKggg1179

P. fragile

P. knowlesi

4 NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346 ARTICLE

as divergent from the remaining ape and human parasites as were other Plasmodium species, and thus likely represents a previously unidentied species. All other ape parasite sequences were closely related to each other and to human P. vivax sequences, and thus appear to represent a single species (Fig. 2). Within this P. vivax clade, chimpanzee- and gorilla-derived sequences were interspersed, but all human-derived sequences formed a single well-supported lineage that fell within the radiation of the ape parasites. Inclusion of previously published non-SGA sequences conrmed this topology, although many of the database sequences exhibited long branches suggestive of PCR errors (Supplementary Fig. 4a). Interestingly, the one P. vivax sequence recently identied in a European traveller who became infected after working in a central African forest13 did not fall within the human P. vivax lineage, but clustered with parasites obtained from wild-living chimpanzees and gorillas (Supplementary Fig. 4a). This conrms the suspicion that this traveller acquired his infection by cross-species transmission from a wild ape.

To examine the robustness of these phylogenetic relationships, we selected additional genomic regions that had previously been used for evolutionary studies of P. vivax7,26. These included portions of the apicoplast caseinolytic protease C (clpC) gene as well as the nuclear genes encoding lactate dehydrogenase (ldh), adenylosuccinate lyase (asl), cell division cycle 2-related kinase (crk2) and b-tubulin (b-tub). Although fewer ape samples yielded amplication products, this was most likely due to lower copy numbers of apicoplast and nuclear genomes compared with mtDNA. Nonetheless, many samples yielded more than one P. vivax haplotype, indicating coinfection with multiple locally circulating variants (Supplementary Table 3). To increase the number of suitable reference sequences, we also amplied these same fragments from P. vivax-positive human samples, and for some gene regions from related macaque parasites. The resulting phylogenies yielded very similar topologies, with human parasites always forming a monophyletic lineage. Moreover, this lineage fell within the radiation of the P. vivax-like ape sequences for ve of the six loci tested (Fig. 3; Supplementary Figs 57). In contrast, chimpanzee and gorilla sequences were again interspersed, suggesting that P. vivax is often transmitted between the two ape hosts.

The apicoplast and nuclear gene sequences also conrmed the existence of the second, closely related, ape Plasmodium species. Sequences obtained for the clpC, ldh and crk2 genes (Fig. 3; Supplementary Fig. 6) from one or the other of the same two chimpanzee samples that yielded divergent mtDNA sequences (DGptt540 and BQptt392) were again clearly distinct fromP. vivax. In the mtDNA, ldh and 50 crk2 phylogenies, this new species was the closest relative of P. vivax (although strong support for this was only found in the latter two trees), while in the other trees relationships were not well resolved. While the new species appears to be rare, it was found in samples from two locations about 110 km apart. Moreover, its detection depended on the cross-reactivity of P. vivax-specic PCR primers, indicating that its prevalence and host association remain to be determined. Nonetheless, existing parasite sequences indicate that African apes harbour not only P. vivax, but also its closest relative.

Relative diversity of ape and human P. vivax. Since the phylogenetic analyses indicated that human P. vivax strains were derived from within the radiation of ape parasites, we expected the human strains to exhibit lower genetic diversity. To test this directly, we calculated the relative nucleotide diversity of SGA-derived ape and human P. vivax sequences. Indeed, values for the average number of nucleotide differences per site (p) were higher

for ape parasites than for human parasites at all loci tested. However, the extent of this increased diversity varied among genes (Table 2). For example, the diversity of mtDNA sequences was only 1.4 times higher among the ape parasites than among the human parasites, while for the apicoplast sequence this value was 6 times higher. For nuclear genes the ape sequences were 9 (asl, ldh) to 50 (crk2) times more diverse (Table 2). If the levels of diversity were in fact similar in the ape and human parasites, it would be most unlikely to observe this difference consistently across this number of loci (for example, MannWhitney U-test applied to the four nuclear loci, P 0.014).

The relative diversity of the various genes also differed between ape and human parasites. For example, for ape parasites, nuclear gene sequences were 1225 times more diverse than mtDNA, whereas this ratio was only 12 for human parasites (Table 2). In the absence of positive selection or demographic changes, relative diversities within species should be similar to those between species. Thus, for comparison, we calculated distances between orthologous sequences from P. vivax and the closely relatedP. cynomolgi9. For four of the ve loci, the relative diversity (scaled to the mtDNA value) among ape P. vivax strains was remarkably similar to the relative interspecic divergence (Table 2); the single exception was crk2, which was unusually diverse among the ape parasites (compared with other nuclear loci), but also the most conserved nuclear gene between species. In contrast, the diversity values for nuclear and apicoplast genes (relative to mtDNA) among human P. vivax strains were strikingly low (Table 2). This reduced diversity among the human P. vivax strains most likely reects a recent bottleneck that, depending on the composition of the founder population, could have affected the relative diversity of organelle and nuclear genomes differently.

DiscussionOur nding that wild-living apes in central Africa show widespread infection with diverse strains of P. vivax and harbour a distinct but related Plasmodium species provides new insight into the evolutionary history of human P. vivax, and potentially solves the paradox that a mutation conferring resistance toP. vivax occurs at high frequency in the very region where this parasite is absent. These results indicate that human P. vivax arose from within a Plasmodium species that infects chimpanzees and gorillas and indicate an origin in Africa rather than, as previously assumed, in Asia68. One interpretation of the phylogenies is that a single host switch from apes gave rise to human P. vivax, analogous to the origin of humanP. falciparum10. However, this seems unlikely in this case since ape P. vivax does not divide into gorilla- and chimpanzee-specic lineages and since humans are susceptible to both natural13 and experimental27 ape P. vivax infections. Thus, a more plausible interpretation is that an ancestral P. vivax stock was able to infect humans, gorillas and chimpanzees in Africa until the Duffy-negative mutation started to spread (perhaps around 30,000 years ago28) and eliminated P. vivax from humans there. Under this scenario, extant human-infecting P. vivax represents a bottlenecked lineage that survived after spreading out of Africa. Much more recently, concomitant with host migrations, humanP. vivax has been reintroduced to Africa29.

Several alternative scenarios for the origins of ape and human

P. vivax have recently been discussed13, but none of these seems plausible in light of the present data. All previous models assumed that P. vivax originated in humans in Asia, following the cross-species transmission of a monkey parasite, and that humans then brought the parasite to African apes. This assumption has been based on the fact that the closest known relatives

NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346

a b

KApts1680

MOptt51114_SGA5.1

SYpte56_SGA20.6

NKggg1167_SGA2.3

KApts1680 KApts1680

SYpte56_SGA20.9

DGptt540

LIpte50577

LIpte50579

SYptt17_SGA5.2

SYptt1

DPptt143_SGA2.10

SYptt43_SGA2.5

SYptt47

DGptt540

SYptt41

DGptt540

LIpte50577

LIpte50579

MBptt323_SGA5.8

MOptt51114

MOptt51114

SYpte56

Salvador I

EKggg1179

DGptt540_SGA20.6

NKggg1167_SGA5.8

Hu_hap4

Hu_hap5

Brazil I

1 DGptt540

LBggg1222_SGA2.9

MOptt51114

SYptt1

SYptt17

EKggg514_SGA5.2

DGptt540

SYptt43

EKggg1179_SGA5.1

NKggg1167_SGA2.7

LIpte50579

SYptt1_SGA5.1

BQptt392_SGA5.3

Mauritania I

India VII

SYptt47_SGA5.4

n = 18

n = 3

n = 73

EKggg1179_SGA5.2

n = 3

EKggg1179_SGA5.8

Hu_hap3

92 1

EKggg1179_SGA5.3

Brazil I

North Korean

Hu_hap1

Mauritania I

EKggg1179

North Korean

SYptt43

SYptt1

DGptt540

LIpte50579

n = 3

n = 2

Salvador I

Hu_hap2

EKggg1179

DSggg83_SGA5.8

MOggg51106

EKggg1179

P. hylobati

EKggg1179

P. inui

EKggg1179

5 nt

1 nt

BQptt392_SGA5.5

P. simiovale

EKggg1179

1 EKggg1179

EKggg1179

P. coatneyi

P. knowlesi

99 1

100 1 EKggg1179

DGptt540

DGptt540

P. cynomolgi

DGptt540

P. fieldi

P. cynomolgi*

LIpte50577

EKggg1179

P. simiovale*

P. fragile

P. fragile*

P. knowlesi

Figure 3 | Evolutionary relationships of ape and human P. vivax parasites in nuclear and apicoplast gene regions. The phylogenetic positions of (a) lactate dehydrogenase (ldh) gene (711 bp) and (b) caseinolytic protease C (clpC) gene (574 bp) sequences from ape and human P. vivax strains are shown in relation to human and macaque parasite reference sequences. All sequences were generated by SGA20, except for human (Salvador I, India VII, Mauritania I, North Korean and Brazil I) and simian reference strains from the database (asterisks indicate SGA-derived ldh sequences for P. simiovale, P. cynomolgi and P. fragile; see Supplementary Tables 68 for GenBank accession numbers). Newly derived ape P. vivax sequences are labelled and colour coded as in Fig. 2. Human and simian reference sequences are shown in black. Human ldh haplotypes are depicted as described in Fig. 2. Related parasite sequences from chimpanzee samples DGptt540 (ldh) and BQptt392 (clpC) likely represent a new Plasmodium species. Trees were inferred using maximum likelihood methods56. Numbers above and below nodes indicate bootstrap values (Z70%) and Bayesian posterior probabilities (Z0.95), respectively (the scale bar represents ve and one nucleotide (nt) substitutions, respectively).

6 NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346 ARTICLE

Table 2 | Nucleotide diversity in human and ape P. vivax lineages.

P. vivax genome

Locus Length (bp) Human P. vivax* Ape P. vivax* Ape/human ratiow Distance toP. cynomolgiz

No. p ( 10 3)y Ratio|| No. p ( 10 3)y Ratio|| Distance Ratio||

Mitochondrial cox1; cytb 2,443 138 0.75 1.0 62 1.02 1.0 1.4 0.012 1.0 Nuclear ldh 679 114 1.58 2.1 42 14.07 13.8 8.9 0.161 13.5 Nuclear asl 838 97 1.45 1.9 21 12.78 12.5 8.8 0.197 16.5 Nuclear crk2 666 134 0.54 0.7 32 25.92 25.4 49.9 0.081 6.8 Nuclear b-tub 684 81 0.81 1.1 12 11.85 11.6 14.6 0.161 13.5

Apicoplast clpC 574 70 0.34 0.5 21 1.99 2.0 5.9 0.024 2.0

*All human parasite sequences were derived using SGA methods from a global sampling of P. vivax strains (Supplementary Table 5), except for ve completely sequenced reference genomes from the database (Salvador I, Brazil I, Mauritania I, North Korea and India VII). Ape P. vivax sequences were also derived by SGA from faecal and blood samples of wild-living and sanctuary apes (Supplementary Table 3).

wRatio of nucleotide diversity (p) values in ape and human parasites.

zAverage distance between P. vivax and the orthologous sequences from P.cynomolgi.yNucleotide diversity (p). Human and ape sequences were compared over the same length of sequence.

||Diversity (or distance) value expressed relative to that for mtDNA.

of P. vivax all seem to infect Asian primates69. However, we now show that chimpanzees harbour a Plasmodium species that is more closely related to P. vivax than are any of the Asian primate parasites. Thus, it is more parsimonious to assume that the common ancestor of these two species existed in Africa. How this lineage was introduced into African apes remains unknown; however, this appears to have occurred a long time before the origin of P. vivax.

To explain current levels of genetic diversity in ape and humanP. vivax strains, previous models invoking an Asian origin either require human P. vivax in Asia to have gone extinct, prior to repopulation from Africa, or necessitate Asian P. vivax to have gone through a bottleneck (of unknown cause)13. In contrast, the African origin model does not require such an ancestral (now extinct) human P. vivax population in Asia, but explains the reduced diversity of human parasites, as resulting from an out-of-Africa bottleneck as seen in P. falciparum30, and in humans themselves31. It has also been suggested that P. vivax is more likely to have spread from Asia to Africa, because human P. vivax strains in Asia seem to be the most diverse13,32 and because phylogeographic analyses indicated elevated migration rates from Asia (especially India) to Africa6. However, this is now explained more simply by the extinction of human P. vivax in Africa, which would have had high diversity, due to the spread of the Duffy-negative mutation. P. vivax strains currently infecting humans in Africa are indeed of Asian origin, but this reects a reintroduction and occurred only recently, perhaps with the peopling of Madagascar from Asia within the last few thousand years33.

If the origin of P. vivax had been due to transmission from macaques in southeast Asia, this would imply a convoluted evolutionary history, given the timescales that have been invoked. Estimates of the time of the last common ancestor of humanP. vivax are generally on the order of hundreds of thousands of years ago. For example, using mtDNA sequences this ancestor was estimated to have existed about 400,000 years ago32, and a more recent comparison of nuclear genome sequences suggested a similar date34. Modern humans are thought to have evolved in Africa and to rst have entered Asia no more than 60,000 years ago35. Thus, if the estimates of the timescale of the coalescence of human P. vivax lineages are correct, the recipient of the transmission from macaques must have been some earlier hominin species, and P. vivax must have diversied for a long time in that host before numerous lineages were transmitted to modern humans after they emerged from Africa.

The existence of a sylvatic P. vivax reservoir has public health implications. First, it solves the mystery of P. vivax infections in travellers returning from regions where 99% of the human

population is Duffy negative36,37. Second, it raises the possibility that humans living in close proximity of chimpanzees and gorillas may become infected by ape P. vivax. A recent study of individuals attending a health clinic in the Republic of Congo revealed that 10% carried antibodies specic for preerythrocytic stages of P. vivax, suggesting continuing exposure to P. vivax sporozoites from an unidentied source38. Since ape P. vivax is highly prevalent, especially in west central Africa, wild-living chimpanzees and gorillas could serve as an infection reservoir, especially in areas where an inux of Duffy-positive humans through commerce and travel coincides with increasing forest encroachment and ape habitat destruction. Although Duffy-negative individuals are generally protected from blood-stage infections, recent studies in Madagascar39 and Ethiopia40 have shown that P. vivax is not absolutely dependent on the Duffy receptor. It will thus be important to assess the potential of ape parasites to acquire this phenotype, once the underlying genetic determinants have been identied in human strains. The possibility that ape P. vivax may spread via international travel to countries where human P. vivax is actively transmitted should also be considered. Since ape P. vivax is much more diverse than human P. vivax (Figs 2 and 3; Table 2), it is potentially more versatile to escape treatment and prevention measures, especially if human and ape parasites were able to recombine. Given the documented propensity of P. vivax for host switching13,27, it seems important to screen Duffy-positive and -negative humans in west central Africa, as well as transmitting mosquito vectors, for the presence of ape P. vivax. Such studies are now possible through the development of molecular tools that distinguish ape from human P. vivax, which also permit the screening of faecal samples for liver-derived parasite DNA in the absence of patent blood infection. These results are necessary to inform malaria control and eradication efforts and to assess future human zoonotic risk.

Methods

Ape samples. A total of 5,469 faecal samples newly (TL site) or previously collected from wild-living chimpanzees (P. troglodytes), western gorillas (G. gorilla), eastern gorillas (G. beringei) and bonobos (P. paniscus) for molecular epidemio-logical studies of simian retroviruses18,19,4144 and Laverania parasites10 were selected for P. vivax screening. All specimens were derived from non-habituated apes, except for 170 samples from chimpanzees in Gombe National Park (GM), Tanzania, and 26 samples from gorillas in the Dzanga-Sangha Reserve (DS), Central African Republic, who were habituated to the presence of human observers. Samples were collected (1:1 vol/vol) in RNAlater (Life Technologies), transported at ambient temperatures and stored at 80 C. Faecal DNA was

extracted using the QIAamp Stool DNA Mini Kit (Qiagen, Valencia, CA) and used to amplify portions of the host mitochondrial genome to conrm species and subspecies origin18,19,4144. A subset was also subjected to microsatellite analysis at

NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346

48 polymorphic loci10,18,4144 to estimate the number of sampled individuals (Supplementary Table 1). In addition to faecal samples from wild populations, we also obtained stool and blood samples from sanctuary chimpanzees and gorillas that were kept in outside enclosures immediately adjacent to the habitat of wild apes. These included 113 faecal and 66 blood samples from chimpanzees housed at the Sanaga Yong Rescue Centre (SY), 2 faecal samples from gorillas and 14 blood samples from chimpanzees housed at the Limbe Wildlife Centre (LI) and 8 blood samples from 6 chimpanzees and 2 gorillas housed at the Mfou National Park Wildlife Rescue Centre (MO), all located in Cameroon. Faecal samples were collected from known individuals under direct observation. Blood samples were obtained by veni-puncture (dried blood spots, whole blood, buffy coats, red blood cells) and represented left-over material from routine health examinations or were collected for specic veterinary (diagnostic) purposes. Blood and faecal collections were approved by the Ministry of Forestry and Wildlife of Cameroon. Two chimpanzees at the SY sanctuary were suspected to suffer from malaria and sampled during severe febrile illnesses. One had positive blood smears on site and was subsequently identied to be PCR positive for P. reichenowi, while the other was PCR positive for ape P. vivax. Both responded to malaria treatment. A few additional chimpanzees who were PCR positive for either Laverania or ape P. vivax exhibited milder symptoms at or near the time of sampling, but the majority of captive apes, including several who were Laverania and/or non-Laverania sequence positive by blood analysis, were asymptomatic at the time of blood collection. These individuals were also blood smear negative. Samples were shipped in compliance with Convention on International Trade in Endangered Species of Wild Fauna and Flora regulations and country specic import and export permits. DNA was extracted from whole blood and dried blood spots using the QIAamp Blood DNA Mini Kit (Qiagen, Valencia, CA).

Monkey samples. To investigate the full host range of P. vivax, we screened blood samples from 998 non-human primates sampled in Cameroon and the Democratic Republic of the Congo (Supplementary Table 4). Collection and molecular characterization of these samples have been described24,25,45, except for samples from the BO, MO, MA and MS eld sites (Supplementary Fig. 3). Briey, samples were collected from primate bushmeat as dried blood spots or whole blood using strategies specically designed not to increase demand24,25,45. DNA was extracted using the QIAamp blood kit (Qiagen, Courtaboeuf, France) and the species origin was determined by amplifying a 386 bp mtDNA fragment spanning the 12S rRNA gene using single round PCR24,25. The non-human primate samples were obtained with approval from the Ministries of Health and Environment and the National Ethics Committees of Cameroon and the Democratic Republic of the Congo.

Human samples. To increase the number of human P. vivax reference sequences and to generate sequences devoid of PCR error, we obtained a global sampling of human P. vivax infections. These included dried blood spots or DNA samples from(i) 61 international travellers diagnosed at the Malaria Reference Laboratory of the London School of Hygiene and Tropical Medicine, London, UK (designated MRL),(ii) 35 international travellers diagnosed at the Department of Infectious Diseases, Karolinska University Hospital, Stockholm, Sweden (designated SW), and (iii) 32 residents of malaria endemic areas in China (designated GX; n 10), Myanmar

(V0; n 10), Thailand (PVAR; n 2) and India (IN; n 10), who sought treat

ment at local health clinics (Supplementary Table 5). In all cases, P. vivax infection was initially identied by microscopy and then conrmed by diagnostic PCR followed by direct amplicon sequencing. DNA was extracted from whole blood or dried blood spots using QIAamp DNA Mini Kit (Qiagen, Valencia, CA). Anon-ymised DNA samples previously collected by the MRL were provided under its remit to undertake epidemiological surveillance relevant to imported malaria in the UK. All other subjects provided written informed consent for the collection and analysis of samples, which were sent without patient identiers and other patient information, except for the country of (known or presumed) P. vivax acquisition. The study was approved by the Institutional Review Boards of the University of Pennsylvania, the Karolinska Institute and the Pennsylvania State University.

Other Plasmodium species. Genomic DNA from, or dried blood spots containing, P. inui (catalogue number MRA-486F), P. simiovale (MRA-488F), P. cynomolgi (MRA-350G) and P. fragile (MRA-352G) were obtained from the Malaria Research and Reference Reagent Resource Center (MR4) of the American Type Culture Collection (ATCC, Manassas, VA) to generate single-genome-derived reference sequences for phylogenetic analyses.

Microsatellite analyses. Ape faecal samples that have previously been subjected to microsatellite analyses have been reported10, except for bonobo samples obtained at the TL eld site, which were obtained more recently. Briey, faecal DNA was extracted and used to amplify eight polymorphic microsatellite loci. Amplication products were analysed on an automated sequencer (Applied Biosystems, CA) and sized using GeneMapper 4.0 (Applied Biosystems). For individual identication, samples were rst grouped by mtDNA haplotype. Within each haplotype, samples were then grouped by microsatellite genotypes. Since samples were partially degraded due to prolonged storage at ambient temperatures, we allowed allelic mismatches at up to six loci to guard against allelic dropout, but

only if they represented a missing allele. This conservative genotyping approach likely underestimates the number of individuals screened and thus reects a minimum estimate. Samples with evidence of DNA admixture (multiple peaks for the same locus) were discarded.

P. vivax diagnostic PCR. Faecal and blood samples were screened by conventional PCR for P. vivax mitochondrial (cox1) sequences as described46. Although originally designed to amplify P. vivax sequences from samples coinfected with Laverania species, the primers were also found to amplify Hepatocystis spp. mtDNA (Supplementary Table 4). Primers Pv2768p (50-GTATGGATCGAATC

TTACTTATTC-30) and Pv3287n (50-AATACCAGATACTAAAAGACCAACA ATGATA-30) were used in the rst round of PCR, while Pv2870p (50-TTGCA ATCATAAAACTTTAGGTC-30) and Pv3185n (50-TCCTCCAAATTCTGCTGCT GTAGATAAAATG-30) were used in the second round, yielding a 297 bp amplicon (Supplementary Fig. 1a). Amplicons were gel puried and sequenced directly to conrm Plasmodium infection.

SGA of P. vivax gene sequences. To derive P. vivax sequences devoid of PCR errors, including Taq polymerase-induced misincorporations and template switching, blood and faecal samples from P. vivax-infected apes and humans were subjected to SGA methods, essentially as described10,20. DNA was extracted and end point diluted such that fewer than 30% of PCR reactions yielded an amplication product. According to a Poisson distribution, a well yielding a PCR product at this dilution will contain only a single ampliable template more than 83% of the time. Multiple different gene regions were amplied, including mitochondrial fragments S, A, B, C and D (Supplementary Fig. 1a), the apicoplast clpC gene (574 bp), and the nuclear genes ldh (713724 bp), asl (838 bp), crk2 (666 bp) and b-tub (684 bp) (Supplementary Table 3). Primers were designed in conserved regions of the P. vivax genome, except for a set of ldh primers intended to amplify more diverse macaque parasites, which were designed using both P. vivax and P. knowlesi consensus sequences (Supplementary Table 9). PCR cycling conditions were as described10, except for varying elongation times depending on the length of the amplicon (1min per 1 kb amplicon length). Amplication products were sequenced directly and analysed using Sequencher (Gene Codes Corporation, Ann Arbor, MI). All sequences with double peaks in the chromatogram were discarded, except for sequences derived from ape faecal samples that contained a single double peak. Single double peaks in SGA sequences either indicate the presence of two P. vivax variants differing by a single nucleotide or a PCR misincorporation within the rst or second round of PCR. We resolved these ambiguous sites by selecting the nucleotide that was identical to the apeP. vivax consensus. This conservative approach avoided the loss of otherwise valuable sequence information. GenBank accession codes of newly derived SGA sequences are listed in Supplementary Tables 6 and 7.

DARC genotyping. To examine whether P. vivax-infected and -uninfected apes exhibited sequence polymorphisms in their DARC promoter and/or adjacent coding sequences, we amplied a 1,286 bp region of the DARC gene in two partially overlapping PCR fragments from 134 ape (faecal and blood) samples (Supplementary Fig. 2a). Fragment A was amplied using DARCpF1 (50-GCTG

TCCCATTGTCCCCTAG-30) and DARCpR8 (50-GGCCCCATACTCACCCT GTGC-30) in the rst round of PCR, and DARCpF3 (50-GCACAATGATACA CAGCAAAC-30) and FYPdn (50-CCATGGCACCGTTTGGTTCAGG-30)47 in the second round. Fragment B was amplied using DARC_F5 (50-AGGCAGTGGG CGTGGGGTAAG-30) and DARC_R5 (50-AGCCATACCAGACACAGTAG CC-30) in the rst round, and DARC_NTF1 (50-TTGGCTCTTATCTTGG AAGCAC-30) and DARC_NTR1 (50-TGGTGAGGATGAAGAAGGGCAGT-30) in the second round. PCR conditions were the same as those used to amplify P. vivax mtDNA. Amplicons were directly sequenced and analysed using Sequencher.

Ape P. vivax infection rates. For sites where the number of sampled chimpanzees was known (Supplementary Table 1), P. vivax infection rates were estimated based on the proportion of positive individuals (but correcting for test sensitivity; see below), with 95% condence limits determined assuming binomial sampling. For eld sites where the number of sampled individuals was not known, infection rates were estimated based on the number of faecal samples, but correcting for specimen degradation, oversampling and the sensitivity of the diagnostic PCR test. The latter was estimated by determining the proportion of PCR-positive specimens fromP. vivax-infected apes that were sampled more than once on the same day, as previously described10. Including data from 14 such apes (6 central chimpanzees, 5 eastern chimpanzees and 3 western gorillas), we estimated the sensitivity to be 32% (17 of 53 samples were positive; condence limits were determined assuming binomial sampling). Using previously reported values for oversampling10 (1.77,1.84, 3.74, 2.06 and 1.84 for central chimpanzees, western gorillas, eastern chimpanzees, eastern gorillas and bonobos, respectively) and sample degradation10(0.13), the proportion of P. vivax-sequence-positive apes was estimated for each eld site (Table 1 and Supplementary Table 1). However, since P. vivax detection in faecal samples is considerably less sensitive than in blood, the resulting prevalence rates should be interpreted as minimum estimates.

8 NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346 ARTICLE

Nucleotide diversity calculations. Newly generated ape and human P. vivax sequences were aligned with sequences from the complete genomes of ve P. vivax reference strains (Salvador I, India VII, Mauritania I, North Korean and Brazil I). Other than these ve reference strains, only SGA-derived sequences were included in the analysis, so as to not inate diversity values by including Taq polymerase-induced errors. Identical sequences from the same sample or individual were excluded; identical sequences from different samples or individuals were retained. For each gene, sequences were aligned using ClustalW48, and both ape and humanP. vivax alignments were trimmed to the same length. The nucleotide diversity (p) of human and ape parasite sequences was calculated for six P. vivax gene fragments (mtDNA fragment D, apicoplast clpC, nuclear ldh, asl, crk2 and b-tub) using

DnaSP version 5.10 (ref. 49).

Distance calculation from P. cynomolgi. For each P. vivax gene, the same set of sequences as used for the diversity calculation was aligned with one P. cynomolgi reference sequence9. The distance from each P. vivax sequence to the P. cynomolgi reference was calculated using the R package version 3.08 with the TamuraNei correction50,51, and the mean distance was calculated for each gene.

Phylogenetic analyses. Sequence alignments were constructed using ClustalW version 2.1 (ref. 48) and manually adjusted using MacClade52. Regions that could not be unambiguously aligned were omitted from subsequent phylogenetic analyses. Nuclear gene sequences were subjected to recombination analysis using GARD53. Evolutionary models for phylogenetic analyses were determined using the Akaike information criterion with Modeltest (version 3.7)54 and PAUP*(ref. 55). Maximum likelihood phylogenies with bootstrap support (100 replicates) were estimated jointly with model parameter values by means of PhyML (version3)56 using both nearest-neighbour interchange and subtree pruning and regrafting with Neighbor Joining and 10 random-addition starting trees57. Posterior probabilities for nodes in phylogenetic trees were calculated using MrBayes (version 3.2)58, using an average standard deviation of partition frequencies o0.01 as a convergence diagnostic. Trees were constructed from mitochondrial fragment

D sequences (2,524 bp, Fig. 2 and Supplementary Fig. 4a), apicoplast clpC sequences (574 bp, Fig. 3b), and nuclear b-tub (664 bp, Supplementary Fig. 5), ldh (711 bp, Fig. 3a and Supplementary Fig. 4b), crk2 (271 and 372 bp, Supplementary Fig. 6) and asl (838 bp, Supplementary Fig. 7) sequences.

References

1. Gething, P. W. et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl. Trop. Dis. 6, e1814 (2012).

2. Mendis, K., Sina, B. J., Marchesini, P. & Carter, R. The neglected burden of Plasmodium vivax malaria. Am. J. Trop. Med. Hyg. 64, 97106 (2001).

3. Miller, L. H., Mason, S. J., Clyde, D. F. & McGinniss, M. H. The resistance factor to Plasmodium vivax in blacks. The Duffy-blood-group genotype, FyFy. New Engl. J. Med. 295, 302304 (1976).

4. Livingstone, F. B. The Duffy blood groups, vivax malaria, and malaria selection in human populations: a review. Hum. Biol. 56, 413425 (1984).

5. Carter, R. & Mendis, K. N. Evolutionary and historical aspects of the burden of malaria. Clin. Microbiol. Rev. 15, 564594 (2002).

6. Mu, J. et al. Host switch leads to emergence of Plasmodium vivax malaria in humans. Mol. Biol. Evol. 22, 16861693 (2005).

7. Escalante, A. A. et al. A monkeys tale: the origin of Plasmodium vivax as a human malaria parasite. Proc. Natl Acad. Sci. USA 102, 19801985 (2005).

8. Carlton, J. M., Das, A. & Escalante, A. A. Genomics, population genetics and evolutionary history of Plasmodium vivax. Adv. Parasitol. 81, 203222 (2013).

9. Tachibana, S. et al. Plasmodium cynomolgi genome sequences provide insight into Plasmodium vivax and the monkey malaria clade. Nat. Genet. 44, 10511055 (2012).

10. Liu, W. et al. Origin of the human malaria parasite Plasmodium falciparum in gorillas. Nature 467, 420425 (2010).

11. Krief, S. et al. On the diversity of malaria parasites in African apes and the origin of Plasmodium falciparum from bonobos. PLoS Pathog. 6, e1000765 (2010).

12. Kaiser, M. et al. Wild chimpanzees infected with ve Plasmodium species. Emerg. Infect. Dis. 16, 19561959 (2010).

13. Prugnolle, F. et al. Diversity, host switching and evolution of Plasmodium vivax infecting African great apes. Proc. Natl Acad. Sci. USA 110, 81238128 (2013).

14. Paupy, C. et al. Anopheles moucheti and Anopheles vinckei are candidate vectors of ape Plasmodium parasites, including Plasmodium praefalciparum in Gabon. PLoS One 8, e57294 (2013).

15. Duval, L. et al. African apes as reservoirs of Plasmodium falciparum and the origin and diversication of the Laverania subgenus. Proc. Natl Acad. Sci. USA 107, 1056110566 (2010).

16. Santiago, M. L. et al. SIVcpz in wild chimpanzees. Science 295, 465 (2002).17. Santiago, M. L. et al. Foci of endemic simian immunodeciency virus infection in wild-living eastern chimpanzees (Pan troglodytes schweinfurthii). J. Virol. 77, 75457562 (2003).

18. Keele, B. F. et al. Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science 313, 523526 (2006).

19. Li, Y. et al. Eastern chimpanzees, but not bonobos, represent a simian immunodeciency virus reservoir. J. Virol. 86, 1077610791 (2012).

20. Liu, W. et al. Single genome amplication and direct amplicon sequencing of Plasmodium spp. DNA from ape fecal samples. Protoc. Exch. doi:http://dx.doi.org/10.1038/nprot.2010.156

Web End =10.1038/ http://dx.doi.org/10.1038/nprot.2010.156

Web End =nprot.2010.156 (2010).

21. Rayner, J. C., Liu, W., Peeters, M., Sharp, P. M. & Hahn, B. H. A plethora of Plasmodium species in wild apes: a source of human infection? Trends Parasitol. 27, 222229 (2011).

22. Tournamille, C., Colin, Y., Cartron, J. P. & Le Van Kim, C. Disruption of a GATA motif in the Duffy gene promoter abolishes erythroid gene expression in Duffy negative individuals. Nat. Genet. 10, 224228 (1995).

23. King, C. L. et al. Fya/Fyb antigen polymorphism in human erythrocyte Duffy antigen affects susceptibility to Plasmodium vivax malaria. Proc. Natl Acad. Sci. USA 108, 2011320118 (2011).

24. Aghokeng, A. F. et al. Extensive survey on the prevalence and genetic diversity of SIVs in primate bushmeat provides insights into risks for potential new cross-species transmissions. Infect. Genet. Evol. 10, 386396 (2010).

25. Ahuka-Mundeke, S. et al. Novel multiplexed HIV/simian immunodeciency virus antibody detection assay. Emerg. Infect. Dis. 17, 22772286 (2011).

26. Kedzierski, L. et al. Phylogenetic analysis of the genus Plasmodium based on the gene encoding adenylosuccinate lyase. Infect. Genet. Evol. 1, 297301 (2002).

27. Contacos, P. G. et al. Transmission of Plasmodium schwetzi from the chimpanzee to man by mosquito bite. Am. J. Trop. Med. Hyg. 19, 190195 (1970).

28. Hamblin, M. T. & Di-Rienzo, A. Detection of the signature of natural selection in humans: evidence from the Duffy blood group locus. Am. J. Hum. Genet. 66, 16691679 (2000).

29. Culleton, R. & Carter, R. African Plasmodium vivax: distribution and origins. Int. J. Parasitol. 42, 10911097 (2012).

30. Tanabe, K. et al. Plasmodium falciparum accompanied the human expansion out of Africa. Curr. Biol. 20, 12831289 (2010).

31. Ramachandran, S. et al. Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa. Proc. Natl Acad. Sci. USA 102, 1594215947 (2005).

32. Cornejo, O. E. & Escalante, A. A. The origin and age of Plasmodium vivax. Trends Parasitol. 22, 558563 (2006).

33. Tofanelli, S. et al. On the origins and admixture of Malagasy: new evidence from high-resolution analyses of paternal and maternal lineages. Mol. Biol. Evol. 26, 21092124 (2009).

34. Neafsey, D. E. et al. The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum. Nat. Genet. 44, 10461050 (2012).

35. Mellars, P. Going east: new genetic and archaeological perspectives on the modern human colonization of Eurasia. Science 313, 796800 (2006).

36. Carter, R. Speculations on the origins of Plasmodium vivax malaria. Trends Parasitol. 13, 214219 (2003).

37. Guerra, C. A. et al. The international limits and population at risk of Plasmodium vivax transmission in 2009. PLoS Negl. Trop. Dis. 4, e774 (2010).

38. Culleton, R. et al. Evidence for the transmission of Plasmodium vivax in the Republic of the Congo, West Central Africa. J. Infect. Dis. 200, 14651469 (2009).

39. Mnard, D. et al. Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people. Proc. Natl Acad. Sci. USA 107, 59675971 (2010).

40. Woldearegai, T. G. et al. Plasmodium vivax malaria in Duffy-negative individuals from Ethiopia. Trans. R Soc. Trop. Med. Hyg. 107, 328331 (2013).

41. Keele, B. F. et al. Increased mortality and AIDS-like immunopathology in wild chimpanzees infected with SIVcpz. Nature 460, 515519 (2009).

42. Rudicell, R. S. et al. High prevalence of simian immunodeciency virus infection in a community of savanna chimpanzees. J. Virol. 85, 99189928 (2010).

43. Neel, C. et al. Molecular epidemiology of simian immunodeociency virus infection in wild-living gorillas. J. Virol. 84, 14641476 (2009).

44. Etienne, L. et al. Non-invasive follow-up of simian immunodeciency virus infection in wild-living nonhabituated western lowland gorillas in Cameroon.J. Virol. 86, 97609772 (2012).45. Aghokeng, A. F. et al. Widely varying SIV prevalence rates in naturally infected primate species from Cameroon. Virology 345, 174189 (2006).

46. Sundararaman, S. A. et al. Plasmodium falciparum-like parasites infecting wild apes in southern Cameroon do not represent a recurrent source of human malaria. Proc. Natl Acad. Sci. USA 110, 70207025 (2013).

47. Zimmerman, P. A. et al. Emergence of FY*Anull in a Plasmodium vivax-endemic region of Papua New Guinea. Proc. Natl Acad. Sci. USA 96, 1397313977 (1999).

NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4346

48. Larkin, M. A. et al. Clustal W and Clustal X version 2.0. Bioinformatics 23, 29472948 (2007).

49. Librado, P. & Rozas, J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 14511452 (2009).

50. Paradis, E., Claude, J. & Strimmer, K. APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics 20, 289290 (2004).

51. Tamura, K. & Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 10, 512526 (1993).

52. Maddison, W. P. & Maddison, D. R. MacCladeAnalysis of Phylogeny and Character EvolutionVersion 4 (Sinauer Associates, Inc., 2001).

53. Kosakovsky Pond, S. L., Posada, D., Gravenor, M. B., Woelk, C. H. & Frost, S.D. Automated phylogenetic detection of recombination using a genetic algorithm. Mol. Biol. Evol. 23, 18911901 (2006).54. Posada, D. & Buckley, T. Model selection and model averaging in phylogenetics: advantages of akaike information criterion and bayesian approaches over likelihood ratio tests. Syst. Biol. 53, 793808 (2004).

55. PAUP* 4.0: Phylogenetic Analysis Using Parsimony (*and Other Methods) v.4.0b2a (Sinauer Associates, Inc., Sunderland, MA, 1999).56. Guindon, S., Delsuc, F., Dufayard, J. F. & Gascuel, O. Estimating maximum likelihood phylogenies with PhyML. Methods Mol. Biol. 537, 113137 (2009).

57. Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307321 (2010).

58. Ronquist, F. & Huelsenbeck, J. P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19, 15721574 (2003).

Acknowledgements

We thank the staff of project PRESICA, the World Wildlife Fund for Nature (WWF/ DRC), the Institut National de Recherches Biomdicales (INRB, Kinshasa, DRC), Global Viral Cameroon, the Bonobo Conservation Initiative and Vie Sauvage, as well as Didier Mazongo, Octavie Lunguya, Muriel Aloni and Valentin Mbenz for eldwork in Cameroon and the DRC; eld assistants in Gombe National Park and Ugalla for sample collection in Tanzania, the staff of the Sanaga Yong, Limbe and Ape Action Africa/Mfou National Park Wildlife Rescue Centres for collecting blood samples from captive apes; the Cameroonian Ministries of Health, Forestry and Wildlife, and Scientic Research and Innovation for permission to collect samples in Cameroon; the Water and Forest Ministry for permission to collect samples in the Central African Republic; the Ministry of Forest Economy and Sustainable Development for permission to collect samples in the Republic of Congo; the Ministry of Scientic Research and Technology, the Department of Ecology and Management of Plant and Animal Resources of the University of Kisangani, the Ministries of Health and Environment and the National Ethics Committee for permission to collect samples in the DRC; the Tanzania Commission for Science and

Technology and the Tanzania Wildlife Research Institute for permission to conduct research in Gombe National Park and Ugalla. This work was supported by grants from the National Institutes of Health (R01 AI091595, R37 AI050529, R01 AI58715, T32 AI007532, P30 AI045008), the Agence Nationale de Recherche sur le Sida (ANRS 12125/ 12182/12255), the Agence Nationale de Recherche (Programme Blanc, Sciences de la Vie, de la Sant et des Ecosystmes and ANR 11 BSV3 021 01, Projet PRIMAL), Harvard University, the Arthur L. Greene Fund, the Jane Goodall Institute, the Wellcome Trust (098051), the Leakey Foundation, Google.org and the Skoll Foundation. This study was also made possible by the generous support of the American people through the United States Agency for International Development (USAID) Emerging Pandemic Threats PREDICT. The contents are the responsibility of the authors and do not necessarily reect the views of USAID or the United States Government.

Author contributions

All authors contributed to the acquisition, analysis and interpretation of the data; W.L., R.C., R.L.C., G.M.S., J.C.R., M.P., B.H.H. and P.M.S. initiated and conceived the study; W.L., Y.L., K.S.S., J.A.M, A.G. and C.F.D. performed non-invasive P. vivax testing and SGA analyses; J.A.M, M.A.R., P.A.C., A.G.S. and F.B.-R. performed genotyping; S.A.S. calculated ape P. vivax infection rates; G.H.L., L.J.P. and P.M.S. performed phylogenetic analyses; A.A., S.L., A.E., F.M., E.G., C.B. and S.A.-M. performed non-invasive P. vivax testing on apes and monkeys; S.L., S.A.-M., B.-I.I., J.-B.N.N., S.S., C.M.S., D.B.M., M.K.G., P.J.K., P.D.W., A.V.G., M.N.M., A.K.P., F.A.S., M.L.W., A.E.P., J.A.H., T.B.H., P.B., M.L., B.T., J.K., B.S.S., N.D.W., E.M.-N. and E.D. conducted or supervised eldwork; L.C., Z.W., A.F., C.J.S. and D.N. characterized human P. vivax infections; and W.L., R.C., R.L.C., G.M.S., J.C.R., M.P., B.H.H. and P.M.S. coordinated the contributions of all authors and wrote the paper.

Additional information

Accession codes: All newly derived ape and human Plasmodium sequences have been deposited in the GenBank Nucleotide Sequence Database under accession codes KF591752KF591851 and KF618374KF618618, respectively. DARC sequences have been deposited in the GenBank Nucleotide Sequence Database under accession codes KF618448KF618495.

Supplementary Information accompanies this paper at http://www.nature.com/nature

Web End =http://www.nature.com/nature communications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsand

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsand

Web End =reprintsand permissions/

How to cite this article: Liu, W. et al. African origin of the malaria parasite Plasmodium vivax. Nat. Commun. 5:3346 doi: 10.1038/ncomms4346 (2014).

10 NATURE COMMUNICATIONS | 5:3346 | DOI: 10.1038/ncomms4346 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.