Long read assemblies of geographically dispersed Plasmodium falciparum isolates reveal highly structured subtelomeres

Annotation and visualization of parasite, fungi and arthropod genomes with Companion

As sequencing genomes has become increasingly popular, the need for annotation of the resulting assemblies is growing. Structural and functional annotation is still challenging as it includes finding the correct gene sequences, annotating other elements such as RNA and being able to submit those data to databases to share it with the community. Compared to de novo assembly where contiguous chromosomes are a sign of high quality, it is difficult to visualize and assess the quality of annotation. We developed the Companion web server to allow non-experts to annotate their genome using a reference-based method, enabling them to assess the output before submitting to public databases. In this update paper, we describe how we have included novel methods for gene finding and made the Companion server more efficient for annotation of genomes of up to 1 Gb in size. The reference set was increased to include genomes of interest for human and animal health from the fungi and arthropod kingdoms. We show that Companion outperforms existing comparable tools where closely related references are available.

Plasmodium falciparum CLAG Paralogs All Traffic to the Host Membrane but Knockouts Have Distinct Phenotypes

Malaria parasites increase their host erythrocyte's permeability to obtain essential nutrients from plasma and facilitate intracellular growth. In the human Plasmodium falciparum pathogen, this increase is mediated by the plasmodial surface anion channel (PSAC) and has been linked to CLAG3, a protein integral to the host erythrocyte membrane and encoded by a member of the conserved clag multigene family. Whether paralogs encoded by other clag genes also insert at the host membrane is unknown; their contributions to PSAC formation and other roles served are also unexplored. Here, we generated transfectant lines carrying epitope-tagged versions of each CLAG. Each paralog is colocalized with CLAG3, with concordant trafficking via merozoite rhoptries to the host erythrocyte membrane of newly invaded erythrocytes. Each also exists within infected cells in at least two forms: an alkaline-extractable soluble form and a form integral to the host membrane. Like CLAG3, CLAG2 has a variant region cleaved by extracellular proteases, but CLAG8 and CLAG9 are protease resistant. Paralog knockout lines, generated through CRISPR/Cas9 transfection, exhibited uncompromised growth in PGIM, a modified medium with higher physiological nutrient levels; this finding is in marked contrast to a recently reported CLAG3 knockout parasite. CLAG2 and CLAG8 knockout lines exhibited compensatory increases in the transcription of the remaining clags and associated rhoph genes, yielding increased PSAC-mediated uptake for specific solutes. We also report on the distinct transport properties of these knockout lines. Similar membrane topologies at the host membrane are consistent with each CLAG paralog contributing to PSAC, but other roles require further examination.

Imported malaria into Australia: surveillance insights and opportunities

BACKGROUND

Malaria continues to pose a significant burden in endemic countries, many of which lack access to molecular surveillance. Insights from malaria cases in travellers returning to non-endemic areas can provide valuable data to inform endemic country programmes. To evaluate the potential for novel global insights into malaria, we examined epidemiological and molecular data from imported malaria cases to Australia.

METHODS

We analysed malaria cases reported in Australia from 2012 to 2022 using National Notifiable Disease Surveillance System data. Molecular data on imported malaria cases were obtained from literature searches.

RESULTS

Between 2012 and 2022, 3204 malaria cases were reported in Australia. Most cases (69%) were male and 44% occurred in young adults aged 20-39 years. Incidence rates initially declined between 2012 and 2015, then increased until 2019. During 2012-2019, the incidence in travellers ranged from 1.34 to 7.71 per 100 000 trips. Cases were primarily acquired in Sub-Saharan Africa (n = 1433; 45%), Oceania (n = 569; 18%) and Southern and Central Asia (n = 367; 12%). The most common countries of acquisition were Papua New Guinea (n = 474) and India (n = 277). Plasmodium falciparum accounted for 58% (1871/3204) of cases and was predominantly acquired in Sub-Saharan Africa, and Plasmodium vivax accounted for 32% (1016/3204), predominantly from Oceania and Asia. Molecular studies of imported malaria cases to Australia identified genetic mutations and deletions associated with drug resistance and false-negative rapid diagnostic test results, and led to the establishment of reference genomes for P. vivax and Plasmodium malariae.

CONCLUSIONS

Our analysis highlights the continuing burden of imported malaria into Australia. Molecular studies have offered valuable insights into drug resistance and diagnostic limitations, and established reference genomes. Integrating molecular data into national surveillance systems could provide important infectious disease intelligence to optimize treatment guidelines for returning travellers and support endemic country surveillance programmes.

Cryptosporidium Genomics - Current Understanding, Advances, and Applications

Purpose of Review

Here we highlight the significant contribution that genomics-based approaches have had on the field of Cryptosporidium research and the insights these approaches have generated into Cryptosporidium biology and transmission.

Recent Findings

There are advances in genomics, genetic manipulation, gene expression, and single-cell technologies. New and better genome sequences have revealed variable sub-telomeric gene families and genes under selection. RNA expression data now include single-cell and post-infection time points. These data have provided insights into the Cryptosporidium life cycle and host-pathogen interactions. Antisense and ncRNA transcripts are abundant. The critical role of the dsRNA virus is becoming apparent.

Summary

The community's ability to identify genomic targets in the abundant, yet still lacking, collection of genomic data, combined with their increased ability to assess function via gene knock-out, is revolutionizing the field. Advances in the detection of virulence genes, surveillance, population genomics, recombination studies, and epigenetics are upon us.

Transcriptional plasticity of virulence genes provides malaria parasites with greater adaptive capacity for avoiding host immunity

Chronic, asymptomatic malaria infections contribute substantially to disease transmission and likely represent the most significant impediment preventing malaria elimination and eradication. Plasmodium falciparum parasites evade antibody recognition through transcriptional switching between members of the var gene family, which encodes the major virulence factor and surface antigen on infected red blood cells. This process can extend infections for up to a year; however, infections have been documented to last for over a decade, constituting an unseen reservoir of parasites that undermine eradication and control efforts. How parasites remain immunologically "invisible" for such lengthy periods is entirely unknown. Here we show that in addition to the accepted paradigm of mono-allelic var gene expression, individual parasites can simultaneously express multiple var genes or enter a state in which little or no var gene expression is detectable. This unappreciated flexibility provides parasites with greater adaptive capacity than previously understood and challenges the dogma of mutually exclusive var gene expression. It also provides an explanation for the antigenically "invisible" parasites observed in chronic asymptomatic infections.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

