VEuPathDB: the eukaryotic pathogen, vector and host bioinformatics resource center in 2023

The Eukaryotic Pathogen, Vector and Host Informatics Resource (VEuPathDB, https://veupathdb.org) is a Bioinformatics Resource Center funded by the National Institutes of Health with additional funding from the Wellcome Trust. VEuPathDB supports >600 organisms that comprise invertebrate vectors, eukaryotic pathogens (protists and fungi) and relevant free-living or non-pathogenic species or hosts. Since 2004, VEuPathDB has analyzed omics data from the public domain using contemporary bioinformatic workflows, including orthology predictions via OrthoMCL, and integrated the analysis results with analysis tools, visualizations, and advanced search capabilities. The unique data mining platform coupled with >3000 pre-analyzed data sets facilitates the exploration of pertinent omics data in support of hypothesis driven research. Comparisons are easily made across data sets, data types and organisms. A Galaxy workspace offers the opportunity for the analysis of private large-scale datasets and for porting to VEuPathDB for comparisons with integrated data. The MapVEu tool provides a platform for exploration of spatially resolved data such as vector surveillance and insecticide resistance monitoring. To address the growing body of omics data and advances in laboratory techniques, VEuPathDB has added several new data types, searches and features, improved the Galaxy workspace environment, redesigned the MapVEu interface and updated the infrastructure to accommodate these changes.

Malaria disrupts the rhesus macaque gut microbiome

Previous studies have suggested that a relationship exists between severity and transmissibility of malaria and variations in the gut microbiome, yet only limited information exists on the temporal dynamics of the gut microbial community during a malarial infection. Here, using a rhesus macaque model of relapsing malaria, we investigate how malaria affects the gut microbiome. In this study, we performed 16S sequencing on DNA isolated from rectal swabs of rhesus macaques over the course of an experimental malarial infection with Plasmodium cynomolgi and analyzed gut bacterial taxa abundance across primary and relapsing infections. We also performed metabolomics on blood plasma from the animals at the same timepoints and investigated changes in metabolic pathways over time. Members of Proteobacteria (family Helicobacteraceae) increased dramatically in relative abundance in the animal's gut microbiome during peak infection while Firmicutes (family Lactobacillaceae and Ruminococcaceae), Bacteroidetes (family Prevotellaceae) and Spirochaetes amongst others decreased compared to baseline levels. Alpha diversity metrics indicated decreased microbiome diversity at the peak of parasitemia, followed by restoration of diversity post-treatment. Comparison with healthy subjects suggested that the rectal microbiome during acute malaria is enriched with commensal bacteria typically found in the healthy animal's mucosa. Significant changes in the tryptophan-kynurenine immunomodulatory pathway were detected at peak infection with P. cynomolgi, a finding that has been described previously in the context of P. vivax infections in humans. During relapses, which have been shown to be associated with less inflammation and clinical severity, we observed minimal disruption to the gut microbiome, despite parasites being present. Altogether, these data suggest that the metabolic shift occurring during acute infection is associated with a concomitant shift in the gut microbiome, which is reversed post-treatment.

MaHPIC malaria systems biology data from Plasmodium cynomolgi sporozoite longitudinal infections in macaques

Plasmodium cynomolgi causes zoonotic malarial infections in Southeast Asia and this parasite species is important as a model for Plasmodium vivax and Plasmodium ovale. Each of these species produces hypnozoites in the liver, which can cause relapsing infections in the blood. Here we present methods and data generated from iterative longitudinal systems biology infection experiments designed and performed by the Malaria Host-Pathogen Interaction Center (MaHPIC) to delve deeper into the biology, pathogenesis, and immune responses of P. cynomolgi in the Macaca mulatta host. Infections were initiated by sporozoite inoculation. Blood and bone marrow samples were collected at defined timepoints for biological and computational experiments and integrative analyses revolving around primary illness, relapse illness, and subsequent disease and immune response patterns. Parasitological, clinical, haematological, immune response, and -omic datasets (transcriptomics, proteomics, metabolomics, and lipidomics) including metadata and computational results have been deposited in public repositories. The scope and depth of these datasets are unprecedented in studies of malaria, and they are projected to be a F.A.I.R., reliable data resource for decades.

Plasmodium knowlesi Cytoadhesion Involves SICA Variant Proteins

Plasmodium knowlesi poses a health threat throughout Southeast Asian communities and currently causes most cases of malaria in Malaysia. This zoonotic parasite species has been studied in Macaca mulatta (rhesus monkeys) as a model for severe malarial infections, chronicity, and antigenic variation. The phenomenon of Plasmodium antigenic variation was first recognized during rhesus monkey infections. Plasmodium-encoded variant proteins were first discovered in this species and found to be expressed at the surface of infected erythrocytes, and then named the Schizont-Infected Cell Agglutination (SICA) antigens. SICA expression was shown to be spleen dependent, as SICA expression is lost after P. knowlesi is passaged in splenectomized rhesus. Here we present data from longitudinal P. knowlesi infections in rhesus with the most comprehensive analysis to date of clinical parameters and infected red blood cell sequestration in the vasculature of tissues from 22 organs. Based on the histopathological analysis of 22 tissue types from 11 rhesus monkeys, we show a comparative distribution of parasitized erythrocytes and the degree of margination of the infected erythrocytes with the endothelium. Interestingly, there was a significantly higher burden of parasites in the gastrointestinal tissues, and extensive margination of the parasites along the endothelium, which may help explain gastrointestinal symptoms frequently reported by patients with P. knowlesi malarial infections. Moreover, this margination was not observed in splenectomized rhesus that were infected with parasites not expressing the SICA proteins. This work provides data that directly supports the view that a subpopulation of P. knowlesi parasites cytoadheres and sequesters, likely via SICA variant antigens acting as ligands. This process is akin to the cytoadhesive function of the related variant antigen proteins, namely Erythrocyte Membrane Protein-1, expressed by Plasmodium falciparum.

Clinical recovery of Macaca fascicularis infected with Plasmodium knowlesi

BACKGROUND

Kra monkeys (Macaca fascicularis), a natural host of Plasmodium knowlesi, control parasitaemia caused by this parasite species and escape death without treatment. Knowledge of the disease progression and resilience in kra monkeys will aid the effective use of this species to study mechanisms of resilience to malaria. This longitudinal study aimed to define clinical, physiological and pathological changes in kra monkeys infected with P. knowlesi, which could explain their resilient phenotype.

METHODS

Kra monkeys (n = 15, male, young adults) were infected intravenously with cryopreserved P. knowlesi sporozoites and the resulting parasitaemias were monitored daily. Complete blood counts, reticulocyte counts, blood chemistry and physiological telemetry data (n = 7) were acquired as described prior to infection to establish baseline values and then daily after inoculation for up to 50 days. Bone marrow aspirates, plasma samples, and 22 tissue samples were collected at specific time points to evaluate longitudinal clinical, physiological and pathological effects of P. knowlesi infections during acute and chronic infections.

RESULTS

As expected, the kra monkeys controlled acute infections and remained with low-level, persistent parasitaemias without anti-malarial intervention. Unexpectedly, early in the infection, fevers developed, which ultimately returned to baseline, as well as mild to moderate thrombocytopenia, and moderate to severe anaemia. Mathematical modelling and the reticulocyte production index indicated that the anaemia was largely due to the removal of uninfected erythrocytes and not impaired production of erythrocytes. Mild tissue damage was observed, and tissue parasite load was associated with tissue damage even though parasite accumulation in the tissues was generally low.

CONCLUSIONS

Kra monkeys experimentally infected with P. knowlesi sporozoites presented with multiple clinical signs of malaria that varied in severity among individuals. Overall, the animals shared common mechanisms of resilience characterized by controlling parasitaemia 3-5 days after patency, and controlling fever, coupled with physiological and bone marrow responses to compensate for anaemia. Together, these responses likely minimized tissue damage while supporting the establishment of chronic infections, which may be important for transmission in natural endemic settings. These results provide new foundational insights into malaria pathogenesis and resilience in kra monkeys, which may improve understanding of human infections.

Distinct amino acid and lipid perturbations characterize acute versus chronic malaria

Chronic malaria is a major public health problem and significant challenge for disease eradication efforts. Despite its importance, the biological factors underpinning chronic malaria are not fully understood. Recent studies have shown that host metabolic state can influence malaria pathogenesis and transmission, but its role in chronicity is not known. Here, with the goal of identifying distinct modifications in the metabolite profiles of acute versus chronic malaria, metabolomics was performed on plasma from Plasmodium-infected humans and nonhuman primates with a range of parasitemias and clinical signs. In rhesus macaques infected with Plasmodium coatneyi, significant alterations in amines, carnitines, and lipids were detected during a high parasitemic acute phase and many of these reverted to baseline levels once a low parasitemic chronic phase was established. Plasmodium gene expression, studied in parallel in the macaques, revealed transcriptional changes in amine, fatty acid, lipid and energy metabolism genes, as well as variant antigen genes. Furthermore, a common set of amines, carnitines, and lipids distinguished acute from chronic malaria in plasma from human Plasmodium falciparum cases. In summary, distinct host-parasite metabolic environments have been uncovered that characterize acute versus chronic malaria, providing insights into the underlying host-parasite biology of malaria disease progression.

Assembly of highly repetitive genomes using short reads: the genome of discrete typing unit III Trypanosoma cruzi strain 231

Next-generation sequencing (NGS) methods are low-cost high-throughput technologies that produce thousands to millions of sequence reads. Despite the high number of raw sequence reads, their short length, relative to Sanger, PacBio or Nanopore reads, complicates the assembly of genomic repeats. Many genome tools are available, but the assembly of highly repetitive genome sequences using only NGS short reads remains challenging. Genome assembly of organisms responsible for important neglected diseases such as Trypanosoma cruzi, the aetiological agent of Chagas disease, is known to be challenging because of their repetitive nature. Only three of six recognized discrete typing units (DTUs) of the parasite have their draft genomes published and therefore genome evolution analyses in the taxon are limited. In this study, we developed a computational workflow to assemble highly repetitive genomes via a combination of de novo and reference-based assembly strategies to better overcome the intrinsic limitations of each, based on Illumina reads. The highly repetitive genome of the human-infecting parasite T. cruzi 231 strain was used as a test subject. The combined-assembly approach shown in this study benefits from the reference-based assembly ability to resolve highly repetitive sequences and from the de novo capacity to assemble genome-specific regions, improving the quality of the assembly. The acceptable confidence obtained by analyzing our results showed that our combined approach is an attractive option to assemble highly repetitive genomes with NGS short reads. Phylogenomic analysis including the 231 strain, the first representative of DTU III whose genome was sequenced, was also performed and provides new insights into T. cruzi genome evolution.

Integrative analysis associates monocytes with insufficient erythropoiesis during acute Plasmodium cynomolgi malaria in rhesus macaques

Background

Mild to severe anaemia is a common complication of malaria that is caused in part by insufficient erythropoiesis in the bone marrow. This study used systems biology to evaluate the transcriptional and alterations in cell populations in the bone marrow during Plasmodium cynomolgi infection of rhesus macaques (a model of Plasmodium vivax malaria) that may affect erythropoiesis.

Results

An appropriate erythropoietic response did not occur to compensate for anaemia during acute cynomolgi malaria despite an increase in erythropoietin levels. During this period, there were significant perturbations in the bone marrow transcriptome. In contrast, relapses did not induce anaemia and minimal changes in the bone marrow transcriptome were detected. The differentially expressed genes during acute infection were primarily related to ongoing inflammatory responses with significant contributions from Type I and Type II Interferon transcriptional signatures. These were associated with increased frequency of intermediate and non-classical monocytes. Recruitment and/or expansion of these populations was correlated with a decrease in the erythroid progenitor population during acute infection, suggesting that monocyte-associated inflammation may have contributed to anaemia. The decrease in erythroid progenitors was associated with downregulation of genes regulated by GATA1 and GATA2, two master regulators of erythropoiesis, providing a potential molecular basis for these findings.

Conclusions

These data suggest the possibility that malarial anaemia may be driven by monocyte-associated disruption of GATA1/GATA2 function in erythroid progenitors resulting in insufficient erythropoiesis during acute infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-2029-z) contains supplementary material, which is available to authorized users.

