New insights into the evolutionary history of Plasmodium falciparum from mitochondrial genome sequence analyses of Indian isolates

Intra-individual variability in ancient plasmodium DNA recovery highlights need for enhanced sampling

Malaria has been a leading cause of death in human populations for centuries and remains a major public health challenge in African countries, especially affecting children. Among the five Plasmodium species infecting humans, Plasmodium falciparum is the most lethal. Ancient DNA research has provided key insights into the origins, evolution, and virulence of pathogens that affect humans. However, extensive screening of ancient skeletal remains for Plasmodium DNA has shown that such genomic material is rare, with no studies so far addressing potential intra-individual variability. Consequently, the pool of ancient mitochondrial DNA (mtDNA) or genomic sequences for P. falciparum is extremely limited, with fewer than 20 ancient sequences available for genetic analysis, and no complete P. falciparum mtDNA from Classical antiquity published to date. To investigate intra-individual diversity and genetic origins of P. falciparum from the Roman period, we generated 39 sequencing libraries from multiple teeth and two from the femur of a Roman malaria-infected individual. The results revealed considerable variability in P. falciparum recovery across different dental samples within the individual, while the femur samples showed no preservation of Plasmodium DNA. The reconstructed 43-fold P. falciparum mtDNA genome supports the hypothesis of an Indian origin for European P. falciparum and suggests mtDNA continuity in Europe over the past 2000 years.

The evolutionary history of Plasmodium falciparum from mitochondrial and apicoplast genomes of China-Myanmar border isolates

BACKGROUND

The frequent communication between African and Southeast Asian (SEA) countries has led to the risk of imported malaria cases in the China-Myanmar border (CMB) region. Therefore, tracing the origins of new malaria infections is important in the maintenance of malaria-free zones in this border region. A new genotyping tool based on a robust mitochondrial (mt) /apicoplast (apico) barcode was developed to estimate genetic diversity and infer the evolutionary history of Plasmodium falciparum across the major distribution ranges. However, the mt/apico genomes of P. falciparum isolates from the CMB region to date are poorly characterized, even though this region is highly endemic to P. falciparum malaria.

METHODS

We have sequenced the whole mt/apico genome of 34 CMB field isolates and utilized a published data set of 147 mt/apico genome sequences to present global genetic diversity and to revisit the evolutionary history of the CMB P. falciparum.

RESULTS

Genetic differentiation based on mt/apico genome of P. falciparum revealed that the CMB (Lazan, Myanmar) isolates presented high genetic diversity with several characteristics of ancestral populations and shared many of the genetic features with West Thailand (Mae Sot; WTH) and to some extent West African (Banjul, Gambia; Navrongo, Ghana; WAF) isolates. The reconstructed haplotype network displayed that the CMB and WTH P. falciparum isolates have the highest representation (five) in the five ancestral (central) haplotypes (H1, H2, H4, H7, and H8), which are comparatively older than isolates from other SEA populations as well as the WAF populations. In addition, the highest estimate of the time to the Most Recent Common Ancestor (TMRCA) of 42,400 (95% CI 18,300-82100) years ago was presented by the CMB P. falciparum compared to the other regional populations. The statistically significant negative values of Fu's Fs with unimodal distribution in pairwise mismatch distribution curves indicate past demographic expansions in CMB P. falciparum with slow population expansion between approximately 12,500-20,000 ybp.

CONCLUSIONS

The results on the complete mt/apico genome sequence analysis of the CMB P. falciparum indicated high genetic diversity with ancient population expansion and TMRCA, and it seems probable that P. falciparum might have existed in CMB, WTH, and WAF for a long time before being introduced into other Southeast Asian countries or regions. To reduce the impact of sample size or geographic bias on the estimate of the evolutionary timeline, future studies need to expand the range of sample collection and ensure the representativeness of samples across geographic distributions. Additionally, by mapping global patterns of mt/apico genome polymorphism, we will gain valuable insights into the evolutionary history of P. falciparum and optimised strategies for controlling P. falciparum malaria at international borders.

Vivax malaria: a possible stumbling block for malaria elimination in India

