Low-complexity regions in Plasmodium falciparum: missing links in the evolution of an extreme genome

Cryo-EM structure of 1-deoxy-D-xylulose 5-phosphate synthase DXPS from Plasmodium falciparum reveals a distinct N-terminal domain

Plasmodium falciparum is the main causative agent of malaria, a deadly disease that mainly affects children under five years old. Artemisinin-based combination therapies have been pivotal in controlling the disease, but resistance has arisen in various regions, increasing the risk of treatment failure. The non-mevalonate pathway is essential for the isoprenoid synthesis in Plasmodium and provides several under-explored targets to be used in the discovery of new antimalarials. 1-deoxy-D-xylulose-5-phosphate synthase (DXPS) is the first and rate-limiting enzyme of the pathway. Despite its importance, there are no structures available for any Plasmodium spp., due to the complex sequence which contains large regions of high disorder, making crystallisation a difficult task. In this manuscript, we use cryo-electron microscopy to solve the P. falciparum DXPS structure at a final resolution of 2.42 Å. Overall, the structure resembles other DXPS enzymes but includes a distinct N-terminal domain exclusive to the Plasmodium genus. Mutational studies show that destabilization of the cap domain interface negatively impacts protein stability and activity. Additionally, a density for the co-factor thiamine diphosphate is found in the active site. Our work highlights the potential of cryo-EM to obtain structures of P. falciparum proteins that are unfeasible by means of crystallography.

A Plasmodium apicoplast-targeted unique exonuclease/FEN exhibits interspecies functional differences attributable to an insertion that alters DNA-binding

The human malaria parasite Plasmodium falciparum genome is among the most A + T rich, with low complexity regions (LCRs) inserted in coding sequences including those for proteins targeted to its essential relict plastid (apicoplast). Replication of the apicoplast genome (plDNA), mediated by the atypical multifunctional DNA polymerase PfPrex, would require additional enzymatic functions for lagging strand processing. We identified an apicoplast-targeted, [4Fe-4S]-containing, FEN/Exo (PfExo) with a long LCR insertion and detected its interaction with PfPrex. Distinct from other known exonucleases across organisms, PfExo recognized a wide substrate range; it hydrolyzed 5'-flaps, processed dsDNA as a 5'-3' exonuclease, and was a bipolar nuclease on ssDNA and RNA-DNA hybrids. Comparison with the rodent P. berghei ortholog PbExo, which lacked the insertion and [4Fe-4S], revealed interspecies functional differences. The insertion-deleted PfExoΔins behaved like PbExo with a limited substrate repertoire because of compromised DNA binding. Introduction of the PfExo insertion into PbExo led to gain of activities that the latter initially lacked. Knockout of PbExo indicated essentiality of the enzyme for survival. Our results demonstrate the presence of a novel apicoplast exonuclease with a functional LCR that diversifies substrate recognition, and identify it as the candidate flap-endonuclease and RNaseH required for plDNA replication and maintenance.

Genetic Sequence Variation in the Plasmodium falciparum Histidine-Rich Protein 2 Gene from Field Isolates in Tanzania: Impact on Malaria Rapid Diagnosis

Malaria rapid diagnosis test (RDT) is crucial for managing the disease, and the effectiveness of detection depends on parameters such as sensitivity and specificity of the RDT. Several factors can affect the performance of RDT. In this study, we focused on the pfhrp2 sequence variation and its impact on RDTs targeted by antigens encoded by Plasmodium falciparum histidine-rich protein 2 (pfhrp2). Field samples collected during cross-sectional surveys in Tanzania were sequenced to investigate the pfhrp2 sequence diversity and evaluate the impact on HRP2-based RDT performance. We observed significant mean differences in amino acid repeats between current and previous studies. Several new amino acid repeats were found to occur at different frequencies, including types AAY, AHHAHHAAN, and AHHAA. Based on the abundance of types 2 and 7 amino acid repeats, the binary predictive model was able to predict RDT insensitivity by about 69% in the study area. About 85% of the major epitopes targeted by monoclonal antibodies (MAbs) in RDT were identified. Our study suggested that the extensive sequence variation in pfhrp2 can contribute to reduced RDT sensitivity. The correlation between the different combinations of amino acid repeats and the performance of RDT in different malaria transmission settings should be investigated further.

Deletion of the Plasmodium falciparum exported protein PTP7 leads to Maurer’s clefts vesiculation, host cell remodeling defects, and loss of surface presentation of EMP1

Presentation of the variant antigen, Plasmodium falciparum erythrocyte membrane protein 1 (EMP1), at knob-like protrusions on the surface of infected red blood cells, underpins the parasite’s pathogenicity. Here we describe a protein PF3D7_0301700 (PTP7), that functions at the nexus between the intermediate trafficking organelle, the Maurer’s cleft, and the infected red blood cell surface. Genetic disruption of PTP7 leads to accumulation of vesicles at the Maurer’s clefts, grossly aberrant knob morphology, and failure to deliver EMP1 to the red blood cell surface. We show that an expanded low complexity sequence in the C-terminal region of PTP7, identified only in the Laverania clade of Plasmodium, is critical for efficient virulence protein trafficking.

We describe a malaria parasite protein, PTP7, involved in virulence factor trafficking that is associated with Maurer’s clefts and other trafficking compartments. Upon disruption of the PTP7 locus, the Maurer’s clefts become decorated with vesicles; the knobby protrusions on the host red blood cell surface are fewer and distorted; and trafficking of the virulence protein, EMP1, to the host red blood cell surface is ablated. We provide evidence that a region of PTP7 with low sequence complexity plays an important role in virulence protein trafficking from the Maurer’s clefts.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

The relationship between protein domains and homopeptides in the Plasmodium falciparum proteome

The proteome of the malaria parasite Plasmodium falciparum is notable for the pervasive occurrence of homopeptides or low-complexity regions (i.e., regions that are made from a small subset of amino-acid residue types). The most prevalent of these are made from residues encoded by adenine/thymidine (AT)-rich codons, in particular asparagine. We examined homopeptide occurrences within protein domains in P. falciparum. Homopeptide enrichments occur for hydrophobic (e.g., valine), or small residues (alanine or glycine) in short spans (<5 residues), but these enrichments disappear for longer lengths. We observe that short asparagine homopeptides (<10 residues long) have a dramatic relative depletion inside protein domains, indicating some selective constraint to keep them from forming. We surmise that this is possibly linked to co-translational protein folding, although there are specific protein domains that are enriched in longer asparagine homopeptides (≥10 residues) indicating a functional linkage for specific poly-asparagine tracts. Top gene ontology functional category enrichments for homopeptides associated with diverse protein domains include “vesicle-mediated transport”, and “DNA-directed 5′-3′ RNA polymerase activity”, with various categories linked to “binding” evidencing significant homopeptide depletions. Also, in general homopeptides are substantially enriched in the parts of protein domains that are near/in IDRs. The implications of these findings are discussed.

Plasmodium falciparum translational machinery condones polyadenosine repeats

Plasmodium falciparum is a causative agent of human malaria. Sixty percent of mRNAs from its extremely AT-rich (81%) genome harbor long polyadenosine (polyA) runs within their ORFs, distinguishing the parasite from its hosts and other sequenced organisms. Recent studies indicate polyA runs cause ribosome stalling and frameshifting, triggering mRNA surveillance pathways and attenuating protein synthesis. Here, we show that P. falciparum is an exception to this rule. We demonstrate that both endogenous genes and reporter sequences containing long polyA runs are efficiently and accurately translated in P. falciparum cells. We show that polyA runs do not elicit any response from No Go Decay (NGD) or result in the production of frameshifted proteins. This is in stark contrast to what we observe in human cells or T. thermophila, an organism with similar AT-content. Finally, using stalling reporters we show that Plasmodium cells evolved not to have a fully functional NGD pathway.

