Characterization of the repertoire diversity of the Plasmodium falciparum stevor multigene family in laboratory and field isolates

Diamonds in the rif: Alignment-free comparative genomics analysis reveals strain-transcendent Plasmodium falciparum antigens amidst extensive genetic diversity

The repetitive interspersed family (rif) and subtelomeric variable open reading frames (stevor) are highly diverse multi-gene families in the malaria parasite Plasmodium falciparum. Embedded on the surface of infected erythrocytes, RIFIN and STEVOR proteins are involved in cytoadherence and immune evasion, but the extent of family-wide sequence diversity across strains has yet to be comprehensively investigated in light of improved resolution of the subtelomeric genome sequences. Using a k-mer frequency approach, we analyzed long-read genomic sequence data from 18 geographically diverse P. falciparum genome assemblies, including lab strains and clinical isolates. We hypothesized that k-mer sequence comparison can identify existing RIFIN and STEVOR subgroups, identify novel subgroups, and generate more robust and reliable estimates of family-wide sequence diversity. Full-length RIFIN and STEVOR proteins shared on average 49.5% and 61.1% amino acid k-mer similarity, respectively, which fell to 25.1% and 20% in the hypervariable regions alone. Despite this diversity, we identified 11 RIFINs and five STEVORs that were conserved across strains above expected thresholds. A subset of these strain-transcendent genes was similar and syntenic to genes in related Plasmodium species, suggesting an ancient origin. Additionally, in silico structural predictions from AlphaFold showed that three-dimensional structures of RIFIN receptor-binding regions were more conserved than their sequences suggested. Evolutionarily constrained RIFINs and STEVORs may have critical functions in parasite survival or pathogenesis. This study provides a framework for investigating diversity in highly variable multi-gene families and highlights the potential of strain-transcendent RIFIN and STEVOR proteins as vaccine candidates.

Antigenic Variation in Plasmodium falciparum

Plasmodium falciparum is the protozoan parasite that causes most malaria-associated morbidity and mortality in humans with over 500,000 deaths annually. The disease symptoms are associated with repeated cycles of invasion and asexual multiplication inside red blood cells of the parasite. Partial, non-sterile immunity to P. falciparum malaria develops only after repeated infections and continuous exposure. The successful evasion of the human immune system relies on the large repertoire of antigenically diverse parasite proteins displayed on the red blood cell surface and on the merozoite membrane where they are exposed to the human immune system. Expression switching of these polymorphic proteins between asexual parasite generations provides an efficient mechanism to adapt to the changing environment in the host and to maintain chronic infection. This chapter discusses antigenic diversity and variation in the malaria parasite and our current understanding of the molecular mechanisms that direct the expression of these proteins.

Assessing the genetic diversity of the vir genes in Indian Plasmodium vivax population

Variant surface antigens (VSAs) present on the surface of parasitized erythrocytes facilitate many Plasmodium spp. to escape the host immune system during infection. Multigene families coding for VSAs exist in several Plasmodium spp. and are located on telomeric and subtelomeric regions of the chromosomes. P. vivax genome also contains a multigene superfamily vir (variant interspersed repeats), present in the subtelomeric region with a possible role in immune evasion like the var gene in P. falciparum. Blood samples from 148 symptomatic malaria cases were collected from five different regions of India, viz. Mangalore, Rourkela, Goa, Delhi and Jabalpur. P. vivax isolates (74 single infections) were sequenced for four vir genes (viz. vir 27, vir 4, vir 12 and vir 21) and analyzed for the genetic variability existing in different populations of India. The results indicate that vir genes in different P. vivax populations in India are highly divergent both within and between the isolates. High levels of single nucleotide polymorphisms (SNPs) were observed attributing to the existing polymorphism for all the four vir genes studied across the population. Detailed knowledge of the genetic variation among the vir genes will help in understanding the evolutionary aspects of vir genes and may also provide basis for understanding the disease chronicity.

Expression and strain variation of the novel "small open reading frame" (smorf) multigene family in Babesia bovis

Small open reading frame (smorf) genes comprise the second largest Babesia bovis multigene family. All known 44 variant smorf genes are located in close chromosomal proximity to ves1 genes, which encode proteins that mediate cytoadhesion and contribute to immune evasion. In this study, we characterised the general topology of smorf genes and investigated the gene repertoire, transcriptional profile and SMORF expression in two distinct strains, T2Bo and Mo7. Sequence analysis using degenerate primers identified additional smorf genes in each strain and demonstrated that the smorf gene repertoire varies between strains, with conserved and unique genes in both. Smorf genes have multiple semi-conserved and variable blocks, and a large hypervariable insertion in 20 of the 44 genes defines two major branches of the family, termed smorf A and smorf B. A total of 32 smorf genes are simultaneously transcribed in T2Bo strain B. bovis merozoites obtained from deep brain tissue of an acutely infected animal. SMORF peptide-specific antiserum bound in immunoblots to multiple proteins with a range of sizes predicted by smorf genes, confirming translation of smorf gene products from these transcripts. These results indicate that the smorf multigene family is larger than previously described and demonstrate that smorf genes are expressed and are undergoing variation, both within strains and in a lineage-specific pattern independent of strain specificity. The function of these novel proteins is unknown.

