Comparative transcriptomics reveal differential gene expression among Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms.

The documentation of Plasmodium vivax malaria across Africa especially in regions where Duffy negatives are dominant suggests possibly alternative erythrocyte invasion mechanisms. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa is unclear. In this study, we examined the expression of 4,404 gene transcripts belong to 12 functional groups and 43 erythrocyte binding gene candidates in Ethiopian isolates and compared them with the Cambodian and Brazilian P. vivax transcriptomes. Overall, there were 10-26% differences in the gene expression profile amongst geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian isolates. Members of the reticulocyte binding protein PvRBP2a and PvRBP3 expressed six-fold higher than Duffy binding protein PvDBP1 and 60-fold higher than PvEBP/DBP2 in the Ethiopian isolates. Other genes including PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression patterns were observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, gametocyte genes including PvAP2-G, PvGAP (female gametocytes), and Pvs47 (male gametocytes) were highly expressed across geographical samples.

Comparative transcriptomics reveal differential gene expression in Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms.

Plasmodium vivax uses Duffy binding protein (PvDBP1) to bind to the Duffy Antigen-Chemokine Receptor (DARC) to invade human erythrocytes. Individuals who lack DARC expression (Duffy-negative) are thought to be resistance to P. vivax. In recent years, P. vivax malaria is becoming more prevalent in Africa with a portion of these cases detected in Duffy-negatives. Apart from DBP1, members of the reticulocyte binding protein (RBP) and tryptophan-rich antigen (TRAg) families may also play a role in erythrocyte invasion. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa and more specifically the expression level of several erythrocyte binding gene candidates as compared to DBP1 are largely unknown. This paper characterized the first P. vivax transcriptome in Africa and compared with those from the Southeast Asian and South American isolates. The expression of 4,404 gene transcripts belong to 12 functional groups including 43 specific erythrocyte binding gene candidates were examined. Overall, there were 10-26% differences in the gene expression profile amongst the geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian P. vivax. Erythrocyte binding genes including PvRBP2a and PvRBP3 expressed six-fold higher than PvDBP1and 60-fold higher than PvEBP/DBP2. Other genes including PvRBP1a, PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression was observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, the standard biomarker for detecting female gametocytes, PvAP2-G (PVP01_1440800), GAP (PVP01_1403000), and Pvs47 (PVP01_1208000) were highly expressed across geographical samples. These findings provide an important baseline for future comparisons of P. vivax transcriptomes from Duffy-negative infections and highlight potential biomarkers for improved gametocyte detection.

Prevalence of G6PD deficiency and distribution of its genetic variants among malaria-suspected patients visiting Metehara health centre, Eastern Ethiopia.

BACKGROUND: Glucose-6-phosphate dehydrogenase (G6PD) is cytosolic enzyme, which has a vital role for the integrity and functioning of red blood cells. Lower activity of this enzyme leads to the occurrence of acute haemolytic anaemia after exposure to oxidative stressors like primaquine. Primaquine is an important drug for the radical cure of Plasmodium vivax and blocking transmission of Plasmodium falciparum, and thereby enhancing malaria elimination. However, there is a need to identify G6PD deficient individuals and administer the drug with caution due to its haemolytic side effects. The main objective of this study is to determine the prevalence of G6PD deficiency among malaria-suspected individuals. METHODS: A facility-based cross-sectional study was conducted from September 2020 to September 2021 in Metehara Health Centre, Eastern Ethiopia. A structured questionnaire was used to collect the socio-demographic and clinical information of the study participants. Capillary and venous blood samples were collected based on standard procedures for onsite screening, dried blood spot preparation, and malaria microscopy. The G6PD enzyme activity was measured by careSTART™ G6PD biosensor analyzer. Data was entered and analysed by SPSS. RESULTS: A total of 498 study participants were included in the study, of which 62% (309) were males. The overall prevalence of G6PD deficiency based on the biosensor screening was 3.6% (18/498), of which 2.9% and 4.8% were males and females, respectively. Eleven of the G6PD deficient samples had mutations confirmed by G6PD gene sequencing analysis. Mutations were detected in G267 + 119C/T, A376T, and ChrX:154535443. A significant association was found in sex and history of previous malaria infection with G6PD deficiency. CONCLUSIONS: The study showed that the G6PD deficient phenotype exists in Metehara even if the prevalence is not very high. G267 + 119C/T mutation is the predominant G6PD variant in this area. Therefore, malaria patient treatment using primaquine should be monitored closely for any adverse effects.