Polymorphic markers for identification of parasite population in Plasmodium malariae

Background

Molecular genotyping in Plasmodium serves many aims including providing tools for studying parasite population genetics and distinguishing recrudescence from reinfection. Microsatellite typing, insertion-deletion (INDEL) and single nucleotide polymorphisms is used for genotyping, but only limited information is available for Plasmodium malariae, an important human malaria species. This study aimed to provide a set of genetic markers to facilitate the study of P. malariae population genetics.

Methods

Markers for microsatellite genotyping and pmmsp1 gene polymorphisms were developed and validated in symptomatic P. malariae field isolates from Myanmar (N = 37). Fragment analysis was used to determine allele sizes at each locus to calculate multiplicity of infections (MOI), linkage disequilibrium, heterozygosity and construct dendrograms. Nucleotide diversity (π), number of haplotypes, and genetic diversity (H_d) were assessed and a phylogenetic tree was constructed. Genome-wide microsatellite maps with annotated regions of newly identified markers were constructed.

Results

Six microsatellite markers were developed and tested in 37 P. malariae isolates which showed sufficient heterozygosity (0.530–0.922), and absence of linkage disequilibrium (I_A^S=0.03, p value > 0.05) (N = 37). In addition, a tandem repeat (VNTR)-based pmmsp1 INDEL polymorphisms marker was developed and assessed in 27 P. malariae isolates showing a nucleotide diversity of 0.0976, haplotype gene diversity of 0.698 and identified 14 unique variants. The size of VNTR consensus repeat unit adopted as allele was 27 base pairs. The markers Pm12_426 and pmmsp1 showed greatest diversity with heterozygosity scores of 0.920 and 0.835, respectively. Using six microsatellites markers, the likelihood that any two parasite strains would have the same microsatellite genotypes was 8.46 × 10⁻⁴ and was further reduced to 1.66 × 10⁻⁴ when pmmsp1 polymorphisms were included.

Conclusions

Six novel microsatellites genotyping markers and a set of pmmsp1 VNTR-based INDEL polymorphisms markers for P. malariae were developed and validated. Each marker could be independently or in combination employed to access genotyping of the parasite. The newly developed markers may serve as a useful tool for investigating parasite diversity, population genetics, molecular epidemiology and for distinguishing recrudescence from reinfection in drug efficacy studies.

Adaptation of Translational Machinery in Malaria Parasites to Accommodate Translation of Poly-Adenosine Stretches Throughout Its Life Cycle

Malaria is caused by unicellular apicomplexan parasites of the genus Plasmodium, which includes the major human parasite Plasmodium falciparum. The complex cycle of the malaria parasite in both mosquito and human hosts has been studied extensively. There is tight control of gene expression in each developmental stage, and at every level of gene synthesis: from RNA transcription, to its subsequent translation, and finally post-translational modifications of the resulting protein. Whole-genome sequencing of P. falciparum has laid the foundation for significant biological advances by revealing surprising genomic information. The P. falciparum genome is extremely AT-rich (∼80%), with a substantial portion of genes encoding intragenic polyadenosine (polyA) tracks being expressed throughout the entire parasite life cycle. In most eukaryotes, intragenic polyA runs act as negative regulators of gene expression. Recent studies have shown that translation of mRNAs containing 12 or more consecutive adenosines results in ribosomal stalling and frameshifting; activating mRNA surveillance mechanisms. In contrast, P. falciparum translational machinery can efficiently and accurately translate polyA tracks without activating mRNA surveillance pathways. This unique feature of P. falciparum raises interesting questions: (1) How is P. falciparum able to efficiently and correctly translate polyA track transcripts, and (2) What are the specifics of the translational machinery and mRNA surveillance mechanisms that separate P. falciparum from other organisms? In this review, we analyze possible evolutionary shifts in P. falciparum protein synthesis machinery that allow efficient translation of an AU rich-transcriptome. We focus on physiological and structural differences of P. falciparum stage specific ribosomes, ribosome-associated proteins, and changes in mRNA surveillance mechanisms throughout the complete parasite life cycle, with an emphasis on the mosquito and liver stages.

Evidence for diversifying selection of genetic regions of encoding putative collagen-like host-adhesive fibers in Pasteuria penetrans

Pasteuria spp. belong to a group of genetically diverse endospore-forming bacteria (phylum: Firmicutes) that are known to parasitize plant-parasitic nematodes and water fleas (Daphnia spp.). Collagen-like fibres form the nap on the surface of endospores and the genes encoding these sequences have been hypothesised to be involved in the adhesion of the endospores of Pasteuria spp. to their hosts. We report a group of 17 unique collagen-like genes putatively encoded by Pasteuria penetrans (strain: Res148) that formed five different phylogenetic clusters and suggest that collagen-like proteins are an important source of genetic diversity in animal pathogenic Firmicutes including Pasteuria. Additionally, and unexpectedly, we identified a putative collagen-like sequence which had a very different sequence structure to the other collagen-like proteins but was similar to the protein sequences in Megaviruses that are involved in host-parasite interactions. We, therefore, suggest that these diverse endospore surface proteins in Pasteuria are involved in biological functions, such as cellular adhesion; however, they are not of monophyletic origin and were possibly obtained de novo by mutation or possibly through selection acting upon several historic horizontal gene transfer events.

De Novo Mutations Resolve Disease Transmission Pathways in Clonal Malaria

Detecting de novo mutations in viral and bacterial pathogens enables researchers to reconstruct detailed networks of disease transmission and is a key technique in genomic epidemiology. However, these techniques have not yet been applied to the malaria parasite, Plasmodium falciparum, in which a larger genome, slower generation times, and a complex life cycle make them difficult to implement. Here, we demonstrate the viability of de novo mutation studies in P. falciparum for the first time. Using a combination of sequencing, library preparation, and genotyping methods that have been optimized for accuracy in low-complexity genomic regions, we have detected de novo mutations that distinguish nominally identical parasites from clonal lineages. Despite its slower evolutionary rate compared with bacterial or viral species, de novo mutation can be detected in P. falciparum across timescales of just 1–2 years and evolutionary rates in low-complexity regions of the genome can be up to twice that detected in the rest of the genome. The increased mutation rate allows the identification of separate clade expansions that cannot be found using previous genomic epidemiology approaches and could be a crucial tool for mapping residual transmission patterns in disease elimination campaigns and reintroduction scenarios.

Tail Wags the Dog? Functional Gene Classes Driving Genome-Wide GC Content in Plasmodium spp

