Antigenic variation in Plasmodium falciparum: gene organization and regulation of the var multigene family

A novel computational pipeline for var gene expression augments the discovery of changes in the Plasmodium falciparum transcriptome during transition from in vivo to short-term in vitro culture

The pathogenesis of severe Plasmodium falciparum malaria involves cytoadhesive microvascular sequestration of infected erythrocytes, mediated by P. falciparum erythrocyte membrane protein 1 (PfEMP1). PfEMP1 variants are encoded by the highly polymorphic family of var genes, the sequences of which are largely unknown in clinical samples. Previously, we published new approaches for var gene profiling and classification of predicted binding phenotypes in clinical P. falciparum isolates (Wichers et al., 2021), which represented a major technical advance. Building on this, we report here a novel method for var gene assembly and multidimensional quantification from RNA-sequencing that outperforms the earlier approach of Wichers et al., 2021, on both laboratory and clinical isolates across a combination of metrics. Importantly, the tool can interrogate the var transcriptome in context with the rest of the transcriptome and can be applied to enhance our understanding of the role of var genes in malaria pathogenesis. We applied this new method to investigate changes in var gene expression through early transition of parasite isolates to in vitro culture, using paired sets of ex vivo samples from our previous study, cultured for up to three generations. In parallel, changes in non-polymorphic core gene expression were investigated. Modest but unpredictable var gene switching and convergence towards var2csa were observed in culture, along with differential expression of 19% of the core transcriptome between paired ex vivo and generation 1 samples. Our results cast doubt on the validity of the common practice of using short-term cultured parasites to make inferences about in vivo phenotype and behaviour.

Effects of the G-quadruplex-binding drugs quarfloxin and CX-5461 on the malaria parasite Plasmodium falciparum

Plasmodium falciparum is the deadliest causative agent of human malaria. This parasite has historically developed resistance to most drugs, including the current frontline treatments, so new therapeutic targets are needed. Our previous work on guanine quadruplexes (G4s) in the parasite's DNA and RNA has highlighted their influence on parasite biology, and revealed G4 stabilising compounds as promising candidates for repositioning. In particular, quarfloxin, a former anticancer agent, kills blood-stage parasites at all developmental stages, with fast rates of kill and nanomolar potency. Here we explored the molecular mechanism of quarfloxin and its related derivative CX-5461. In vitro, both compounds bound to P. falciparum-encoded G4 sequences. In cellulo, quarfloxin was more potent than CX-5461, and could prevent establishment of blood-stage malaria in vivo in a murine model. CX-5461 showed clear DNA damaging activity, as reported in human cells, while quarfloxin caused weaker signatures of DNA damage. Both compounds caused transcriptional dysregulation in the parasite, but the affected genes were largely different, again suggesting different modes of action. Therefore, CX-5461 may act primarily as a DNA damaging agent in both Plasmodium parasites and mammalian cells, whereas the complete antimalarial mode of action of quarfloxin may be parasite-specific and remains somewhat elusive.

A paradoxical population structure of var DBLα types in Africa

The var multigene family encodes the P. falciparum erythrocyte membrane protein 1 (PfEMP1), which is important in host-parasite interaction as a virulence factor and major surface antigen of the blood stages of the parasite, responsible for maintaining chronic infection. Whilst important in the biology of P. falciparum, these genes (50 to 60 genes per parasite genome) are routinely excluded from whole genome analyses due to their hyper-diversity, achieved primarily through recombination. The PfEMP1 head structure almost always consists of a DBLα-CIDR tandem. Categorised into different groups (upsA, upsB, upsC), different head structures have been associated with different ligand-binding affinities and disease severities. We study how conserved individual DBLα types are at the country, regional, and local scales in Sub-Saharan Africa. Using publicly-available sequence datasets and a novel ups classification algorithm, cUps, we performed an in silico exploration of DBLα conservation through time and space in Africa. In all three ups groups, the population structure of DBLα types in Africa consists of variants occurring at rare, low, moderate, and high frequencies. Non-rare variants were found to be temporally stable in a local area in endemic Ghana. When inspected across different geographical scales, we report different levels of conservation; while some DBLα types were consistently found in high frequencies in multiple African countries, others were conserved only locally, signifying local preservation of specific types. Underlying this population pattern is the composition of DBLα types within each isolate DBLα repertoire, revealed to also consist of a mix of types found at rare, low, moderate, and high frequencies in the population. We further discuss the adaptive forces and balancing selection, including host genetic factors, potentially shaping the evolution and diversity of DBLα types in Africa.

Characterization of Aspartic Proteases from Paracoccidioides brasiliensis and Their Role in Fungal Thermo-Dimorphism

Paracoccidioidomycosis (PCM) is the most prevalent systemic mycosis in Latin America and is caused by fungi from the Paracoccidioides genus. The infection begins after inhalation of the fungal propagules and their thermo-dimorphic shift to yeast form. Proteases play an important role in the host invasion process and immune modulation in many pathogenic microorganisms. Aspartyl proteases are virulence factors in many human fungal pathogens that play an important role in the host invasion process morphogenesis, cellular function, immunity, and nutrition. In the present study, we characterized the modulation of acid proteases from Paracoccidioides brasiliensis. We detected four aspartyl proteases in P. brasiliensis with high homology to aspartic protease from Saccharomyces cerevisiae Pep4. Furthermore, we demonstrated that Pepstatin A can inhibit dimorphic switching (mycelium→yeast) in P. brasiliensis. In addition, these genes were modulated during thermo-dimorphism (M→Y transition) in the presence or absence of carbon and nitrogen sources and during growth at pH 4 during 24 and 48 h. We also observed that P. brasiliensis increase the secretion of aspartic proteases when cultivated at pH 4, and these acid proteases cleave BSA, collagen, and hemoglobin. These data suggest that aspartyl proteases are modulated by environmental conditions and during fungal thermo-dimorphism. Thus, this work brings new possibilities for studying the role of aspartyl proteases in the host-pathogen relationship and P. brasiliensis biology.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

CD36-A Host Receptor Necessary for Malaria Parasites to Establish and Maintain Infection

Plasmodium falciparum-infected erythrocytes (PfIEs) present P. falciparum erythrocyte membrane protein 1 proteins (PfEMP1s) on the cell surface, via which they cytoadhere to various endothelial cell receptors (ECRs) on the walls of human blood vessels. This prevents the parasite from passing through the spleen, which would lead to its elimination. Each P. falciparum isolate has about 60 different PfEMP1s acting as ligands, and at least 24 ECRs have been identified as interaction partners. Interestingly, in every parasite genome sequenced to date, at least 75% of the encoded PfEMP1s have a binding domain for the scavenger receptor CD36 widely distributed on host endothelial cells and many other cell types. Here, we discuss why the interaction between PfIEs and CD36 is optimal to maintain a finely regulated equilibrium that allows the parasite to multiply and spread while causing minimal harm to the host in most infections.

Innate immunity to malaria: The good, the bad and the unknown