An evolutionary approach to identify potentially protective B cell epitopes involved in naturally acquired immunity to malaria and the role of EBA-175 in protection amongst denizens of Bolifamba, Cameroon

Background

The search for a vaccine against malaria caused by Plasmodium falciparum has lasted for more than 100 years, with considerable progress in the identification of a number of vaccine candidates. The post-genomic era offers new opportunities for an expedited search using rational vaccine design and prioritization of key B-cell epitopes involved in natural acquired immunity.

Methods

Malaria vaccine candidate genes that have reached clinical trial were searched on an evolutionary relationship tree, to determine their level of lineage-specificity. Ten other genes with similar protein features and level of lineage specificity to the vaccine candidates were randomly selected, and computationally evaluated for the presence of B-cell epitopes. The protein fragment with maximum probability of putative epitopes were synthesized and used in an ELISA experiment to determine the presence of antibodies to these peptides, in the serum of malaria patients and healthy malaria uninfected inhabitants from a malaria endemic region (Bolifamba), alongside with a vaccine candidate EBA-175.

Results

Two peptide fragments of 25 and 30 amino acid length from PF3D7_1233400 and PF3D7_1437500 respectively, coded as PF4-123 and PF4-143 were shown to contain B-cell epitope(s). Total IgG antibodies to these peptides were not significantly different between sick and healthy participants, but cytophilic antibodies to these peptides were significantly higher in healthy participants (p < 0.03). Total IgG to the vaccine candidate EBA-175 was significantly higher in sick participants than in healthy participants, likewise cytophilic antibodies (p < 0.04). Antibodies to the peptides PF4-123 and PF4-143 correlated negatively (p = 0.025 and 0.008 and r = −0.291 and −0.345, respectively) to parasite load. Total IgG antibodies to EBA-175 showed a negative correlation to parasite load (r = −0.144), which was not significant (p = 0.276). Duration of stay in Bolifamba also negatively correlated with parasite load (p = 0.026, r = −0.419) and total IgG to PF4-143 was significantly associated with prolonged duration of stay in the locality of Bolifamba, Cameroon (p = 0.006, r = 0.361).