Over recent years, progress in molecular markers for genotyping malaria parasites has enabled informative studies of epidemiology and transmission dynamics. Results have highlighted the value of these tools for surveillance to support malaria control and elimination strategies. There are many different types and panels of markers available for malaria parasite genotyping, and for end users, the nuances of these markers with respect to 'use case', resolution, and accuracy, are not well defined. This review clarifies issues surrounding different molecular markers and their application to malaria control and elimination. We describe available marker panels, use cases, implications for different transmission settings, limitations, access, cost, and data accuracy. The information provided can be used as a guide for molecular epidemiology and surveillance of malaria.

Selective whole-genome sequencing of Plasmodium parasites directly from blood samples by nanopore adaptive sampling

IMPORTANCE

Malaria is caused by parasites of the genus Plasmodium, and reached a global disease burden of 247 million cases in 2021. To study drug resistance mutations and parasite population dynamics, whole-genome sequencing of patient blood samples is commonly performed. However, the predominance of human DNA in these samples imposes the need for time-consuming laboratory procedures to enrich Plasmodium DNA. We used the Oxford Nanopore Technologies' adaptive sampling feature to circumvent this problem and enrich Plasmodium reads directly during the sequencing run. We demonstrate that adaptive nanopore sequencing efficiently enriches Plasmodium reads, which simplifies and shortens the timeline from blood collection to parasite sequencing. In addition, we show that the obtained data can be used for monitoring genetic markers, or to generate nearly complete genomes. Finally, owing to its inherent mobility, this technology can be easily applied on-site in endemic areas where patients would benefit the most from genomic surveillance.

Plasmodium falciparum merozoite surface protein 1 as asexual blood stage malaria vaccine candidate

INTRODUCTION

Malaria represents a public health challenge in tropical and subtropical regions, and currently deployed control strategies are likely insufficient to drive elimination of malaria. Development and improvement of malaria vaccines might be key to reduce disease burden. Vaccines targeting asexual blood stages of the parasite have shown limited efficacy when studied in human trials conducted over the past decades.

AREAS COVERED

Vaccine candidates based on the merozoite surface protein 1 (MSP1) were initially envisioned as one of the most promising approaches to provide immune protection against asexual blood-stage malaria. Successful immunization studies in monkey involved the use of the full-length MSP1 (MSP1FL) as vaccine construct. Vaccines using MSP1FL for immunization have the potential benefit of including numerous conserved B-cell and T-cell epitopes. This could result in improved parasite strain-transcending, protective immunity in the field. We review outcomes of clinical trials that utilized a variety of MSP1 constructs and formulations, including MSP1FL, either alone or in combination with other antigens, in both animal models and humans.

EXPERT OPINION

Novel approaches to analyze breadth and magnitude of effector functions of MSP1-targeting antibodies in volunteers undergoing experimental vaccination and controlled human malaria infection will help to define correlates of protective immunity.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

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.

Whole-genome surveillance identifies markers of Plasmodium falciparum drug resistance and novel genomic regions under selection in Mozambique

IMPORTANCE

Malaria is a devastating disease caused by Plasmodium parasites. The evolution of parasite drug resistance continues to hamper progress toward malaria elimination, and despite extensive efforts to control malaria, it remains a leading cause of death in Mozambique and other countries in the region. The development of successful vaccines and identification of molecular markers to track drug efficacy are essential for managing the disease burden. We present an analysis of the parasite genome in Mozambique, a country with one of the highest malaria burdens globally and limited available genomic data, revealing current selection pressure. We contribute additional evidence to limited prior studies supporting the effectiveness of SWGA in producing reliable genomic data from complex clinical samples. Our results provide the identity of genomic loci that may be associated with current antimalarial drug use, including artemisinin and lumefantrine, and reveal selection pressure predicted to compromise the efficacy of current vaccine candidates.

A new Plasmodium vivax reference genome for South American isolates

BACKGROUND

Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering.

RESULTS

This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions.

CONCLUSIONS

Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

Unique Properties of Nutrient Channels on Plasmodium-Infected Erythrocytes

Intracellular malaria parasites activate an ion and organic solute channel on their host erythrocyte membrane to acquire a broad range of essential nutrients. This plasmodial surface anion channel (PSAC) facilitates the uptake of sugars, amino acids, purines, some vitamins, and organic cations, but remarkably, it must exclude the small Na+ ion to preserve infected erythrocyte osmotic stability in plasma. Although molecular, biochemical, and structural studies have provided fundamental mechanistic insights about PSAC and advanced potent inhibitors as exciting antimalarial leads, important questions remain about how nutrients and ions are transported. Here, I review PSAC's unusual selectivity and conductance properties, which should guide future research into this important microbial ion channel.

Telomerase-independent survival leads to a mosaic of complex subtelomere rearrangements in Chlamydomonas reinhardtii

Telomeres and subtelomeres, the genomic regions located at chromosome extremities, are essential for genome stability in eukaryotes. In the absence of the canonical maintenance mechanism provided by telomerase, telomere shortening induces genome instability. The landscape of the ensuing genome rearrangements is not accessible by short-read sequencing. Here, we leverage Oxford Nanopore Technologies long-read sequencing to survey the extensive repertoire of genome rearrangements in telomerase mutants of the model green microalga Chlamydomonas reinhardtii In telomerase-mutant strains grown for hundreds of generations, most chromosome extremities were capped by short telomere sequences that were either recruited de novo from other loci or maintained in a telomerase-independent manner. Other extremities did not end with telomeres but only with repeated subtelomeric sequences. The subtelomeric elements, including rDNA, were massively rearranged and involved in breakage-fusion-bridge cycles, translocations, recombinations, and chromosome circularization. These events were established progressively over time and displayed heterogeneity at the subpopulation level. New telomere-capped extremities composed of sequences originating from more internal genomic regions were associated with high DNA methylation, suggesting that de novo heterochromatin formation contributes to the restoration of chromosome end stability in C. reinhardtii The diversity of alternative strategies present in the same organism to maintain chromosome integrity and the variety of rearrangements found in telomerase mutants are remarkable, and illustrate genome plasticity at short timescales.

Malaria Genomics, Vaccine Development, and Microbiome