Malaria is the cause of 600.000 deaths annually. However, these deaths represent only a tiny fraction of total malaria cases. Repeated natural infections with the causative agent, Plasmodium sp. parasites, induce protection from severe disease but not sterile immunity. Thus, immunity to Plasmodium is incomplete. Conversely, immunization with attenuated sporozoite stage parasites can induce sterile immunity albeit after multiple vaccinations. These different outcomes are likely to be influenced strongly by the innate immune response to different stages of the parasite lifecycle. Even small numbers of sporozoites can induce a robust proinflammatory type I interferon response, which is believed to be driven by the sensing of parasite RNA. Moreover, induction of innate like gamma-delta cells contributes to the development of adaptive immune responses. Conversely, while blood stage parasites can induce a strong proinflammatory response, regulatory mechanisms are also triggered. In agreement with this, intact parasites are relatively weakly sensed by innate immune cells, but isolated parasite molecules, notably DNA and RNA can induce strong responses. Thus, the innate response to Plasmodium parasite likely represents a trade-off between strong pro-inflammatory responses that may potentiate immunity and regulatory processes that protect the host from cytokine storms that can induce life threatening illness.

PfARID Regulates P. falciparum Malaria Parasite Male Gametogenesis and Female Fertility and Is Critical for Parasite Transmission to the Mosquito Vector

Sexual reproduction of Plasmodium falciparum parasites is critical to the spread of malaria in the human population. The factors that regulate gene expression underlying formation of fertilization-competent gametes, however, remain unknown. Here, we report that P. falciparum expresses a protein with an AT-rich interaction domain (ARID) which, in other organisms, is part of chromatin remodeling complexes. P. falciparum ARID (PfARID) localized to the parasite nucleus and is critical for the formation of male gametes and fertility of female gametes. PfARID gene deletion (Pfarid^–) gametocytes showed downregulation of gene expression important for gametogenesis, antigenic variation, and cell signaling and for parasite development in the mosquito. Our study identifies PfARID as a critical nuclear protein involved in regulating the gene expression landscape of mature gametocytes. This establishes fertility and also prepares the parasite for postfertilization events that are essential for infection of the mosquito vector.

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).

"Immunoinformatic Identification of T-Cell and B-Cell Epitopes From Giardia lamblia Immunogenic Proteins as Candidates to Develop Peptide-Based Vaccines Against Giardiasis"

Giardiasis is one of the most common gastrointestinal infections worldwide, mainly in developing countries. The etiological agent is the Giardia lamblia parasite. Giardiasis mainly affects children and immunocompromised people, causing symptoms such as diarrhea, dehydration, abdominal cramps, nausea, and malnutrition. In order to develop an effective vaccine against giardiasis, it is necessary to understand the host-Giardia interactions, the immunological mechanisms involved in protection against infection, and to characterize the parasite antigens that activate the host immune system. In this study, we identify and characterize potential T-cell and B-cell epitopes of Giardia immunogenic proteins by immunoinformatic approaches, and we discuss the potential role of those epitopes to stimulate the host´s immune system. We selected the main immunogenic and protective proteins of Giardia experimentally investigated. We predicted T-cell and B-cell epitopes using immunoinformatic tools (NetMHCII and BCPREDS). Variable surface proteins (VSPs), structural (giardins), metabolic, and cyst wall proteins were identified as the more relevant immunogens of G. lamblia. We described the protein sequences with the highest affinity to bind MHC class II molecules from mouse (I-A^k and I-A^d) and human (DRB1*03:01 and DRB1*13:01) alleles, as well as we selected promiscuous epitopes, which bind to the most common range of MHC class II molecules in human population. In addition, we identified the presence of conserved epitopes within the main protein families (giardins, VSP, CWP) of Giardia. To our knowledge, this is the first in silico study that analyze immunogenic proteins of G. lamblia by combining bioinformatics strategies to identify potential T-cell and B-cell epitopes, which can be potential candidates in the development of peptide-based vaccines. The bioinformatics analysis demonstrated in this study provides a deeper understanding of the Giardia immunogens that bind to critical molecules of the host immune system, such as MHC class II and antibodies, as well as strategies to rational design of peptide-based vaccine against giardiasis.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

The impact of indoor residual spraying on Plasmodium falciparum microsatellite variation in an area of high seasonal malaria transmission in Ghana, West Africa

Here, we report the first population genetic study to examine the impact of indoor residual spraying (IRS) on Plasmodium falciparum in humans. This study was conducted in an area of high seasonal malaria transmission in Bongo District, Ghana. IRS was implemented during the dry season (November-May) in three consecutive years between 2013 and 2015 to reduce transmission and attempt to bottleneck the parasite population in humans towards lower diversity with greater linkage disequilibrium. The study was done against a background of widespread use of long-lasting insecticidal nets, typical for contemporary malaria control in West Africa. Microsatellite genotyping with 10 loci was used to construct 392 P. falciparum multilocus infection haplotypes collected from two age-stratified cross-sectional surveys at the end of the wet seasons pre- and post-IRS. Three-rounds of IRS, under operational conditions, led to a >90% reduction in transmission intensity and a 35.7% reduction in the P. falciparum prevalence (p < .001). Despite these declines, population genetic analysis of the infection haplotypes revealed no dramatic changes with only a slight, but significant increase in genetic diversity (He : pre-IRS = 0.79 vs. post-IRS = 0.81, p = .048). Reduced relatedness of the parasite population (p < .001) was observed post-IRS, probably due to decreased opportunities for outcrossing. Spatiotemporal genetic differentiation between the pre- and post-IRS surveys (D = 0.0329 [95% CI: 0.0209 - 0.0473], p = .034) was identified. These data provide a genetic explanation for the resilience of P. falciparum to short-term IRS programmes in high-transmission settings in sub-Saharan Africa.

Mechanistic within-host models of the asexual Plasmodium falciparum infection: a review and analytical assessment

Background

Malaria blood-stage infection length and intensity are important drivers of disease and transmission; however, the underlying mechanisms of parasite growth and the host’s immune response during infection remain largely unknown. Over the last 30 years, several mechanistic mathematical models of malaria parasite within-host dynamics have been published and used in malaria transmission models.

Methods

Mechanistic within-host models of parasite dynamics were identified through a review of published literature. For a subset of these, model code was reproduced and descriptive statistics compared between the models using fitted data. Through simulation and model analysis, key features of the models were compared, including assumptions on growth, immune response components, variant switching mechanisms, and inter-individual variability.

Results

The assessed within-host malaria models generally replicate infection dynamics in malaria-naïve individuals. However, there are substantial differences between the model dynamics after disease onset, and models do not always reproduce late infection parasitaemia data used for calibration of the within host infections. Models have attempted to capture the considerable variability in parasite dynamics between individuals by including stochastic parasite multiplication rates; variant switching dynamics leading to immune escape; variable effects of the host immune responses; or via probabilistic events. For models that capture realistic length of infections, model representations of innate immunity explain early peaks in infection density that cause clinical symptoms, and model representations of antibody immune responses control the length of infection. Models differed in their assumptions concerning variant switching dynamics, reflecting uncertainty in the underlying mechanisms of variant switching revealed by recent clinical data during early infection. Overall, given the scarce availability of the biological evidence there is limited support for complex models.

Conclusions

