The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees

An essential vesicular-trafficking phospholipase mediates neutral lipid synthesis and contributes to hemozoin formation in Plasmodium falciparum

Background

Plasmodium falciparum is the pathogen responsible for the most devastating form of human malaria. As it replicates asexually in the erythrocytes of its human host, the parasite feeds on haemoglobin uptaken from these cells. Heme, a toxic by-product of haemoglobin utilization by the parasite, is neutralized into inert hemozoin in the food vacuole of the parasite. Lipid homeostasis and phospholipid metabolism are crucial for this process, as well as for the parasite’s survival and propagation within the host. P. falciparum harbours a uniquely large family of phospholipases, which are suggested to play key roles in lipid metabolism and utilization.

Results

Here, we show that one of the parasite phospholipase (P. falciparum lysophospholipase, PfLPL1) plays an essential role in lipid homeostasis linked with the haemoglobin degradation and heme conversion pathway. Fluorescence tagging showed that the PfLPL1 in infected blood cells localizes to dynamic vesicular structures that traffic from the host-parasite interface at the parasite periphery, through the cytosol, to get incorporated into a large vesicular lipid rich body next to the food-vacuole. PfLPL1 is shown to harbour enzymatic activity to catabolize phospholipids, and its transient downregulation in the parasite caused a significant reduction of neutral lipids in the food vacuole-associated lipid bodies. This hindered the conversion of heme, originating from host haemoglobin, into the hemozoin, and disrupted the parasite development cycle and parasite growth. Detailed lipidomic analyses of inducible knock-down parasites deciphered the functional role of PfLPL1 in generation of neutral lipid through recycling of phospholipids. Further, exogenous fatty-acids were able to complement downregulation of PfLPL1 to rescue the parasite growth as well as restore hemozoin levels.

Conclusions

We found that the transient downregulation of PfLPL1 in the parasite disrupted lipid homeostasis and caused a reduction in neutral lipids essentially required for heme to hemozoin conversion. Our study suggests a crucial link between phospholipid catabolism and generation of neutral lipids (TAGs) with the host haemoglobin degradation pathway.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01042-z.

Gene tree species tree reconciliation with gene conversion

Gene tree/species tree reconciliation is a recent decisive progress in phylogenetic methods, accounting for the possible differences between gene histories and species histories. Reconciliation consists in explaining these differences by gene-scale events such as duplication, loss, transfer, which translates mathematically into a mapping between gene tree nodes and species tree nodes or branches. Gene conversion is a frequent and important evolutionary event, which results in the replacement of a gene by a copy of another from the same species and in the same gene tree. Including this event in reconciliation models has never been attempted because it introduces a dependency between lineages, and standard algorithms based on dynamic programming become ineffective. We propose here a novel mathematical framework including gene conversion as an evolutionary event in gene tree/species tree reconciliation. We describe a randomized algorithm that finds, in polynomial running time, a reconciliation minimizing the number of duplications, losses and conversions in the case when their weights are equal. We show that the space of optimal reconciliations includes an analog of the last common ancestor reconciliation, but is not limited to it. Our algorithm outputs any optimal reconciliation with a non-null probability. We argue that this study opens a research avenue on including gene conversion in reconciliation, and discuss its possible importance in biology.

Position of human blood group O(H) and phenotype-determining enzymes in growth and infectious disease

The human ABO(H) blood group phenotypes arise from the evolutionarily oldest genetic system found in primate populations. While the blood group antigen A is considered the ancestral primordial structure, under the selective pressure of life‐threatening diseases blood group O(H) came to dominate as the most frequently occurring blood group worldwide. Non‐O(H) phenotypes demonstrate impaired formation of adaptive and innate immunoglobulin specificities due to clonal selection and phenotype formation in plasma proteins. Compared with individuals with blood group O(H), blood group A individuals not only have a significantly higher risk of developing certain types of cancer but also exhibit high susceptibility to malaria tropica or infection by Plasmodium falciparum. The phenotype‐determining blood group A glycotransferase(s), which affect the levels of anti‐A/Tn cross‐reactive immunoglobulins in phenotypic glycosidic accommodation, might also mediate adhesion and entry of the parasite to host cells via trans‐species O‐GalNAc glycosylation of abundantly expressed serine residues that arise throughout the parasite's life cycle, while excluding the possibility of antibody formation against the resulting hybrid Tn antigen. In contrast, human blood group O(H), lacking this enzyme, is indicated to confer a survival advantage regarding the overall risk of developing cancer, and individuals with this blood group rarely develop life‐threatening infections involving evolutionarily selective malaria strains.