Plasmodium parasites are valuable models to understand how nucleotide composition affects mutation, diversification, and adaptation. No other observed eukaryotes have undergone such large changes in genomic Guanine-Cytosine (GC) content as seen in the genus Plasmodium (∼30% within 35-40 Myr). Although mutational biases are known to influence GC content in the human-infective Plasmodium vivax and Plasmodium falciparum; no study has addressed how different gene functional classes contribute to genus-wide compositional changes, or if Plasmodium GC content variation is driven by natural selection. Here, we tested the hypothesis that certain gene processes and functions drive variation in global GC content between Plasmodium species. We performed a large-scale comparative genomic analysis using the genomes and predicted genes of 17 Plasmodium species encompassing a wide genomic GC content range. Genic GC content was sorted and divided into ten equally sized quantiles that were then assessed for functional enrichment classes. In agreement that selection on gene classes may drive genomic GC content, trans-membrane proteins were enriched within extreme GC content quantiles (Q1 and Q10). Specifically, variant surface antigens, which primarily interact with vertebrate immune systems, showed skewed GC content distributions compared with other trans-membrane proteins. Although a definitive causation linking GC content, expression, and positive selection within variant surface antigens from Plasmodium vivax, Plasmodium berghei, and Plasmodium falciparum could not be established, we found that regardless of genomic nucleotide composition, genic GC content and expression were positively correlated during trophozoite stages. Overall, these data suggest that, alongside mutational biases, functional protein classes drive Plasmodium GC content change.

Disordered epitopes as peptide vaccines

The development of clinically useful peptide-based vaccines remains a long-standing goal. This review highlights that intrinsically disordered protein antigens, which lack an ordered three-dimensional structure, represent excellent starting points for the development of such vaccines. Disordered proteins represent an important class of antigen in a wide range of human pathogens, and, contrary to widespread belief, they are frequently targets of protective antibody responses. Importantly, disordered epitopes appear invariably to be linear epitopes, rendering them ideally suited to incorporation into a peptide vaccine. Nonetheless, the conformational properties of disordered antigens, and hence their recognition by antibodies, frequently depend on the interactions they make and the context in which they are presented to the immune system. These effects must be considered in the design of an effective vaccine. Here we discuss these issues and propose design principles that may facilitate the development of peptide vaccines targeting disordered antigens.

Repetitive sequences in malaria parasite proteins

Five species of parasite cause malaria in humans with the most severe disease caused by Plasmodium falciparum. Many of the proteins encoded in the P. falciparum genome are unusually enriched in repetitive low-complexity sequences containing a limited repertoire of amino acids. These repetitive sequences expand and contract dynamically and are among the most rapidly changing sequences in the genome. The simplest repetitive sequences consist of single amino acid repeats such as poly-asparagine tracts that are found in approximately 25% of P. falciparum proteins. More complex repeats of two or more amino acids are also common in diverse parasite protein families. There is no universal explanation for the occurrence of repetitive sequences and it is possible that many confer no function to the encoded protein and no selective advantage or disadvantage to the parasite. However, there are increasing numbers of examples where repetitive sequences are important for parasite protein function. We discuss the diverse roles of low-complexity repetitive sequences throughout the parasite life cycle, from mediating protein-protein interactions to enabling the parasite to evade the host immune system.

Diverse single-amino-acid repeat profiles in the genus Cryptosporidium

Genome sequencing has greatly contributed to our understanding of parasitic protozoa. This is particularly the case for Cryptosporidium species (phylum Apicomplexa) which are difficult to propagate. Because of their polymorphic nature, simple sequence repeats have been used extensively as genotypic markers to differentiate between isolates, but no global analysis of amino acid repeats in Cryptosporidium genomes has been reported. Taking advantage of several newly sequenced Cryptosporidium genomes, a comparative analysis of single-amino-acid repeats (SAARs) in seven species was undertaken. This analysis revealed a striking difference between the SAAR profile of the gastric and intestinal species which infect mammals and one species which infects birds. In average, total SAAR length in gastric species is only 25% of the cumulative SAAR length in the genome of Cryptosporidium parvum, Cryptosporidium hominis and Cryptosporidium meleagridis, species infectious to humans. The SAAR profile in the avian parasite Cryptosporidium baileyi stands out due to the presence of long asparagine repeats. Cryptosporidium baileyi proteins with repeats ⩾20 residues are significantly enriched in regulatory functions. As postulated for the related apicomplexan species Plasmodium falciparum, these observations suggest that Cryptosporidium SAARs evolve in response to selective pressure. The putative selective mechanisms driving SAAR evolution in Cryptosporidium species are unknown.

Plasmodium parasites of birds have the most AT-rich genes of eukaryotes

The genomic architecture of organisms, including nucleotide composition, can be highly variable, even among closely-related species. To better understand the causes leading to structural variation in genomes, information on distinct and diverse genomic features is needed. Malaria parasites are known for encompassing a wide range of genomic GC-content and it has long been thought that Plasmodium falciparum, the virulent malaria parasite of humans, has the most AT-biased eukaryotic genome. Here, I perform comparative genomic analyses of the most AT-rich eukaryotes sequenced to date, and show that the avian malaria parasites Plasmodium gallinaceum, P. ashfordi, and P. relictum have the most extreme coding sequences in terms of AT-bias. Their mean GC-content is 21.21, 21.22 and 21.60 %, respectively, which is considerably lower than the transcriptome of P. falciparum (23.79 %) and other eukaryotes. This information enables a better understanding of genome evolution and raises the question of how certain organisms are able to prosper despite severe compositional constraints.

Comparative analysis of low complexity regions in Plasmodia

Low complexity regions (LCRs) are a common feature shared by many genomes, but their evolutionary and functional significance remains mostly unknown. At the core of the uncertainty is a poor understanding of the mechanisms that regulate their retention in genomes, whether driven by natural selection or neutral evolution. Applying a comparative approach of LCRs to multiple strains and species is a powerful approach to identify patterns of conservation in these regions. Using this method, we investigate the evolutionary history of LCRs in the genus Plasmodium based on orthologous protein coding genes shared by 11 species and strains from primate and rodent-infecting pathogens. We find multiple lines of evidence in support of natural selection as a major evolutionary force shaping the composition and conservation of LCRs through time and signatures that their evolutionary paths are species specific. Our findings add a comparative analysis perspective to the debate on the evolution of LCRs and harness the power of sequence comparisons to identify potential functionally important LCR candidates.

Extreme mutation bias and high AT content in Plasmodium falciparum

For reasons that remain unknown, the Plasmodium falciparum genome has an exceptionally high AT content compared to other Plasmodium species and eukaryotes in general - nearly 80% in coding regions and approaching 90% in non-coding regions. Here, we examine how this phenomenon relates to genome-wide patterns of de novo mutation. Mutation accumulation experiments were performed by sequential cloning of six P. falciparum isolates growing in human erythrocytes in vitro for 4 years, with 279 clones sampled for whole genome sequencing at different time points. Genome sequence analysis of these samples revealed a significant excess of G:C to A:T transitions compared to other types of nucleotide substitution, which would naturally cause AT content to equilibrate close to the level seen across the P. falciparum reference genome (80.6% AT). These data also uncover an extremely high rate of small indel mutation relative to other species, primarily associated with repetitive AT-rich sequences, in addition to larger-scale structural rearrangements focused in antigen-coding var genes. In conclusion, high AT content in P. falciparum is driven by a systematic mutational bias and ultimately leads to an unusual level of microstructural plasticity, raising the question of whether this contributes to adaptive evolution.

Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum

The malaria parasite Plasmodium falciparum has a great capacity for evolutionary adaptation to evade host immunity and develop drug resistance. Current understanding of parasite evolution is impeded by the fact that a large fraction of the genome is either highly repetitive or highly variable and thus difficult to analyze using short-read sequencing technologies. Here, we describe a resource of deep sequencing data on parents and progeny from genetic crosses, which has enabled us to perform the first genome-wide, integrated analysis of SNP, indel and complex polymorphisms, using Mendelian error rates as an indicator of genotypic accuracy. These data reveal that indels are exceptionally abundant, being more common than SNPs and thus the dominant mode of polymorphism within the core genome. We use the high density of SNP and indel markers to analyze patterns of meiotic recombination, confirming a high rate of crossover events and providing the first estimates for the rate of non-crossover events and the length of conversion tracts. We observe several instances of meiotic recombination within copy number variants associated with drug resistance, demonstrating a mechanism whereby fitness costs associated with resistance mutations could be compensated and greater phenotypic plasticity could be acquired.

Profiles of low complexity regions in Apicomplexa

BACKGROUND

Low complexity regions (LCRs) are a ubiquitous feature in genomes and yet their evolutionary history and functional roles are unclear. Previous studies have shown contrasting evidence in favor of both neutral and selective mechanisms of evolution for different sets of LCRs suggesting that modes of identification of these regions may play a role in our ability to discern their evolutionary history. To further investigate this issue, we used a multiple threshold approach to identify species-specific profiles of proteome complexity and, by comparing properties of these sets, determine the influence that starting parameters have on evolutionary inferences.

RESULTS

We find that, although qualitatively similar, quantitatively each species has a unique LCR profile which represents the frequency of these regions within each genome. Inferences based on these profiles are more accurate in comparative analyses of genome complexity as they allow to determine the relative complexity of multiple genomes as well as the type of repetitiveness that is most common in each. Based on the multiple threshold LCR sets obtained, we identified predominant evolutionary mechanisms at different complexity levels, which show neutral mechanisms acting on highly repetitive LCRs (e.g., homopolymers) and selective forces becoming more important as heterogeneity of the LCRs increases.

CONCLUSIONS

Our results show how inferences based on LCRs are influenced by the parameters used to identify these regions. Sets of LCRs are heterogeneous aggregates of regions that include homo- and heteropolymers and, as such, evolve according to different mechanisms. LCR profiles provide a new way to investigate genome complexity across species and to determine the driving mechanism of their evolution.

Antibody Recognition of Disordered Antigens

Disordered proteins are important antigens in a range of infectious diseases. Little is known, however, about the molecular details of recognition of disordered antigens by their cognate antibodies. Using a large dataset of protein antigens, we show that disordered epitopes are as likely to be recognized by antibodies as ordered epitopes. Moreover, the affinity with which antigens are recognized is, unexpectedly, only weakly dependent on the degree of disorder within the epitope. Structurally defined complexes of ordered and disordered protein antigens with their cognate antibodies reveal that disordered epitopes are smaller than their ordered counterparts, but are more efficient in their interactions with antibody. Our results demonstrate that disordered antigens are bona fide targets of antibody recognition, and that recognition of disordered epitopes is particularly sensitive to epitope variation, a finding with implications for the effects of disorder on the specificity of molecular recognition more generally.

Characterization of the differences in the cyclopiazonic acid binding mode to mammalian and P. Falciparum Ca2+ pumps: a computational study

Despite the investments in malaria research, an effective vaccine has not yet been developed and the causative parasites are becoming increasingly resistant to most of the available drugs. PfATP6, the sarco/endoplasmic reticulum Ca2+ pump (SERCA) of P. falciparum, has been recently genetically validated as a potential antimalarial target and cyclopiazonic acid (CPA) has been found to be a potent inhibitor of SERCAs in several organisms, including P. falciparum. In position 263, PfATP6 displays a leucine residue, whilst the corresponding position in the mammalian SERCA is occupied by a glutamic acid. The PfATP6 L263E mutation has been studied in relation to the artemisinin inhibitory effect on P. falciparum and recent studies have provided evidence that the parasite with this mutation is more susceptible to CPA. Here, we characterized, for the first time, the interaction of CPA with PfATP6 and its mammalian counterpart to understand similarities and differences in the mode of binding of the inhibitor to the two Ca2+ pumps. We found that, even though CPA does not directly interact with the residue in position 263, the presence of a hydrophobic residue in this position in PfATP6 rather than a negatively charged one, as in the mammalian SERCA, entails a conformational arrangement of the binding pocket which, in turn, determines a relaxation of CPA leading to a different binding mode of the compound. Our findings highlight differences between the plasmodial and human SERCA CPA-binding pockets that may be exploited to design CPA derivatives more selective toward PfATP6. Proteins 2015; 83:564–574. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Asparagine requirement in Plasmodium berghei as a target to prevent malaria transmission and liver infections

The proteins of Plasmodium, the malaria parasite, are strikingly rich in asparagine. Plasmodium depends primarily on host haemoglobin degradation for amino acids and has a rudimentary pathway for amino acid biosynthesis, but retains a gene encoding asparagine synthetase (AS). Here we show that deletion of AS in Plasmodium berghei (Pb) delays the asexual- and liver-stage development with substantial reduction in the formation of ookinetes, oocysts and sporozoites in mosquitoes. In the absence of asparagine synthesis, extracellular asparagine supports suboptimal survival of PbAS knockout (KO) parasites. Depletion of blood asparagine levels by treating PbASKO-infected mice with asparaginase completely prevents the development of liver stages, exflagellation of male gametocytes and the subsequent formation of sexual stages. In vivo supplementation of asparagine in mice restores the exflagellation of PbASKO parasites. Thus, the parasite life cycle has an absolute requirement for asparagine, which we propose could be targeted to prevent malaria transmission and liver infections.

Malaria parasites obtain amino acids primarily from the host, but possess a gene encoding a putative asparagine synthetase. Here, the authors show that this enzyme is functional and that asparagine is crucial for the development of the parasite's sexual stages in mosquitoes and liver stages in mice.

Structural basis for epitope masking and strain specificity of a conserved epitope in an intrinsically disordered malaria vaccine candidate

Merozoite surface protein 2 (MSP2) is an intrinsically disordered, membrane-anchored antigen of the malaria parasite Plasmodium falciparum. MSP2 can elicit a protective, albeit strain-specific, antibody response in humans. Antibodies are generated to the conserved N- and C-terminal regions but many of these react poorly with the native antigen on the parasite surface. Here we demonstrate that recognition of a conserved N-terminal epitope by mAb 6D8 is incompatible with the membrane-bound conformation of that region, suggesting a mechanism by which native MSP2 escapes antibody recognition. Furthermore, crystal structures and NMR spectroscopy identify transient, strain-specific interactions between the 6D8 antibody and regions of MSP2 beyond the conserved epitope. These interactions account for the differential affinity of 6D8 for the two allelic families of MSP2, even though 6D8 binds to a fully conserved epitope. These results highlight unappreciated mechanisms that may modulate the specificity and efficacy of immune responses towards disordered antigens.

Structure of Plasmodium falciparum orotate phosphoribosyltransferase with autologous inhibitory protein-protein interactions

The most severe form of malaria is caused by the obligate parasite Plasmodium falciparum. Orotate phosphoribosyltransferase (OPRTase) is the fifth enzyme in the de novo pyrimidine-synthesis pathway in the parasite, which lacks salvage pathways. Among all of the malaria de novo pyrimidine-biosynthesis enzymes, the structure of P. falciparum OPRTase (PfOPRTase) was the only one unavailable until now. PfOPRTase that could be crystallized was obtained after some low-complexity sequences were removed. Four catalytic dimers were seen in the asymmetic unit (a total of eight polypeptides). In addition to revealing unique amino acids in the PfOPRTase active sites, asymmetric dimers in the larger structure pointed to novel parasite-specific protein-protein interactions that occlude the catalytic active sites. The latter could potentially modulate PfOPRTase activity in parasites and possibly provide new insights for blocking PfOPRTase functions.