This study suggests that much of the inter-individual variability observed in clinical malaria infections has traditionally been attributed in models to random variability, rather than mechanistic disease dynamics. Thus, it is proposed that newly developed models should assume simple immune dynamics that minimally capture mechanistic understandings and avoid over-parameterization and large stochasticity which inaccurately represent unknown disease mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03813-z.

The mechanisms of action of Plasmodium infection against cancer

Our murine cancer model studies have demonstrated that Plasmodium infection activates the immune system that has been inhibited by cancer cells, counteracts tumor immunosuppressive microenvironment, inhibits tumor angiogenesis, inhibits tumor growth and metastasis, and prolongs the survival time of tumor-bearing mice. Based on these studies, three clinical trials of Plasmodium immunotherapy for advanced cancers have been approved and are ongoing in China. After comparing the mechanisms of action of Plasmodium immunotherapy with those of immune checkpoint blockade therapy, we propose the notion that cancer is an ecological disease and that Plasmodium immunotherapy is a systemic ecological counterattack therapy for this ecological disease, with limited side effects and without danger to public health based on the use of artesunate and other measures. Recent reports of tolerance to treatment and limitations in majority of patients associated with the use of checkpoint blockers further support this notion. We advocate further studies on the mechanisms of action of Plasmodium infection against cancer and investigations on Plasmodium-based combination therapy in the coming future.

Common virulence gene expression in adult first-time infected malaria patients and severe cases

Sequestration of Plasmodium falciparum(P. falciparum)-infected erythrocytes to host endothelium through the parasite-derived P. falciparum erythrocyte membrane protein 1 (PfEMP1) adhesion proteins is central to the development of malaria pathogenesis. PfEMP1 proteins have diversified and expanded to encompass many sequence variants, conferring each parasite a similar array of human endothelial receptor-binding phenotypes. Here, we analyzed RNA-seq profiles of parasites isolated from 32 P. falciparum-infected adult travellers returning to Germany. Patients were categorized into either malaria naive (n = 15) or pre-exposed (n = 17), and into severe (n = 8) or non-severe (n = 24) cases. For differential expression analysis, PfEMP1-encoding var gene transcripts were de novo assembled from RNA-seq data and, in parallel, var-expressed sequence tags were analyzed and used to predict the encoded domain composition of the transcripts. Both approaches showed in concordance that severe malaria was associated with PfEMP1 containing the endothelial protein C receptor (EPCR)-binding CIDRα1 domain, whereas CD36-binding PfEMP1 was linked to non-severe malaria outcomes. First-time infected adults were more likely to develop severe symptoms and tended to be infected for a longer period. Thus, parasites with more pathogenic PfEMP1 variants are more common in patients with a naive immune status, and/or adverse inflammatory host responses to first infections favor the growth of EPCR-binding parasites.

DNA Double-Strand Breaks: A Double-Edged Sword for Trypanosomatids

For nearly all eukaryotic cells, stochastic DNA double-strand breaks (DSBs) are one of the most deleterious types of DNA lesions. DSB processing and repair can cause sequence deletions, loss of heterozygosity, and chromosome rearrangements resulting in cell death or carcinogenesis. However, trypanosomatids (single-celled eukaryotes parasites) do not seem to follow this premise strictly. Several studies have shown that trypanosomatids depend on DSBs to perform several events of paramount importance during their life cycle. For Trypanosoma brucei, DSBs formation is associated with host immune evasion via antigenic variation. In Trypanosoma cruzi, DSBs play a crucial role in the genetic exchange, a mechanism that is still little explored but appear to be of fundamental importance for generating variability. In Leishmania spp., DSBs are necessary to generate genomic changes by gene copy number variation (CNVs), events that are essential for these organisms to overcome inhospitable conditions. As DSB repair in trypanosomatids is primarily conducted via homologous recombination (HR), most of the events associated with DSBs are HR-dependent. This review will discuss the latest findings on how trypanosomatids balance the benefits and inexorable challenges caused by DSBs.

Setting the stage: The initial immune response to blood-stage parasites

Individuals growing up in malaria endemic areas gradually develop protection against clinical malaria and passive transfer experiments in humans have demonstrated that this protection is mediated in part by protective antibodies. However, neither the target antigens, specific effector mechanisms, nor the role of continual parasite exposure have been elucidated, which complicates vaccine development. Progress has been made in defining the innate signaling pathways activated by parasite components, including DNA, RNA, hemozoin, and phospholipids, which initiate the immune response and will be the focus of this review. The challenge that remains within the field is to understand the role of these early responses in the development of protective adaptive responses that clear iRBC and block merozoite invasion so that optimal vaccines and therapeutics may be produced.

The Impact of Malaria Parasites on Dendritic Cell-T Cell Interaction

Malaria is caused by apicomplexan parasites of the genus Plasmodium. While infection continues to pose a risk for the majority of the global population, the burden of disease mainly resides in Sub-Saharan Africa. Although immunity develops against disease, this requires years of persistent exposure and is not associated with protection against infection. Repeat infections occur due to the parasite's ability to disrupt or evade the host immune responses. However, despite many years of study, the mechanisms of this disruption remain unclear. Previous studies have demonstrated a parasite-induced failure in dendritic cell (DCs) function affecting the generation of helper T cell responses. These T cells fail to help B cell responses, reducing the production of antibodies that are necessary to control malaria infection. This review focuses on our current understanding of the effect of Plasmodium parasite on DC function, DC-T cell interaction, and T cell activation. A better understanding of how parasites disrupt DC-T cell interactions will lead to new targets and approaches to reinstate adaptive immune responses and enhance parasite immunity.

Identification of a conserved var gene in different Plasmodium falciparum strains

BACKGROUND

The multicopy var gene family of Plasmodium falciparum is of crucial importance for pathogenesis and antigenic variation. So far only var2csa, the var gene responsible for placental malaria, was found to be highly conserved among all P. falciparum strains. Here, a new conserved 3D7 var gene (PF3D7_0617400) is identified in several field isolates.

METHODS

DNA sequencing, transcriptional analysis, Cluster of Differentiation (CD) 36-receptor binding, indirect immunofluorescence with PF3D7_0617400-antibodies and quantification of surface reactivity against semi-immune sera were used to characterize an NF54 clone and a Gabonese field isolate clone (MOA C3) transcribing the gene. A population of 714 whole genome sequenced parasites was analysed to characterize the conservation of the locus in African and Asian isolates. The genetic diversity of two var2csa fragments was compared with the genetic diversity of 57 microsatellites fragments in field isolates.

RESULTS

PFGA01_060022400 was identified in a Gabonese parasite isolate (MOA) from a chronic infection and found to be 99% identical with PF3D7_0617400 of the 3D7 genome strain. Transcriptional analysis and immunofluorescence showed expression of the gene in an NF54 and a MOA clone but CD36 binding assays and surface reactivity to semi-immune sera differed markedly in the two clones. Long-read Pacific bioscience whole genome sequencing showed that PFGA01_060022400 is located in the internal cluster of chromosome 6. The full length PFGA01_060022400 was detected in 36 of 714 P. falciparum isolates and 500 bp fragments were identified in more than 100 isolates. var2csa was in parts highly conserved (H_e = 0) but in other parts as variable (H_e = 0.86) as the 57 microsatellites markers (H_e = 0.8).