Conclusions

The present study has identified two genes PF3D7_1233400 and PF3D7_1437500 containing peptide fragment (PF4-123 and PF4-143) with B-cell epitopes that are correlated with naturally acquired immunity to malaria. A pipeline has been developed for rapid identification of other B-cell epitopes involved in naturally acquired immunity.

The Genome of Haemoproteus tartakovskyi and Its Relationship to Human Malaria Parasites

The phylogenetic relationships among hemosporidian parasites, including the origin of Plasmodium falciparum, the most virulent malaria parasite of humans, have been heavily debated for decades. Studies based on multiple-gene sequences have helped settle many of these controversial phylogenetic issues. However, denser taxon sampling and genome-wide analyses are needed to confidently resolve the evolutionay relationships among hemosporidian parasites. Genome sequences of several Plasmodium parasites are available but only for species infecting primates and rodents. To root the phylogenetic tree of Plasmodium, genomic data from related parasites of birds or reptiles are required. Here, we use a novel approach to isolate parasite DNA from microgametes and describe the first genome of a bird parasite in the sister genus to Plasmodium, Haemoproteus tartakovskyi. Similar to Plasmodium parasites, H. tartakovskyi has a small genome (23.2 Mb, 5,990 genes) and a GC content (25.4%) closer to P. falciparum (19.3%) than to Plasmodium vivax (42.3%). Combined with novel transcriptome sequences of the bird parasite Plasmodium ashfordi, our phylogenomic analyses of 1,302 orthologous genes demonstrate that mammalian-infecting malaria parasites are monophyletic, thus rejecting the repeatedly proposed hypothesis that the ancestor of Laverania parasites originated from a secondary host shift from birds to humans. Genes and genomic features previously found to be shared between P. falciparum and bird malaria parasites, but absent in other mammal malaria parasites, are therefore signatures of maintained ancestral states. We foresee that the genome of H. tartakovskyi will open new directions for comparative evolutionary analyses of malarial adaptive traits.

