Population genetic and biophysical evidences reveal that purifying selection shapes the genetic landscape of Plasmodium falciparum RH ligands in Chhattisgarh and West Bengal, India.

BACKGROUND: Reticulocyte binding protein-like homologs (RHs) are currently being evaluated as anti-erythrocytic stage vaccine targets against Plasmodium falciparum malaria. Present study explores the possible evolutionary drivers shaping the genetic organization of Pfrhs in Indian parasite population. It simultaneously evaluates a putative gain-of-function variant of PfRH5, a keystone member of PfRH family. METHODS: Receptor binding regions of Pfrh1, Pfrh2a/b, Pfrh4 and whole Pfrh5 were amplified using blood samples of P. falciparum malaria patients from Chhattisgarh and West Bengal and sequenced. Assembled sequences were analysed using MEGA7 and DnaSPv6. Binding affinities of recombinant PfRH5 proteins with basigin (BSG) were compared using in silico (CHARMM and AUTODOCK) and in vitro (Circular dichroism, fluorescence spectroscopy and isothermal titration calorimetry) methods. RESULTS: Pfrh1 (0.5), Pfrh2a/b (0.875), Pfrh4 (0.667) and Pfrh5 (0.778) sequence changes corresponded to low frequency (< 0.05) variants which resulted in an overall negative Tajima's D. Since mismatch distribution of none of the Pfrh loci corroborated with the model of demographic expansion, a possible role of natural selection formulating Pfrh sequence diversity was investigated. Among the 5 members, Pfrh5 displayed very high dN/dS (5.7) ratio. Nevertheless, the model of selective sweep due to presence of any advantageous substitutions could not be invoked as polymorphic nonsynonymous sites (17/18) for Pfrh5 exceeded significantly over the divergent (62/86) ones (p = 0.0436). The majority of extant PfRH5 sequences (52/83) differed from the reference Pf3D7 allele by a single amino acid mismatch (C203Y). This non-conservative alteration was predicted to lower the total interaction energy of that PfRH5variant with BSG, compared to PfRH53D7. Biophysical evidences validated the proposition that PfRH5variant formed a more stable complex with BSG. Thermodynamic association constant for interaction of BSG with PfRH5variant was also found to be higher (Kavariant = 3.63E6 ± 2.02E6 M-1 and Ka3D7 = 1.31E6 ± 1.21E6 M-1). CONCLUSIONS: Together, the study indicates that the genetic architecture of Pfrhs is principally shaped by purifying selection. The most abundant and ubiquitous PfRH5 variant harbouring 203Y, exhibits a greater affinity for BSG compared to PfRH53D7 possessing 203C allele. The study underscores the importance of selecting the functional allele that best represents circulating strains in natural parasite populations as vaccine targets.

Diversity analysis of MSP1 identifies conserved epitope organization in block 2 amidst high sequence variability in Indian Plasmodium falciparum isolates.

BACKGROUND: Despite its immunogenicity, the polymorphic nature of merozoite surface protein 1, an important vaccine candidate for Plasmodium falciparum malaria, remains a concern. This study analyses the impact of genetic variability and parasite population structure on epitope organization of different MSP1 segments. METHODS: Altogether 98 blood samples collected from P. falciparum infected mild and severe malaria patients of Chhattisgarh and West Bengal were used to sequence regions encoding block 2 and MSP1-19 of msp1. Sequences were analysed using MEGA7, DnaSPv5, Arlequin3.5 and BepiPred. RESULTS: All three major MSP1 block 2 allele families namely K1, MAD20 and RO33 were detected in the samples and they together resulted in 41 indel variants. Chhattisgarh samples displayed an average MOI of 2.07 ± 1.59 which was higher in mild malaria and in age group < 18 years. Ultra-structure of block 2 alleles revealed that mutation and repeat expansion were two major mechanisms responsible for allelic variability of K1 and MAD20. Regions flanking block 2 were highly variable in Chhattisgarh with average mismatch differences (k) ranging from 1.198 to 5.156 for three families. In contrast, region encompassing MSP1-19 exhibited limited heterogeneity (kChhattisgarh = 1.45, kWest Bengal = 1.363). Of the 50 possible B cell linear epitopes predicted from block 2 variants, 94.9% (131 of 138) of the parasites could be represented by three conserved antigens. CONCLUSIONS: Present data indicates that natural selection and transmission intensity jointly play a role in controlling allelic diversity of MSP1 in Indian parasite isolates. Despite remarkable genetic variability, a limited number of predominant and conserved epitopes are present in Indian parasite isolates reinstating the importance of MSP1 as a promising malaria vaccine candidate.

Deciphering genetic regulation of CD14 by SP1 through characterization of peripheral blood mononuclear transcriptome of P. faiciparum and P. vivax infected malaria patients.

BACKGROUND: Plasmodium falciparum and Plasmodium vivax are two major parasites responsible for malaria which remains a threat to almost 50% of world's population despite decade-long eradication program. One possible reason behind this conundrum is that the bases of clinical variability in malaria caused by either species are complex and poorly understood. METHODS: Whole-genome transcriptome was analyzed to identify the active and predominant pathways in the PBMC of P. falciparum and P. vivax infected malaria patients. Deregulated genes were identified and annotated using R Bioconductor and DAVID/KEGG respectively. Genetic and functional regulation of CD14, a prioritized candidate, were established by quantitative RT-PCR, genotyping using RFLP and resequencing, mapping of transcription factor binding using CONSITE and TFBIND, dual luciferase assay, western blot analysis, RNAi- mediated gene knockdown and chromatin-immunoprecipation. FINDINGS: The study highlighted that deregulation of host immune and inflammatory genes particularly CD14 as a key event in P. falciparum malaria. An abundance of allele-C of rs5744454, located in CD14 promoter, in severe malaria motivated us to establish an allele-specific regulation of CD14 by SP1. An enhancement of SP1 and CD14 expression was observed in artemisinin treated human monocyte cell line. INTERPRETATION: Our data not only reinstates that CD14 of TLR pathway plays a predominant role in P. falciparum malaria, it establishes a functional basis for genetic association of rs5744454 with P. falciparum severe malaria by demonstrating a cis-regulatory role of this promoter polymorphism. Moreover, the study points towards a novel pharmacogenetic aspect of artemisinin-based anti-malarial therapy. FUND: DST-SERB, Govt. of India, SR/SO/HS-0056/2013.