Plasmodium falciparum SERA5 plays a non-enzymatic role in the malarial asexual blood-stage lifecycle

The malaria parasite Plasmodium falciparum replicates in an intraerythrocytic parasitophorous vacuole (PV). The most abundant P. falciparum PV protein, called SERA5, is essential in blood stages and possesses a papain-like domain, prompting speculation that it functions as a proteolytic enzyme. Unusually however, SERA5 possesses a Ser residue (Ser596) at the position of the canonical catalytic Cys of papain-like proteases, and the function of SERA5 or whether it performs an enzymatic role is unknown. In this study, we failed to detect proteolytic activity associated with the Ser596-containing parasite-derived or recombinant protein. However, substitution of Ser596 with a Cys residue produced an active recombinant enzyme with characteristics of a cysteine protease, demonstrating that SERA5 can bind peptides. Using targeted homologous recombination in P. falciparum, we substituted Ser596 with Ala with no phenotypic consequences, proving that SERA5 does not perform an essential enzymatic role in the parasite. We could also replace an internal segment of SERA5 with an affinity-purification tag. In contrast, using almost identical targeting constructs, we could not truncate or C-terminally tag the SERA5 gene, or replace Ser596 with a bulky Arg residue. Our findings show that SERA5 plays an indispensable but non-enzymatic role in the P. falciparum blood-stage life cycle.

A comparative proteomic analysis of the simple amino acid repeat distributions in Plasmodia reveals lineage specific amino acid selection

Background

Microsatellites have been used extensively in the field of comparative genomics. By studying microsatellites in coding regions we have a simple model of how genotypic changes undergo selection as they are directly expressed in the phenotype as altered proteins. The simplest of these tandem repeats in coding regions are the tri-nucleotide repeats which produce a repeat of a single amino acid when translated into proteins. Tri-nucleotide repeats are often disease associated, and are also known to be unstable to both expansion and contraction. This makes them sensitive markers for studying proteome evolution, in closely related species.

Results

The evolutionary history of the family of malarial causing parasites Plasmodia is complex because of the life-cycle of the organism, where it interacts with a number of different hosts and goes through a series of tissue specific stages. This study shows that the divergence between the primate and rodent malarial parasites has resulted in a lineage specific change in the simple amino acid repeat distribution that is correlated to A–T content. The paper also shows that this altered use of amino acids in SAARs is consistent with the repeat distributions being under selective pressure.

Conclusions

The study shows that simple amino acid repeat distributions can be used to group related species and to examine their phylogenetic relationships. This study also shows that an outgroup species with a similar A–T content can be distinguished based only on the amino acid usage in repeats, and suggest that this might be a useful feature for proteome clustering. The lineage specific use of amino acids in repeat regions suggests that comparative studies of SAAR distributions between proteomes gives an insight into the mechanisms of expansion and the selective pressures acting on the organism.

Structural insights into the protease-like antigen Plasmodium falciparum SERA5 and its noncanonical active-site serine

The sera genes of the malaria-causing parasite Plasmodium encode a family of unique proteins that are maximally expressed at the time of egress of parasites from infected red blood cells. These multi-domain proteins are unique, containing a central papain-like cysteine-protease fragment enclosed between the disulfide-linked N- and C-terminal domains. However, the central fragment of several members of this family, including serine repeat antigen 5 (SERA5), contains a serine (S596) in place of the active-site cysteine. Here we report the crystal structure of the central protease-like domain of Plasmodium falciparum SERA5, revealing a number of anomalies in addition to the putative nucleophilic serine: (1) the structure of the putative active site is not conducive to binding substrate in the canonical cysteine-protease manner; (2) the side chain of D594 restricts access of substrate to the putative active site; and (3) the S(2) specificity pocket is occupied by the side chain of Y735, reducing this site to a small depression on the protein surface. Attempts to determine the structure in complex with known inhibitors were not successful. Thus, despite having revealed its structure, the function of the catalytic domain of SERA5 remains an enigma.

