Comparative genome-wide analysis and evolutionary history of haemoglobin-processing and haem detoxification enzymes in malarial parasites

Patrath Ponsuwanna, Theerarat Kochakarn, Duangkamon Bunditvorapoom, Krittikorn Kümpornsin, Thomas D. Otto, Chase Ridenour, Kesinee Chotivanich, Prapon Wilairat, Nicholas J. White, Olivo Miotto, Thanat Chookajorn

Malaria parasites have evolved a series of intricate mechanisms to survive and propagate within host red blood cells. Intra-erythrocytic parasitism requires these organisms to digest haemoglobin and detoxify iron-bound haem. These tasks are executed by haemoglobin-specific proteases and haem biocrystallization factors that are components of a large multi-subunit complex. Since haemoglobin processing machineries are functionally and genetically linked to the modes of action and resistance mechanisms of several anti-malarial drugs, an understanding of their evolutionary history is important for drug development and drug resistance prevention. Maximum likelihood trees of genetic repertoires encoding haemoglobin processing machineries within Plasmodium species, and with the representatives of Apicomplexan species with various host tropisms, were created. Genetic variants were mapped onto existing three-dimensional structures. Genome-wide single nucleotide polymorphism data were used to analyse the selective pressure and the effect of these mutations at the structural level. Recent expansions in the falcipain and plasmepsin repertoires are unique to human malaria parasites especially in the Plasmodium falciparum and P. reichenowi lineage. Expansion of haemoglobin-specific plasmepsins occurred after the separation event of Plasmodium species, but the other members of the plasmepsin family were evolutionarily conserved with one copy for each sub-group in every Apicomplexan species. Haemoglobin-specific falcipains are separated from invasion-related falcipain, and their expansions within one specific locus arose independently in both P. falciparum and P. vivax lineages. Gene conversion between P. falciparum falcipain 2A and 2B was observed in artemisinin-resistant strains. Comparison between the numbers of non-synonymous and synonymous mutations suggests a strong selective pressure at falcipain and plasmepsin genes. The locations of amino acid changes from non-synonymous mutations mapped onto protein structures revealed clusters of amino acid residues in close proximity or near the active sites of proteases. A high degree of polymorphism at the haemoglobin processing genes implicates an imposition of selective pressure. The identification in recent years of functional redundancy of haemoglobin-specific proteases makes them less appealing as potential drug targets, but their expansions, especially in the human malaria parasite lineages, unequivocally point toward their functional significance during the independent and repetitive adaptation events in malaria parasite evolutionary history.