Natural selection and population genetic structure of domain-I of Plasmodium falciparum apical membrane antigen-1 in India.

Development of a vaccine against Plasmodium falciparum infection is an urgent priority particularly because of widespread resistance to most traditionally used drugs. Multiple evidences point to apical membrane antigen-1(AMA-1) as a prime vaccine candidate directed against P. falciparum asexual blood-stages. To gain understanding of the genetic and demographic forces shaping the parasite sequence diversity in Kolkata, a part of Pfama-1 gene covering domain-I was sequenced from 100 blood samples of malaria patients. Statistical and phylogenetic analyses of the sequences were performed using DnaSP and MEGA. Very high haplotype diversity was detected both at nucleotide (0.998±0.002) and amino-acid (0.996±0.001) levels. An abundance of low frequency polymorphisms (Tajima's D=-1.190, Fu & Li's D(∗) and F(∗)=-3.068 and -2.722), unimodal mismatch distribution and a star-like median-joining network of ama-1 haplotypes indicated a recent population expansion among Kolkata parasites. The high minimum number of recombination events (Rm=26) and a significantly high dN/dS of 3.705 (P<0.0001) in Kolkata suggested recombination and positive selection as major forces in the generation and maintenance of ama-1 allelic diversity. To evaluate the impact of observed non-synonymous substitutions in the context of AMA-1 functionality, PatchDock and FireDock protein-protein interaction solutions were mapped between PfAMA-1-PfRON2 and PfAMA-1-host IgNAR. Alterations in the desolvation and global energies of PfAMA-1-PfRON2 interaction complexes at the hotspot contact residues were observed together with redistribution of surface electrostatic potentials at the variant alleles with respect to referent Pf3D7 sequence. Finally, a comparison of P. falciparum subpopulations in five Indian regional isolates retrieved from GenBank revealed a significant level of genetic differentiation (FST=0.084-0.129) with respect to Kolkata sequences. Collectively, our results indicated a very high allelic and haplotype diversity, a high recombination rate and a signature of natural selection favoring accumulation of non-synonymous substitutions that facilitated PfAMA-1-PfRON2 interaction and hence parasite growth in Kolkata clinical isolates.

Nutritional status of Kamar tribal children in Chhattisgarh.

OBJECTIVE: Assessment of the nutritional status among 309 Kamar children (161 boys and 148 girls) aged 4-12 years. METHODS: In this cross-sectional study, 24-hours dietary recall method was used to assess dietary intakes of children. Height and weight were recorded and children were classified by WHO criterion (Z-score) using nutritional indices i.e. weight for age, height for age and weight for height. NCHS data were used as reference. Mean energy and protein intake per day were measured and compared with Recommended Dietary Allowances (RDA) of Indian standards. RESULTS: More than 90 percent of children (both boys and girls) in the age group of 4-6 years suffered by underweight (<-2 SD weight for age), which was comparatively lower in 7-9 and 10-12 years age group children. 84.51 percent of boys suffered by stunting (<-2 SD height for age), which was much higher than girls (47.54%) in 4-6 years age group. Similarly, 80 percent of 4-6 years age group children were affected by wasting (<-2 SD weight for height). The consumption of energy and protein were also much lower among Kamar children than the RDA of India throughout the ages. CONCLUSION: Undernutrition in the form of underweight, stunting and wasting and low consumption of dietary intake (energy and protein) was found to be widely prevalent among Kamar tribal children. Therefore, an urgent dietary intervention programme is necessary. Further studies are required to investigate into problem and to supplement the key nutrients which are required to ensure a good nutritional status in children.

Polarity and temporality of high-resolution y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists.

Although considerable cultural impact on social hierarchy and language in South Asia is attributable to the arrival of nomadic Central Asian pastoralists, genetic data (mitochondrial and Y chromosomal) have yielded dramatically conflicting inferences on the genetic origins of tribes and castes of South Asia. We sought to resolve this conflict, using high-resolution data on 69 informative Y-chromosome binary markers and 10 microsatellite markers from a large set of geographically, socially, and linguistically representative ethnic groups of South Asia. We found that the influence of Central Asia on the pre-existing gene pool was minor. The ages of accumulated microsatellite variation in the majority of Indian haplogroups exceed 10,000-15,000 years, which attests to the antiquity of regional differentiation. Therefore, our data do not support models that invoke a pronounced recent genetic input from Central Asia to explain the observed genetic variation in South Asia. R1a1 and R2 haplogroups indicate demographic complexity that is inconsistent with a recent single history. Associated microsatellite analyses of the high-frequency R1a1 haplogroup chromosomes indicate independent recent histories of the Indus Valley and the peninsular Indian region. Our data are also more consistent with a peninsular origin of Dravidian speakers than a source with proximity to the Indus and with significant genetic input resulting from demic diffusion associated with agriculture. Our results underscore the importance of marker ascertainment for distinguishing phylogenetic terminal branches from basal nodes when attributing ancestral composition and temporality to either indigenous or exogenous sources. Our reappraisal indicates that pre-Holocene and Holocene-era--not Indo-European--expansions have shaped the distinctive South Asian Y-chromosome landscape.