Recent advances in malaria genetics and genomics have transformed many aspects of malaria research in areas of molecular evolution, epidemiology, transmission, host-parasite interaction, drug resistance, pathogenicity, and vaccine development. Here, in addition to introducing some background information on malaria parasite biology, parasite genetics/genomics, and genotyping methods, we discuss some applications of genetic and genomic approaches in vaccine development and in studying interactions with microbiota. Genetic and genomic data can be used to search for novel vaccine targets, design an effective vaccine strategy, identify protective antigens in a whole-organism vaccine, and evaluate the efficacy of a vaccine. Microbiota has been shown to influence disease outcomes and vaccine efficacy; studying the effects of microbiota in pathogenicity and immunity may provide information for disease control. Malaria genetics and genomics will continue to contribute greatly to many fields of malaria research.

Measuring changes in Plasmodium falciparum census population size in response to sequential malaria control interventions

Here we introduce a new endpoint "census population size" to evaluate the epidemiology and control of Plasmodium falciparum infections, where the parasite, rather than the infected human host, is the unit of measurement. To calculate census population size, we rely on a definition of parasite variation known as multiplicity of infection (M O I v a r ), based on the hyper-diversity of the v a r multigene family. We present a Bayesian approach to estimate M O I v a r from sequencing and counting the number of unique DBLα tags (or DBLα types) of v a r genes, and derive from it census population size by summation of M O I v a r in the human population. We track changes in this parasite population size and structure through sequential malaria interventions by indoor residual spraying (IRS) and seasonal malaria chemoprevention (SMC) from 2012 to 2017 in an area of high-seasonal malaria transmission in northern Ghana. Following IRS, which reduced transmission intensity by > 90% and decreased parasite prevalence by ~40-50%, significant reductions in v a r diversity, M O I v a r , and population size were observed in ~2,000 humans across all ages. These changes, consistent with the loss of diverse parasite genomes, were short lived and 32-months after IRS was discontinued and SMC was introduced, v a r diversity and population size rebounded in all age groups except for the younger children (1-5 years) targeted by SMC. Despite major perturbations from IRS and SMC interventions, the parasite population remained very large and retained the v a r population genetic characteristics of a high-transmission system (high v a r diversity; low v a r repertoire similarity) demonstrating the resilience of P. falciparum to short-term interventions in high-burden countries of sub-Saharan Africa.

Decoding the impact of nuclear organization on antigenic variation in parasites

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA)

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

An optimized GATK4 pipeline for Plasmodium falciparum whole genome sequencing variant calling and analysis

BACKGROUND

Accurate variant calls from whole genome sequencing (WGS) of Plasmodium falciparum infections are crucial in malaria population genomics. Here a falciparum variant calling pipeline based on GATK version 4 (GATK4) was optimized and applied to 6626 public Illumina WGS samples.

METHODS

Control WGS and accurate PacBio assemblies of 10 laboratory strains were leveraged to optimize parameters that control the heterozygosity, local assembly region size, ploidy, mapping and base quality in both GATK HaplotypeCaller and GenotypeGVCFs. From these controls, a high-quality training dataset was generated to recalibrate the raw variant data.

RESULTS

On current high-quality samples (read length = 250 bp, insert size = 405-524 bp), the optimized pipeline shows improved sensitivity (86.6 ± 1.7% for SNPs and 82.2 ± 5.9% for indels) compared to the default GATK4 pipeline (77.7 ± 1.3% for SNPs; and 73.1 ± 5.1% for indels, adjusted P < 0.001) and previous variant calling with GATK version 3 (GATK3, 70.3 ± 3.0% for SNPs and 59.7 ± 5.8% for indels, adjusted P < 0.001). Its sensitivity on simulated mixed infection samples (80.8 ± 6.1% for SNPs and 78.3 ± 5.1% for indels) was again improved relative to default GATK4 (68.8 ± 6.0% for SNPs and 38.9 ± 0.7% for indels, adjusted, adjusted P < 0.001). Precision was high and comparable across all pipelines on each type of data tested. The resulting combination of high-quality SNPs and indels increases the resolution of local population population structure detection in sub-Saharan Africa. Finally, increasing ploidy improves the detection of drug resistance mutations and estimation of complexity of infection.

CONCLUSIONS

Overall, this study provides an optimized falciparum GATK4 pipeline resource for variant calling which should help improve genomic studies of malaria.

Knobs, Adhesion, and Severe Falciparum Malaria

Plasmodium falciparum can cause a severe disease with high mortality. A major factor contributing to the increased virulence of P. falciparum, as compared to other human malarial parasites, is the sequestration of infected erythrocytes in the capillary beds of organs and tissues. This sequestration is due to the cytoadherence of infected erythrocytes to endothelial cells. Cytoadherence is primarily mediated by a parasite protein expressed on the surface of the infected erythrocyte called P. falciparum erythrocyte membrane protein-1 (PfEMP1). PfEMP1 is embedded in electron-dense protuberances on the surface of the infected erythrocytes called knobs. These knobs are assembled on the erythrocyte membrane via exported parasite proteins, and the knobs function as focal points for the cytoadherence of infected erythrocytes to endothelial cells. PfEMP1 is a member of the var gene family, and there are approximately 60 antigenically distinct PfEMP1 alleles per parasite genome. Var gene expression exhibits allelic exclusion, with only a single allele being expressed by an individual parasite. This results in sequential waves of antigenically distinct infected erythrocytes and this antigenic variation allows the parasite to establish long-term chronic infections. A wide range of endothelial cell receptors can bind to the various PfEMP1 alleles, and thus, antigenic variation also results in a change in the cytoadherence phenotype. The cytoadherence phenotype may result in infected erythrocytes sequestering in different tissues and this difference in sequestration may explain the wide range of possible clinical manifestations associated with severe falciparum malaria.

The exception that proves the rule: Virulence gene expression at the onset of Plasmodium falciparum blood stage infections

