Plasmodium vivax lineages: geographical distribution, tandem repeat polymorphism, and phylogenetic relationship

Lizard host abundances and climatic factors explain phylogenetic diversity and prevalence of blood parasites on an oceanic island

Host abundance might favour the maintenance of a high phylogenetic diversity of some parasites via rapid transmission rates. Blood parasites of insular lizards represent a good model to test this hypothesis because these parasites can be particularly prevalent in islands and host lizards highly abundant. We applied deep amplicon sequencing and analysed environmental predictors of blood parasite prevalence and phylogenetic diversity in the endemic lizard Gallotia galloti across 24 localities on Tenerife, an island in the Canary archipelago that has experienced increasing warming and drought in recent years. Parasite prevalence assessed by microscopy was over 94%, and a higher proportion of infected lizards was found in warmer and drier locations. A total of 33 different 18s rRNA parasite haplotypes were identified, and the phylogenetic analyses indicated that they belong to two genera of Adeleorina (Apicomplexa: Coccidia), with Karyolysus as the dominant genus. The most important predictor of between-locality variation in parasite phylogenetic diversity was the abundance of lizard hosts. We conclude that a combination of climatic and host demographic factors associated with an insular syndrome may be favouring a rapid transmission of blood parasites among lizards on Tenerife, which may favour the maintenance of a high phylogenetic diversity of parasites.

Compendium-Wide Analysis of Pseudomonas aeruginosa Core and Accessory Genes Reveals Transcriptional Patterns across Strains PAO1 and PA14

Pseudomonas aeruginosa is an opportunistic pathogen that causes difficult-to-treat infections. Two well-studied divergent P. aeruginosa strain types, PAO1 and PA14, have significant genomic heterogeneity, including diverse accessory genes present in only some strains. Genome content comparisons find core genes that are conserved across both PAO1 and PA14 strains and accessory genes that are present in only a subset of PAO1 and PA14 strains. Here, we use recently assembled transcriptome compendia of publicly available P. aeruginosa RNA sequencing (RNA-seq) samples to create two smaller compendia consisting of only strain PAO1 or strain PA14 samples with each aligned to their cognate reference genome. We confirmed strain annotations and identified other samples for inclusion by assessing each sample's median expression of PAO1-only or PA14-only accessory genes. We then compared the patterns of core gene expression in each strain. To do so, we developed a method by which we analyzed genes in terms of which genes showed similar expression patterns across strain types. We found that some core genes had consistent correlated expression patterns across both compendia, while others were less stable in an interstrain comparison. For each accessory gene, we also determined core genes with correlated expression patterns. We found that stable core genes had fewer coexpressed neighbors that were accessory genes. Overall, this approach for analyzing expression patterns across strain types can be extended to other groups of genes, like phage genes, or applied for analyzing patterns beyond groups of strains, such as samples with different traits, to reveal a deeper understanding of regulation. IMPORTANCE Pseudomonas aeruginosa is a ubiquitous pathogen. There is much diversity among P. aeruginosa strains, including two divergent but well-studied strains, PAO1 and PA14. Understanding how these different strain-level traits manifest is important for identifying targets that regulate different traits of interest. With the availability of thousands of PAO1 and PA14 samples, we created two strain-specific RNA-seq compendia where each one contains hundreds of samples from PAO1 or PA14 strains and used them to compare the expression patterns of core genes that are conserved in both strain types and to determine which core genes have expression patterns that are similar to those of accessory genes that are unique to one strain or the other using an approach that we developed. We found a subset of core genes with different transcriptional patterns across PAO1 and PA14 strains and identified those core genes with expression patterns similar to those of strain-specific accessory genes.

Multi-Omics Advancements towards Plasmodium vivax Malaria Diagnosis

Plasmodium vivax malaria is one of the most lethal infectious diseases, with 7 million infections annually. One of the roadblocks to global malaria elimination is the lack of highly sensitive, specific, and accurate diagnostic tools. The absence of diagnostic tools in particular has led to poor differentiation among parasite species, poor prognosis, and delayed treatment. The improvement necessary in diagnostic tools can be broadly grouped into two categories: technologies-driven and omics-driven progress over time. This article discusses the recent advancement in omics-based malaria for identifying the next generation biomarkers for a highly sensitive and specific assay with a rapid and antecedent prognosis of the disease. We summarize the state-of-the-art diagnostic technologies, the key challenges, opportunities, and emerging prospects of multi-omics-based sensors.

Evidence of a Recent Bottleneck in Plasmodium falciparum Populations on the Honduran-Nicaraguan Border