Reconciliation with non-binary species trees

Reconciliation extracts information from the topological incongruence between gene and species trees to infer duplications and losses in the history of a gene family. The inferred duplication-loss histories provide valuable information for a broad range of biological applications, including ortholog identification, estimating gene duplication times, and rooting and correcting gene trees. While reconciliation for binary trees is a tractable and well studied problem, there are no algorithms for reconciliation with non-binary species trees. Yet a striking proportion of species trees are non-binary. For example, 64% of branch points in the NCBI taxonomy have three or more children. When applied to non-binary species trees, current algorithms overestimate the number of duplications because they cannot distinguish between duplication and incomplete lineage sorting. We present the first algorithms for reconciling binary gene trees with non-binary species trees under a duplication-loss parsimony model. Our algorithms utilize an efficient mapping from gene to species trees to infer the minimum number of duplications in O(|V(G) | x (k(S) + h(S))) time, where |V(G)| is the number of nodes in the gene tree, h(S) is the height of the species tree and k(S) is the size of its largest polytomy. We present a dynamic programming algorithm which also minimizes the total number of losses. Although this algorithm is exponential in the size of the largest polytomy, it performs well in practice for polytomies with outdegree of 12 or less. We also present a heuristic which estimates the minimal number of losses in polynomial time. In empirical tests, this algorithm finds an optimal loss history 99% of the time. Our algorithms have been implemented in NOTUNG, a robust, production quality, tree-fitting program, which provides a graphical user interface for exploratory analysis and also supports automated, high-throughput analysis of large data sets.

Inhibition of malaria parasite development by a cyclic peptide that targets the vital parasite protein SERA5

The serine repeat antigen (SERA) proteins of the malaria parasites Plasmodium spp. contain a putative enzyme domain similar to that of papain family cysteine proteases. In Plasmodium falciparum parasites, more than half of the SERA family proteins, including the most abundantly expressed form, SERA5, have a cysteine-to-serine substitution within the putative catalytic triad of the active site. Although SERA5 is required for blood-stage parasite survival, the occurrence of a noncanonical catalytic triad casts doubt on the importance of the enzyme domain in this function. We used phage display to identify a small (14-residue) disulfide-bonded cyclic peptide (SBP1) that targets the enzyme domain of SERA5. Biochemical characterization of the interaction shows that it is dependent on the conformation of both the peptide and protein. Addition of this peptide to parasite cultures compromised development of late-stage parasites compared to that of control parasites or those incubated with equivalent amounts of the carboxymethylated peptide. This effect was similar in two different strains of P. falciparum as well as in a transgenic strain where the gene encoding the related serine-type parasitophorous vacuole protein SERA4 was deleted. In compromised parasites, the SBP1 peptide crosses both the erythrocyte and parasitophorous vacuole membranes and accumulates within the parasitophorous vacuole. In addition, both SBP1 and SERA5 were identified in the parasite cytosol, indicating that the plasma membrane of the parasite was compromised as a result of SBP1 treatment. These data implicate an important role for SERA5 in the regulation of the intraerythrocytic development of late-stage parasites and as a target for drug development.

Malarial proteases and host cell egress: an 'emerging' cascade

Malaria is a scourge of large swathes of the globe, stressing the need for a continuing effort to better understand the biology of its aetiological agent. Like all pathogens of the phylum Apicomplexa, the malaria parasite spends part of its life inside a host cell or cyst. It eventually needs to escape (egress) from this protective environment to progress through its life cycle. Egress of Plasmodium blood-stage merozoites, liver-stage merozoites and mosquito midgut sporozoites relies on protease activity, so the enzymes involved have potential as antimalarial drug targets. This review examines the role of parasite proteases in egress, in the light of current knowledge of the mechanics of the process. Proteases implicated in egress include the cytoskeleton-degrading malarial proteases falcipain-2 and plasmepsin II, plus a family of putative papain-like proteases called SERA. Recent revelations have shown that activation of the SERA proteases may be triggered by regulated secretion of a subtilisin-like serine protease called SUB1. These findings are discussed in the context of the potential for development of new chemotherapeutics targeting this stage in the parasite's life cycle.

Expression and processing of Plasmodium berghei SERA3 during liver stages