Controlled human malaria infections (CHMI) are a valuable tool to study parasite gene expression in vivo under defined conditions. In previous studies, virulence gene expression was analyzed in samples from volunteers infected with the Plasmodium falciparum (Pf) NF54 isolate, which is of African origin. Here, we provide an in-depth investigation of parasite virulence gene expression in malaria-naïve European volunteers undergoing CHMI with the genetically distinct Pf 7G8 clone, originating in Brazil. Differential expression of var genes, encoding major virulence factors of Pf, PfEMP1s, was assessed in ex vivo parasite samples as well as in parasites from the in vitro cell bank culture that was used to generate the sporozoites (SPZ) for CHMI (Sanaria PfSPZ Challenge (7G8)). We report broad activation of mainly B-type subtelomeric located var genes at the onset of a 7G8 blood stage infection in naïve volunteers, mirroring the NF54 expression study and suggesting that the expression of virulence-associated genes is generally reset during transmission from the mosquito to the human host. However, in 7G8 parasites, we additionally detected a continuously expressed single C-type variant, Pf7G8_040025600, that was most highly expressed in both pre-mosquito cell bank and volunteer samples, suggesting that 7G8, unlike NF54, maintains expression of some previously expressed var variants during transmission. This suggests that in a new host, the parasite may preferentially express the variants that previously allowed successful infection and transmission. Trial registration: ClinicalTrials.gov - NCT02704533; 2018-004523-36.

Molecular epidemiology of continued Plasmodium falciparum disease transmission after an outbreak in Ecuador

To better understand the factors underlying the continued incidence of clinical episodes of falciparum malaria in E-2025 countries targeting elimination, we characterized the molecular epidemiology of Plasmodium falciparum disease transmission after a clonal outbreak in Ecuador. Here we study disease transmission by documenting the diversity and population structure of the major variant surface antigen of the blood stages of P. falciparum encoded by the var multigene family. We used a high-resolution genotyping method, “varcoding”, involving targeted amplicon sequencing to fingerprint the DBLα encoding region of var genes to describe both antigenic var diversity and var repertoire similarity or relatedness in parasite isolates from clinical cases. We identified nine genetic varcodes in 58 P. falciparum isolates causing clinical disease in 2013-2015. Network analyses revealed that four of the varcodes were highly related to the outbreak varcode, with identification of possible diversification of the outbreak parasites by recombination as seen in three of those varcodes. The majority of clinical cases in Ecuador were associated with parasites with highly related or recombinant varcodes to the outbreak clone and due to local transmission rather than recent importation of parasites from other endemic countries. Sharing of types in Ecuadorian varcodes to those sampled in South American varcodes reflects historical parasite importation of some varcodes, especially from Colombia and Peru. Our findings highlight the translational application of varcoding for outbreak surveillance in epidemic/unstable malaria transmission, such as in E-2025 countries, and point to the need for surveillance of local reservoirs of infection in Ecuador to achieve the malaria elimination goal by 2025.

An Optimized GATK4 Pipeline for Plasmodium falciparum Whole Genome Sequencing Variant Calling and Analysis

Background Accurate variant calls from whole genome sequencing (WGS) of Plasmodium falciparum infections are crucial in malaria population genomics. Here we optimized a falciparum variant calling pipeline based on GATK version 4 (GATK4) and applied it to 6,626 public Illumina WGS samples. Methods We optimized parameters that control the heterozygosity, local assembly region size, ploidy, mapping and base quality in both GATK HaplotypeCaller and GenotypeGVCFs leveraging control WGS and accurate PacBio assemblies of 10 laboratory strains. From these controls we generated a high-quality training dataset to recalibrate the raw variant data. Results On current high-quality samples (read length = 250bp, insert size = 405 - 524 bp ), we show improved sensitivity (86.6 ± 1.7% for SNPs and 82.2 ± 5.9% for indels) compared to the default GATK4 pipeline (77.7 ± 1.3% for SNPs; and 73.1 ± 5.1% for indels, adjusted P < 0.001) and previous variant calling with GATK version 3 (GATK3, 70.3 ± 3.0% for SNPs and 59.7 ± 5.8% for indels, adjusted P < 0.001). The sensitivity of our pipeline on simulated mixed infection samples (80.8 ± 6.1% for SNPs and 78.3 ± 5.1% for indels) was again improved relative to default GATK4 (68.8 ± 6.0% for SNPs and 38.9 ± 0.7% for indels, adjusted P < 0.001). Precision was high and comparable across all pipelines on each type of data tested. We further show that using the combination of high-quality SNPs and indels increases the resolution of local population population structure detection in sub-Saharan Africa. We finally demonstrate that increasing ploidy improves the detection of drug resistance mutations and estimation of complexity of infection. Conclusions Overall, we provide an optimized GATK4 pipeline and resource for falciparum variant calling which should help improve genomic studies of malaria.

Antigenic diversity in malaria parasites is maintained on extrachromosomal DNA

Sequence variation among antigenic var genes enables Plasmodium falciparum malaria parasites to evade host immunity. Using long sequence reads from haploid clones from a mutation accumulation experiment, we detect var diversity inconsistent with simple chromosomal inheritance. We discover putatively circular DNA that is strongly enriched for var genes, which exist in multiple alleles per locus separated by recombination and indel events. Extrachromosomal DNA likely contributes to rapid antigenic diversification in P. falciparum.

Disulfide bond and crosslinking analyses reveal inter-domain interactions that contribute to the rigidity of placental malaria VAR2CSA structure and formation of CSA binding channel

VAR2CSA, a multidomain Plasmodium falciparum protein, mediates the adherence of parasite-infected red blood cells to chondroitin 4-sulfate (C4S) in the placenta, contributing to placental malaria. Therefore, detailed understanding of VAR2CSA structure likely help developing strategies to treat placental malaria. The VAR2CSA ectodomain consists of an N-terminal segment (NTS), six Duffy binding-like (DBL) domains, and three interdomains (IDs) present in sequence NTS-DBL1x-ID1-DBL2x-ID2-DBL3x-DBL4ε-ID3-DBL5ε-DBL6ε. Recent electron microscopy studies showed that VAR2CSA is compactly organized into a globular structure containing C4S-binding channel, and that DBL5ε-DBL6ε arm is attached to the NTS-ID3 core structure. However, the structural elements involved in inter-domain interactions that stabilize the VAR2CSA structure remain largely not understood. Here, limited proteolysis and peptide mapping by mass spectrometry showed that VAR2CSA contains several inter-domain disulfide bonds that stabilize its compact structure. Chemical crosslinking-mass spectrometry showed that all IDs interact with DBL4ε; additionally, IDs interact with other DBL domains, demonstrating that IDs are the key structural scaffolds that shape the functional NTS-ID3 core. Ligand binding analysis suggested that NTS considerably restricts the C4S binding. Overall, our study revealed that inter-domain disulfide bonds and interactions between IDs and DBL domains contribute to the stability of VAR2CSA structural architecture and formation of C4S-binding channel.