Comparative transcriptomics and metabolomics in a rhesus macaque drug administration study

We describe a multi-omic approach to understanding the effects that the anti-malarial drug pyrimethamine has on immune physiology in rhesus macaques (Macaca mulatta). Whole blood and bone marrow (BM) RNA-Seq and plasma metabolome profiles (each with over 15,000 features) have been generated for five naïve individuals at up to seven timepoints before, during and after three rounds of drug administration. Linear modeling and Bayesian network analyses are both considered, alongside investigations of the impact of statistical modeling strategies on biological inference. Individual macaques were found to be a major source of variance for both omic data types, and factoring individuals into subsequent modeling increases power to detect temporal effects. A major component of the whole blood transcriptome follows the BM with a time-delay, while other components of variation are unique to each compartment. We demonstrate that pyrimethamine administration does impact both compartments throughout the experiment, but very limited perturbation of transcript or metabolite abundance was observed following each round of drug exposure. New insights into the mode of action of the drug are presented in the context of pyrimethamine's predicted effect on suppression of cell division and metabolism in the immune system.

The dynamics of LTR retrotransposon accumulation across 25 million years of panicoid grass evolution

Sample sequence analysis was employed to investigate the repetitive DNAs that were most responsible for the evolved variation in genome content across seven panicoid grasses with >5-fold variation in genome size and different histories of polyploidy. In all cases, the most abundant repeats were LTR retrotransposons, but the particular families that had become dominant were found to be different in the Pennisetum, Saccharum, Sorghum and Zea lineages. One element family, Huck, has been very active in all of the studied species over the last few million years. This suggests the transmittal of an active or quiescent autonomous set of Huck elements to this lineage at the founding of the panicoids. Similarly, independent recent activity of Ji and Opie elements in Zea and of Leviathan elements in Sorghum and Saccharum species suggests that members of these families with exceptional activation potential were present in the genome(s) of the founders of these lineages. In a detailed analysis of the Zea lineage, the combined action of several families of LTR retrotransposons were observed to have approximately doubled the genome size of Zea luxurians relative to Zea mays and Zea diploperennis in just the last few million years. One of the LTR retrotransposon amplification bursts in Zea may have been initiated by polyploidy, but the great majority of transposable element activations are not. Instead, the results suggest random activation of a few or many LTR retrotransposons families in particular lineages over evolutionary time, with some families especially prone to future activation and hyper-amplification.