The countries of Central America and the island of Hispaniola have set the goal of eliminating malaria in less than a decade. Although efforts to reduce the malaria burden in the region have been successful, there has been an alarming increase in cases in the Nicaraguan Moskitia since 2014. The continuous decrease in cases between 2000 and 2014, followed by a rapid expansion from 2015 to the present, has generated a potential bottleneck effect in the populations of Plasmodium spp. Consequently, this study aimed to evaluate the genetic diversity of P. falciparum and the decrease in allelic richness in this population. The polymorphic regions of the pfmsp-1 and pfmsp-2 genes of patients with falciparum malaria from Honduras and Nicaragua were analyzed using nested PCR and sequencing. Most of the samples were classified into the K1 allelic subfamily of the pfmsp-1 gene and into the 3D7 subfamily of the pfmsp-2 gene. Despite the low genetic diversity found, more than half of the samples presented a polyclonal K1/RO33 haplotype. No sequence polymorphisms were found within each allelic subfamily. This study describes a notable decrease in the genetic diversity of P. falciparum in the Moskitia region after a bottleneck phenomenon. These results will be useful for future epidemiological investigations and the monitoring of malaria transmission in Central America.

Ookinete-Specific Genes and 18S SSU rRNA Evidenced in Plasmodium vivax Selection and Adaptation by Sympatric Vectors

In the southern Pacific coast of Chiapas, Mexico (SM), the two most abundant vector species, Nyssorhynchus albimanus and Anopheles pseudopunctipennis, were susceptible to different Plasmodium vivax Pvs25/28 haplotypes. To broaden our understanding of the existing P. vivax in the area, genes encoding proteins relevant for ookinete development and the 18S rRNA were studied. P. vivax infectivity (percentage of infected mosquitoes and oocyst numbers) was evaluated by simultaneously feeding infected blood samples from patients to Ny. albimanus and An. pseudopunctipennis female mosquitoes. Three infectivity patterns were identified: one group of parasites were more infective to An. pseudopunctipennis than to Ny. albimanus, another group was more infective to Ny. albimanus, while a third group infected both vectors similarly. In 29 parasite isolates, the molecular variations of ookinete-specific genes and the 18S rRNA-type S were analyzed. Using concatenated sequences, phylogenetic trees, and Structure analysis, parasite clustering within SM isolates and between these and those from other geographical origins were investigated. A ML phylogenetic tree resolved two parasite lineages: PvSM-A and PvSM-B. They were associated to a different 18S rRNA variant. PvSM-A parasites had 18S rRNA variant rV2 and correspond to parasites causing high oocyst infection in Ny. albimanus. A new ML tree and Structure analysis, both comprising global sequences, showed PvSM-A clustered with Latin American parasites. Meanwhile, all isolates of PvSM-B had 18S rRNA variant rV1 and remained as unique genetic cluster comprising two subgroups: PvSM-Ba, producing high infection in An. pseudopunctipennis, and PvSM-Bb, causing similar oocyst infection in both vector species. PvSM-A parasites were genetically similar to parasites from South America. Meanwhile, PvSM-B were exclusive to southern Mexico and share ancestry with Asian parasites. The results suggest that these lineages evolved separately, likely by geographic and vector restriction.

Genetic profiling of the Plasmodium falciparum parasite population in uncomplicated malaria from India

Background

The genetic complexity and the existence of several polymorphisms in parasites are the major hindrances for the malaria control programmes of the country. The genetic profiling in the parasite populations in India will provide useful baseline data for future studies elucidating the parasite structure and distribution of drug resistance genotypes in different regions.

Methods

The blood samples of symptomatic patients were collected and analysed for drug resistance genes (Pfcrt, Pfmdr-1, dhfr, dhps and k13) and gametocyte genes (Pfs25, Pfg377); in vitro drug sensitivity assay by schizont maturation inhibition (SMI) was also performed in adapted field isolates.

Results

Of the 122 field isolates analysed; 65.5% showed Pfcrt K76T mutant alleles, 61.4% Pfmdr-1 N86Y mutants, 59.5% dhfr mutants, 59.8% dhps mutants was observed, but no polymorphism was seen for k13. The sequence analysis of Pfg377 gene revealed five types of populations in the field isolates. The inhibitory concentrations (IC₅₀) for anti-malarial drugs viz chloroquine (CQ), artesunate (AS), were in the range of 10.11–113.2 nM and 2.26–4.08 nM, respectively, in the field isolates evaluate by in vitro assay. The IC₅₀ values for CQ have shown a remarkable reduction on comparison with the previous available data, whereas a slight increase in the IC₅₀ values for AS was observed in the study.

Conclusions

The increase in mutation rate in drug resistance allelic loci with inhibitory concentration of CQ and AS drugs was observed in the field isolates and high diversity in Pfg377 gametocyte gene indicate towards parasite multifactorial behaviour. The knowledge of the prevalent drug resistance genes is important for intervention measures to be successful and efforts should also be made to prevent transmission of P. falciparum.