Unravelling var complexity: Relationship between DBLα types and var genes in Plasmodium falciparum

The enormous diversity and complexity of var genes that diversify rapidly by recombination has led to the exclusion of assembly of these genes from major genome initiatives (e.g., Pf6). A scalable solution in epidemiological surveillance of var genes is to use a small ‘tag’ region encoding the immunogenic DBLα domain as a marker to estimate var diversity. As var genes diversify by recombination, it is not clear the extent to which the same tag can appear in multiple var genes. This relationship between marker and gene has not been investigated in natural populations. Analyses of in vitro recombination within and between var genes have suggested that this relationship would not be exclusive. Using a dataset of publicly-available assembled var sequences, we test this hypothesis by studying DBLα-var relationships for four study sites in four countries: Pursat (Cambodia) and Mae Sot (Thailand), representing low malaria transmission, and Navrongo (Ghana) and Chikwawa (Malawi), representing high malaria transmission. In all study sites, DBLα-var relationships were shown to be predominantly 1-to-1, followed by a second largest proportion of 1-to-2 DBLα-var relationships. This finding indicates that DBLα tags can be used to estimate not just DBLα diversity but var gene diversity when applied in a local endemic area. Epidemiological applications of this result are discussed.

The mystery of persistent, asymptomatic Plasmodium falciparum infections

Plasmodium falciparum causes millions of malaria infections and hundreds of thousands of deaths annually. These parasites avoid the adaptive immune response by systematically cycling through a limited repertoire of variant surface antigens after which the number of circulating parasites drops to extremely low levels, coinciding with a loss of symptoms and eventual clearance of the infection. However, in regions with extended dry seasons or in individuals who no longer reside in endemic areas, asymptomatic infections have been observed to persist for many months or years, potentially serving as reservoirs for transmission. Recent work suggests the possibility that parasites can assume a state in which no variant surface antigens are expressed, thus rendering them virtually invisible to the immune system and enabling them to persist at low levels indefinitely.

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.

Fitness Costs of pfhrp2 and pfhrp3 Deletions Underlying Diagnostic Evasion in Malaria Parasites

BACKGROUND

Rapid diagnostic tests based on detection of histidine-rich proteins (HRPs) are widely used for malaria diagnosis, but parasites carrying pfhrp deletions can evade detection and are increasing in frequency in some countries. Models aim to predict conditions under which pfhrp2 and/or pfhrp3 deletions will increase, but a key parameter-the fitness cost of deletions-is unknown.

METHODS

We removed pfhrp2 and/or pfhrp3 from a Malawian parasite clone using gene editing approaches) and measured fitness costs by conducting pairwise competition experiments.

RESULTS

We observed significant fitness costs of 0.087 ± 0.008 (1 standard error) per asexual cycle for pfhrp2 deletion and 0.113 ± 0.008 for the pfhrp2/3 double deletion, relative to the unedited progenitor parasite. Selection against deletions is strong and comparable to that resulting from drug resistance mutations.

CONCLUSIONS

Prior modeling suggested that diagnostic selection may drive increased frequency of pfhrp deletions only when fitness costs are mild. Our experiments show that costs of pfhrp deletions are higher than these thresholds, but modeling and empirical results can be reconciled if the duration of infection is short. These results may inform future modeling to understand why pfhrp2/3 deletions are increasing in some locations (Ethiopia and Eritrea) but not in others (Mekong region).

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

A nuclear redox sensor modulates gene activation and var switching in Plasmodium falciparum

The virulence of Plasmodium falciparum, which causes the deadliest form of human malaria, is attributed to its ability to evade the human immune response. These parasites "choose" to express a single variant from a repertoire of surface antigens called PfEMP1, which are placed on the surface of the infected red cell. Immune evasion is achieved by switches in expression between var genes, each encoding a different PfEMP1 variant. While the mechanisms that regulate mutually exclusive expression of var genes are still elusive, antisense long-noncoding RNAs (lncRNAs) transcribed from the intron of the active var gene were implicated in the "choice" of the single active var gene. Here, we show that this lncRNA colocalizes with the site of var mRNA transcription and is anchored to the var locus via DNA:RNA interactions. We define the var lncRNA interactome and identify a redox sensor, P. falciparum thioredoxin peroxidase I (PfTPx-1), as one of the proteins associated with the var antisense lncRNA. We show that PfTPx-1 localizes to a nuclear subcompartment associated with active transcription on the nuclear periphery, in ring-stage parasite, when var transcription occurs. In addition, PfTPx-1 colocalizes with S-adenosylmethionine synthetase (PfSAMS) in the nucleus, and its overexpression leads to activation of var2csa, similar to overexpression of PfSAMS. Furthermore, we show that PfTPx-1 knockdown alters the var switch rate as well as activation of additional gene subsets. Taken together, our data indicate that nuclear PfTPx-1 plays a role in gene activation possibly by providing a redox-controlled nuclear microenvironment ideal for active transcription.

A GPU-accelerated compute framework for pathogen genomic variant identification to aid genomic epidemiology of infectious disease: a malaria case study

As recently demonstrated by the COVID-19 pandemic, large-scale pathogen genomic data are crucial to characterize transmission patterns of human infectious diseases. Yet, current methods to process raw sequence data into analysis-ready variants remain slow to scale, hampering rapid surveillance efforts and epidemiological investigations for disease control. Here, we introduce an accelerated, scalable, reproducible, and cost-effective framework for pathogen genomic variant identification and present an evaluation of its performance and accuracy across benchmark datasets of Plasmodium falciparum malaria genomes. We demonstrate superior performance of the GPU framework relative to standard pipelines with mean execution time and computational costs reduced by 27× and 4.6×, respectively, while delivering 99.9% accuracy at enhanced reproducibility.

Screening strategies and laboratory assays to support Plasmodium falciparum histidine-rich protein deletion surveillance: where we are and what is needed

