Skip to content
1887

Abstract

As the importance of transcriptional variation and regulation for becomes more apparent, advances for non- species are hindered by our reliance on natural infections to study parasite biology. Untargeted transcriptomic research is also complicated by low-parasite densities and high proportions of human genetic material, highlighting the need for optimized sample processing protocols. In this study, we used a culture diluted in whole blood as a mock natural infection to compare white blood cell (WBC), rRNA and globin depletion methods and RNA-seq library preparation kits to create an optimized protocol for low-volume sample processing. Depletion techniques and library preparation kit selection are both crucial to enrich parasite RNA in low-parasite density and low-volume samples. WBC depletion methods can increase the percentage of mapped paired-end (PE) reads from 5.4 to 37.2–54.5%, depending on the method. Additionally, globin RNA and rRNA depletion significantly decreased the human globin RNA and rRNA combined gene counts up to 90%. When comparing the library preparation kits, the transcription profiles of the mRNA kits were highly correlated and enriched for protein-coding genes. The total RNA library prep kit was dominated by rRNA, despite having the highest percentage of mapped PE reads. Parasite sequencing output can be optimized when WBC and globin depletion are used in combination with an mRNA library preparation kit bringing untargeted transcriptomics to resource-limited settings which will lead to new insights in parasite biology.

Funding
This study was supported by the:
  • Fonds Wetenschappelijk Onderzoek (Award PhD-SB Fellowship (1SC5522N))
    • Principal Award Recipient: ErinSauve
  • Departement Economie, Wetenschap en Innovatie (Award SOFI)
    • Principal Award Recipient: AnnaRosanas-Urgell
  • This is an open-access article distributed under the terms of the Creative Commons Attribution License.
Loading

Article metrics loading...

/content/journal/mgen/10.1099/mgen.0.001546
2025-11-20
2025-12-16

Metrics

Loading full text...

Full text loading...

/deliver/fulltext/mgen/11/11/mgen001546.html?itemId=/content/journal/mgen/10.1099/mgen.0.001546&mimeType=html&fmt=ahah