Low Density Plasmodium Infections and G6PD Deficiency Among Malaria Suspected Febrile Individuals in Ethiopia.

The identification and management of low parasitemia infections have become increasingly challenging for malaria control and elimination. Submicroscopic Plasmodium infections and G6PD deficiency among febrile patients require more sensitive diagnostic methods to improve detection and careful treatment regime of these infections. In Ethiopia, information on the low density submicroscopic malarial infections and frequency of G6PD deficiency (G6PDd) is scarce. In this study, 297 malaria suspected febrile patient samples were collected from health facilities of Bonga town in southwestern Ethiopia. The positivity rates of Plasmodium infection were determined by microscopy and quantitative PCR. G6PD activity level was determined by careSTART™ G6PD biosensor and the frequency of three common variants: G6PD*A (A376G), G6PD*A- (G202A) and Mediterranean (C563T) were investigated. G6PD gene sequencing was performed to detect mutations in exons 2-11 for both G6PD normal and deficient samples based on the phenotypic assay. More than twice Plasmodium infected samples was detected by qPCR (52/297; 17.4%) than microscopy (21/297; 7.0%). About 31 (10%) of the infections were submicroscopic. Bednet usage and age had a significant association with Plasmodium infection. Of the 271 participants who were tested for G6PD phenotype, 19 (7.0%) had low G6PD level. No mutations were observed in A376G, G202A, and C563T in the G6PDd samples, but three novel non-synonymous mutations in exon 2 including a C to T transition at position ChrX:6504 (Arg to Thr), G to T at ChrX:6369 (Ser to IIe), and G to C at ChrX:6664 (Gln to His) were detected. A high number of submicroscopic Plasmodium infections observed in this study pose a challenge for accurate and timely diagnosis, which could hinder malaria control efforts. G6PD deficiency in malaria patients pose danger when treating patients with primaquine. The three novel mutations detected in exon 2 of the G6PD gene merit further investigation on the hemolytic risk when exposed to oxidative antimalarials, their prevalence, and clinical significance.

Contrasting epidemiology and genetic variation of Plasmodium vivax infecting Duffy-negative individuals across Africa.

OBJECTIVES: Plasmodium vivax malaria was thought to be rare in Africans who lack the Duffy blood group antigen expression. However, recent studies indicate that P. vivax can infect Duffy-negative individuals and has spread into areas of high Duffy negativity across Africa. Our study compared epidemiological and genetic features of P. vivax between African regions. METHODS: A standardized approach was used to identify and quantify P. vivax from Botswana, Ethiopia, and Sudan, where Duffy-positive and Duffy-negative individuals coexist. The study involved sequencing the Duffy binding protein (DBP) gene and inferring genetic relationships among P. vivax populations across Africa. RESULTS: Among 1215 febrile patients, the proportions of Duffy negativity ranged from 20-36% in East Africa to 84% in southern Africa. Average P. vivax prevalence among Duffy-negative populations ranged from 9.2% in Sudan to 86% in Botswana. Parasite density in Duffy-negative infections was significantly lower than in Duffy-positive infections. P. vivax in Duffy-negative populations were not monophyletic, with P. vivax in Duffy-negative and Duffy-positive populations sharing similar DBP haplotypes and occurring in multiple, well-supported clades. CONCLUSIONS: Duffy-negative Africans are not resistant to P. vivax, and the public health significance of this should not be neglected. Our study highlights the need for a standardized approach and more resources/training directed towards the diagnosis of vivax malaria in Africa.

