A Plasmodium knowlesi A1-H.1 transcriptome time course focusing on the late asexual blood stages

This study presents the first transcriptome time-course of the *Plasmodium knowlesi* A1-H.1 strain grown in human erythrocytes, revealing strong genome-wide conservation of temporal gene expression with *P. vivax* and providing an interactive webtool to facilitate comparative functional studies using *P. knowlesi* as a model for *P. vivax*.

The Gist

The Big Picture: A "Model Train" for Malaria Research

Imagine you are trying to study a very shy, elusive animal that lives in a dense jungle. It's hard to catch, hard to keep in a cage, and you can't easily ask it questions. This animal is Plasmodium vivax, a parasite that causes a massive amount of malaria in humans. Because it's so hard to study in a lab, scientists are stuck guessing how it works.

Enter the Plasmodium knowlesi parasite. It's a close cousin of P. vivax (like a tiger and a lion). The best part? P. knowlesi is easy to keep in a lab, grows fast, and can even infect human blood cells in a dish. Scientists have long hoped to use P. knowlesi as a "model train" to understand how the real "train" (P. vivax) works.

But there was a problem: We didn't have a complete instruction manual (a transcriptome) for this specific lab version of P. knowlesi when it's growing in human blood. Without the manual, we didn't know exactly which "parts" (genes) were turning on and off at which times.

This paper is that new instruction manual.

---

What Did the Scientists Do?

Think of the life cycle of a malaria parasite inside a red blood cell as a 27-hour movie. The movie has distinct scenes:

1. The Arrival (Ring stage): The parasite enters the cell.
2. The Growth (Trophozoite stage): It eats and grows.
3. The Reproduction (Schizont stage): It splits into many new babies.
4. The Escape: The babies burst out to infect new cells.

The scientists took a "time-lapse video" of this movie. They grew the parasite in human blood and took snapshots of its genetic activity at five specific moments (5, 14, 20, 24, and 27 hours after it entered the cell).

They used a technology called RNA-seq (think of it as a high-speed camera that takes a picture of every single gene that is "active" or "speaking" at that exact moment).

The Key Findings (The "Plot Twists")

1. The "Just-in-Time" Factory

The study confirmed that the parasite is a master of efficiency. It doesn't keep all its tools out at once. Instead, it uses a "Just-in-Time" delivery system.

• Analogy: Imagine a construction site. You don't bring the paint, the bricks, and the roof all at the start of the day. You bring the bricks when you need to build the wall, and the paint only when the wall is done.
• The Result: The scientists saw that genes for "building the wall" (DNA replication) turned on early, while genes for "painting the door" (invasion tools) turned on only right before the parasite was ready to burst out.

2. The "Cheat Sheet" for Future Experiments

The scientists created two very useful lists:

• The "Steady Eddies": A list of 8 genes that are always "on" at the same volume, no matter what stage the parasite is in.
  • Use: These are like the "volume knob" on a stereo. If you want to measure how loud a specific gene is singing, you compare it to these steady genes to make sure your measurement is accurate.
• The "Stage Markers": A list of 20 genes that scream "I am here!" only at specific times.
  • Use: If you look at a single parasite under a microscope and can't tell if it's a baby or a teenager, you can check if these "marker genes" are active to tell you exactly how old it is. This helps clean up messy data from other studies.

3. The "Twin Test" (Comparing to P. vivax)

This is the most important part. The scientists compared their new P. knowlesi manual with the existing (but incomplete) manual for P. vivax.

• The Verdict: 75% of the genes are twins. They turn on and off at the exact same relative times in both parasites.
• The Metaphor: Imagine two identical twins. One is wearing a blue shirt (P. knowlesi in the lab), and the other is wearing a red shirt (P. vivax in the wild). Even though they are in different places, they are doing the exact same dance steps at the same time.
• Why it matters: This proves that if you study a gene in the lab (P. knowlesi), you can be 95% sure it works the same way in the wild (P. vivax). This is a huge win because it means we can finally test drugs and vaccines in the lab that will actually work on the real human malaria parasite.

4. The "Special Forces" (Invasion Genes)

The scientists looked closely at the genes responsible for the parasite breaking into new cells (the "invasion" genes).