Rapid diagnostic tests (RDTs) detecting Plasmodium falciparum histidine-rich protein 2 (HRP2) have been an important tool for malaria diagnosis, especially in resource-limited settings lacking quality microscopy. Plasmodium falciparum parasites with deletion of the pfhrp2 gene encoding this antigen have now been identified in dozens of countries across Asia, Africa, and South America, with new reports revealing a high prevalence of deletions in some selected regions. To determine whether HRP2-based RDTs are appropriate for continued use in a locality, focused surveys and/or surveillance activities of the endemic P. falciparum population are needed. Various survey and laboratory methods have been used to determine parasite HRP2 phenotype and pfhrp2 genotype, and the data collected by these different methods need to be interpreted in the appropriate context of survey and assay utilized. Expression of the HRP2 antigen can be evaluated using point-of-care RDTs or laboratory-based immunoassays, but confirmation of a deletion (or mutation) of pfhrp2 requires more intensive laboratory molecular assays, and new tools and strategies for rigorous but practical data collection are particularly needed for large surveys. Because malaria diagnostic strategies are typically developed at the national level, nationally representative surveys and/or surveillance that encompass broad geographical areas and large populations may be required. Here is discussed contemporary assays for the phenotypic and genotypic evaluation of P. falciparum HRP2 status, consider their strengths and weaknesses, and highlight key concepts relevant to timely and resource-conscious workflows required for efficient diagnostic policy decision making.

Characterising genome architectures using genome decomposition analysis

Genome architecture describes how genes and other features are arranged in genomes. These arrangements reflect the evolutionary pressures on genomes and underlie biological processes such as chromosomal segregation and the regulation of gene expression. We present a new tool called Genome Decomposition Analysis (GDA) that characterises genome architectures and acts as an accessible approach for discovering hidden features of a genome assembly. With the imminent deluge of high-quality genome assemblies from projects such as the Darwin Tree of Life and the Earth BioGenome Project, GDA has been designed to facilitate their exploration and the discovery of novel genome biology. We highlight the effectiveness of our approach in characterising the genome architectures of single-celled eukaryotic parasites from the phylum Apicomplexa and show that it scales well to large genomes.

De Novo Assembly of Plasmodium knowlesi Genomes From Clinical Samples Explains the Counterintuitive Intrachromosomal Organization of Variant SICAvar and kir Multiple Gene Family Members

Plasmodium knowlesi, a malaria parasite of Old World macaque monkeys, is used extensively to model Plasmodium biology. Recently, P. knowlesi was found in the human population of Southeast Asia, particularly Malaysia. P. knowlesi causes uncomplicated to severe and fatal malaria in the human host with features in common with the more prevalent and virulent malaria caused by Plasmodium falciparum. As such, P. knowlesi presents a unique opportunity to develop experimental translational model systems for malaria pathophysiology informed by clinical data from same-species human infections. Experimental lines of P. knowlesi represent well-characterized genetically stable parasites, and to maximize their utility as a backdrop for understanding malaria pathophysiology, genetically diverse contemporary clinical isolates, essentially wild-type, require comparable characterization. The Oxford Nanopore PCR-free long-read sequencing platform was used to sequence and de novo assemble P. knowlesi genomes from frozen clinical samples. The sequencing platform and assembly pipelines were designed to facilitate capturing data and describing, for the first time, P. knowlesi schizont-infected cell agglutination (SICA) var and Knowlesi-Interspersed Repeats (kir) multiple gene families in parasites acquired from nature. The SICAvar gene family members code for antigenically variant proteins analogous to the virulence-associated P. falciparum erythrocyte membrane protein (PfEMP1) multiple var gene family. Evidence presented here suggests that the SICAvar family members have arisen through a process of gene duplication, selection pressure, and variation. Highly evolving genes including PfEMP1family members tend to be restricted to relatively unstable sub-telomeric regions that drive change with core genes protected in genetically stable intrachromosomal locations. The comparable SICAvar and kir gene family members are counter-intuitively located across chromosomes. Here, we demonstrate that, in contrast to conserved core genes, SICAvar and kir genes occupy otherwise gene-sparse chromosomal locations that accommodate rapid evolution and change. The novel methods presented here offer the malaria research community not only new tools to generate comprehensive genome sequence data from small clinical samples but also new insight into the complexity of clinically important real-world parasites.

Epigenetics of malaria parasite nutrient uptake, but why?

The conserved plasmodial surface anion channel (PSAC) mediates nutrient uptake by bloodstream malaria parasites and is an antimalarial target. This pathogen-associated channel is linked to the clag multigene family, which is variably expanded in Plasmodium spp. Member genes are under complex epigenetic regulation, with the clag3 genes of the human P. falciparum pathogen exhibiting monoallelic transcription and mutually exclusive surface exposure on infected erythrocytes. While other multigene families use monoallelic expression to evade host immunity, the reasons of epigenetic control of clag genes are unclear. I consider existing models and their implications for nutrient acquisition and immune evasion. Understanding the reasons for epigenetic regulation of PSAC-mediated nutrient uptake will help clarify host-pathogen interactions and guide development of therapies resistant to allele switching.

Shared Mechanisms for Mutually Exclusive Expression and Antigenic Variation by Protozoan Parasites

Cellular decision-making at the level of gene expression is a key process in the development and evolution of every organism. Variations in gene expression can lead to phenotypic diversity and the development of subpopulations with adaptive advantages. A prime example is the mutually exclusive activation of a single gene from within a multicopy gene family. In mammals, this ranges from the activation of one of the two immunoglobulin (Ig) alleles to the choice in olfactory sensory neurons of a single odorant receptor (OR) gene from a family of more than 1,000. Similarly, in parasites like Trypanosoma brucei, Giardia lamblia or Plasmodium falciparum, the process of antigenic variation required to escape recognition by the host immune system involves the monoallelic expression of vsg, vsp or var genes, respectively. Despite the importance of this process, understanding how this choice is made remains an enigma. The development of powerful techniques such as single cell RNA-seq and Hi-C has provided new insights into the mechanisms these different systems employ to achieve monoallelic gene expression. Studies utilizing these techniques have shown how the complex interplay between nuclear architecture, physical interactions between chromosomes and different chromatin states lead to single allele expression. Additionally, in several instances it has been observed that high-level expression of a single gene is preceded by a transient state where multiple genes are expressed at a low level. In this review, we will describe and compare the different strategies that organisms have evolved to choose one gene from within a large family and how parasites employ this strategy to ensure survival within their hosts.