Conformational dynamics and antigenicity in the disordered malaria antigen merozoite surface protein 2

Merozoite surface protein 2 (MSP2) of Plasmodium falciparum is an abundant, intrinsically disordered protein that is GPI-anchored to the surface of the invasive blood stage of the malaria parasite. Recombinant MSP2 has been trialled as a component of a malaria vaccine, and is one of several disordered proteins that are candidates for inclusion in vaccines for malaria and other diseases. Nonetheless, little is known about the implications of protein disorder for the development of an effective antibody response. We have therefore undertaken a detailed analysis of the conformational dynamics of the two allelic forms of MSP2 (3D7 and FC27) using NMR spectroscopy. Chemical shifts and NMR relaxation data indicate that conformational and dynamic properties of the N- and C-terminal conserved regions in the two forms of MSP2 are essentially identical, but significant variation exists between and within the central variable regions. We observe a strong relationship between the conformational dynamics and the antigenicity of MSP2, as assessed with antisera to recombinant MSP2. Regions of increased conformational order in MSP2, including those in the conserved regions, are more strongly antigenic, while the most flexible regions are minimally antigenic. This suggests that modifications that increase conformational order may offer a means to tune the antigenicity of MSP2 and other disordered antigens, with implications for vaccine design.

Genome-wide patterns of genetic polymorphism and signatures of selection in Plasmodium vivax

Plasmodium vivax is the most prevalent human malaria parasite outside of Africa. Yet, studies aimed to identify genes with signatures consistent with natural selection are rare. Here, we present a comparative analysis of the pattern of genetic variation of five sequenced isolates of P. vivax and its divergence with two closely related species, Plasmodium cynomolgi and Plasmodium knowlesi, using a set of orthologous genes. In contrast to Plasmodium falciparum, the parasite that causes the most lethal form of human malaria, we did not find significant constraints on the evolution of synonymous sites genome wide in P. vivax. The comparative analysis of polymorphism and divergence across loci allowed us to identify 87 genes with patterns consistent with positive selection, including genes involved in the “exportome” of P. vivax, which are potentially involved in evasion of the host immune system. Nevertheless, we have found a pattern of polymorphism genome wide that is consistent with a significant amount of constraint on the replacement changes and prevalent negative selection. Our analyses also show that silent polymorphism tends to be larger toward the ends of the chromosomes, where many genes involved in antigenicity are located, suggesting that natural selection acts not only by shaping the patterns of variation within the genes but it also affects genome organization.

Increased substitution rates surrounding low-complexity regions within primate proteins

Previous studies have found that DNA-flanking low-complexity regions (LCRs) have an increased substitution rate. Here, the substitution rate was confirmed to increase in the vicinity of LCRs in several primate species, including humans. This effect was also found among human sequences from the 1000 Genomes Project. A strong correlation was found between average substitution rate per site and distance from the LCR, as well as the proportion of genes with gaps in the alignment at each site and distance from the LCR. Along with substitution rates, d_N/d_S ratios were also determined for each site, and the proportion of sites undergoing negative selection was found to have a negative relationship with distance from the LCR.

Homopolymer tract organization in the human malarial parasite Plasmodium falciparum and related Apicomplexan parasites

BACKGROUND

Homopolymeric tracts, particularly poly dA.dT, are enriched within the intergenic sequences of eukaryotic genomes where they appear to act as intrinsic regulators of nucleosome positioning. A previous study of the incomplete genome of the human malarial parasite Plasmodium falciparum reports a higher than expected enrichment of poly dA.dT tracts, far above that anticipated even in this highly AT rich genome. Here we report an analysis of the relative frequency, length and spatial arrangement of homopolymer tracts for the complete P. falciparum genome, extending this analysis to twelve additional genomes of Apicomplexan parasites important to human and animal health. In addition, using nucleosome-positioning data available for P. falciparum, we explore the correlation of poly dA.dT tracts with nucleosome-positioning data over key expression landmarks within intergenic regions.

RESULTS

We describe three apparent lineage-specific patterns of homopolymeric tract organization within the intergenic regions of these Apicomplexan parasites. Moreover, a striking pattern of enrichment of overly long poly dA.dT tracts in the intergenic regions of Plasmodium spp. uniquely extends into protein coding sequences. There is a conserved spatial arrangement of poly dA.dT immediately flanking open reading frames and over predicted core promoter sites. These key landmarks are all relatively depleted in nucleosomes in P. falciparum, as would be expected for poly dA.dT acting as nucleosome exclusion sequences.

CONCLUSIONS

Previous comparative studies of homopolymer tract organization emphasize evolutionary diversity; this is the first report of such an analysis within a single phylum. Our data provide insights into the evolution of homopolymeric tracts and the selective pressures at play in their maintenance and expansion.

Protective epitopes of the Plasmodium falciparum SERA5 malaria vaccine reside in intrinsically unstructured N-terminal repetitive sequences

The malaria vaccine candidate antigen, SE36, is based on the N-terminal 47 kDa domain of Plasmodium falciparum serine repeat antigen 5 (SERA5). In epidemiological studies, we have previously shown the inhibitory effects of SE36 specific antibodies on in vitro parasite growth and the negative correlation between antibody level and malaria symptoms. A phase 1 b trial of the BK-SE36 vaccine in Uganda elicited 72% protective efficacy against symptomatic malaria in children aged 6–20 years during the follow-up period 130–365 days post–second vaccination. Here, we performed epitope mapping with synthetic peptides covering the whole sequence of SE36 to identify and map dominant epitopes in Ugandan adult serum presumed to have clinical immunity to P. falciparum malaria. High titer sera from the Ugandan adults predominantly reacted with peptides corresponding to two successive N-terminal regions of SERA5 containing octamer repeats and serine rich sequences, regions of SERA5 that were previously reported to have limited polymorphism. Affinity purified antibodies specifically recognizing the octamer repeats and serine rich sequences exhibited a high antibody-dependent cellular inhibition (ADCI) activity that inhibited parasite growth. Furthermore, protein structure predictions and structural analysis of SE36 using spectroscopic methods indicated that N-terminal regions possessing inhibitory epitopes are intrinsically unstructured. Collectively, these results suggest that strict tertiary structure of SE36 epitopes is not required to elicit protective antibodies in naturally immune Ugandan adults.

The origin and diversification of the merozoite surface protein 3 (msp3) multi-gene family in Plasmodium vivax and related parasites

The genus Plasmodium is a diversified group of parasites with more than 200 known species that includes those causing malaria in humans. These parasites use numerous proteins in a complex process that allows them to invade the red blood cells of their vertebrate hosts. Many of those proteins are part of multi-gene families; one of which is the merozoite surface protein-3 (msp3) family. The msp3 multi-gene family is considered important in the two main human parasites, Plasmodium vivax and Plasmodium falciparum, as its paralogs are simultaneously expressed in the blood stage (merozoite) and are immunogenic. There are large differences among Plasmodium species in the number of paralogs in this family. Such differences have been previously explained, in part, as adaptations that allow the different Plasmodium species to invade their hosts. To investigate this, we characterized the array containing msp3 genes among several Plasmodium species, including P. falciparum and P. vivax. We first found no evidence indicating that the msp3 family of P. falciparum was homologous to that of P. vivax. Subsequently, by focusing on the diverse clade of nonhuman primate parasites to which P. vivax is closely related, where homology was evident, we found no evidence indicating that the interspecies variation in the number of paralogs was an adaptation related to changes in host range or host switches. Overall, we hypothesize that the evolution of the msp3 family in P. vivax is consistent with a model of multi-allelic diversifying selection where the paralogs may have functionally redundant roles in terms of increasing antigenic diversity. Thus, we suggest that the expressed MSP3 proteins could serve as "decoys", via antigenic diversity, during the critical process of invading the host red blood cells.