Cysteine proteases mediate liberation of Plasmodium berghei merozoites from infected hepatocytes. In an attempt to identify the responsible parasite proteases, we screened the genome of P. berghei for cysteine protease-encoding genes. RT-PCR analyses revealed that transcription of four out of five P. berghei serine repeat antigen (PbSERA) genes was strongly upregulated in late liver stages briefly before the parasitophorous vacuole membrane ruptured to release merozoites into the host cell cytoplasm, suggesting a role of PbSERA proteases in these processes. In order to characterize PbSERA3 processing, we raised an antiserum against a non-conserved region of the protein and generated a transgenic P. berghei strain expressing a TAP-tagged PbSERA3 under the control of the endogenous promoter. Immunofluorescence assays revealed that PbSERA3 leaks into the host cell cytoplasm during merozoite development, where it might contribute to host cell death or activate host cell proteases that execute cell death. Importantly, processed PbSERA3 has been detected by Western blot analysis in cell extracts of schizont-infected cells and merozoite-infected detached hepatic cells.

Evidence for a common role for the serine-type Plasmodium falciparum serine repeat antigen proteases: implications for vaccine and drug design

Serine repeat antigens (SERAs) are a family of secreted "cysteine-like" proteases of Plasmodium parasites. Several SERAs possess an atypical active-site serine residue in place of the canonical cysteine. The human malaria parasite Plasmodium falciparum possesses six "serine-type" (SERA1 to SERA5 and SERA9) and three "cysteine-type" (SERA6 to SERA8) SERAs. Here, we investigate the importance of the serine-type SERAs to blood-stage parasite development and examine the extent of functional redundancy among this group. We attempted to knock out the four P. falciparum serine-type SERA genes that have not been disrupted previously. SERA1, SERA4, and SERA9 knockout lines were generated, while only SERA5, the most strongly expressed member of the SERA family, remained refractory to genetic deletion. Interestingly, we discovered that while SERA4-null parasites completed the blood-stage cycle normally, they exhibited a twofold increase in the level of SERA5 mRNA. The inability to disrupt SERA5 and the apparent compensatory increase in SERA5 expression in response to the deletion of SERA4 provides evidence for an important blood-stage function for the serine-type SERAs and supports the notion of functional redundancy among this group. Such redundancy is consistent with our phylogenetic analysis, which reveals a monophyletic grouping of the serine-type SERAs across the genus Plasmodium and a predominance of postspeciation expansion. While SERA5 is to some extent further validated as a target for vaccine and drug development, our data suggest that the expression level of other serine-type SERAs is the only barrier to escape from anti-SERA5-specific interventions.

Structural evolution of nrDNA ITS in Pinaceae and its phylogenetic implications

Nuclear ribosomal DNA (nrDNA) has been considered as an important tool for inferring phylogenetic relationships at many taxonomic levels. In comparison with its fast concerted evolution in angiosperms, nrDNA is symbolized by slow concerted evolution and substantial ITS region length variation in gymnosperms, particularly in Pinaceae. Here we studied structure characteristics, including subrepeat composition, size, GC content and secondary structure, of nrDNA ITS regions of all Pinaceae genera. The results showed that the ITS regions of all taxa studied contained subrepeat units, ranging from 2 to 9 in number, and these units could be divided into two types, longer subrepeat (LSR) without the motif (5'-GGCCACCCTAGTC) and shorter subrepeat (SSR) with the motif. Phylogenetic analyses indicate that the homology of some SSRs still can be recognized, providing important informations for the evolutionary history of nrDNA ITS and phylogeny of Pinaceae. In particular, the adjacent tandem SSRs are not more closely related to one another than they are to remote SSRs in some genera, which may imply that multiple structure variations such as recombination have occurred in the ITS1 region of these groups. This study also found that GC content in the ITS1 region is relevant to its sequence length and subrepeat number, and could provide some phylogenetic information, especially supporting the close relationships among Picea, Pinus, and Cathaya. Moreover, several characteristics of the secondary structure of Pinaceae ITS1 were found as follows: (1) the structure is dominated by several extended hairpins; (2) the configuration complexity is positively correlated with subrepeat number; (3) paired subrepeats often partially overlap at the conserved motif (5'-GGCCACCCTAGTC), and form a long stem, while other subrepeats fold onto itself, leaving part of the conserved motif exposed in hairpin loops.