Varia: a tool for prediction, analysis and visualisation of variable genes

Background

Parasites use polymorphic gene families to evade the immune system or interact with the host. Assessing the diversity and expression of such gene families in pathogens can inform on the repertoire or host interaction phenotypes of clinical relevance. However, obtaining the sequences and quantifying their expression is a challenge. In Plasmodium falciparum, the highly polymorphic var genes encode the major virulence protein, PfEMP1, which bind a range of human receptors through varying combinations of DBL and CIDR domains. Here we present a tool, Varia, to predict near full-length gene sequences and domain compositions of query genes from database genes sharing short sequence tags. Varia generates output through two complementary pipelines. Varia_VIP returns all putative gene sequences and domain compositions of the query gene from any partial sequence provided, thereby enabling experimental validation of specific genes of interest and detailed assessment of their putative domain structure. Varia_GEM accommodates rapid profiling of var gene expression in complex patient samples from DBLα expression sequence tags (EST), by computing a sample overall transcript profile stratified by PfEMP1 domain types.

Results

Varia_VIP was tested querying sequence tags from all DBL domain types using different search criteria. On average 92% of query tags had one or more 99% identical database hits, resulting in the full-length query gene sequence being identified (> 99% identical DNA > 80% of query gene) among the five most prominent database hits, for ~ 33% of the query genes. Optimized Varia_GEM settings allowed correct prediction of > 90% of domains placed among the four most N-terminal domains, including the DBLα domain, and > 70% of C-terminal domains. With this accuracy, N-terminal domains could be predicted for > 80% of queries, whereas prediction rates of C-terminal domains dropped with the distance from the DBLα from 70 to 40%.

Conclusion

Prediction of var sequence and domain composition is possible from short sequence tags. Varia can be used to guide experimental validation of PfEMP1 sequences of interest and conduct high-throughput analysis of var type expression in patient samples.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04573-6.

STRIDE: a command-line HMM-based identifier and sub-classifier of Plasmodium falciparum RIFIN and STEVOR variant surface antigen families

BACKGROUND

RIFINs and STEVORs are variant surface antigens expressed by P. falciparum that play roles in severe malaria pathogenesis and immune evasion. These two highly diverse multigene families feature multiple paralogs, making their classification challenging using traditional bioinformatic methods.

RESULTS

STRIDE (STevor and RIfin iDEntifier) is an HMM-based, command-line program that automates the identification and classification of RIFIN and STEVOR protein sequences in the malaria parasite Plasmodium falciparum. STRIDE is more sensitive in detecting RIFINs and STEVORs than available PFAM and TIGRFAM tools and reports RIFIN subtypes and the number of sequences with a FHEYDER amino acid motif, which has been associated with severe malaria pathogenesis.

CONCLUSIONS

STRIDE will be beneficial to malaria research groups analyzing genome sequences and transcripts of clinical field isolates, providing insight into parasite biology and virulence.

A coordinated transcriptional switching network mediates antigenic variation of human malaria parasites

Malaria parasites avoid immune clearance through their ability to systematically alter antigens exposed on the surface of infected red blood cells. This is accomplished by tightly regulated transcriptional control of individual members of a large, multicopy gene family called var and is the key to both the virulence and chronic nature of malaria infections. Expression of var genes is mutually exclusive and controlled epigenetically, however how large populations of parasites coordinate var gene switching to avoid premature exposure of the antigenic repertoire is unknown. Here, we provide evidence for a transcriptional network anchored by a universally conserved gene called var2csa that coordinates the switching process. We describe a structured switching bias that shifts overtime and could shape the pattern of var expression over the course of a lengthy infection. Our results provide an explanation for a previously mysterious aspect of malaria infections and shed light on how parasites possessing a relatively small repertoire of variant antigen-encoding genes can coordinate switching events to limit antigen exposure, thereby maintaining chronic infections.

Resetting var Gene Transcription in Plasmodium falciparum

One of the key mechanisms contributing to the virulence of Plasmodium falciparum is its ability to undergo antigenic switching among antigenically distinct variants of the PfEMP1 adhesive proteins, encoded by the var gene family. To avoid premature exposure of its antigenic repertoire, the parasite transcribes its var genes in a mutually exclusive manner, and switch expression at a very slow rate. This process is epigenetically regulated and it relies on "epigenetic memory," which imprints the single active var gene to remain active for multiple replication cycles. Erasing this epigenetic memory in parasites grown in culture resembles parasites, which egress from the liver. It could therefore be of interest for investigating var switching patterns at the onset of malaria infections. In addition, this procedure could be used for creating heterogeneity of var expression among parasite populations. The methodology described here for resetting of var gene expression is based on promoter titration, also known as molecular sponging.

Successful Profiling of Plasmodium falciparum var Gene Expression in Clinical Samples via a Custom Capture Array

var genes encode Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens. These highly diverse antigens are displayed on the surface of infected erythrocytes and play a critical role in immune evasion and sequestration of infected erythrocytes. Studies of var expression using non-leukocyte-depleted blood are challenging because of the predominance of host genetic material and lack of conserved var segments. Our goal was to enrich for parasite RNA, allowing de novo assembly of var genes and detection of expressed novel variants. We used two overall approaches: (i) enriching for total mRNA in the sequencing library preparations and (ii) enriching for parasite RNA with a custom capture array based on Roche’s SeqCap EZ enrichment system. The capture array was designed with probes based on the whole 3D7 reference genome and an additional >4,000 full-length var gene sequences from other P. falciparum strains. We tested each method on RNA samples from Malian children with severe or uncomplicated malaria infections. All reads mapping to the human genome were removed, the remaining reads were assembled de novo into transcripts, and from these, var-like transcripts were identified and annotated. The capture array produced the longest maximum length and largest numbers of var gene transcripts in each sample, particularly in samples with low parasitemia. Identifying the most-expressed var gene sequences in whole-blood clinical samples without the need for extensive processing or generating sample-specific reference genome data is critical for understanding the role of PfEMP1s in malaria pathogenesis.