Small molecule screen for candidate antimalarials targeting Plasmodium Kinesin-5

Plasmodium falciparum and vivax are responsible for the majority of malaria infections worldwide, resulting in over a million deaths annually. Malaria parasites now show measured resistance to all currently utilized drugs. Novel antimalarial drugs are urgently needed. The Plasmodium Kinesin-5 mechanoenzyme is a suitable "next generation" target. Discovered via small molecule screen experiments, the human Kinesin-5 has multiple allosteric sites that are "druggable." One site in particular, unique in its sequence divergence across all homologs in the superfamily and even within the same family, exhibits exquisite drug specificity. We propose that Plasmodium Kinesin-5 shares this allosteric site and likewise can be targeted to uncover inhibitors with high specificity. To test this idea, we performed a screen for inhibitors selective for Plasmodium Kinesin-5 ATPase activity in parallel with human Kinesin-5. Our screen of nearly 2000 compounds successfully identified compounds that selectively inhibit both P. vivax and falciparum Kinesin-5 motor domains but, as anticipated, do not impact human Kinesin-5 activity. Of note is a candidate drug that did not biochemically compete with the ATP substrate for the conserved active site or disrupt the microtubule-binding site. Together, our experiments identified MMV666693 as a selective allosteric inhibitor of Plasmodium Kinesin-5; this is the first identified protein target for the Medicines of Malaria Venture validated collection of parasite proliferation inhibitors. This work demonstrates that chemical screens against human kinesins are adaptable to homologs in disease organisms and, as such, extendable to strategies to combat infectious disease.

Structural polymorphism in the promoter of pfmrp2 confers Plasmodium falciparum tolerance to quinoline drugs

Drug resistance in Plasmodium falciparum remains a challenge for the malaria eradication programmes around the world. With the emergence of artemisinin resistance, the efficacy of the partner drugs in the artemisinin combination therapies (ACT) that include quinoline-based drugs is becoming critical. So far only few resistance markers have been identified from which only two transmembrane transporters namely PfMDR1 (an ATP-binding cassette transporter) and PfCRT (a drug-metabolite transporter) have been experimentally verified. Another P. falciparum transporter, the ATP-binding cassette containing multidrug resistance-associated protein (PfMRP2) represents an additional possible factor of drug resistance in P. falciparum. In this study, we identified a parasite clone that is derived from the 3D7 P. falciparum strain and shows increased resistance to chloroquine, mefloquine and quinine through the trophozoite and schizont stages. We demonstrate that the resistance phenotype is caused by a 4.1 kb deletion in the 5′ upstream region of the pfmrp2 gene that leads to an alteration in the pfmrp2 transcription and thus increased level of PfMRP2 protein. These results also suggest the importance of putative promoter elements in regulation of gene expression during the P. falciparum intra-erythrocytic developmental cycle and the potential of genetic polymorphisms within these regions to underlie drug resistance.

Repeat polymorphisms in the low-complexity regions of Plasmodium falciparum ABC transporters and associations with in vitro antimalarial responses

The Plasmodium falciparum genome is rich in regions of low amino acid complexity which evolve with few constraints on size. To explore the extent of diversity in these loci, we sequenced repeat regions in pfmdr1, pfmdr5, pfmdr6, pfmrp2, and the antigenic locus pfmsp8 in laboratory and cultured-adapted clinical isolates. We further assessed associations between the repeats and parasite in vitro responses to 7 antimalarials to determine possible adaptive roles of these repeats in drug tolerance. Our results show extensive repeat variations in the reference and clinical isolates in all loci. We also observed a modest increase in dihydroartemisinin activity in parasites harboring the pfmdr1 sequence profile 7-2-10 (reflecting the number of asparagine repeats, number of aspartate repeats, and number of asparagine repeats in the final series of the gene product) (P = 0.0321) and reduced sensitivity to chloroquine, mefloquine, quinine, and dihydroartemisinin in those with the 7-2-11 profile (P = 0.0051, 0.0068, 0.0011, and 0.0052, respectively). Interestingly, we noted an inverse association between two drugs whereby isolates with 6 asparagine repeats encoded by pfmdr6 were significantly more susceptible to piperaquine than those with 8 (P = 0.0057). Against lumefantrine, those with 8 repeats were, however, more sensitive (P = 0.0144). In pfmrp2, the 7-DNNNTS/NNNNTS (number of DNNNTS or NNNNTS motifs; underlining indicates dimorphism) repeat group was significantly associated with a higher lumefantrine 50% inhibitory concentration (IC50) (P = 0.008) than in those without. No associations were observed with pfmsp8. These results hint at the probable utility of some repeat conformations as markers of in vitro antimalarial response; hence, biochemical functional studies to ascertain their role in P. falciparum are required.

Malarial parasite diversity in chimpanzees: the value of comparative approaches to ascertain the evolution of Plasmodium falciparum antigens

BACKGROUND

Plasmodium falciparum shares its most recent common ancestor with parasites found in African apes; these species constitute the so-called Laverania clade. In this investigation, the evolutionary history of Plasmodium lineages found in chimpanzees (Pan troglodytes) was explored.

METHODS

Here, the remainders of 74 blood samples collected as part of the chimpanzees' routine health examinations were studied. For all positive samples with parasite lineages belonging to the Laverania clade, the complete mitochondrial genome (mtDNA), the gene encoding dihydrofolate reductase-thymidylate synthase (dhfr-ts), the chloroquine resistance transporter (Pfcrt), the circumsporozoite protein (csp), merozoite surface protein 2 (msp2), and the DBL-1 domain from var2CSA were amplified, cloned, and sequenced. Other Plasmodium species were included in the mtDNA, dhfr-ts, and csp analyses. Phylogenetic and evolutionary genetic analyses were performed, including molecular clock analyses on the mtDNA.

RESULTS/CONCLUSIONS

Nine chimpanzees were malaria positive (12.2%); four of those infections were identified as P. falciparum, two as a Plasmodium reichenowi-like parasite or Plasmodium sp., one as Plasmodium gaboni, and two as Plasmodium malariae. All P. falciparum isolates were resistant to chloroquine indicating that the chimpanzees acquired such infections from humans in recent times. Such findings, however, are not sufficient for implicating chimpanzees as an animal reservoir for P. falciparum.Timing estimates support that the Laverania clade has co-existed with hominids for a long-period of time. The proposed species P. gaboni, Plasmodium billbrayi, and Plasmodium billcollinsi are monophyletic groups supporting that they are indeed different species.An expanded CSP phylogeny is presented, including all the Laverania species and other malarial parasites. Contrasting with other Plasmodium, the Laverania csp exhibits great conservation at the central tandem repeat region. Msp2 and var2CSA, however, show extended recent polymorphism in P. falciparum that likely originated after the P. reichenowi-P. falciparum split. The accumulation of such diversity may indicate adaptation to the human host. These examples support the notion that comparative approaches among P. falciparum and its related species will be of great value in understanding the evolution of proteins that are important in parasite invasion of the human red blood cell, as well as those involved in malaria pathogenesis.