Plasmodium vivax is geographically the most widely dispersed human malaria parasite species. It has shown resilience and a great deal of adaptability. Genomic studies suggest that P. vivax originated from Asia or Africa and moved to the rest of the world. Although P. vivax is evolutionarily an older species than Plasmodium falciparum, its biology, transmission, pathology, and control still require better elucidation. P. vivax poses problems for malaria elimination because of the ability of a single primary infection to produce multiple relapses over months and years. P. vivax malaria elimination program needs early diagnosis, and prompt and complete radical treatment, which is challenging, to simultaneously exterminate the circulating parasites and dormant hypnozoites lodged in the hepatocytes of the host liver. As prompt surveillance and effective treatments are rolled out, preventing primaquine toxicity in the patients having glucose-6-phosphate dehydrogenase (G6PD) deficiency should be a priority for the vivax elimination program. This review sheds light on the burden of P. vivax, changing epidemiological patterns, the hurdles in elimination efforts, and the essential tools needed not just in India but globally. These tools encompass innovative treatments for eliminating dormant parasites, coping with evolving drug resistance, and the development of potential vaccines against the parasite.

How can the complex epidemiology of malaria in India impact its elimination?

Malaria is a human health hazard in the tropical and subtropical zones of the globe and is poised to be eliminated by the year 2030. Despite a decrease in incidence in the past two decades, many endemic countries, including India, report cases regularly. The epidemiology of malaria in India is unique owing to several features of the Plasmodium parasites, Anopheles vectors, ecoepidemiological situations conducive to disease transmission, and susceptible humans living in rural and forested areas. Limitations in public health reach, and poor health-seeking behaviour of vulnerable populations living in hard-to-reach areas, add to the problem. We bring all of these factors together in a comprehensive framework and opine that, in spite of complexities, targeted elimination of malaria in India is achievable with planned programmatic approaches.

Genetic affinities of an eradicated European Plasmodium falciparum strain

Malaria was present in most of Europe until the second half of the 20th century, when it was eradicated through a combination of increased surveillance and mosquito control strategies, together with cross-border and political collaboration. Despite the severe burden of malaria on human populations, it remains contentious how the disease arrived and spread in Europe. Here, we report a partial Plasmodium falciparum nuclear genome derived from a set of antique medical slides stained with the blood of malaria-infected patients from Spain’s Ebro Delta, dating to the 1940s. Our analyses of the genome of this now eradicated European P. falciparum strain confirms stronger phylogeographical affinity to present-day strains in circulation in central south Asia, rather than to those in Africa. This points to a longitudinal, rather than a latitudinal, spread of malaria into Europe. In addition, this genome displays two derived alleles in the pfmrp1 gene that have been associated with drug resistance. Whilst this could represent standing variation in the ancestral P. falciparum population, these mutations may also have arisen due to the selective pressure of quinine treatment, which was an anti-malarial drug already in use by the time the sample we sequenced was mounted on a slide.

Using the Plasmodium mitochondrial genome for classifying mixed-species infections and inferring the geographical origin of P. falciparum parasites imported to the U.S

The ability to identify mixed-species infections and track the origin of Plasmodium parasites can further enhance the development of treatment and prevention recommendations as well as outbreak investigations. Here, we explore the utility of using the full Plasmodium mitochondrial genome to classify Plasmodium species, detect mixed infections, and infer the geographical origin of imported P. falciparum parasites to the United States (U.S.). Using the recently developed standardized, high-throughput Malaria Resistance Surveillance (MaRS) protocol, the full Plasmodium mitochondrial genomes of 265 malaria cases imported to the U.S. from 2014-2017 were sequenced and analyzed. P. falciparum infections were found in 94.7% (251/265) of samples. Five percent (14/265) of samples were identified as mixed- Plasmodium species or non-P. falciparum, including P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. P. falciparum mitochondrial haplotypes analysis revealed greater than eighteen percent of samples to have at least two P. falciparum mitochondrial genome haplotypes, indicating either heteroplasmy or multi-clonal infections. Maximum-likelihood phylogenies of 912 P. falciparum mitochondrial genomes with known country origin were used to infer the geographical origin of thirteen samples from persons with unknown travel histories as: Africa (country unspecified) (n = 10), Ghana (n = 1), Southeast Asia (n = 1), and the Philippines (n = 1). We demonstrate the utility and current limitations of using the Plasmodium mitochondrial genome to classify samples with mixed-infections and infer the geographical origin of imported P. falciparum malaria cases to the U.S. with unknown travel history.

K. S, Chandramohan B, Shanmuganatham V, Karanth KP. Reinvestigating the status of malaria parasite (Plasmodium sp.) in Indian non-human primates