Design of a variant surface antigen-supplemented microarray chip for whole transcriptome analysis of multiple Plasmodium falciparum cytoadherent strains, and identification of strain-transcendent rif and stevor genes

Background

The cytoadherence of Plasmodium falciparum is thought to be mediated by variant surface antigens (VSA), encoded by var, rif, stevor and pfmc-2tm genes. The last three families have rarely been studied in the context of cytoadherence. As most VSA genes are unique, the variability among sequences has impeded the functional study of VSA across different P. falciparum strains. However, many P. falciparum genomes have recently been sequenced, allowing the development of specific microarray probes for each VSA gene.

Methods

All VSA sequences from the HB3, Dd2 and IT/FCR3 genomes were extracted using HMMer software. Oligonucleotide probes were designed with OligoRankPick and added to the 3D7-based microarray chip. As a proof of concept, IT/R29 parasites were selected for and against rosette formation and the transcriptomes of isogenic rosetting and non-rosetting parasites were compared by microarray.

Results

From each parasite strain 50-56 var genes, 125-132 rif genes, 26-33 stevor genes and 3-8 pfmc-2tm genes were identified. Bioinformatic analysis of the new VSA sequences showed that 13 rif genes and five stevor genes were well-conserved across at least three strains (83-100% amino acid identity). The ability of the VSA-supplemented microarray chip to detect cytoadherence-related genes was assessed using P. falciparum clone IT/R29, in which rosetting is known to be mediated by PfEMP1 encoded by ITvar9. Whole transcriptome analysis showed that the most highly up-regulated gene in rosetting parasites was ITvar9 (19 to 429-fold up-regulated over six time points). Only one rif gene (IT4rifA_042) was up-regulated by more than four fold (five fold at 12 hours post-invasion), and no stevor or pfmc-2tm genes were up-regulated by more than two fold. 377 non-VSA genes were differentially expressed by three fold or more in rosetting parasites, although none was as markedly or consistently up-regulated as ITvar9.

Conclusions

Probes for the VSA of newly sequenced P. falciparum strains can be added to the 3D7-based microarray chip, allowing the analysis of the entire transcriptome of multiple strains. For the rosetting clone IT/R29, the striking transcriptional upregulation of ITvar9 was confirmed, and the data did not support the involvement of other VSA families in rosette formation.

Cloning of the repertoire of individual Plasmodium falciparum var genes using transformation associated recombination (TAR)

One of the major virulence factors of the malaria causing parasite is the Plasmodium falciparum encoded erythrocyte membrane protein 1 (PfEMP1). It is translocated to It the membrane of infected erythrocytes and expressed from approximately 60 var genes in a mutually exclusive manner. Switching of var genes allows the parasite to alter functional and antigenic properties of infected erythrocytes, to escape the immune defense and to establish chronic infections. We have developed an efficient method for isolating VAR genes from telomeric and other genome locations by adapting transformation-associated recombination (TAR) cloning, which can then be analyzed and sequenced. For this purpose, three plasmids each containing a homologous sequence representing the upstream regions of the group A, B, and C var genes and a sequence homologous to the conserved acidic terminal segment (ATS) of var genes were generated. Co-transfection with P. falciparum strain ITG2F6 genomic DNA in yeast cells yielded 200 TAR clones. The relative frequencies of clones from each group were not biased. Clones were screened by PCR, as well as Southern blotting, which revealed clones missed by PCR due to sequence mismatches with the primers. Selected clones were transformed into E. coli and further analyzed by RFLP and end sequencing. Physical analysis of 36 clones revealed 27 distinct types potentially representing 50% of the var gene repertoire. Three clones were selected for sequencing and assembled into single var gene containing contigs. This study demonstrates that it is possible to rapidly obtain the repertoire of var genes from P. falciparum within a single set of cloning experiments. This technique can be applied to individual isolates which will provide a detailed picture of the diversity of var genes in the field. This is a powerful tool to overcome the obstacles with cloning and assembly of multi-gene families by simultaneously cloning each member.

Small variant surface antigens and Plasmodium evasion of immunity

Antigenic variation at the Plasmodium-infected erythrocyte surface plays a critical role in malaria disease severity and host immune evasion. Our current understanding of the role of Plasmodium variant surface antigens in antigenic variation and immune evasion is largely limited to the extensive work carried out on the Plasmodium falciparum var gene family. Although homologues of var genes are not present in other malaria species, small variant gene families comprising the rif and stevor genes in P. falciparum and the pir genes in Plasmodium vivax, Plasmodium knowlesi and the rodent malaria Plasmodium chabaudi, Plasmodium berghei and Plasmodium yoelii also show features suggesting a role in antigenic variation and immune evasion. In this article, we highlight our current understanding of these variant antigens and provide insights on the mechanisms developed by malaria parasites to effectively avoid the host immune response and establish chronic infection.