Plasmodium falciparum heat shock protein 110 stabilizes the asparagine repeat-rich parasite proteome during malarial fevers

One-fourth of Plasmodium falciparum proteins have asparagine repeats that increase the propensity for aggregation, especially at elevated temperatures that occur routinely in malaria-infected patients. Here we report that a Plasmodium Asn repeat-containing protein (PFI1155w) formed aggregates in mammalian cells at febrile temperatures, as did a yeast Asn/Gln-rich protein (Sup35). Co-expression of the cytoplasmic P. falciparum heat shock protein 110 (PfHsp110c) prevented aggregation. Human or yeast orthologs were much less effective. All-Asn and all-Gln versions of Sup35 were protected from aggregation by PfHsp110c, suggesting that this chaperone is not limited to handling runs of asparagine. PfHsp110c gene-knockout parasites were not viable and conditional knockdown parasites died slowly in the absence of protein-stabilizing ligand. When exposed to brief heat shock, these knockdowns were unable to prevent aggregation of PFI1155w or Sup35 and died rapidly. We conclude that PfHsp110c protects the parasite from harmful effects of its asparagine repeat-rich proteome during febrile episodes.

Hypervariable antigen genes in malaria have ancient roots

Background

The var genes of the human malaria parasite Plasmodium falciparum are highly polymorphic loci coding for the erythrocyte membrane proteins 1 (PfEMP1), which are responsible for the cytoaherence of P. falciparum infected red blood cells to the human vasculature. Cytoadhesion, coupled with differential expression of var genes, contributes to virulence and allows the parasite to establish chronic infections by evading detection from the host’s immune system. Although studying genetic diversity is a major focus of recent work on the var genes, little is known about the gene family's origin and evolutionary history.

Results

Using a novel hidden Markov model-based approach and var sequences assembled from additional isolates and species, we are able to reveal elements of both the early evolution of the var genes as well as recent diversifying events. We compare sequences of the var gene DBLα domains from divergent isolates of P. falciparum (3D7 and HB3), and a closely-related species, Plasmodium reichenowi. We find that the gene family is equally large in P. reichenowi and P. falciparum -- with a minimum of 51 var genes in the P. reichenowi genome (compared to 61 in 3D7 and a minimum of 48 in HB3). In addition, we are able to define large, continuous blocks of homologous sequence among P. falciparum and P. reichenowi var gene DBLα domains. These results reveal that the contemporary structure of the var gene family was present before the divergence of P. falciparum and P. reichenowi, estimated to be between 2.5 to 6 million years ago. We also reveal that recombination has played an important and traceable role in both the establishment, and the maintenance, of diversity in the sequences.

Conclusions

Despite the remarkable diversity and rapid evolution found in these loci within and among P. falciparum populations, the basic structure of these domains and the gene family is surprisingly old and stable. Revealing a common structure as well as conserved sequence among two species also has implications for developing new primate-parasite models for studying the pathology and immunology of falciparum malaria, and for studying the population genetics of var genes and associated virulence phenotypes.

Plasmodium falciparum variability and immune evasion proceed from antigenicity of consensus sequences from DBL6ε; generalization to all DBL from VAR2CSA

We studied all consensus sequences within the four least ‘variable blocks’ (VB) present in the DBL6ε domain of VAR2CSA, the protein involved in the adhesion of infected red blood cells by Plasmodium falciparum that causes the Pregnancy-Associated Malaria (PAM). Characterising consensus sequences with respect to recognition of antibodies and percentage of responders among pregnant women living in areas where P. falciparum is endemic allows the identification of the most antigenic sequences within each VB. When combining these consensus sequences among four serotypes from VB1 or VB5, the most often recognized ones are expected to induce pan-reactive antibodies recognizing VAR2CSA from all plasmodial strains. These sequences are of main interest in the design of an immunogenic molecule. Using a similar approach than for DBL6ε, we studied the five other DBL and the CIDRpam from VAR2CSA, and again identified VB segments with highly conserved consensus sequences. In addition, we identified consensus sequences in other var genes expressed by non-PAM parasites. This finding paves the way for vaccine design against other pathologies caused by P. falciparum.

An evolutionary perspective on the kinome of malaria parasites

Malaria parasites belong to an ancient lineage that diverged very early from the main branch of eukaryotes. The approximately 90-member plasmodial kinome includes a majority of eukaryotic protein kinases that clearly cluster within the AGC, CMGC, TKL, CaMK and CK1 groups found in yeast, plants and mammals, testifying to the ancient ancestry of these families. However, several hundred millions years of independent evolution, and the specific pressures brought about by first a photosynthetic and then a parasitic lifestyle, led to the emergence of unique features in the plasmodial kinome. These include taxon-restricted kinase families, and unique peculiarities of individual enzymes even when they have homologues in other eukaryotes. Here, we merge essential aspects of all three malaria-related communications that were presented at the Evolution of Protein Phosphorylation meeting, and propose an integrated discussion of the specific features of the parasite's kinome and phosphoproteome.

Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

BACKGROUND

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few amino acids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution.

RESULTS

We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance.

CONCLUSION

We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Comprehensive analysis of tandem amino acid repeats from ten angiosperm genomes

BACKGROUND

The presence of tandem amino acid repeats (AARs) is one of the signatures of eukaryotic proteins. AARs were thought to be frequently involved in bio-molecular interactions. Comprehensive studies that primarily focused on metazoan AARs have suggested that AARs are evolving rapidly and are highly variable among species. However, there is still controversy over causal factors of this inter-species variation. In this work, we attempted to investigate this topic mainly by comparing AARs in orthologous proteins from ten angiosperm genomes.

RESULTS

Angiosperm AAR content is positively correlated with the GC content of the protein coding sequence. However, based on observations from fungal AARs and insect AARs, we argue that the applicability of this kind of correlation is limited by AAR residue composition and species' life history traits. Angiosperm AARs also tend to be fast evolving and structurally disordered, supporting the results of comprehensive analyses of metazoans. The functions of conserved long AARs are summarized. Finally, we propose that the rapid mRNA decay rate, alternative splicing and tissue specificity are regulatory processes that are associated with angiosperm proteins harboring AARs.

CONCLUSIONS

Our investigation suggests that GC content is a predictor of AAR content in the protein coding sequence under certain conditions. Although angiosperm AARs lack conservation and 3D structure, a fraction of the proteins that contain AARs may be functionally important and are under extensive regulation in plant cells.

Protein-based signatures of functional evolution in Plasmodium falciparum

Background

It has been known for over a decade that Plasmodium falciparum proteins are enriched in non-globular domains of unknown function. The potential for these regions of protein sequence to undergo high levels of genetic drift provides a fundamental challenge to attempts to identify the molecular basis of adaptive change in malaria parasites.

Results