IMPORTANCE Malaria parasites display antigens on the surface of infected red blood cells in the human host that facilitate attachment to blood vessels, contributing to the severity of infection. These antigens are highly variable, allowing the parasite to evade the immune system. Identifying these expressed antigens is critical to understanding the development of severe malarial disease. However, clinical samples contain limited amounts of parasite genetic material, a challenge for sequencing efforts further compounded by the extreme diversity of the parasite surface antigens. We present a method that enriches for these antigen sequences in clinical samples using a custom capture array, requiring minimal processing in the field. While our results are focused on the malaria parasite Plasmodium falciparum, this approach has broad applicability to other highly diverse antigens from other parasites and pathogens such as those that cause giardiasis and leishmaniasis.

Malaria protection due to sickle haemoglobin depends on parasite genotype

Host genetic factors can confer resistance against malaria¹, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member^2–4PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

A strong association has been found between three regions of the Plasmodium falciparum genome and sickle haemoglobin in children with severe malaria, suggesting parasites have adapted to overcome natural host immunity.

Characterization of Apicomplexan Amino Acid Transporters (ApiATs) in the Malaria Parasite Plasmodium falciparum

During the symptomatic human blood phase, malaria parasites replicate within red blood cells. Parasite proliferation relies on the uptake of nutrients, such as amino acids, from the host cell and blood plasma, requiring transport across multiple membranes. Amino acids are delivered to the parasite through the parasite-surrounding vacuolar compartment by specialized nutrient-permeable channels of the erythrocyte membrane and the parasitophorous vacuole membrane (PVM). However, further transport of amino acids across the parasite plasma membrane (PPM) is currently not well characterized. In this study, we focused on a family of Apicomplexan amino acid transporters (ApiATs) that comprises five members in Plasmodium falciparum. First, we localized four of the P. falciparum ApiATs (PfApiATs) at the PPM using endogenous green fluorescent protein (GFP) tagging. Next, we applied reverse genetic approaches to probe into their essentiality during asexual replication and gametocytogenesis. Upon inducible knockdown and targeted gene disruption, a reduced asexual parasite proliferation was detected for PfApiAT2 and PfApiAT4. Functional inactivation of individual PfApiATs targeted in this study had no effect on gametocyte development. Our data suggest that individual PfApiATs are partially redundant during asexual in vitro proliferation and fully redundant during gametocytogenesis of P. falciparum parasites.

IMPORTANCE Malaria parasites live and multiply inside cells. To facilitate their extremely fast intracellular proliferation, they hijack and transform their host cells. This also requires the active uptake of nutrients, such as amino acids, from the host cell and the surrounding environment through various membranes that are the consequence of the parasite’s intracellular lifestyle. In this paper, we focus on a family of putative amino acid transporters termed ApiAT. We show expression and localization of four transporters in the parasite plasma membrane of Plasmodium falciparum-infected erythrocytes that represent one interface of the pathogen to its host cell. We probed into the impact of functional inactivation of individual transporters on parasite growth in asexual and sexual blood stages of P. falciparum and reveal that only two of them show a modest but significant reduction in parasite proliferation but no impact on gametocytogenesis, pointing toward dispensability within this transporter family.

Gramtools enables multiscale variation analysis with genome graphs

Genome graphs allow very general representations of genetic variation; depending on the model and implementation, variation at different length-scales (single nucleotide polymorphisms (SNPs), structural variants) and on different sequence backgrounds can be incorporated with different levels of transparency. We implement a model which handles this multiscale variation and develop a JSON extension of VCF (jVCF) allowing for variant calls on multiple references, both implemented in our software gramtools. We find gramtools outperforms existing methods for genotyping SNPs overlapping large deletions in M. tuberculosis and is able to genotype on multiple alternate backgrounds in P. falciparum, revealing previously hidden recombination.

A nuclear enterprise: zooming in on nuclear organization and gene expression control in the African trypanosome

African trypanosomes are early divergent protozoan parasites responsible for high mortality and morbidity as well as a great economic burden among the world's poorest populations. Trypanosomes undergo antigenic variation in their mammalian hosts, a highly sophisticated immune evasion mechanism. Their nuclear organization and mechanisms for gene expression control present several conventional features but also a number of striking differences to the mammalian counterparts. Some of these unorthodox characteristics, such as lack of controlled transcription initiation or enhancer sequences, render their monogenic antigen transcription, which is critical for successful antigenic variation, even more enigmatic. Recent technological developments have advanced our understanding of nuclear organization and gene expression control in trypanosomes, opening novel research avenues. This review is focused on Trypanosoma brucei nuclear organization and how it impacts gene expression, with an emphasis on antigen expression. It highlights several dedicated sub-nuclear bodies that compartmentalize specific functions, whilst outlining similarities and differences to more complex eukaryotes. Notably, understanding the mechanisms underpinning antigen as well as general gene expression control is of great importance, as it might help designing effective control strategies against these organisms.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Evolution of transcriptional control of antigenic variation and virulence in human and ape malaria parasites

Background

The most severe form of human malaria is caused by the protozoan parasite Plasmodium falciparum. This unicellular organism is a member of a subgenus of Plasmodium called the Laverania that infects apes, with P. falciparum being the only member that infects humans. The exceptional virulence of this species to humans can be largely attributed to a family of variant surface antigens placed by the parasites onto the surface of infected red blood cells that mediate adherence to the vascular endothelium. These proteins are encoded by a large, multicopy gene family called var, with each var gene encoding a different form of the protein. By changing which var gene is expressed, parasites avoid immune recognition, a process called antigenic variation that underlies the chronic nature of malaria infections.

Results

Here we show that the common ancestor of the branch of the Laverania lineage that includes the human parasite underwent a remarkable change in the organization and structure of elements linked to the complex transcriptional regulation displayed by the var gene family. Unlike the other members of the Laverania, the clade that gave rise to P. falciparum evolved distinct subsets of var genes distinguishable by different upstream transcriptional regulatory regions that have been associated with different expression profiles and virulence properties. In addition, two uniquely conserved var genes that have been proposed to play a role in coordinating transcriptional switching similarly arose uniquely within this clade. We hypothesize that these changes originated at a time of dramatic climatic change on the African continent that is predicted to have led to significant changes in transmission dynamics, thus selecting for patterns of antigenic variation that enabled lengthier, more chronic infections.

Conclusions

These observations suggest that changes in transmission dynamics selected for significant alterations in the transcriptional regulatory mechanisms that mediate antigenic variation in the parasite lineage that includes P. falciparum. These changes likely underlie the chronic nature of these infections as well as their exceptional virulence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01872-z.