Substantial population structure of Plasmodium vivax in Thailand facilitates identification of the sources of residual transmission

Background

Plasmodium vivax transmission in Thailand has been substantially reduced over the past 10 years, yet it remains highly endemic along international borders. Understanding the genetic relationship of residual parasite populations can help track the origins of the parasites that are reintroduced into malaria-free regions within the country.

Methodology/Results

A total of 127 P. vivax isolates were genotyped from two western provinces (Tak and Kanchanaburi) and one eastern province (Ubon Ratchathani) of Thailand using 10 microsatellite markers. Genetic diversity was high, but recent clonal expansion was detected in all three provinces. Substantial population structure and genetic differentiation of parasites among provinces suggest limited gene flow among these sites. There was no haplotype sharing among the three sites, and a reduced panel of four microsatellite markers was sufficient to assign the parasites to their provincial origins.

Conclusion/Significance

Significant parasite genetic differentiation between provinces shows successful interruption of parasite spread within Thailand, but high diversity along international borders implies a substantial parasite population size in these regions. The provincial origin of P. vivax cases can be reliably determined by genotyping four microsatellite markers, which should be useful for monitoring parasite reintroduction after malaria elimination.

This study presents an updated view of the P. vivax populations along the Thai-Myanmar and the Thai-Cambodian borders. Genotyping of parasite samples collected after intensified malaria control demonstrated that despite the decline in overall transmission intensity, the genetic diversity of the P. vivax parasites remained high. Parasite populations from three border provinces showed clear genetic separation. This indicates successful interruption of parasite gene flow within Thailand, but suggests frequent parasite migration across international borders. From the analysis of 10 microsatellite markers, we further refined a set of four that are sufficient for distinguishing the provincial origins of these parasites, which should allow tracking of parasite introduction among these provinces.

Source identification of autochthonous-introduced Plasmodium vivax Malaria, Spain

In 2014, an autochthonous case of introduced malaria caused by Plasmodium vivax was identified in Spain. The strain that infected this patient was identical to that of a prior imported case from Pakistan. This is the first case where the source of infection could be identified since elimination in Spain.

Genotyping of Plasmodium vivax by minisatellite marker and its application in differentiating relapse and new infection

BACKGROUND

Plasmodium vivax malaria is a major public health problem in India. Control of vivax malaria is challenging due to various factors including relapse which increase the burden significantly. There is no well studied marker to differentiate relapse from reinfection. This creates hindrance in search for anti-relapse medicines. The genomic study of minisatellite can help in characterization of relapse and new infection of vivax malaria.

METHODS

Eighty-eight samples of P. vivax were collected from malaria clinic. All the 14 chromosomes of P. vivax were scanned for minisatellite marker by Tandem Repeat Finder software Version 4.07b. Minisatellite marker CH1T1M13779 from chromosome one was applied for genotyping in 88 samples of P. vivax including 2 recurrence cases.

RESULTS

Whole genome of P. vivax was scanned and found to have one hundred minisatellite markers. CH1T1M13779 minisatellite marker from chromosome-1 was used for amplification in 88 samples of P. vivax. Of 66 amplified samples, 14 alleles were found with varied allele frequency. The base size of 280 (13.63 %) 320 bp (13.63 %) and 300 bp (16.66 %) showed the predominant allele in the P. vivax population. Genotyping of two paired samples (day 0 and day relapse) could demonstrate the presence of relapse and reinfection.

CONCLUSION

The CH1T1M13779 can be potential minisatellite marker which can be used to differentiate between relapse and new infection of P. vivax strain.

Low genetic diversity and functional constraint in loci encoding Plasmodium vivax P12 and P38 proteins in the Colombian population

Background

Plasmodium vivax is one of the five species causing malaria in human beings, affecting around 391 million people annually. The development of an anti-malarial vaccine has been proposed as an alternative for controlling this disease. However, its development has been hampered by allele-specific responses produced by the high genetic diversity shown by some parasite antigens. Evaluating these antigens’ genetic diversity is thus essential when designing a completely effective vaccine.

Methods

The gene sequences of Plasmodium vivax p12 (pv12) and p38 (pv38), obtained from field isolates in Colombia, were used for evaluating haplotype polymorphism and distribution by population genetics analysis. The evolutionary forces generating the variation pattern so observed were also determined.

Results

Both pv12 and pv38 were shown to have low genetic diversity. The neutral model for pv12 could not be discarded, whilst polymorphism in pv38 was maintained by balanced selection restricted to the gene’s 5′ region. Both encoded proteins seemed to have functional/structural constraints due to the presence of s48/45 domains, which were seen to be highly conserved.

Conclusions