Evolutionary comparisons were undertaken using a set of forty P. falciparum metabolic enzyme genes, both within the hominid malaria clade (P. reichenowi) and across the genus (P. chabaudi). All genes contained coding elements highly conserved across the genus, but there were also a large number of regions of weakly or non-aligning coding sequence. These displayed remarkable levels of non-synonymous fixed differences within the hominid malaria clade indicating near complete release from purifying selection (dN/dS ratio at residues non-aligning across genus: 0.64, dN/dS ratio at residues identical across genus: 0.03). Regions of low conservation also possessed high levels of hydrophilicity, a marker of non-globularity. The propensity for such regions to act as potent sources of non-synonymous genetic drift within extant P. falciparum isolates was confirmed at chromosomal regions containing genes known to mediate drug resistance in field isolates, where 150 of 153 amino acid variants were located in poorly conserved regions. In contrast, all 22 amino acid variants associated with drug resistance were restricted to highly conserved regions. Additional mutations associated with laboratory-selected drug resistance, such as those in PfATPase4 selected by spiroindolone, were similarly restricted while mutations in another calcium ATPase (PfSERCA, a gene proposed to mediate artemisinin resistance) that reach significant frequencies in field isolates were located exclusively in poorly conserved regions consistent with genetic drift.

Conclusion

Coding sequences of malaria parasites contain prospectively definable domains subject to neutral or nearly neutral evolution on a scale that appears unrivalled in biology. This distinct evolutionary landscape has potential to confound analytical methods developed for other genera. Against this tide of genetic drift, polymorphisms mediating functional change stand out to such an extent that evolutionary context provides a useful signal for identifying the molecular basis of drug resistance in malaria parasites, a finding that is of relevance to both genome-wide and candidate gene studies in this genus.

Rapid whole genome optical mapping of Plasmodium falciparum

Background

Immune evasion and drug resistance in malaria have been linked to chromosomal recombination and gene copy number variation (CNV). These events are ideally studied using comparative genomic analyses; however in malaria these analyses are not as common or thorough as in other infectious diseases, partly due to the difficulty in sequencing and assembling complete genome drafts. Recently, whole genome optical mapping has gained wide use in support of genomic sequence assembly and comparison. Here, a rapid technique for producing whole genome optical maps of Plasmodium falciparum is described and the results of mapping four genomes are presented.

Methods

Four laboratory strains of P. falciparum were analysed using the Argus™ optical mapping system to produce ordered restriction fragment maps of all 14 chromosomes in each genome. Plasmodium falciparum DNA was isolated directly from blood culture, visualized using the Argus™ system and assembled in a manner analogous to next generation sequence assembly into maps (AssemblyViewer™, OpGen Inc.^®). Full coverage maps were generated for P. falciparum strains 3D7, FVO, D6 and C235. A reference P. falciparum in silico map was created by the digestion of the genomic sequence of P. falciparum with the restriction enzyme AflII, for comparisons to genomic optical maps. Maps were then compared using the MapSolver™ software.

Results

Genomic variation was observed among the mapped strains, as well as between the map of the reference strain and the map derived from the putative sequence of that same strain. Duplications, deletions, insertions, inversions and misassemblies of sizes ranging from 3,500 base pairs up to 78,000 base pairs were observed. Many genomic events occurred in areas of known repetitive sequence or high copy number genes, including var gene clusters and rifin complexes.

Conclusions

This technique for optical mapping of multiple malaria genomes allows for whole genome comparison of multiple strains and can assist in identifying genetic variation and sequence contig assembly. New protocols and technology allowed us to produce high quality contigs spanning four P. falciparum genomes in six weeks for less than $1,000.00 per genome. This relatively low cost and quick turnaround makes the technique valuable compared to other genomic sequencing technologies for studying genetic variation in malaria.

Genome cartography: charting the apicomplexan genome

Genes reside in particular genomic contexts that can be mapped at many levels. Historically, 'genetic maps' were used primarily to locate genes. Recent technological advances in the determination of genome sequences have made the analysis and comparison of whole genomes possible and increasingly tractable. What do we see if we shift our focus from gene content (the 'inventory' of genes contained within a genome) to the composition and organization of a genome? This review examines what has been learned about the evolution of the apicomplexan genome as well as the significance and impact of genomic location on our understanding of the eukaryotic genome and parasite biology.

Increased polymorphism near low-complexity sequences across the genomes of Plasmodium falciparum isolates

Low-complexity regions (LCRs) within proteins sequences are often considered to evolve neutrally even though recent studies reported evidence for selection acting on some of them. Because of their widespread distribution among eukaryotes genomes and the potential deleterious effect of expansion/contraction of some of them in humans, low-complexity sequences are of major interest and numerous studies have attempted to describe their dynamic between genomes as well as the factors correlated to their variation and to assess their selective value. However, due to the scarcity of individual genomes within a species, most of the analyses so far have been performed at the species level with the implicit assumption that the variation both in composition and size within species is too small relative to the between-species divergence to affect the conclusions of the analysis. Here we used the available genomes of 14 Plasmodium falciparum isolates to assess the relationship between low-complexity sequence variation and factors such as nucleotide polymorphism across strains, sequence composition, and protein expression. We report that more than half of the 7,711 low-complexity sequences found within aligned coding sequences are variable in size among strains. Across strains, we observed an increasing density of polymorphic sites toward the LCR boundaries. This observation strongly suggests the joint effects of lowered selective constraints on low-complexity sequences and a mutagenic effect of these simple sequences.

High recombination rates and hotspots in a Plasmodium falciparum genetic cross

Background

The human malaria parasite Plasmodium falciparum survives pressures from the host immune system and antimalarial drugs by modifying its genome. Genetic recombination and nucleotide substitution are the two major mechanisms that the parasite employs to generate genome diversity. A better understanding of these mechanisms may provide important information for studying parasite evolution, immune evasion and drug resistance.

Results

Here, we used a high-density tiling array to estimate the genetic recombination rate among 32 progeny of a P. falciparum genetic cross (7G8 × GB4). We detected 638 recombination events and constructed a high-resolution genetic map. Comparing genetic and physical maps, we obtained an overall recombination rate of 9.6 kb per centimorgan and identified 54 candidate recombination hotspots. Similar to centromeres in other organisms, the sequences of P. falciparum centromeres are found in chromosome regions largely devoid of recombination activity. Motifs enriched in hotspots were also identified, including a 12-bp G/C-rich motif with 3-bp periodicity that may interact with a protein containing 11 predicted zinc finger arrays.

Conclusions

These results show that the P. falciparum genome has a high recombination rate, although it also follows the overall rule of meiosis in eukaryotes with an average of approximately one crossover per chromosome per meiosis. GC-rich repetitive motifs identified in the hotspot sequences may play a role in the high recombination rate observed. The lack of recombination activity in centromeric regions is consistent with the observations of reduced recombination near the centromeres of other organisms.

Asparagine repeat function in a Plasmodium falciparum protein assessed via a regulatable fluorescent affinity tag

One in four proteins in Plasmodium falciparum contains asparagine repeats. We probed the function of one such 28-residue asparagine repeat present in the P. falciparum proteasome lid subunit 6, Rpn6. To aid our efforts, we developed a regulatable, fluorescent affinity (RFA) tag that allows cellular localization, manipulation of cellular levels, and affinity isolation of a chosen protein in P. falciparum. The tag comprises a degradation domain derived from Escherichia coli dihydrofolate reductase together with GFP. The expression of RFA-tagged proteins is regulated by the simple folate analog trimethoprim (TMP). Parasite lines were generated in which full-length Rpn6 and an asparagine repeat-deletion mutant of Rpn6 were fused to the RFA tag. The knockdown of Rpn6 upon removal of TMP revealed that this protein is essential for ubiquitinated protein degradation and for parasite survival, but the asparagine repeat is dispensable for protein expression, stability, and function. The data point to a genomic mechanism for repeat perpetuation rather than a positive cellular role. The RFA tag should facilitate study of the role of essential genes in parasite biology.