• They found that the "weapons" (proteins) the parasite uses to break into human cells are highly active and turn on right before the parasite bursts out.
• Interestingly, even though this lab strain was adapted to grow in human blood, it still kept some of its "old weapons" designed for monkey blood. This suggests the parasite is holding onto its past, which is fascinating for understanding how it evolves.

5. The "Secret Switch" (Sexual Development)

There is a master switch gene called AP2-G that tells a parasite to stop growing and start making babies (gametocytes) to infect mosquitoes.

• In the lab strain, this switch was only half-turned on. The scientists found that only the end of the gene was being read, not the beginning.
• The Metaphor: It's like trying to read a book but only the last page is printed. The parasite knows it should make babies, but the instructions are incomplete, so it stays in "growth mode" and never makes the sexual stage. This explains why this lab strain is hard to turn into a mosquito-infecting form.

The "Magic Tool" They Built

To make this data useful for everyone, the scientists built a free, interactive website.

• The Metaphor: Imagine a massive library where you can type in the name of any gene from P. vivax or P. knowlesi, and the website instantly shows you a graph of when that gene is active.
• The Benefit: Any scientist in the world can now quickly check: "If I study this gene in the lab, will it tell me anything about the real disease?" They can answer this in seconds instead of months.

The Bottom Line

This paper is like filling in the missing pages of a map.
Before, scientists were trying to navigate the jungle of P. vivax malaria with a blurry, incomplete map. Now, thanks to this study, they have a high-definition, detailed map of the "model" parasite (P. knowlesi) that matches the real one almost perfectly.

This allows researchers to:

1. Test drugs faster in the lab.
2. Design better vaccines by knowing exactly when the parasite's "weak spots" are exposed.
3. Understand the enemy by knowing exactly which genes are doing what, and when.

It's a major step forward in the fight to eliminate malaria.

Technical Summary

1. Problem Statement

• The Research Gap: Plasmodium vivax is the second most prevalent cause of human malaria but lacks a continuous in vitro culture system, severely limiting functional genomics, drug testing, and vaccine development.
• The Proposed Model: Plasmodium knowlesi (specifically the A1-H.1 line) is phylogenetically close to P. vivax and can be cultured in human erythrocytes, making it an ideal surrogate model. However, previous transcriptome data for P. knowlesi were either generated from macaque-adapted strains (Pk1(A+)) using microarrays (older annotation, lower resolution) or via single-cell RNA-seq (scRNA-seq), which suffers from lower sequencing depth and sensitivity for lowly expressed genes.
• The Specific Need: There was a lack of a high-resolution, bulk RNA-sequencing time course for the P. knowlesi A1-H.1 strain grown in human red blood cells (RBCs). Such a dataset is essential to:
  1. Validate the temporal expression conservation between P. knowlesi and P. vivax orthologues.
  2. Identify stage-specific biomarkers to refine scRNA-seq annotations.
  3. Characterize the expression of multigene families involved in invasion and immune evasion.

2. Methodology

• Parasite Culture: The P. knowlesi A1-H.1 line was cultured in human normocytes. Cultures were tightly synchronized to a 1-hour window.
• Sampling Strategy: Bulk RNA was collected at five specific time points during the 27-hour intraerythrocytic developmental cycle (IDC):
  • 5 hours post-invasion (hpi): Rings
  • 14 hpi: Mid-trophozoites
  • 20 hpi: Late trophozoites
  • 24 hpi: Mid-schizonts
  • 27 hpi: Late schizonts (purified via Nycodenz to remove rings).
• Sequencing: mRNA was extracted, globin-depleted, and sequenced using Illumina NovaSeq 6000 (150bp paired-end).
• Bioinformatics Pipeline:
  • Mapping: Reads aligned to the P. knowlesi strain H reference genome (PlasmoDB v2 r66) using STAR.
  • Normalization: Transcripts Per Million (TPM) normalization followed by log2 transformation and Z-scoring.
  • Comparative Analysis: P. vivax (field isolate smru1) and P. falciparum (strain II3) time course data were downloaded, processed similarly, and aligned to a 0–1 relative time scale to account for different IDC lengths (27h vs. 48h vs. 48h).
  • Orthologue Similarity: Dynamic Time Warping (DTW) was used to quantify the similarity of temporal expression patterns between orthologous gene pairs.
  • Annotation: Gene Ontology (GO) enrichment analysis was performed on genes peaking at specific time points.