Real-time loop-mediated isothermal amplification (RealAmp) for the species-specific identification of Plasmodium vivax

Plasmodium vivax infections remain a major source of malaria-related morbidity and mortality. Early and accurate diagnosis is an integral component of effective malaria control programs. Conventional molecular diagnostic methods provide accurate results but are often resource-intensive, expensive, have a long turnaround time and are beyond the capacity of most malaria-endemic countries. Our laboratory has recently developed a new platform called RealAmp, which combines loop-mediated isothermal amplification (LAMP) with a portable tube scanner real-time isothermal instrument for the rapid detection of malaria parasites. Here we describe new primers for the detection of P. vivax using the RealAmp method. Three pairs of amplification primers required for this method were derived from a conserved DNA sequence unique to the P. vivax genome. The amplification was carried out at 64°C using SYBR Green or SYTO-9 intercalating dyes for 90 minutes with the tube scanner set to collect fluorescence signals at 1-minute intervals. Clinical samples of P. vivax and other human-infecting malaria parasite species were used to determine the sensitivity and specificity of the primers by comparing with an 18S ribosomal RNA-based nested PCR as the gold standard. The new set of primers consistently detected laboratory-maintained isolates of P. vivax from different parts of the world. The primers detected P. vivax in the clinical samples with 94.59% sensitivity (95% CI: 87.48-98.26%) and 100% specificity (95% CI: 90.40-100%) compared to the gold standard nested-PCR method. The new primers also proved to be more sensitive than the published species-specific primers specifically developed for the LAMP method in detecting P. vivax.

Jumbled genomes: missing Apicomplexan synteny