Many human parasites and pathogens have closely related counterparts among non-human primates. For example, non-human primates harbour several species of malaria causing parasites of the genus Plasmodium. Studies suggest that for a better understanding of the origin and evolution of human malaria parasites it is important to know the diversity and evolutionary relationships of these parasites in non-human primates. Much work has been undertaken on malaria parasites in wild great Apes of Africa as well as wild monkeys of Southeast Asia however studies are lacking from South Asia, particularly India. India is one of the major malaria prone regions in the world and exhibits high primate diversity which in turn provides ideal setting for both zoonoses and anthropozoonoses. In this study we report the molecular data for malaria parasites from wild populations of Indian non-human primates. We surveyed 349 fecal samples from five different Indian non-human primates, while 94 blood and tissue samples from one of the Indian non-human primate species (Macaca radiata) and one blood sample from M. mulatta. Our results confirm the presence of P. fragile, P. inui and P. cynomolgi in Macaca radiata. Additionally, we report for the first time the presence of human malarial parasite, P. falciparum, in M. mulatta and M. radiata. Additionally, our results indicate that M. radiata does not exhibit population structure probably due to human mediated translocation of problem monkeys. Human mediated transport of macaques adds an additional level of complexity to tacking malaria in human. This issue has implications for both the spread of primate as well as human specific malarias.

Mitochondrial DNA from the eradicated European Plasmodium vivax and P. falciparum from 70-year-old slides from the Ebro Delta in Spain

Phylogenetic analysis of Plasmodium parasites has indicated that their modern-day distribution is a result of a series of human-mediated dispersals involving transport between Africa, Europe, America, and Asia. A major outstanding question is the phylogenetic affinity of the malaria causing parasites Plasmodium vivax and falciparum in historic southern Europe-where it was endemic until the mid-20th century, after which it was eradicated across the region. Resolving the identity of these parasites will be critical for answering several hypotheses on the malaria dispersal. Recently, a set of slides with blood stains of malaria-affected people from the Ebro Delta (Spain), dated between 1942 and 1944, have been found in a local medical collection. We extracted DNA from three slides, two of them stained with Giemsa (on which Plasmodium parasites could still be seen under the microscope) and another one consisting of dried blood spots. We generated the data using Illumina sequencing after using several strategies aimed at increasing the Plasmodium DNA yield: depletion of the human genomic (g)DNA content through hybridization with human gDNA baits, and capture-enrichment using gDNA derived from P. falciparum Plasmodium mitochondrial genome sequences were subsequently reconstructed from the resulting data. Phylogenetic analysis of the eradicated European P. vivax mtDNA genome indicates that the European isolate is closely related to the most common present-day American haplotype and likely entered the American continent post-Columbian contact. Furthermore, the European P. falciparum mtDNA indicates a link with current Indian strains that is in agreement with historical accounts.

Multi-method assessment of patients with febrile illness reveals over-diagnosis of malaria in rural Uganda

Background

Health clinics in rural Africa are typically resource-limited. As a result, many patients presenting with fever are treated with anti-malarial drugs based only on clinical presentation. This is a considerable issue in Uganda, where malaria is routinely over-diagnosed and over-treated, constituting a wastage of resources and an elevated risk of mortality in wrongly diagnosed patients. However, rapid diagnostic tests (RDTs) for malaria are increasingly being used in health facilities. Being fast, easy and inexpensive, RDTs offer the opportunity for feasible diagnostic capacity in resource-limited areas. This study evaluated the rate of malaria misdiagnosis and the accuracy of RDTs in rural Uganda, where presumptive diagnosis still predominates. Specifically, the diagnostic accuracy of “gold standard” methods, microscopy and PCR, were compared to the most feasible method, RDTs.

Methods

Patients presenting with fever at one of two health clinics in the Kabarole District of Uganda were enrolled in this study. Blood was collected by finger prick and used to administer RDTs, make blood smears for microscopy, and blot Whatman FTA cards for DNA extraction, polymerase chain reaction (PCR) amplification, and sequencing. The accuracy of RDTs and microscopy were assessed relative to PCR, considered the new standard of malaria diagnosis.

Results

A total of 78 patients were enrolled, and 31 were diagnosed with Plasmodium infection by at least one method. Comparing diagnostic pairs determined that RDTs and microscopy performed similarly, being 92.6 and 92.0 % sensitive and 95.5 and 94.4 % specific, respectively. Combining both methods resulted in a sensitivity of 96.0 % and specificity of 100 %. However, both RDTs and microscopy missed one case of non-falciparum malaria (Plasmodium malariae) that was identified and characterized by PCR and sequencing. In total, based on PCR, 62.0 % of patients would have been misdiagnosed with malaria if symptomatic diagnosis was used.