• Tool Development: An interactive web tool was developed to visualize and compare expression profiles of P. knowlesi and P. vivax orthologues.

3. Key Contributions

1. First High-Resolution Bulk Transcriptome: Generated the first bulk RNA-seq time course for P. knowlesi A1-H.1 in human RBCs with high temporal resolution, particularly in late stages (trophozoite/schizont).
2. Identification of Biomarkers:
   • Identified 8 constitutively expressed genes suitable for qPCR normalization.
   • Identified 20 stage-specific biomarker genes per time point (100 total) to annotate and calibrate existing and future scRNA-seq datasets.
3. Comparative Framework: Established a genome-wide comparison of temporal expression between P. knowlesi and P. vivax, demonstrating that ~75% of orthologues share similar expression dynamics.
4. Interactive Web Tool: Created a public resource (https://interactive.itg.be/app/mal-pk-pv-expression-viewer) allowing researchers to query orthologue expression conservation, facilitating the selection of P. knowlesi targets for P. vivax studies.
5. Multigene Family Characterization: Detailed the expression dynamics of large gene families (SICAvar, kir, PHIST, TRAg) and invasion-associated genes (DBP, NBPX, MSP, RON, etc.).

4. Key Results

• Temporal Expression Patterns: The dataset confirmed a "just-in-time" transcriptional program.
  • Early stages (5–14 hpi): Enriched for transcription, translation, and ribosomal proteins.
  • Late Trophozoite (20 hpi): Enriched for DNA replication.
  • Schizont stages (24–27 hpi): Enriched for invasion-related functions (rhoptry and microneme proteins).
• Conservation with P. vivax:
  • 75% (2,920/3,877) of orthologous gene pairs showed similar temporal expression patterns, validating P. knowlesi A1-H.1 as a robust model for P. vivax biology.
  • 25% were classified as "dissimilar," often linked to species-specific functions (e.g., endoplasmic reticulum trafficking in P. vivax rings).
• Invasion Genes: Most invasion-associated genes (DBP, NBPX, TRAg, MSP) showed conserved expression peaks during the schizont stage.
  • Note: Some P. knowlesi invasion genes showed ring-stage expression not seen in P. vivax, potentially due to technical artifacts (dead schizonts) or biological differences.
  • DBP/NBPX Specifics: Despite adaptation to human RBCs, the strain expressed ligands for macaque receptors (DBPγ, NBPXb) at high levels, while DBPβ was absent due to a genomic deletion.
• ApiAP2 Transcription Factors:
  • Most ApiAP2 factors showed conserved timing.
  • PkAP2-G (Sexual Conversion): Showed a unique expression pattern restricted to the 3' end of the gene (late trophozoite/schizont) with no ring-stage expression, differing from P. falciparum but resembling P. vivax. This suggests a potential mechanism for the lack of gametocyte production in this specific in vitro line.
• Multigene Families:
  • SICAvar: Generally low expression (consistent with in vitro culture lacking immune pressure), with one highly expressed member (PKNH_0620500).
  • kir (pir): Mostly low expression, with a few highly expressed genes peaking in rings or schizonts.
  • PHIST & TRAg: High expression levels with distinct temporal patterns; TRAg genes showed dual expression in schizonts and rings.

5. Significance

• Advancing P. vivax Research: This study provides the critical reference framework needed to use P. knowlesi A1-H.1 as a proxy for P. vivax. By confirming that 75% of genes share temporal expression profiles, researchers can confidently study P. vivax gene function using the tractable P. knowlesi system.
• Refining Single-Cell Data: The identified stage-specific biomarkers allow for the precise annotation of P. knowlesi and P. vivax scRNA-seq data, correcting previous annotations based on P. berghei (a rodent malaria) which may not perfectly align with human malaria stages.
• Drug and Vaccine Targets: The detailed characterization of invasion ligands and their expression timing highlights potential targets for transmission-blocking vaccines and therapeutics, specifically those active during the schizont stage.
• Resource Accessibility: The open-access web tool democratizes access to comparative transcriptomics, enabling rapid hypothesis generation for neglected tropical disease research.

In conclusion, this paper delivers a high-quality, human-adapted P. knowlesi transcriptome resource that bridges the gap between P. vivax biology and experimental tractability, offering a powerful tool for the global malaria research community.