Plasmodium vivax From Duffy-Negative and Duffy-Positive Individuals Share Similar Gene Pools in East Africa.

Plasmodium vivax malaria was thought to be rare in Africa, but an increasing number of P. vivax cases reported across Africa and in Duffy-negative individuals challenges this dogma. The genetic characteristics of P. vivax in Duffy-negative infections, the transmission of P. vivax in East Africa, and the impact of environments on transmission remain largely unknown. This study examined genetic and transmission features of P. vivax from 107 Duffy-negative and 305 Duffy-positive individuals in Ethiopia and Sudan. No clear genetic differentiation was found in P. vivax between the 2 Duffy groups, indicating between-host transmission. P. vivax from Ethiopia and Sudan showed similar genetic clusters, except samples from Khartoum, possibly due to distance and road density that inhibited parasite gene flow. This study is the first to show that P. vivax can transmit to and from Duffy-negative individuals and provides critical insights into the spread of P. vivax in sub-Saharan Africa.

Gene Polymorphisms Among Plasmodium vivax Geographical Isolates and the Potential as New Biomarkers for Gametocyte Detection.

The unique biological features of Plasmodium vivax not only make it difficult to control but also to eliminate. For the transmission of the malaria parasite from infected human to the vector, gametocytes play a major role. The transmission potential of a malarial infection is inferred based on microscopic detection of gametocytes and molecular screening of genes in the female gametocytes. Microscopy-based detection methods could grossly underestimate the reservoirs of infection as gametocytes may occur as submicroscopic or as micro- or macro-gametocytes. The identification of genes that are highly expressed and polymorphic in male and female gametocytes is critical for monitoring changes not only in their relative proportions but also the composition of gametocyte clones contributing to transmission over time. Recent transcriptomic study revealed two distinct clusters of highly correlated genes expressed in the P. vivax gametocytes, indicating that the male and female terminal gametocytogeneses are independently regulated. However, the detective power of these genes is unclear. In this study, we compared genetic variations of 15 and 11 genes expressed, respectively, in the female and male gametocytes among P. vivax isolates from Southeast Asia, Africa, and South America. Further, we constructed phylogenetic trees to determine the resolution power and clustering patterns of gametocyte clones. As expected, Pvs25 (PVP01_0616100) and Pvs16 (PVP01_0305600) expressed in the female gametocytes were highly conserved in all geographical isolates. In contrast, genes including 6-cysteine protein Pvs230 (PVP01_0415800) and upregulated in late gametocytes ULG8 (PVP01_1452800) expressed in the female gametocytes, as well as two CPW-WPC family proteins (PVP01_1215900 and PVP01_1320100) expressed in the male gametocytes indicated considerably high nucleotide and haplotype diversity among isolates. Parasite samples expressed in male and female gametocyte genes were observed in separate phylogenetic clusters and likely represented distinct gametocyte clones. Compared to Pvs25, Pvs230 (PVP01_0415800) and a CPW-WPC family protein (PVP01_0904300) showed higher expression in a subset of Ethiopian P. vivax samples. Thus, Pvs230, ULG8, and CPW-WPC family proteins including PVP01_0904300, PVP01_1215900, and PVP01_1320100 could potentially be used as novel biomarkers for detecting both sexes of P. vivax gametocytes in low-density infections and estimating transmission reservoirs.

Alternative Invasion Mechanisms and Host Immune Response to Plasmodium vivax Malaria: Trends and Future Directions.

Plasmodium vivax malaria is a neglected tropical disease, despite being more geographically widespread than any other form of malaria. The documentation of P. vivax infections in different parts of Africa where Duffy-negative individuals are predominant suggested that there are alternative pathways for P. vivax to invade human erythrocytes. Duffy-negative individuals may be just as fit as Duffy-positive individuals and are no longer resistant to P.vivax malaria. In this review, we describe the complexity of P. vivax malaria, characterize pathogenesis and candidate invasion genes of P. vivax, and host immune responses to P. vivax infections. We provide a comprehensive review on parasite ligands in several Plasmodium species that further justify candidate genes in P. vivax. We also summarize previous genomic and transcriptomic studies related to the identification of ligand and receptor proteins in P. vivax erythrocyte invasion. Finally, we identify topics that remain unclear and propose future studies that will greatly contribute to our knowledge of P. vivax.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes.

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale.