Conclusions

Results suggest that diagnosis of malaria based on symptoms alone appears to be highly inaccurate in this setting. Furthermore, RDTs were very effective at diagnosing malaria, performing as well or better than microscopy. However, only PCR and DNA sequencing detected non-P. falciparum species, which highlights an important limitation of this test and a treatment concern for non-falciparum malaria patients. Nevertheless, RDTs appear the only feasible method in rural or resource-limited areas, and therefore offer the best way forward in malaria management in endemic countries.

Sequence analysis of pfcrt and pfmdr1 genes and its association with chloroquine resistance in Southeast Indian Plasmodium falciparum isolates

BACKGROUND

Due to the widespread resistance of Plasmodium falciparum to chloroquine drug, artemisinin-based combination therapy (ACT) has been recommended as the first-line treatment. This study aims to evaluate the extent of chloroquine resistance in P. falciparum infection after the introduction of ACT. This study was carried out based on the mutation analysis in P. falciparum chloroquine resistant transporter (pfcrt) and P. falciparum multidrug resistance 1 (pfmdr1) genes. Identification of these molecular markers plays a significant role in monitoring and assessment of drug resistance as well as in designing an effective antimalarial drug policy in India.

METHODS

Sixty blood samples were collected from patients infected with P. falciparum from JIPMER, Puducherry and MKCG Medical College, Odisha. Polymerase chain reaction-restriction fragment length polymorphism was performed, targeting the point mutation of K76T in pfcrt and N86Y in pfmdr1 gene. The PCR products were sequenced, genotyped and further analysed for amino acid changes in these codons.

RESULTS

The frequency of pfcrt mutation at 76th position was dominant for mutant T allele with 56.7% and wild type K, 43.3%. Majority of pfmdr1 86 allele were wild type, with N (90%) and mutant, Y (10%). Additionally, we found three haplotypes for CQ resistance, SVMNT, CVIET and CVIKT in association with the pfcrt gene. However, a poorly studied SNP in pfmdr1 gene (Y184F) associated with CQ resistance showed high frequency (70%) in P. falciparum isolates.

CONCLUSIONS

The point mutation K76T of pfcrt is high in P. falciparum suggesting a sustained high CQ resistance even after five years of the introduction of ACTs for antimalarial therapy. The present study suggests a strong association of CQ resistance with pfcrt T76, but not with the pfmdr1 Y86 mutation. However, sequence analysis showed that Y184F mutation on pfmdr1 gene was found to be associated with high resistance. Also, a new pfcrt haplotype 'CVIKT' associated with CQ resistance was found to be present in Indian strains of P. falciparum. The data obtained from this study helps in continuous monitoring of drug resistance in malaria and also suggests the need for careful usage of CQ in Plasmodium vivax malarial treatment.

Mitochondrial population genomic analyses reveal population structure and demography of Indian Plasmodium falciparum

Inference on the genetic diversity of Plasmodium falciparum populations could help in better management of malaria. A very recent study with mitochondrial (mt) genomes in global P. falciparum had revealed interesting evolutionary genetic patterns of Indian isolates in comparison to global ones. However, no population genetic study using the whole mt genome sequences of P. falciparum isolates collected in the entire distribution range in India has yet been performed. We herewith have analyzed 85 whole mt genomes (48 already published and 37 entirely new) sampled from eight differentially endemic Indian locations to estimate genetic diversity and infer population structure and historical demography of Indian P. falciparum. We found 19 novel Indian-specific Single Nucleotide Polymorphisms (SNPs) and 22 novel haplotypes segregating in Indian P. falciparum. Accordingly, high haplotype and nucleotide diversities were detected in Indian P. falciparum in comparison to many other global isolates. Indian P. falciparum populations were found to be moderately sub-structured with four different genetic clusters. Interestingly, group of local populations aggregate to form each cluster; while samples from Jharkhand and Odisha formed a single cluster, P. falciparum isolates from Asom formed an independent one. Similarly, Surat, Bilaspur and Betul formed a single cluster and Goa and Mangalore formed another. Interestingly, P. falciparum isolates from the two later populations were significantly genetically differentiated from isolates collected in other six Indian locations. Signature of historical population expansion was evident in five population samples, and the onset of expansion event was found to be very similar to African P. falciparum. In agreement with the previous finding, the estimated Time to Most Recent Common Ancestor (TMRCA) and the effective population size were high in Indian P. falciparum. All these genetic features of Indian P. falciparum with high mt genome diversity are somehow similar to Africa, but quite different from other Asian population samples.