References

  1. Zhu L, Van Der Pluijm RW, Kucharski M, Nayak S, Tripathi J et al. Artemisinin resistance in the malaria parasite, Plasmodium falciparum, originates from its initial transcriptional response. Commun Biol 2022; 5:274 [View Article] [PubMed]
    [Google Scholar]
  2. Nayak S, Peto TJ, Kucharski M, Tripura R, Callery JJ et al. Population genomics and transcriptomics of Plasmodium falciparum in Cambodia and Vietnam uncover key components of the artemisinin resistance genetic background. Nat Commun 2024; 15:10625 [View Article]
    [Google Scholar]
  3. Llorà-Batlle O, Michel-Todó L, Witmer K, Toda H, Fernández-Becerra C et al. Conditional expression of PfAP2-G for controlled massive sexual conversion in Plasmodium falciparum. Sci Adv 2020; 6:eaaz5057 [View Article] [PubMed]
    [Google Scholar]
  4. Prajapati SK, Ayanful-Torgby R, Pava Z, Barbeau MC, Acquah FK et al. The transcriptome of circulating sexually committed Plasmodium falciparum ring stage parasites forecasts malaria transmission potential. Nat Commun 2020; 11:6159 [View Article] [PubMed]
    [Google Scholar]
  5. Thomson-Luque R, Votborg-Novél L, Ndovie W, Andrade CM, Niangaly M et al. Plasmodium falciparum transcription in different clinical presentations of malaria associates with circulation time of infected erythrocytes. Nat Commun 2021; 12:4711 [View Article] [PubMed]
    [Google Scholar]
  6. Tebben K, Yirampo S, Coulibaly D, Koné AK, Laurens MB et al. Immune gene expression changes more during a malaria transmission season than between consecutive seasons. Microbiol Spectr 2024; 12:e0096024 [View Article] [PubMed]
    [Google Scholar]
  7. Sriprawat K, Kaewpongsri S, Suwanarusk R, Leimanis ML, Lek-Uthai U et al. Effective and cheap removal of leukocytes and platelets from Plasmodium vivax infected blood. Malar J 2009; 8:115 [View Article] [PubMed]
    [Google Scholar]
  8. Venkatesan M, Amaratunga C, Campino S, Auburn S, Koch O et al. Using CF11 cellulose columns to inexpensively and effectively remove human DNA from Plasmodium falciparum-infected whole blood samples. Malar J 2012; 11:41 [View Article] [PubMed]
    [Google Scholar]
  9. Shaw-Saliba K, Thomson-Luque R, Obaldía N, Nuñez M, Dutary S et al. Insights into an optimization of Plasmodium vivax Sal-1 in vitro culture: the Aotus primate model. PLoS Negl Trop Dis 2016; 10:e0004870 [View Article] [PubMed]
    [Google Scholar]
  10. Li J, Tao Z, Li Q, Brashear A, Wang Y et al. Further evaluation of the NWF filter for the purification of Plasmodium vivax-infected erythrocytes. Malar J 2017; 16:201 [View Article] [PubMed]
    [Google Scholar]
  11. Tao Z-Y, Xia H, Cao J, Gao Q. Development and evaluation of a prototype non-woven fabric filter for purification of malaria-infected blood. Malar J 2011; 10:251 [View Article] [PubMed]
    [Google Scholar]
  12. JM, Cannon MV, Caleon RL, Wellems TE, Serre D. Single-cell transcription analysis of Plasmodium vivax blood-stage parasites identifies stage- and species-specific profiles of expression. PLoS Biol 2020; 18:e3000711 [View Article] [PubMed]
    [Google Scholar]
  13. Hazzard B, JM, Ellis AC, Pascini TV, Amin S et al. Long read single cell RNA sequencing reveals the isoform diversity of Plasmodium vivax transcripts. PLoS Negl Trop Dis 2022; 16:e0010991 [View Article] [PubMed]
    [Google Scholar]
  14. Hazzard B, JM, Bogale HN, Pascini TV, Ellis AC et al. Single-cell analyses of polyclonal Plasmodium vivax infections and their consequences on parasite transmission. Nat Commun 2024; 15:7625 [View Article] [PubMed]
    [Google Scholar]
  15. LeRoux M, Lakshmanan V, Daily JP. Plasmodium falciparum biology: analysis of in vitro versus in vivo growth conditions. Trends Parasitol 2009; 25:474–481 [View Article] [PubMed]
    [Google Scholar]
  16. Lapp SA, Mok S, Zhu L, Wu H, Preiser PR et al. Plasmodium knowlesi gene expression differs in ex vivo compared to in vitro blood-stage cultures. Malar J 2015; 14:110 [View Article] [PubMed]
    [Google Scholar]
  17. Andradi-Brown C, Wichers-Misterek JS, von Thien H, Höppner YD, Scholz JAM et al. A novel computational pipeline for var gene expression augments the discovery of changes in the Plasmodium falciparum transcriptome during transition from in vivo to short-term in vitro culture. Elife 2024; 12:RP87726 [View Article] [PubMed]
    [Google Scholar]
  18. Simantov K, Goyal M, Dzikowski R. Emerging biology of noncoding RNAs in malaria parasites. PLoS Pathog 2022; 18:e1010600 [View Article] [PubMed]
    [Google Scholar]
  19. Wahl A, Huptas C, Neuhaus K. Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing. Sci Rep 2022; 12:5765 [View Article] [PubMed]
    [Google Scholar]
  20. Adiconis X, Borges-Rivera D, Satija R, DeLuca DS, Busby MA et al. Comparative analysis of RNA sequencing methods for degraded or low-input samples. Nat Methods 2013; 10:623–629 [View Article] [PubMed]
    [Google Scholar]
  21. Jang JS, Berg B, Holicky E, Eckloff B, Mutawe M et al. Comparative evaluation for the globin gene depletion methods for mRNA sequencing using the whole blood-derived total RNAs. BMC Genomics 2020; 21:890 [View Article] [PubMed]
    [Google Scholar]
  22. Haile S, Corbett RD, Bilobram S, Mungall K, Grande BM et al. Evaluation of protocols for rRNA depletion-based RNA sequencing of nanogram inputs of mammalian total RNA. PLoS One 2019; 14:e0224578 [View Article] [PubMed]
    [Google Scholar]
  23. Zhao S, Zhang Y, Gamini R, Zhang B, von Schack D. Evaluation of two main RNA-seq approaches for gene quantification in clinical RNA sequencing: polyA+ selection versus rRNA depletion. Sci Rep 2018; 8:4781 [View Article]
    [Google Scholar]
  24. Moon RW, Hall J, Rangkuti F, Ho YS, Almond N et al. Adaptation of the genetically tractable malaria pathogen Plasmodium knowlesi to continuous culture in human erythrocytes. Proc Natl Acad Sci USA 2013; 110:531–536 [View Article] [PubMed]
    [Google Scholar]
  25. Wingett SW, Andrews S. FastQ screen: a tool for multi-genome mapping and quality control. F1000Res 2018; 7:1338 [View Article] [PubMed]
    [Google Scholar]
  26. Patel H, Ewels P, Peltzer A, Manning J, Botvinnik O et al.nf-core/rnaseq: nf-core/rnaseq v3.14.0 - Hassium Honey Badger (3.14.0) Zenodo; 2024
  27. Di Tommaso P, Chatzou M, Floden EW, Barja PP, Palumbo E et al. Nextflow enables reproducible computational workflows. Nat Biotechnol 2017; 35:316–319 [View Article] [PubMed]
    [Google Scholar]
  28. Ewels PA, Peltzer A, Fillinger S, Patel H, Alneberg J et al. The nf-core framework for community-curated bioinformatics pipelines. Nat Biotechnol 2020; 38:276–278 [View Article] [PubMed]
    [Google Scholar]
  29. Andrews S. FastQC: a quality control tool for high throughput sequence data. n.d
  30. Krueger F, James F, Ewels P, Afyounian E, Schuster-Boeckler B. FelixKrueger/TrimGalore: v0.6.7 - DOI via Zenodo (0.6.7). Zenodo; 2021
  31. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J 2011; 17:10 [View Article]
    [Google Scholar]
  32. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013; 29:15–21 [View Article] [PubMed]
    [Google Scholar]
  33. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods 2017; 14:417–419 [View Article] [PubMed]
    [Google Scholar]
  34. Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V et al. Twelve years of SAMtools and BCFtools. Gigascience 2021; 10:giab008 [View Article] [PubMed]
    [Google Scholar]
  35. R Core Team R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna; 2021 https://www.R-project.org
  36. Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res 2015; 4:1521 [View Article] [PubMed]
    [Google Scholar]
  37. Wickham H. ggplot2: Elegant Graphics for Data Analysis, 2nd ed Cham: Springer international publishing; 2016
    [Google Scholar]
  38. Garnier S, Ross N, Rudis R, Camargo PA, Sciaini M et al. viridis(Lite) - colorblind-friendly color maps for R; 2024 https://sjmgarnier.github.io/viridis/
  39. Schloerke B, Cook D, Larmarange J, Briatte F, Marbach M et al. GGally: extension to 'ggplot2'; 2024 https://CRAN.R-project.org/package=GGally
  40. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 2014; 15:550 [View Article] [PubMed]
    [Google Scholar]
  41. Kolde R. pheatmap: pretty heatmaps; 2019 https://CRAN.R-project.org/package=pheatmap
  42. Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol 2019; 37:773–782 [View Article] [PubMed]
    [Google Scholar]
  43. Tebben K, Dia A, Serre D. Determination of the stage composition of Plasmodium infections from bulk gene expression data. mSystems 2022; 7:e0025822 [View Article] [PubMed]
    [Google Scholar]
  44. Kim A, Popovici J, Menard D, Serre D. Plasmodium vivax transcriptomes reveal stage-specific chloroquine response and differential regulation of male and female gametocytes. Nat Commun 2019; 10:371 [View Article]
    [Google Scholar]
  45. Ribaut C, Berry A, Chevalley S, Reybier K, Morlais I et al. Concentration and purification by magnetic separation of the erythrocytic stages of all human Plasmodium species. Malar J 2008; 7:45 [View Article] [PubMed]
    [Google Scholar]
  46. Janse CJ, Camargo A, Portillo HAD, Waters AP, Herrera S et al. Removal of leucocytes from Plasmodium vivax-infected blood. Ann Trop Med Parasitol 1994; 88:213–216 [View Article]
    [Google Scholar]
  47. Kim A, Popovici J, Vantaux A, Samreth R, Bin S et al. Characterization of P. vivax blood stage transcriptomes from field isolates reveals similarities among infections and complex gene isoforms. Sci Rep 2017; 7:7761 [View Article] [PubMed]
    [Google Scholar]
  48. Kucharski M, Tripathi J, Nayak S, Zhu L, Wirjanata G et al. A comprehensive RNA handling and transcriptomics guide for high-throughput processing of Plasmodium blood-stage samples. Malar J 2020; 19:363 [View Article] [PubMed]
    [Google Scholar]
  49. Bourgard C, Lopes SCP, Lacerda MVG, Albrecht L, Costa FTM. A suitable RNA preparation methodology for whole transcriptome shotgun sequencing harvested from Plasmodium vivax-infected patients. Sci Rep 2021; 11:5089 [View Article] [PubMed]
    [Google Scholar]
  50. World Health Organization.. World malaria report 2024: addressing inequity in the global malaria response. Geneva: World Health Organization. Licence: CC BY-NC-SA 3.0 IGO 2024
    [Google Scholar]
  51. Bozdech Z, Mok S, Hu G, Imwong M, Jaidee A et al. The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites. Proc Natl Acad Sci USA 2008; 105:16290–16295 [View Article] [PubMed]
    [Google Scholar]
  52. Zhu L, Mok S, Imwong M, Jaidee A, Russell B et al. New insights into the Plasmodium vivax transcriptome using RNA-Seq. Sci Rep 2016; 6:20498 [View Article]
    [Google Scholar]
  53. Siegel SV, Chappell L, Hostetler JB, Amaratunga C, Suon S et al. Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions. Sci Rep 2020; 10:16667 [View Article] [PubMed]
    [Google Scholar]
  54. Rangel GW, Clark MA, Kanjee U, Goldberg JM, MacInnis B et al. Plasmodium vivax transcriptional profiling of low input cryopreserved isolates through the intraerythrocytic development cycle. PLoS Negl Trop Dis 2020; 14:e0008104 [View Article] [PubMed]
    [Google Scholar]
  55. De Meulenaere K, Prajapati SK, Villasis E, Cuypers B, Kattenberg JH et al. Band 3–mediated Plasmodium vivax invasion is associated with transcriptional variation in PvTRAg genes. Front Cell Infect Microbiol 2022; 12:1011692 [View Article]
    [Google Scholar]
  56. Kepple D, Ford CT, Williams J, Abagero B, Li S et al. Comparative transcriptomics reveal differential gene expression among Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms. PLoS Negl Trop Dis 2024; 18:e0011926 [View Article] [PubMed]
    [Google Scholar]
  57. Shapiro HM, Apte SH, Chojnowski GM, Hänscheid T, Rebelo M et al. Cytometry in malaria--a practical replacement for microscopy?. Curr Protoc Cytom 2013; Chapter 11:11 [View Article] [PubMed]
    [Google Scholar]
  58. Reuling IJ, Stone WJR, van de Vegte-Bolmer M, van Gemert G-J, Siebelink-Stoter R et al. Concentration of Plasmodium falciparum gametocytes in whole blood samples by magnetic cell sorting enhances parasite infection rates in mosquito feeding assays. Malar J 2017; 16:315 [View Article] [PubMed]
    [Google Scholar]
  59. Omorou R, Bin Sa’id I, Delves M, Severini C, Kouakou YI et al. Protocols for Plasmodium gametocyte production in vitro: an integrative review and analysis. Parasites Vectors 2022; 15:451 [View Article]
    [Google Scholar]
  60. Stewart LB, Freville A, Voss TS, Baker DA, Awandare GA et al. Plasmodium falciparum sexual commitment rate variation among clinical isolates and diverse laboratory-adapted lines. Microbiol Spectr 2022; 10:e0223422 [View Article] [PubMed]
    [Google Scholar]
/content/journal/mgen/10.1099/mgen.0.001546
Loading
/content/journal/mgen/10.1099/mgen.0.001546
Loading

Data & Media loading...

Supplements

Supplementary material 1

PDF
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An error occurred
Approval was partially successful, following selected items could not be processed due to error