Frequent expansion of Plasmodium vivax Duffy Binding Protein in Ethiopia and its epidemiological significance.

Plasmodium vivax invasion of human erythrocytes depends on the Duffy Binding Protein (PvDBP) which interacts with the Duffy antigen. PvDBP copy number has been recently shown to vary between P. vivax isolates in Sub-Saharan Africa. However, the extent of PvDBP copy number variation, the type of PvDBP multiplications, as well as its significance across broad samples are still unclear. We determined the prevalence and type of PvDBP duplications, as well as PvDBP copy number variation among 178 Ethiopian P. vivax isolates using a PCR-based diagnostic method, a novel quantitative real-time PCR assay and whole genome sequencing. For the 145 symptomatic samples, PvDBP duplications were detected in 95 isolates, of which 81 had the Cambodian and 14 Malagasy-type PvDBP duplications. PvDBP varied from 1 to >4 copies. Isolates with multiple PvDBP copies were found to be higher in symptomatic than asymptomatic infections. For the 33 asymptomatic samples, PvDBP was detected with two copies in two of the isolates, and both were the Cambodian-type PvDBP duplication. PvDBP copy number in Duffy-negative heterozygotes was not significantly different from that in Duffy-positives, providing no support for the hypothesis that increased copy number is a specific association with Duffy-negativity, although the number of Duffy-negatives was small and further sampling is required to test this association thoroughly.

Mosaic representations of odors in the input and output layers of the mouse olfactory bulb.

The elementary stimulus features encoded by the olfactory system remain poorly understood. We examined the relationship between 1,666 physical-chemical descriptors of odors and the activity of olfactory bulb inputs and outputs in awake mice. Glomerular and mitral and tufted cell responses were sparse and locally heterogeneous, with only a weak dependence of their positions on physical-chemical properties. Odor features represented by ensembles of mitral and tufted cells were overlapping but distinct from those represented in glomeruli, which is consistent with an extensive interplay between feedforward and feedback inputs to the bulb. This reformatting was well described as a rotation in odor space. The physical-chemical descriptors accounted for a small fraction in response variance, and the similarity of odors in the physical-chemical space was a poor predictor of similarity in neuronal representations. Our results suggest that commonly used physical-chemical properties are not systematically represented in bulbar activity and encourage further searches for better descriptors of odor space.

Deconstructing Odorant Identity via Primacy in Dual Networks.

In the olfactory system, odor percepts retain their identity despite substantial variations in concentration, timing, and background. We study a novel strategy for encoding intensity-invariant stimulus identity that is based on representing relative rather than absolute values of stimulus features. For example, in what is known as the primacy coding model, odorant identities are represented by the conditions that some odorant receptors are activated more strongly than others. Because, in this scheme, odorant identity depends only on the relative amplitudes of olfactory receptor responses, identity is invariant to changes in both intensity and monotonic nonlinear transformations of its neuronal responses. Here we show that sparse vectors representing odorant mixtures can be recovered in a compressed sensing framework via elastic net loss minimization. In the primacy model, this minimization is performed under the constraint that some receptors respond to a given odorant more strongly than others. Using duality transformation, we show that this constrained optimization problem can be solved by a neural network whose Lyapunov function represents the dual Lagrangian and whose neural responses represent the Lagrange coefficients of primacy and other constraints. The connectivity in such a dual network resembles known features of connectivity in olfactory circuits. We thus propose that networks in the piriform cortex implement dual computations to compute odorant identity with the sparse activities of individual neurons representing Lagrange coefficients. More generally, we propose that sparse neuronal firing rates may represent Lagrange multipliers, which we call the dual brain hypothesis. We show such a formulation is well suited to solve problems with multiple interacting relative constraints.