Genetic similarities between Cyclospora cayetanensis and cecum-infecting avian Eimeria spp. in apicoplast and mitochondrial genomes

Background

Cyclospora cayetanensis is an important cause for diarrhea in children in developing countries and foodborne outbreaks of cyclosporiasis in industrialized nations. To improve understanding of the basic biology of Cyclospora spp. and development of molecular diagnostic tools and therapeutics, we sequenced the complete apicoplast and mitochondrial genomes of C. cayetanensis.

Methods

The genome of one Chinese C. cayetanensis isolate was sequenced using Roche 454 and Illumina technologies. The assembled genomes of the apicoplast and mitochondrion were retrieved, annotated, and compared with reference genomes for other apicomplexans to infer genome organizations and phylogenetic relationships. Sequence variations in the mitochondrial genome were identified by comparison of two C. cayetanensis nucleotide sequences from this study and a recent publication.

Results

The apicoplast and mitochondrial genomes of C. cayetanensis are 34,155 and 6,229 bp in size and code for 65 and 5 genes, respectively. Comparative genomic analysis showed high similarities between C. cayetanensis and Eimeria tenella in both genomes; they have 85.6 % and 90.4 % nucleotide sequence similarities, respectively, and complete synteny in gene organization. Phylogenetic analysis of the genomic sequences confirmed the genetic similarities between cecum-infecting avian Eimeria spp. and C. cayetanensis. Like in other coccidia, both genomes of C. cayetanensis are transcribed bi-directionally. The apicoplast genome is circular, codes for the complete machinery for protein biosynthesis, and contains two inverted repeats that differ slightly in LSU rRNA gene sequences. In contrast, the mitochondrial genome has a linear concatemer or circular mapping topology. Eight single-nucleotide and one 7-bp multiple-nucleotide variants were detected between the mitochondrial genomes of C. cayetanensis from this and recent studies.

Conclusions

The apicoplast and mitochondrial genomes of C. cayetanensis are highly similar to those of cecum-infecting avian Eimeria spp. in both genome organization and sequences. The availability of sequence data beyond rRNA and heat shock protein genes could facilitate studies of C. cayetanensis biology and development of genotyping tools for investigations of cyclosporiasis outbreaks.

The distinctive features of Indian malaria parasites

Malaria and factors driving malaria are heterogeneous in India, unlike in other countries, and the epidemiology of malaria therefore is considered 'highly complex'. This complexity is primarily attributed to several unique features of the malaria parasites, mosquito vectors, malaria-susceptible populations, and ecoclimatic variables in India. Recent research on the genetic epidemiology of Indian malaria parasites has been successful in partly unraveling the mysteries underlying these complexities.

DNA sequence polymorphisms of the pfmdr1 gene and association of mutations with the pfcrt gene in Indian Plasmodium falciparum isolates

Mutations in the Plasmodium falciparum multidrug resistance (pfmdr1) gene are known to provide compensatory fitness benefits to the chloroquine (CQ)-resistant malaria parasites and are often associated with specific mutations in the P. falciparum CQ resistant transporter (pfcrt) gene. Prevalence of the specific mutations in these two genes across different malaria endemic regions was mostly studies. However, reports on mutations in the pfmdr1 gene and their genetic associations with mutations in the pfcrt gene in Indian P. falciparum field isolates are scarce. We have sequenced a 560 bp region of pfmdr1 coding sequence in 64 P. falciparum isolates collected from different malaria endemic populations in India. Twenty out of these 64 isolates were laboratory cultured with known in vitro CQ sensitiveness (10 sensitive and 10 resistant). Three low frequency mutations (two non-synonymous and one synonymous) in the pfmdr1 gene were segregating in Indian isolates in addition to the predominant Y₈₆ and Y₁₈₄ ones, with high haplotype and nucleotide diversity in the field isolates in comparison to the cultured ones. No statistically significant genetic association between the mutations in the pfmdr1 and pfcrt gene could be detected; almost all observed associations were intragenic in nature. The results on the genetic diversity of the pfmdr1 gene were discussed in term of evolutionary perspectives in Indian P. falciparum, with possible future potential of gaining further insights on this gene in view of evolving malaria parasites resistant to artemisinin partner drugs.