Whole-genome comparisons provide insight into genome evolution by informing on gene repertoires, gene gains/losses, and genome organization. Most of our knowledge about eukaryotic genome evolution is derived from studies of multicellular model organisms. The eukaryotic phylum Apicomplexa contains obligate intracellular protist parasites responsible for a wide range of human and veterinary diseases (e.g., malaria, toxoplasmosis, and theileriosis). We have developed an in silico protein-encoding gene based pipeline to investigate synteny across 12 apicomplexan species from six genera. Genome rearrangement between lineages is extensive. Syntenic regions (conserved gene content and order) are rare between lineages and appear to be totally absent across the phylum, with no group of three genes found on the same chromosome and in the same order within 25 kb up- and downstream of any orthologous genes. Conserved synteny between major lineages is limited to small regions in Plasmodium and Theileria/Babesia species, and within these conserved regions, there are a number of proteins putatively targeted to organelles. The observed overall lack of synteny is surprising considering the divergence times and the apparent absence of transposable elements (TEs) within any of the species examined. TEs are ubiquitous in all other groups of eukaryotes studied to date and have been shown to be involved in genomic rearrangements. It appears that there are different criteria governing genome evolution within the Apicomplexa relative to other well-studied unicellular and multicellular eukaryotes.

Discovery and assembly of repeat family pseudomolecules from sparse genomic sequence data using the Assisted Automated Assembler of Repeat Families (AAARF) algorithm

BACKGROUND

Higher eukaryotic genomes are typically large, complex and filled with both genes and multiple classes of repetitive DNA. The repetitive DNAs, primarily transposable elements, are a rapidly evolving genome component that can provide the raw material for novel selected functions and also indicate the mechanisms and history of genome evolution in any ancestral lineage. Despite their abundance, universality and significance, studies of genomic repeat content have been largely limited to analyses of the repeats in fully sequenced genomes.

RESULTS

In order to facilitate a broader range of repeat analyses, the Assisted Automated Assembler of Repeat Families algorithm has been developed. This program, written in PERL and with numerous adjustable parameters, identifies sequence overlaps in small shotgun sequence datasets and walks them out to create long pseudomolecules representing the most abundant repeats in any genome. Testing of this program in maize indicated that it found and assembled all of the major repeats in one or more pseudomolecules, including coverage of the major Long Terminal Repeat retrotransposon families. Both Sanger sequence and 454 datasets were appropriate.

CONCLUSION

These results now indicate that hundreds of higher eukaryotic genomes can be efficiently characterized for the nature, abundance and evolution of their major repetitive DNA components.

Discovery and assembly of repeat family pseudomolecules from sparse genomic sequence data using the Assisted Automated Assembler of Repeat Families (AAARF) algorithm

Background

Higher eukaryotic genomes are typically large, complex and filled with both genes and multiple classes of repetitive DNA. The repetitive DNAs, primarily transposable elements, are a rapidly evolving genome component that can provide the raw material for novel selected functions and also indicate the mechanisms and history of genome evolution in any ancestral lineage. Despite their abundance, universality and significance, studies of genomic repeat content have been largely limited to analyses of the repeats in fully sequenced genomes.

Results

In order to facilitate a broader range of repeat analyses, the Assisted Automated Assembler of Repeat Families algorithm has been developed. This program, written in PERL and with numerous adjustable parameters, identifies sequence overlaps in small shotgun sequence datasets and walks them out to create long pseudomolecules representing the most abundant repeats in any genome. Testing of this program in maize indicated that it found and assembled all of the major repeats in one or more pseudomolecules, including coverage of the major Long Terminal Repeat retrotransposon families. Both Sanger sequence and 454 datasets were appropriate.

Conclusion

These results now indicate that hundreds of higher eukaryotic genomes can be efficiently characterized for the nature, abundance and evolution of their major repetitive DNA components.

The Contribution of LTR Retrotransposon Sequences to Gene Evolution in Mus musculus

Approximately 1.5% of mouse genes (Mus musculus) contain long terminal repeat retrotransposon sequences (LRS). Consistent with earlier findings in Caenorhabditis elegans, Drosophila melanogaster, and Homo sapiens, LRS are more likely to be associated with newly evolved genes. Evidence is presented that LRS are often recruited as novel exons or as spliced additions to existing exons. These novel gene configurations may be expressed initially as alternative transcripts providing an opportunity for the evolution of new gene function.