Phylogeny and evolution of the SERA multigene family in the genus Plasmodium

The serine repeat antigen gene family of Plasmodium falciparum (Pf-SERA) consists of nine gene members. By sequence similarity search, 45 genes were identified to be homologous to the Pf-SERA genes in the ongoing seven Plasmodium genome sequencing project databases for the species: P. reichenowi, P. vivax, P. knowlesi, P. yoelii, P. berghei, P. chabaudi, and P. gallinaceum. In combination with additional PCR-based sequencing, we found that almost all SERA genes in each species were aligned in a tandem cluster and sandwiched between two conserved hypothetical protein genes, except for P. reichenowi, which could not be confirmed. The minimum and maximum numbers of clustered genes were 2 and 12 for P. gallinaceum and P. vivax, respectively. The best tree of the maximum likelihood analysis demonstrated that all Plasmodium SERA homologues, except for SERA1 of P. gallinaceum (Pg-SERA1), can be classified into four groups, represented by Pf-SERA5, Pf-SERA6, Pf-SERA7, and Pf-SERA8. Genes in the Pf-SERA8 group, although highly divergent and distantly related to the sequences of other groups, were not pseudogenes. P. berghei SERA5, the counterpart of Pf-SERA8, was expressed in the mosquito stage. P. gallinaceum lacks the orthologues to Pf-SERA5, Pf-SERA6, and Pf-SERA7, suggesting that P. gallinaceum diverged from a common ancestor of all eight Plasmodium species examined before gene duplication(s) occurred to generate these paralogous groups. Here, we reveal an evolutionary trail of SERA gene cluster in the genus Plasmodium and discuss a phylogeny of Plasmodium species from the viewpoint of the evolution of a multigene family.

Loss of the flagellum happened only once in the fungal lineage: phylogenetic structure of kingdom Fungi inferred from RNA polymerase II subunit genes

Background

At present, there is not a widely accepted consensus view regarding the phylogenetic structure of kingdom Fungi although two major phyla, Ascomycota and Basidiomycota, are clearly delineated. Regarding the lower fungi, Zygomycota and Chytridiomycota, a variety of proposals have been advanced. Microsporidia may or may not be fungi; the Glomales (vesicular-arbuscular mycorrhizal fungi) may or may not constitute a fifth fungal phylum, and the loss of the flagellum may have occurred either once or multiple times during fungal evolution. All of these issues are capable of being resolved by a molecular phylogenetic analysis which achieves strong statistical support for major branches. To date, no fungal phylogeny based upon molecular characters has satisfied this criterion.

Results

Using the translated amino acid sequences of the RPB1 and RPB2 genes, we have inferred a fungal phylogeny that consists largely of well-supported monophyletic phyla. Our major results, each with significant statistical support, are: (1) Microsporidia are sister to kingdom Fungi and are not members of Zygomycota; that is, Microsporidia and fungi originated from a common ancestor. (2) Chytridiomycota, the only fungal phylum having a developmental stage with a flagellum, is paraphyletic and is the basal lineage. (3) Zygomycota is monophyletic based upon sampling of Trichomycetes, Zygomycetes, and Glomales. (4) Zygomycota, Basidiomycota, and Ascomycota form a monophyletic group separate from Chytridiomycota. (5) Basidiomycota and Ascomycota are monophyletic sister groups.

Conclusion

In general, this paper highlights the evolutionary position and significance of the lower fungi (Zygomycota and Chytridiomycota). Our results suggest that loss of the flagellum happened only once during early stages of fungal evolution; consequently, the majority of fungi, unlike plants and animals, are nonflagellated. The phylogeny we infer from gene sequences is the first one that is congruent with the widely accepted morphology-based classification of Fungi. We find that, contrary to what has been published elsewhere, the four morphologically defined phyla (Ascomycota, Basidiomycota, Zygomycota and Chytridiomycota) do not overlap with one another. Microsporidia are not included within kingdom Fungi; rather they are a sister-group to the Fungi. Our study demonstrates the applicability of protein sequences from large, slowly-evolving genes to the derivation of well-resolved and highly supported phylogenies across long evolutionary distances.