CONCLUSIONS

Individual var gene sequences exhibit conservation in the global parasite population suggesting that purifying selection may limit overall genetic diversity of some var genes. Notably, field and laboratory isolates expressing the same var gene exhibit markedly different phenotypes.

VAPPER: High-throughput variant antigen profiling in African trypanosomes of livestock

BACKGROUND

Analysing variant antigen gene families on a population scale is a difficult challenge for conventional methods of read mapping and variant calling due to the great variability in sequence, copy number, and genomic loci. In African trypanosomes, hemoparasites of humans and animals, this is complicated by variant antigen repertoires containing hundreds of genes subject to various degrees of sequence recombination.

FINDINGS

We introduce Variant Antigen Profiler (VAPPER), a tool that allows automated analysis of the variant surface glycoprotein repertoires of the most prevalent livestock African trypanosomes. VAPPER produces variant antigen profiles for any isolate of the veterinary pathogens Trypanosoma congolense and Trypanosoma vivax from genomic and transcriptomic sequencing data and delivers publication-ready figures that show how the queried isolate compares with a database of existing strains. VAPPER is implemented in Python. It can be installed to a local Galaxy instance from the ToolShed (https://toolshed.g2.bx.psu.edu/) or locally on a Linux platform via the command line (https://github.com/PGB-LIV/VAPPER). The documentation, requirements, examples, and test data are provided in the Github repository.

CONCLUSION

By establishing two different, yet comparable methodologies, our approach is the first to allow large-scale analysis of African trypanosome variant antigens, large multi-copy gene families that are otherwise refractory to high-throughput analysis.

Immunity against sexual stage Plasmodium falciparum and Plasmodium vivax parasites

The efficient spread of malaria from infected humans to mosquitoes is a major challenge for malaria elimination initiatives. Gametocytes are the only Plasmodium life stage infectious to mosquitoes. Here, we summarize evidence for naturally acquired anti-gametocyte immunity and the current state of transmission blocking vaccines (TBV). Although gametocytes are intra-erythrocytic when present in infected humans, developing Plasmodium falciparum gametocytes may express proteins on the surface of red blood cells that elicit immune responses in naturally exposed individuals. This immune response may reduce the burden of circulating gametocytes. For both P. falciparum and Plasmodium vivax, there is a solid evidence that antibodies against antigens present on the gametocyte surface, when co-ingested with gametocytes, can influence transmission to mosquitoes. Transmission reducing immunity, reducing the burden of infection in mosquitoes, is a well-acknowledged but poorly quantified phenomenon that forms the basis for the development of TBV. Transmission enhancing immunity, increasing the likelihood or intensity of transmission to mosquitoes, is more speculative in nature but is convincingly demonstrated for P. vivax. With the increased interest in malaria elimination, TBV and monoclonal antibodies have moved to the center stage of malaria vaccine development. Methodologies to prioritize and evaluate products are urgently needed.

A unique aspartyl protease gene expansion in Talaromyces marneffei plays a role in growth inside host phagocytes

Aspartyl proteases are a widely represented class of proteolytic enzymes found in eukaryotes and retroviruses. They have been associated with pathogenicity in a range of disease-causing microorganisms. The dimorphic human-pathogenic fungus Talaromyces marneffei has a large expansion of these proteases identified through genomic analyses. Here we characterize the expansion of these genes (pop - paralogue of pep) and their role in T. marneffei using computational and molecular approaches. Many of the genes in this monophyletic family show copy number variation and positive selection despite the preservation of functional regions and possible redundancy. We show that the expression profile of these genes differs and some are expressed during intracellular growth in the host. Several of these proteins have distinctive localization as well as both additive and epistatic effects on the formation of yeast cells during macrophage infections. The data suggest that this is a recently evolved aspartyl protease gene family which affects intracellular growth and contributes to the pathogenicity of T. marneffei.

Is reliance on an inaccurate genome sequence sabotaging your experiments?

Advances in genomics have made whole genome studies increasingly feasible across the life sciences. However, new technologies and algorithmic advances do not guarantee flawless genomic sequences or annotation. Bias, errors, and artifacts can enter at any stage of the process from library preparation to annotation. When planning an experiment that utilizes a genome sequence as the basis for the design, there are a few basic checks that, if performed, may better inform the experimental design and ideally help avoid a failed experiment or inconclusive result.

Controlled human malaria infection with Plasmodium falciparum demonstrates impact of naturally acquired immunity on virulence gene expression

The pathogenesis of Plasmodium falciparum malaria is linked to the variant surface antigen PfEMP1, which mediates tethering of infected erythrocytes to the host endothelium and is encoded by approximately 60 var genes per parasite genome. Repeated episodes of malaria infection result in the gradual acquisition of protective antibodies against PfEMP1 variants. The antibody repertoire is believed to provide a selective pressure driving the clonal expansion of parasites expressing unrecognized PfEMP1 variants, however, due to the lack of experimental in vivo models there is only limited experimental evidence in support of this concept. To get insight into the impact of naturally acquired immunity on the expressed var gene repertoire early during infection we performed controlled human malaria infections of 20 adult African volunteers with life-long malaria exposure using aseptic, purified, cryopreserved P. falciparum sporozoites (Sanaria PfSPZ Challenge) and correlated serological data with var gene expression patterns from ex vivo parasites. Among the 10 African volunteers who developed patent infections, individuals with low antibody levels showed a steep rise in parasitemia accompanied by broad activation of multiple, predominantly subtelomeric var genes, similar to what we previously observed in naïve volunteers. In contrast, individuals with intermediate antibody levels developed asymptomatic infections and the ex vivo parasite populations expressed only few var gene variants, indicative of clonal selection. Importantly, in contrast to parasites from naïve volunteers, expression of var genes coding for endothelial protein C receptor (EPCR)-binding PfEMP1 that are associated with severe childhood malaria was rarely detected in semi-immune adult African volunteers. Moreover, we followed var gene expression for up to six parasite replication cycles and demonstrated for the first time in vivo a shift in the dominant var gene variant. In conclusion, our data suggest that P. falciparum activates multiple subtelomeric var genes at the onset of blood stage infection facilitating rapid expansion of parasite clones which express PfEMP1 variants unrecognized by the host's immune system, thus promoting overall parasite survival in the face of host immunity.

Within-host modeling of blood-stage malaria

Malaria infection continues to be a major health problem worldwide and drug resistance in the major human parasite species, Plasmodium falciparum, is increasing in South East Asia. Control measures including novel drugs and vaccines are in development, and contributions to the rational design and optimal usage of these interventions are urgently needed. Infection involves the complex interaction of parasite dynamics, host immunity, and drug effects. The long life cycle (48 hours in the common human species) and synchronized replication cycle of the parasite population present significant challenges to modeling the dynamics of Plasmodium infection. Coupled with these, variation in immune recognition and drug action at different life cycle stages leads to further complexity. We review the development and progress of "within-host" models of Plasmodium infection, and how these have been applied to understanding and interpreting human infection and animal models of infection.

Epigenetic Players of Chromatin Structure Regulation in Plasmodium falciparum

The protozoan parasite Plasmodium has evolved to survive in different hosts and environments. The diverse strategies of adaptation to different niches involve differential gene expression mechanisms mediated by chromatin plasticity that are poorly characterized in Plasmodium. The parasite employs a wide variety of regulatory mechanisms to complete their life cycle and survive inside hosts. Among them, epigenetic-mediated mechanisms have been implicated for controlling chromatin organization, gene regulation, morphological differentiation, and antigenic variation. The differential gene expression in parasite is largely dependent on the nature of the chromatin structure. The histone core methylation marks and methyl mark readers contribute to chromatin dynamics. Here, we review the recent developments on various epigenetic marks and its enzymes in the Plasmodium falciparum, how these marks play a key role in the regulation of transcriptional activity of variable genes and coordinate the differential gene expression. We also discuss the possible roles of these epigenetic marks in chromatin structure regulation and plasticity at various stages of its development.

Evidence of the Red-Queen Hypothesis from Accelerated Rates of Evolution of Genes Involved in Biotic Interactions in Pneumocystis

Pneumocystis species are ascomycete fungi adapted to live inside the lungs of mammals. These ascomycetes show extensive stenoxenism, meaning that each species of Pneumocystis infects a single species of host. Here, we study the effect exerted by natural selection on gene evolution in the genomes of three Pneumocystis species. We show that genes involved in host interaction evolve under positive selection. In the first place, we found strong evidence of episodic diversifying selection in Major surface glycoproteins (Msg). These proteins are located on the surface of Pneumocystis and are used for host attachment and probably for immune system evasion. Consistent with their function as antigens, most sites under diversifying selection in Msg code for residues with large relative surface accessibility areas. We also found evidence of positive selection in part of the cell machinery used to export Msg to the cell surface. Specifically, we found that genes participating in glycosylphosphatidylinositol (GPI) biosynthesis show an increased rate of nonsynonymous substitutions (dN) versus synonymous substitutions (dS). GPI is a molecule synthesized in the endoplasmic reticulum that is used to anchor proteins to membranes. We interpret the aforementioned findings as evidence of selective pressure exerted by the host immune system on Pneumocystis species, shaping the evolution of Msg and several proteins involved in GPI biosynthesis. We suggest that genome evolution in Pneumocystis is well described by the Red-Queen hypothesis whereby genes relevant for biotic interactions show accelerated rates of evolution.

Identifying functional groups among the diverse, recombining antigenic var genes of the malaria parasite Plasmodium falciparum from a local community in Ghana

A challenge in studying diverse multi-copy gene families is deciphering distinct functional types within immense sequence variation. Functional changes can in some cases be tracked through the evolutionary history of a gene family; however phylogenetic approaches are not possible in cases where gene families diversify primarily by recombination. We take a network theoretical approach to functionally classify the highly recombining var antigenic gene family of the malaria parasite Plasmodium falciparum. We sample var DBLα sequence types from a local population in Ghana, and classify 9,276 of these variants into just 48 functional types. Our approach is to first decompose each sequence type into its constituent, recombining parts; we then use a stochastic block model to identify functional groups among the parts; finally, we classify the sequence types based on which functional groups they contain. This method for functional classification does not rely on an inferred phylogenetic history, nor does it rely on inferring function based on conserved sequence features. Instead, it infers functional similarity among recombining parts based on the sharing of similar co-occurrence interactions with other parts. This method can therefore group sequences that have undetectable sequence homology or even distinct origination. Describing these 48 var functional types allows us to simplify the antigenic diversity within our dataset by over two orders of magnitude. We consider how the var functional types are distributed in isolates, and find a nonrandom pattern reflecting that common var functional types are non-randomly distinct from one another in terms of their functional composition. The coarse-graining of var gene diversity into biologically meaningful functional groups has important implications for understanding the disease ecology and evolution of this system, as well as for designing effective epidemiological monitoring and intervention.

Functional diversification accompanies gene family expansion of MED2 homologs in Candida albicans

Gene duplication facilitates functional diversification and provides greater phenotypic flexibility to an organism. Expanded gene families arise through repeated gene duplication but the extent of functional divergence that accompanies each paralogous gene is generally unexplored because of the difficulty in isolating the effects of single family members. The telomere-associated (TLO) gene family is a remarkable example of gene family expansion, with 14 members in the more pathogenic Candida albicans relative to two TLO genes in the closely-related species C. dubliniensis. TLO genes encode interchangeable Med2 subunits of the major transcriptional regulatory complex Mediator. To identify biological functions associated with each C. albicans TLO, expression of individual family members was regulated using a Tet-ON system and the strains were assessed across a range of phenotypes involved in growth and virulence traits. All TLOs affected multiple phenotypes and a single phenotype was often affected by multiple TLOs, including simple phenotypes such as cell aggregation and complex phenotypes such as virulence in a Galleria mellonella model of infection. No phenotype was regulated by all TLOs, suggesting neofunctionalization or subfunctionalization of ancestral properties among different family members. Importantly, regulation of three phenotypes could be mapped to individual polymorphic sites among the TLO genes, including an indel correlated with two phenotypes, growth in sucrose and macrophage killing. Different selective pressures have operated on the TLO sequence, with the 5’ conserved Med2 domain experiencing purifying selection and the gene/clade-specific 3’ end undergoing extensive positive selection that may contribute to the impact of individual TLOs on phenotypic variability. Therefore, expansion of the TLO gene family has conferred unique regulatory properties to each paralog such that it influences a range of phenotypes. We posit that the genetic diversity associated with this expansion contributed to C. albicans success as a commensal and opportunistic pathogen.

Gene duplication is a rapid mechanism to generate additional sequences for natural selection to act upon and confer greater organismal fitness. If additional copies of the gene are beneficial, this process may be repeated to produce an expanded gene family containing many copies of related sequences. Following duplication, individual gene family members may retain functions of the ancestral gene or acquire new functions through mutation. How functional diversification accompanies expansion into large gene families remains largely unexplored due to the difficulty in assessing individual genes in the presence of the remaining family members. Here, we addressed this question using an inducible promoter to regulate expression of individual genes of the TLO gene family in the commensal yeast and opportunistic pathogen Candida albicans, which encode components of a major transcriptional regulator. Induced expression of individual TLOs affected a wide range of phenotypes such that significant functional overlap occurred among TLO genes and most phenotypes were affected by more than one TLO. Induced expression of individual TLOs did not produce massive phenotypic effects in most cases, suggesting that functional overlap among TLO genes may buffer new mutations that arise. Specific sequence variants among the TLO genes correlated with certain phenotypes and these sequence variants did not necessarily correlate with sequence similarity across the entire gene. Therefore, individual TLO family members evolved specific functional roles following duplication that likely reflect a combination of inherited function and new mutation.

Signatures of competition and strain structure within the major blood-stage antigen of Plasmodium falciparum in a local community in Ghana

The concept of niche partitioning has received considerable theoretical attention at the interface of ecology and evolution of infectious diseases. Strain theory postulates that pathogen populations can be structured into distinct nonoverlapping strains by frequency-dependent selection in response to intraspecific competition for host immune space. The malaria parasite Plasmodium falciparum presents an opportunity to investigate this phenomenon in nature, under conditions of high recombination rate and extensive antigenic diversity. The parasite's major blood-stage antigen, Pf EMP1, is encoded by the hyperdiverse var genes. With a dataset that includes thousands of var DBLα sequence types sampled from asymptomatic cases within an area of high endemicity in Ghana, we address how var diversity is distributed within isolates and compare this to the distribution of microsatellite allelic diversity within isolates to test whether antigenic and neutral regions of the genome are structured differently. With respect to var DBLα sequence types, we find that on average isolates exhibit significantly lower overlap than expected randomly, but that there also exists frequent pairs of isolates that are highly related. Furthermore, the linkage network of var DBLα sequence types reveals a pattern of nonrandom modularity unique to these antigenic genes, and we find that modules of highly linked DBLα types are not explainable by neutral forces related to var recombination constraints, microsatellite diversity, sampling location, host age, or multiplicity of infection. These findings of reduced overlap and modularity among the var antigenic genes are consistent with a role for immune selection as proposed by strain theory. Identifying the evolutionary and ecological dynamics that are responsible for the nonrandom structure in P. falciparum antigenic diversity is important for designing effective intervention in endemic areas.

PhyloChromoMap, a Tool for Mapping Phylogenomic History along Chromosomes, Reveals the Dynamic Nature of Karyotype Evolution in Plasmodium falciparum

The genome of Plasmodium falciparum, the causative agent of malaria in Africa, has been extensively studied since it was first fully sequenced in 2002. However, many open questions remain, including understanding the chromosomal context of molecular evolutionary changes (e.g., relationship between chromosome map and phylogenetic conservation, patterns of gene duplication, and patterns of selection). Here, we present PhyloChromoMap, a method that generates a phylogenomic map of chromosomes from a custom-built bioinformatics pipeline. Using P. falciparum 3D7 as a model, we analyze 2,116 genes with homologs in up to 941 diverse eukaryotic, bacterial and archaeal lineages. We estimate the level of conservation along chromosomes based on conservation across clades, and identify “young” regions (i.e., those with recent or fast evolving genes) that are enriched in subtelomeric regions as compared with internal regions. We also demonstrate that patterns of molecular evolution for paralogous genes differ significantly depending on their location as younger paralogs tend to be found in subtelomeric regions whereas older paralogs are enriched in internal regions. Combining these observations with analyses of synteny, we demonstrate that subtelomeric regions are actively shuffled among chromosome ends, which is consistent with the hypothesis that these regions are prone to ectopic recombination. We also assess patterns of selection by comparing dN/dS ratios of gene family members in subtelomeric versus internal regions, and we include the important antigenic gene family var. These analyses illustrate the highly dynamic nature of the karyotype of P. falciparum, and provide a method for exploring genome dynamics in other lineages.

Characterizing Cell Heterogeneity Using PCR Fingerprinting of Surface Multigene Families in Protozoan Parasites

Parasites counteract the action of the immune system and other environmental pressures by modulating and changing the composition of their cell surfaces. Surface multigene protein families are defined not only by highly variable regions in length and/or sequence exposed to the outer space but also by conserved sequences codifying for the signal peptide, hydrophobic C-terminal regions necessary for GPI modifications, as well as conserved UTR regions for mRNA regulation. The method here presented exploits these conserved signatures for characterizing variations in the mRNA expression of clonal cell populations of protozoan parasites using a combination of nested PCR amplification and capillary electrophoresis. With this workflow, in silico gels from isolated cell clones can be generated, thus providing an excellent tool for analyzing cellular heterogeneity in protozoan parasites.

Variant Gene Expression and Antigenic Variation by Malaria Parasites

Malaria is a significant threat throughout the developing world. Among the most fascinating aspects of the protozoan parasites responsible for this disease are the methods they employ to avoid the immune system and perpetuate chronic infections. Key among these is antigenic variation: By systematically altering antigens that are displayed to the host's immune system, the parasite renders the adaptive immune response ineffective. For Plasmodium falciparum, the species responsible for the most severe form of human malaria, this process involves a complicated molecular mechanism that results in continuously changing patterns of variant-antigen-encoding gene expression. Although many features of this process remain obscure, significant progress has been made in recent years to decipher various molecular aspects of the regulatory cascade that causes chronic infection.

Pathogenicity Determinants of the Human Malaria Parasite Plasmodium falciparum Have Ancient Origins

Cytoadhesion of P. falciparum-infected erythrocytes in the microcirculation is a major virulence determinant. P. falciparum is descended from a subgenus of parasites that also infect chimpanzees and gorillas and exhibits strict host species specificity. Despite their high genetic similarity to P. falciparum, it is unknown whether ape parasites encode adhesion properties similar to those of P. falciparum or are as virulent in their natural hosts. Consequently, it has been unclear when virulent adhesion traits arose in P. falciparum and how long they have been present in the parasite population. It is also unknown whether cytoadhesive interactions pose a barrier to cross-species transmission. We show that parasite domains from the chimpanzee malaria parasite P. reichenowi bind human receptors with specificity similar to that of P. falciparum. Our findings suggest that parasite adhesion traits associated with both mild and severe malaria have much earlier origins than previously appreciated and have important implications for virulence evolution in a major human pathogen.

Plasmodium falciparum, the most deadly of the human malaria parasites, is a member of the Laverania subgenus that also infects African Great Apes. The virulence of P. falciparum is related to cytoadhesion of infected erythrocytes in microvasculature, but the origin of dangerous parasite adhesion traits is poorly understood. To investigate the evolutionary history of the P. falciparum cytoadhesion pathogenicity determinant, we studied adhesion domains from the chimpanzee malaria parasite P. reichenowi. We demonstrate that the P. reichenowi var gene repertoire encodes cysteine-rich interdomain region (CIDR) domains which bind human CD36 and endothelial protein C receptor (EPCR) with the same levels of affinity and at binding sites similar to those bound by P. falciparum. Moreover, P. reichenowi domains interfere with the protective function of the activated protein C-EPCR pathway on endothelial cells, a presumptive virulence trait in humans. These findings provide evidence for ancient evolutionary origins of two key cytoadhesion properties of P. falciparum that contribute to human infection and pathogenicity.

IMPORTANCE Cytoadhesion of P. falciparum-infected erythrocytes in the microcirculation is a major virulence determinant. P. falciparum is descended from a subgenus of parasites that also infect chimpanzees and gorillas and exhibits strict host species specificity. Despite their high genetic similarity to P. falciparum, it is unknown whether ape parasites encode adhesion properties similar to those of P. falciparum or are as virulent in their natural hosts. Consequently, it has been unclear when virulent adhesion traits arose in P. falciparum and how long they have been present in the parasite population. It is also unknown whether cytoadhesive interactions pose a barrier to cross-species transmission. We show that parasite domains from the chimpanzee malaria parasite P. reichenowi bind human receptors with specificity similar to that of P. falciparum. Our findings suggest that parasite adhesion traits associated with both mild and severe malaria have much earlier origins than previously appreciated and have important implications for virulence evolution in a major human pathogen.

Maintenance of phenotypic diversity within a set of virulence encoding genes of the malaria parasite Plasmodium falciparum

Infection by the human malaria parasite Plasmodium falciparum results in a broad spectrum of clinical outcomes, ranging from severe and potentially life-threatening malaria to asymptomatic carriage. In a process of naturally acquired immunity, individuals living in malaria-endemic regions build up a level of clinical protection, which attenuates infection severity in an exposure-dependent manner. Underlying this shift in the immunoepidemiology as well as the observed range in malaria pathogenesis is the var multigene family and the phenotypic diversity embedded within. The var gene-encoded surface proteins Plasmodium falciparum erythrocyte membrane protein 1 mediate variant-specific binding of infected red blood cells to a diverse set of host receptors that has been linked to specific disease manifestations, including cerebral and pregnancy-associated malaria. Here, we show that cross-reactive immune responses, which minimize the within-host benefit of each additionally expressed gene during infection, can cause selection for maximum phenotypic diversity at the genome level. We further show that differential functional constraints on protein diversification stably maintain uneven ratios between phenotypic groups, in line with empirical observation. Our results thus suggest that the maintenance of phenotypic diversity within P. falciparum is driven by an evolutionary trade-off that optimizes between within-host parasite fitness and between-host selection pressure.

Recent advances in malaria genomics and epigenomics

Malaria continues to impose a significant disease burden on low- and middle-income countries in the tropics. However, revolutionary progress over the last 3 years in nucleic acid sequencing, reverse genetics, and post-genome analyses has generated step changes in our understanding of malaria parasite (Plasmodium spp.) biology and its interactions with its host and vector. Driven by the availability of vast amounts of genome sequence data from Plasmodium species strains, relevant human populations of different ethnicities, and mosquito vectors, researchers can consider any biological component of the malarial process in isolation or in the interactive setting that is infection. In particular, considerable progress has been made in the area of population genomics, with Plasmodium falciparum serving as a highly relevant model. Such studies have demonstrated that genome evolution under strong selective pressure can be detected. These data, combined with reverse genetics, have enabled the identification of the region of the P. falciparum genome that is under selective pressure and the confirmation of the functionality of the mutations in the kelch13 gene that accompany resistance to the major frontline antimalarial, artemisinin. Furthermore, the central role of epigenetic regulation of gene expression and antigenic variation and developmental fate in P. falciparum is becoming ever clearer. This review summarizes recent exciting discoveries that genome technologies have enabled in malaria research and highlights some of their applications to healthcare. The knowledge gained will help to develop surveillance approaches for the emergence or spread of drug resistance and to identify new targets for the development of antimalarial drugs and perhaps vaccines.

Recombination Rate Heterogeneity within Arabidopsis Disease Resistance Genes

Meiotic crossover frequency varies extensively along chromosomes and is typically concentrated in hotspots. As recombination increases genetic diversity, hotspots are predicted to occur at immunity genes, where variation may be beneficial. A major component of plant immunity is recognition of pathogen Avirulence (Avr) effectors by resistance (R) genes that encode NBS-LRR domain proteins. Therefore, we sought to test whether NBS-LRR genes would overlap with meiotic crossover hotspots using experimental genetics in Arabidopsis thaliana. NBS-LRR genes tend to physically cluster in plant genomes; for example, in Arabidopsis most are located in large clusters on the south arms of chromosomes 1 and 5. We experimentally mapped 1,439 crossovers within these clusters and observed NBS-LRR gene associated hotspots, which were also detected as historical hotspots via analysis of linkage disequilibrium. However, we also observed NBS-LRR gene coldspots, which in some cases correlate with structural heterozygosity. To study recombination at the fine-scale we used high-throughput sequencing to analyze ~1,000 crossovers within the RESISTANCE TO ALBUGO CANDIDA1 (RAC1) R gene hotspot. This revealed elevated intragenic crossovers, overlapping nucleosome-occupied exons that encode the TIR, NBS and LRR domains. The highest RAC1 recombination frequency was promoter-proximal and overlapped CTT-repeat DNA sequence motifs, which have previously been associated with plant crossover hotspots. Additionally, we show a significant influence of natural genetic variation on NBS-LRR cluster recombination rates, using crosses between Arabidopsis ecotypes. In conclusion, we show that a subset of NBS-LRR genes are strong hotspots, whereas others are coldspots. This reveals a complex recombination landscape in Arabidopsis NBS-LRR genes, which we propose results from varying coevolutionary pressures exerted by host-pathogen relationships, and is influenced by structural heterozygosity.

Polymorphic Mucin-Like Proteins in Schistosoma mansoni, a Variable Antigen and a Key Component of the Compatibility Between the Schistosome and Its Snail Host

The arms race between vertebrate hosts and parasites has led to diversification systems able to generate huge repertoires of immune recognition receptors and antigenic variants. Until recently, the invertebrate immunity was considered to be poorly specific, and consequently, antigenic variability was not expected to be high for their respective parasites. In the present chapter, we show how the study of the interaction between the snail Biomphalaria glabrata and its parasite Schistosome mansoni has shaken this paradigm. We show that the fate of the interaction between the snail and its parasite is at least partly the result of the concordance of highly variable repertoires of immune recognition receptors in the snail and corresponding antigenic variants in the parasite. We call these antigenic variants of the schistosome Schistosoma mansoni polymorphic mucins (SmPoMucs). We show that their high level of diversification is the result of a complex cascade of mechanisms, thus presenting evidence for antigenic variation in a parasite infecting an invertebrate species.

The Alba protein family: Structure and function

Alba family proteins are small, basic, dimeric nucleic acid-binding proteins, which are widely distributed in archaea and a number of eukaryotes. This family of proteins bears the distinct features of regulation through acetylation/deacetylation, hence named as acetylation lowers binding affinity (Alba). Alba family proteins bind DNA cooperatively with no apparent sequence specificity. Besides DNA, Alba proteins also interact with diverse RNA species and associate with ribonucleo-protein complexes. Initially, Alba proteins were recognized as chromosomal proteins and supposed to be involved in the maintenance of chromatin architecture and transcription repression. However, recent studies have shown increasing evidence of functional plasticity among Alba family of proteins that widely range from genome packaging and organization, transcriptional and translational regulation, RNA metabolism, and development and differentiation processes. In recent years, Alba family proteins have attracted growing interest due to their widespread occurrence in large number of organisms. Presence in multiple copies, functional crosstalk, differential binding affinity, and posttranslational modifications are some of the key factors that might regulate the biological functions of Alba family proteins. In this review article, we present an overview of the Alba family proteins, their salient features and emphasize their functional role in different organisms reported so far.

Population Structure Shapes Copy Number Variation in Malaria Parasites

If copy number variants (CNVs) are predominantly deleterious, we would expect them to be more efficiently purged from populations with a large effective population size (Ne) than from populations with a small Ne. Malaria parasites (Plasmodium falciparum) provide an excellent organism to examine this prediction, because this protozoan shows a broad spectrum of population structures within a single species, with large, stable, outbred populations in Africa, small unstable inbred populations in South America and with intermediate population characteristics in South East Asia. We characterized 122 single-clone parasites, without prior laboratory culture, from malaria-infected patients in seven countries in Africa, South East Asia and South America using a high-density single-nucleotide polymorphism/CNV microarray. We scored 134 high-confidence CNVs across the parasite exome, including 33 deletions and 102 amplifications, which ranged in size from <500 bp to 59 kb, as well as 10,107 flanking, biallelic single-nucleotide polymorphisms. Overall, CNVs were rare, small, and skewed toward low frequency variants, consistent with the deleterious model. Relative to African and South East Asian populations, CNVs were significantly more common in South America, showed significantly less skew in allele frequencies, and were significantly larger. On this background of low frequency CNV, we also identified several high-frequency CNVs under putative positive selection using an FST outlier analysis. These included known adaptive CNVs containing rh2b and pfmdr1, and several other CNVs (e.g., DNA helicase and three conserved proteins) that require further investigation. Our data are consistent with a significant impact of genetic structure on CNV burden in an important human pathogen.

Phage-fused epitopes from Leishmania infantum used as immunogenic vaccines confer partial protection against Leishmania amazonensis infection

Two mimotopes ofLeishmania infantumidentified by phage display were evaluated as vaccine candidates in BALB/c mice againstLeishmania amazonensisinfection. The epitope-based immunogens, namely B10 and C01, presented as phage-fused peptides; were used without association of a Th1 adjuvant, and they were administered isolated or in combination into animals. Both clones showed a specific production of interferon-gamma (IFN-γ), interleukin-12 (IL-12) and granulocyte/macrophage colony-stimulating factor (GM-CSF) afterin vitrospleen cells stimulation, and they were able to induce a partial protection against infection. Significant reductions of parasite load in the infected footpads, liver, spleen, bone marrow and paws’ draining lymph nodes were observed in the immunized mice, in comparison with the control groups (saline, saponin, wild-type and non-relevant clones). Protection was associated with an IL-12-dependent production of IFN-γ, mediated mainly by CD8+T cells, against parasite proteins. Protected mice also presented low levels of IL-4 and IL-10, as well as increased levels of parasite-specific IgG2a antibodies. The association of both clones resulted in an improved protection in relation to their individual use. More importantly, the absence of adjuvant did not diminish the cross-protective efficacy againstLeishmaniaspp. infection. This study describes for the first time two epitope-based immunogens selected by phage display technology againstL. infantuminfected dogs sera, which induced a partial protection in BALB/c mice infected withL. amazonensis.

Malaria medicines: a glass half full?

Despite substantial scientific progress over the past two decades, malaria remains a worldwide burden that causes hundreds of thousands of deaths every year. New, affordable and safe drugs are required to overcome increasing resistance against artemisinin-based treatments, treat vulnerable populations, interrupt the parasite life cycle by blocking transmission to the vectors, prevent infection and target malaria species that transiently remain dormant in the liver. In this Review, we discuss how the antimalarial drug discovery pipeline has changed over the past 10 years, grouped by the various target compound or product profiles, to assess progress and gaps, and to recommend priorities.

Real-Time Evolution of a Subtelomeric Gene Family in Candida albicans

Subtelomeric regions of the genome are notable for high rates of sequence evolution and rapid gene turnover. Evidence of subtelomeric evolution has relied heavily on comparisons of historical evolutionary patterns to infer trends and frequencies of these events. Here, we describe evolution of the subtelomeric TLO gene family in Candida albicans during laboratory passaging for over 4000 generations. C. albicans is a commensal and opportunistic pathogen of humans and the TLO gene family encodes a subunit of the Mediator complex that regulates transcription and affects a range of virulence factors. We identified 16 distinct subtelomeric recombination events that altered the TLO repertoire. Ectopic recombination between subtelomeres on different chromosome ends occurred approximately once per 5000 generations and was often followed by loss of heterozygosity, resulting in the complete loss of one TLO gene sequence with expansion of another. In one case, recombination within TLO genes produced a novel TLO gene sequence. TLO copy number changes were biased, with some TLOs preferentially being copied to novel chromosome arms and other TLO genes being frequently lost. The majority of these nonreciprocal recombination events occurred either within the 3′ end of the TLO coding sequence or within a conserved 50-bp sequence element centromere-proximal to TLO coding sequence. Thus, subtelomeric recombination is a rapid mechanism of generating genotypic diversity through alterations in the number and sequence of related gene family members.

Phylogeography of var gene repertoires reveals fine-scale geospatial clustering of Plasmodium falciparum populations in a highly endemic area

Plasmodium falciparum malaria is a major global health problem that is being targeted for progressive elimination. Knowledge of local disease transmission patterns in endemic countries is critical to these elimination efforts. To investigate fine-scale patterns of malaria transmission, we have compared repertoires of rapidly evolving var genes in a highly endemic area. A total of 3680 high-quality DBLα-sequences were obtained from 68 P. falciparum isolates from ten villages spread over two distinct catchment areas on the north coast of Papua New Guinea (PNG). Modelling of the extent of var gene diversity in the two parasite populations predicts more than twice as many var gene alleles circulating within each catchment (Mugil = 906; Wosera = 1094) than previously recognized in PNG (Amele = 369). In addition, there were limited levels of var gene sharing between populations, consistent with local parasite population structure. Phylogeographic analyses demonstrate that while neutrally evolving microsatellite markers identified population structure only at the catchment level, var gene repertoires reveal further fine-scale geospatial clustering of parasite isolates. The clustering of parasite isolates by village in Mugil, but not in Wosera was consistent with the physical and cultural isolation of the human populations in the two catchments. The study highlights the microheterogeneity of P. falciparum transmission in highly endemic areas and demonstrates the potential of var genes as markers of local patterns of parasite population structure.

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.

The revolution of whole genome sequencing to study parasites

Genome sequencing has revolutionized the way in which we approach biological research from fundamental molecular biology to ecology and epidemiology. In the last 10 years the field of genomics has changed enormously as technology has improved and the tools for genomic sequencing have moved out of a few dedicated centers and now can be performed on bench-top instruments. In this review we will cover some of the key discoveries that were catalyzed by some of the first genome projects and discuss how this field is developing, what the new challenges are and how this may impact on research in the near future.

The evolutionary dynamics of variant antigen genes in Babesia reveal a history of genomic innovation underlying host-parasite interaction

Babesia spp. are tick-borne, intraerythrocytic hemoparasites that use antigenic variation to resist host immunity, through sequential modification of the parasite-derived variant erythrocyte surface antigen (VESA) expressed on the infected red blood cell surface. We identified the genomic processes driving antigenic diversity in genes encoding VESA (ves1) through comparative analysis within and between three Babesia species, (B. bigemina, B. divergens and B. bovis). Ves1 structure diverges rapidly after speciation, notably through the evolution of shortened forms (ves2) from 5′ ends of canonical ves1 genes. Phylogenetic analyses show that ves1 genes are transposed between loci routinely, whereas ves2 genes are not. Similarly, analysis of sequence mosaicism shows that recombination drives variation in ves1 sequences, but less so for ves2, indicating the adoption of different mechanisms for variation of the two families. Proteomic analysis of the B. bigemina PR isolate shows that two dominant VESA1 proteins are expressed in the population, whereas numerous VESA2 proteins are co-expressed, consistent with differential transcriptional regulation of each family. Hence, VESA2 proteins are abundant and previously unrecognized elements of Babesia biology, with evolutionary dynamics consistently different to those of VESA1, suggesting that their functions are distinct.