Due to the role that malaria parasite P12 and P38 proteins seem to play during invasion in Plasmodium species, added to the Pv12 and Pv38 antigenic characteristics and the low genetic diversity observed, these proteins might be good candidates to be evaluated in the design of a multistage/multi-antigen vaccine.

Neutral polymorphisms in putative housekeeping genes and tandem repeats unravels the population genetics and evolutionary history of Plasmodium vivax in India

The evolutionary history and age of Plasmodium vivax has been inferred as both recent and ancient by several studies, mainly using mitochondrial genome diversity. Here we address the age of P. vivax on the Indian subcontinent using selectively neutral housekeeping genes and tandem repeat loci. Analysis of ten housekeeping genes revealed a substantial number of SNPs (n = 75) from 100 P. vivax isolates collected from five geographical regions of India. Neutrality tests showed a majority of the housekeeping genes were selectively neutral, confirming the suitability of housekeeping genes for inferring the evolutionary history of P. vivax. In addition, a genetic differentiation test using housekeeping gene polymorphism data showed a lack of geographical structuring between the five regions of India. The coalescence analysis of the time to the most recent common ancestor estimate yielded an ancient TMRCA (232,228 to 303,030 years) and long-term population history (79,235 to 104,008) of extant P. vivax on the Indian subcontinent. Analysis of 18 tandem repeat loci polymorphisms showed substantial allelic diversity and heterozygosity per locus, and analysis of potential bottlenecks revealed the signature of a stable P. vivax population, further corroborating our ancient age estimates. For the first time we report a comparable evolutionary history of P. vivax inferred by nuclear genetic markers (putative housekeeping genes) to that inferred from mitochondrial genome diversity.

Molecular epidemiology of Plasmodium vivax in Latin America: polymorphism and evolutionary relationships of the circumsporozoite gene

Background

The origins and dispersal of Plasmodium vivax to its current worldwide distribution remains controversial. Although progress on P. vivax genetics and genomics has been achieved worldwide, information concerning New World parasites remains fragmented and largely incomplete. More information on the genetic diversity in Latin America (LA) is needed to better explain current patterns of parasite dispersion and evolution.

Methods

Plasmodium vivax circumsporozoite protein gene polymorphism was investigated using polymerase chain reaction amplification and restriction fragment length polymorphism (PCR-RFLP), and Sanger sequencing in isolates from the Pacific Ocean coast of Mexico, Nicaragua, and Peru. In conjunction with worldwide sequences retrieved from the Genbank, mismatch distribution analysis of central repeat region (CRR), frequency estimation of unique repeat types and phylogenetic analysis of the 3′ terminal region, were performed to obtain an integrative view of the genetic relationships between regional and worldwide isolates.

Results

Four RFLP subtypes, vk210a, b, c and d were identified in Southern Mexico and three subtypes vk210a, e and f in Nicaragua. The nucleotide sequences showed that Mexican vk210a and all Nicaraguan isolates were similar to other American parasites. In contrast, vk210b, c and d were less frequent, had a domain ANKKAEDA in their carboxyl end and clustered with Asian isolates. All vk247 isolates from Mexico and Peru had identical RFLP pattern. Their nucleotide sequences showed two copies of GGQAAGGNAANKKAGDAGA at the carboxyl end. Differences in mismatch distribution parameters of the CRR separate vk247 from most vk210 isolates. While vk247 isolates display a homogeneous pattern with no geographical clustering, vk210 isolates display a heterogeneous geographically clustered pattern which clearly separates LA from non-American isolates, except vk210b, c and d from Southern Mexico.

Conclusions

The presence of vk210a in Mexico and vk210e, f and g in Nicaragua are consistent with other previously reported LA isolates and reflect their circulation throughout the continent. The vk210b, c and d are novel genotypes in LA. Their genetic relationships and low variability within these vk210 and/or within the vk247 parasites in Southern Mexico suggest its recent introduction and/or recent expansion to this region. The global analysis of P. vivax csp suggests this parasite introduction to the region and likely LA by different independent events.

Quantifying effect of geographic location on epidemiology of Plasmodium vivax malaria

Recent autochthonous transmission of Plasmodium vivax malaria in previously malaria-free temperate regions has generated renewed interest in the epidemiology of this disease. Accurate estimates of the incubation period and time to relapse are required for effective malaria surveillance; however, this information is currently lacking. By using historical data from experimental human infections with diverse P. vivax strains, survival analysis models were used to obtain quantitative estimates of the incubation period and time to first relapse for P. vivax malaria in broad geographic regions. Results show that Eurasian strains from temperate regions have longer incubation periods, and Western Hemisphere strains from tropical and temperate regions have longer times to relapse compared with Eastern Hemisphere strains. The diversity in these estimates of key epidemiologic parameters for P. vivax supports the need for elucidating local epidemiology to inform clinical follow-up and to build an evidence base toward global elimination of malaria.