# Molecular detection of human Plasmodium species using a multiplex real time PCR

Yassamine Lazrek, Célia Florimond, Béatrice Volney, Manon Discours, Emilie Mosnier, Sandrine Houzé, Stephane Pelleau, Lise Musset

► To cite this version:

Yassamine Lazrek, Célia Florimond, Béatrice Volney, Manon Discours, Emilie Mosnier, et al.. Molecular detection of human Plasmodium species using a multiplex real time PCR. Scientific Reports, 2023, 13 (1), pp.11388. 10.1038/s41598-023-38621-9 . pasteur-04185534

HAL Id: pasteur-04185534

https://pasteur.hal.science/pasteur-04185534

Submitted on 22 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Molecular detection of human
*Plasmodium* species using
a multiplex real time PCR

Yassamine Lazrek<fnref n="1" />✉, Célia Florimond<fnref n="1" />, Béatrice Volney<fnref n="1" />, Manon Discours<fnref n="1" />, Emilie Mosnier², Sandrine Houzé³,⁴, Stéphane Pelleau<fnref n="1" />,⁵ & Lise Musset<fnref n="1" />

Molecular detection methods have revealed higher sensitivity and specificity than conventional microscopy or rapid diagnostic tests for malaria diagnosis. In this study, we implemented, evaluated and validated according to the ISO 15,189 requirements, a multiplex real-time PCR assay to detect and identify the five human malaria parasites. DNA samples were extracted from whole blood or dried blood spots drawn from patients. Based on the External Quality Assessment (whole blood), this method shows 100% sensitivity and specificity. This PCR detected *P. vivax* up to 0.25 p/µl, *P. falciparum* and *P. knowlesi* up to 0.5 p/µl, *P. ovale* up to 1 p/µl and *P. malariae* up to 5 p/µl of blood. From blood spots (extraction from four punches), it detected *P. vivax* at 5 p/µl, *P. falciparum*, *P. ovale* and *P. knowlesi* at 20 p/µl and *P. malariae* at 125 p/µl. In conclusion, this quantitative PCR shows excellent performance, is easy to use and DNA saver. It is especially useful to actively screen large population groups and identify the five human malaria parasites in a context of low malaria transmission.

Malaria is one of the deadliest diseases and has already claim millions of life throughout the world, 619,000 in 2021¹. This vector-borne disease is caused in Humans by five different species of *Plasmodium* genus, *P. falciparum*, *P. vivax*, *P. malariae*, *P. ovale*, and *P. knowlesi*². The last one has been recently transmitted to Humans from primates in Malaysia³. *P. falciparum* generates the highest mortality and morbidity rate. However, *P. vivax* is the most widely distributed across all endemic areas and the second largest contributor to clinical malaria worldwide<fnref n="1" />. This species was originally considered benign but is now recognized as a cause of serious morbidity and mortality⁴,⁵. Accurate malaria diagnosis is crucial to reduce presumptive treatment. Microscopy is the main tool to diagnose malaria together with rapid diagnosis tests (RDTs). Their performances allow the diagnosis of patient with a parasitemia up to 10 parasites/µl of blood. Owing to their sensitivities, they are limited to identify asymptomatic carriers presenting low parasitemia. However, addressing these parasite reservoirs is crucial for national malaria programs moving forward malaria elimination. These programs call for novel diagnosis methods, highly sensitive but costly effective to actively screen large population groups⁶.

The molecular amplification of the small 18S subunit of ribosomal RNA (18S rRNA), first implemented by Snounou et al. using a nested PCR technique, is the most widely used molecular diagnostic tools in medical laboratories or in research programs⁷,⁸. The 18S rRNA gene is commonly the target because of its 5 to 7 copies per genome, which increases sensitivity⁹. It is also conserved in all *Plasmodium* species with specific parts for each of them.

Since then, molecular techniques have evolved and now include, real-time PCR, Reverse Transcriptase PCR, Nucleic Acid Sequence-Based Amplification (NASBA) and LAMP (Loop-mediated isothermal amplification). These PCR can detect a wide range of gene targets and are 10 to 100 times more sensitive than microscopy¹⁰,¹¹. Compare to the original nested-PCR, they also give results more rapidly (within an hour compared to 4.5–10 h).

In this study, we developed and evaluated a sensitive and specific multiplex real-time PCR procedure to detect and identify the five malaria parasites transmitted to Humans (*P. falciparum*, *P. vivax*, *P. malaria*, *P. ovale* and *P. knowlesi*). Adapted from Shokopoulos et al.¹² and de Canale et al.¹³, two real-time multiplex PCR have been

<footnote n="1">¹Laboratoire de Parasitologie, World Health Organization Collaborating Center for Surveillance of Antimalarial Drug Resistance, Centre Nationale de Référence du Paludisme, Institut Pasteur de la Guyane, Cayenne, French Guiana. ²Aix Marseille Univ, INSERM, IRD, SESSTIM, Sciences Economiques & Sociales de la Santé & Traitement de l'Information Médicale, Aix Marseille Institute of Public Health ISSPAM, 13385 Marseille, France. ³IRD, MERIT, Université Paris Cité, 75006 Paris, France. ⁴Centre National de Référence du Paludisme, AP-HP, Hôpital Bichat-Claude Bernard, 75018 Paris, France. ⁵Infectious Disease Epidemiology and Analytics Unit, Institut Pasteur, Université Paris Cité, Paris, France. ✉email: ylazrek@pasteur-cayenne.fr</footnote>
implemented, one targeting *P. falciparum* and *P. vivax* and another targeting *P. malariae*, *P. ovale* and *P. knowlesi*. The multiplex approach has been chosen in order to minimize DNA, time and money consumption. DNA extraction was standardized and controlled using a third simplex PCR targeting the human macroglobulin gene. If a quick and large screening is required during active case detection, we also implemented a *Plasmodium* spp. PCR method adapted from Hassanpour et al.<sup>14</sup>. Those methods have been accredited according to the ISO 15,189 medical biology norm to analyze DNA extracted from venous or capillary blood.

## Results

**Repeatability, reproducibility and internal quality controls of 18S screening real time PCR.** The coefficients of variation (CVs) for intra-assay repeatability and inter-assay reproducibility were showed in Tables 1 and 2.

Internal quality controls (IQC) are included in each sample series. They were first analyzed according to the Levy-Jennings rules based on the coefficient of variation established with repeatability (Table 1). After the first months of implementation, reproducibility values have been determined (Table 2), then simplify and harmonized at 2.5% for all species.

To validate the DNA extraction, a threshold was defined for the internal control human β2Megaglobuline. Based on the results obtained on 58 samples of the “Palustop project”<sup>15</sup> from whole blood (WB) and 31 from dried blood spot (DBS) of the “NRC collection”, Ct ranged between 21 and 24 (mean 22.5 ± 0.6) and between 25.2 and 32.3 (29.1 ± 1.7, Fig. 1) for WB and DBS, respectively. Therefore, the thresholds to validate DNA extraction were fixed, Ct ≤ 24 for WB DNA and, Ct ≤ 33 for DBS DNA.

**Analytical sensitivity and specificity.** The performance of this species-specific real-time PCR was compared to our External Quality Assessment (EQA) to validate the ability of the method to diagnose malaria during patient care. This included samples collected in tubes or dried blood spots infected by one of the five species at different parasitemia, from 0.018 parasites per µl (p/µl) up to 1.1 × 10⁶ p/µl. We also tested several negative EQA. The sensitivity and specificity of the technique for all species in whole blood were 100% (Table 3).

For DNA extracted on dried blood spot (Table 4), the sensitivity was 100% for *P. vivax*, *P. ovale* and *P. knowlesi*. However, for *P. falciparum* and *P. malariae*, the sensitivity and the negative predictive value (NPV) were 73% and 100%, 43% and 100%, respectively. The specificities were at 100%.

**Real-time PCR limits of detection (LoD).** In diagnosis, when the parameter is close to the limit of detection, a sample analyzed several times could be either positive or negative. When a result is based on quantitative

<table><caption>**Table 1.** Repeatability of the real-time PCR approach to detect *Plasmodium* species. Detailed values for high (*) and low (°) parasitemia samples are listed in Table 9.</caption><thead><tr><th></th><th>P. falciparum</th><th>P. vivax</th><th>P. malariae</th><th>P. ovale</th><th>P. knowlesi</th></tr></thead><tbody><tr><td>Coefficient of variation at high-parasitemia*</td><td>0.3%</td><td>2.4%</td><td>1.0%</td><td>0.5%</td><td>0.6%</td></tr><tr><td>Coefficient of variation at low-parasitemia°</td><td>0.6%</td><td>3.1%</td><td>0.5%</td><td>2.2%</td><td>1.0%</td></tr><tr><td>Number of replicates (n)</td><td>n = 3</td><td>n = 3</td><td>n = 3</td><td>n = 3</td><td>n = 3</td></tr></tbody></table>

<table><caption>**Table 2.** Reproducibility of the real-time PCR approach to detect *Plasmodium* species.</caption><thead><tr><th></th><th>P. falciparum</th><th>P. vivax</th><th>P. malariae</th><th>P. ovale</th><th>P. knowlesi</th></tr></thead><tbody><tr><td>Internal control (copy/µl)</td><td>3200</td><td>3300</td><td>100</td><td>3100</td><td>1000</td></tr><tr><td>Coefficient of variation</td><td>1.2%</td><td>1.0%</td><td>3.6%</td><td>2.5%</td><td>2.5%</td></tr><tr><td>Number of replicates (n)</td><td>n = 78</td><td>n = 80</td><td>n = 90</td><td>n = 95</td><td>n = 58</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Human β2Megaglobuline Ct values.</figcaption></figure>
<table><caption>**Table 3.** Sensitivity and specificity regarding the external quality assessment on 200 µl of blood. EQA External Quality Assessment, NPV negative predictive value, PPV positive predictive value. Parasitemia of the positive samples: *P. falciparum* from 0.1 to 1.10⁶ p/µl, *P. vivax* from 0.018 to 2.10⁵ p/µl, *P. malariae* from 1 to 125 p/µl, *P. ovale* from 1 to 200 p/µl and *P. knowlesi* from 1 to 1.10⁴ p/µl.</caption><thead><tr><th></th><th></th><th>Methods</th><th>Positive EQA</th><th>Negative EQA</th><th>Sensitivity (%)</th><th>Specificity (%)</th><th>PPV (%)</th><th>NPV (%)</th></tr></thead><tbody><tr><td rowspan="3">P. falciparum</td><td rowspan="3"></td><td>Positive</td><td>22</td><td>0</td><td rowspan="3">100 ± 13.6</td><td rowspan="3">100 ± 15.8</td><td rowspan="3">100 ± 13.6</td><td rowspan="3">100 ± 15.8</td></tr><tr><td>Negative</td><td>0</td><td>19</td></tr><tr><td>Total</td><td>22</td><td>19</td></tr><tr><td rowspan="3">P. vivax</td><td rowspan="3"></td><td>Positive</td><td>16</td><td>0</td><td rowspan="3">100 ± 18.8</td><td rowspan="3">100 ± 15.8</td><td rowspan="3">100 ± 18.8</td><td rowspan="3">100 ± 15.8</td></tr><tr><td>Negative</td><td>0</td><td>19</td></tr><tr><td>Total</td><td>16</td><td>19</td></tr><tr><td rowspan="3">P. malariae</td><td rowspan="3">Real time PCR</td><td>Positive</td><td>10</td><td>0</td><td rowspan="3">100 ± 30.0</td><td rowspan="3">100 ± 15.8</td><td rowspan="3">100 ± 30.0</td><td rowspan="3">100 ± 15.8</td></tr><tr><td>Negative</td><td>0</td><td>19</td></tr><tr><td>Total</td><td>10</td><td>19</td></tr><tr><td rowspan="3">P. ovale</td><td rowspan="3"></td><td>Positive</td><td>9</td><td>0</td><td rowspan="3">100 ± 33.3</td><td rowspan="3">100 ± 15.8</td><td rowspan="3">100 ± 33.3</td><td rowspan="3">100 ± 15.8</td></tr><tr><td>Negative</td><td>0</td><td>19</td></tr><tr><td>Total</td><td>9</td><td>19</td></tr><tr><td rowspan="3">P. knowlesi</td><td rowspan="3"></td><td>Positive</td><td>9</td><td>0</td><td rowspan="3">100 ± 33.3</td><td rowspan="3">100 ± 15.8</td><td rowspan="3">100 ± 33.3</td><td rowspan="3">100 ± 15.8</td></tr><tr><td>Negative</td><td>0</td><td>19</td></tr><tr><td>Total</td><td>9</td><td>19</td></tr></tbody></table>

<table><caption>**Table 4.** Sensitivity and specificity regarding the External quality assessment in dried blood spot. EQA External Quality Assessment, NPV negative predictive value, PPV positive predictive value. Parasitemia of the positive samples: *P. falciparum* from 0.2 to 1.10⁶ p/µl, *P. vivax* from 5 to 400 p/µl, *P. malariae* from 2.3 to 800 p/µl, *P. ovale* from 20 to 280 p/µl and *P. knowlesi* from 2 to 20 p/µl.</caption><thead><tr><th></th><th></th><th>Methods</th><th>Positive EQA</th><th>Negative EQA</th><th>Sensitivity (%)</th><th>Specificity (%)</th><th>PPV (%)</th><th>NPV (%)</th></tr></thead><tbody><tr><td rowspan="3">P. falciparum</td><td rowspan="3"></td><td>Positive</td><td>8</td><td>0</td><td rowspan="3">73 ± 26</td><td rowspan="3">100 ± 25</td><td rowspan="3">73 ± 26</td><td rowspan="3">100 ± 25</td></tr><tr><td>Negative</td><td>3</td><td>12</td></tr><tr><td>Total</td><td>11</td><td>12</td></tr><tr><td rowspan="3">P. vivax</td><td rowspan="3"></td><td>Positive</td><td>6</td><td>0</td><td rowspan="3">100 ± 50</td><td rowspan="3">100 ± 25</td><td rowspan="3">100 ± 50</td><td rowspan="3">100 ± 25</td></tr><tr><td>Negative</td><td>0</td><td>12</td></tr><tr><td>Total</td><td>6</td><td>12</td></tr><tr><td rowspan="3">P. malariae</td><td rowspan="3">Real time PCR</td><td>Positive</td><td>3</td><td>0</td><td rowspan="3">43 ± 37</td><td rowspan="3">100 ± 25</td><td rowspan="3">43 ± 37</td><td rowspan="3">100 ± 25</td></tr><tr><td>Negative</td><td>4</td><td>12</td></tr><tr><td>Total</td><td>7</td><td>12</td></tr><tr><td rowspan="3">P. ovale</td><td rowspan="3"></td><td>Positive</td><td>4</td><td>0</td><td rowspan="3">100 ± 75</td><td rowspan="3">100 ± 25</td><td rowspan="3">100 ± 75</td><td rowspan="3">100 ± 25</td></tr><tr><td>Negative</td><td>0</td><td>12</td></tr><tr><td>Total</td><td>4</td><td>12</td></tr><tr><td rowspan="3">P. knowlesi</td><td rowspan="3"></td><td>Positive</td><td>3</td><td>0</td><td rowspan="3">100 ± 100</td><td rowspan="3">100 ± 25</td><td rowspan="3">100 ± 100</td><td rowspan="3">100 ± 25</td></tr><tr><td>Negative</td><td>0</td><td>12</td></tr><tr><td>Total</td><td>3</td><td>12</td></tr></tbody></table>

data, it is crucial to determine the limits of this zone, named the grey zone. Outside this zone, the samples should be 100% positive among replicates while within the grey zone, a percentage could be applied. For these malaria PCRs, the grey zone included Ct > 36. To increase sensitivities, samples were analyzed in triplicate. Within the grey zone and according to the specificities of the method, as soon as a sample was positive among the triplicate, the result was positive for malaria infection.

Based on these rules, the detection limits of these real-time PCR have been compared to the historical nested PCR. After DNA extraction of 200 µl of venous blood, the method detected *P. vivax* up to 0.25 p/µl, *P. falciparum* and *P. knowlesi* up to 0.5 p/µl, *P. ovale* up to 1 p/µl (Table 5) and *P. malariae* up to 5 p/µl (Ct = 39.4 ± 1.5, positive wells 3/3). If we include the grey zone results, this detection limits increase up to 0.06 p/µl for *P. knowlesi*, 0.12 p/µl for *P. falciparum* and *P. vivax* and 1 p/µl for *P. malariae*. When DNA was extracted from capillary blood and collected on filter papers, the method allowed the detection of *P. vivax* infection up to 10 p/µl and *P. falciparum* up to 5 p/µl (Supplementary Table S1).

The real time PCR targeting the genus *Plasmodium spp.* demonstrated similar limits of detection than the PCR-species specific (Table 6).
<table><caption>**Table 5.** Comparison of the detection limits of nested and real time PCR for *P. falciparum*, *P. vivax* *P. malariae*, *P. ovale* and *P. knowlesi* based on DNA extracted from 200 µl of whole blood. Results are presented in mean Ct (cycle threshold) ± standard deviation according to the level of parasitemia (p/µl, parasite/µl). In bold and italics, the LoD chosen. NA non-applicable. *Number of positive well in each real time triplicate.</caption><thead><tr><th>Parasitemia (p/µl)</th><th></th><th>100</th><th>20</th><th>10</th><th>1</th><th>0.5</th><th>0.25</th><th>0.12</th><th>0.06</th><th>0.03</th></tr></thead><tbody><tr><td rowspan="3">P. falciparum</td><td>Nested PCR</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Real Time PCR</td><td>28.6 ± 0.4</td><td>32.1 ± 0.3</td><td>32.3 ± 0.1</td><td>36.0 ± 0.2</td><td><strong>37.3 ± 0.8</strong></td><td>36.2 ± 0.7</td><td>38.0 ± NA</td><td>-</td><td>-</td></tr><tr><td>Positive wells*</td><td>3/3</td><td>3/3</td><td>3/3</td><td>3/3</td><td><strong>3/3</strong></td><td>2/3</td><td>1/3</td><td>0/3</td><td>0/3</td></tr><tr><td rowspan="3">P. vivax</td><td>Nested PCR</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Real Time PCR</td><td>28.8 ± 0.0</td><td>31.0 ± 0.1</td><td>31.6 ± 0.1</td><td>36.0 ± 0.9</td><td>36.6 ± 1.0</td><td><strong>37.6 ± 1.3</strong></td><td>38.0 ± NA</td><td>-</td><td>-</td></tr><tr><td>Positive wells*</td><td>3/3</td><td>3/3</td><td>3/3</td><td>3/3</td><td><strong>3/3</strong></td><td>1/3</td><td>0/3</td><td>0/3</td><td>0/3</td></tr><tr><td rowspan="3">P. malariae</td><td>Nested PCR</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Real Time PCR</td><td>35.4 ± 0.4</td><td>38.5 ± 1.0</td><td><strong>38.1 ± 1.0</strong></td><td>41.0 ± NA</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Positive wells*</td><td>3/3</td><td>3/3</td><td><strong>3/3</strong></td><td>1/3</td><td>0/3</td><td>0/3</td><td>0/3</td><td>0/3</td><td>0/3</td></tr><tr><td rowspan="3">P. ovale</td><td>Nested PCR</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Real Time PCR</td><td>31.0 ± 0.2</td><td>33.0 ± 0.2</td><td>35.1 ± 0.8</td><td><strong>36.3 ± 1.2</strong></td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Positive wells*</td><td>3/3</td><td>3/3</td><td>3/3</td><td><strong>3/3</strong></td><td>0/3</td><td>0/3</td><td>0/3</td><td>0/3</td><td>0/3</td></tr><tr><td rowspan="3">P. knowlesi</td><td>Nested PCR</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Real Time PCR</td><td>29.0 ± 0.2</td><td>31.4 ± 0.2</td><td>34.0 ± 0.2</td><td>37.6 ± 0.3</td><td><strong>37.3 ± 0.5</strong></td><td>38.0 ± 0.5</td><td>38.6 ± NA</td><td>40.3 ± NA</td><td>-</td></tr><tr><td>Positive wells*</td><td>3/3</td><td>3/3</td><td>3/3</td><td>3/3</td><td><strong>3/3</strong></td><td>2/3</td><td>1/3</td><td>1/3</td><td>0/3</td></tr></tbody></table>

<table><caption>**Table 6.** Determination of the limits of detection for *Plasmodium* spp. real time PCR targeting 18S genes in whole blood. Results are presented in mean Ct (cycle threshold) ± standard deviation according to the level of parasitemia (p/µl, parasite/µl). In bold and italics, the LoD chosen. NA non-applicable. *Number of positive well in each real time triplicate.</caption><thead><tr><th>Parasitemia (p/µl)</th><th></th><th>100</th><th>20</th><th>10</th><th>1</th><th>0.5</th><th>0.25</th><th>0.12</th><th>0.06</th><th>0.03</th><th>0.015</th></tr></thead><tbody><tr><td rowspan="10">Plasmodium spp.</td><td>P. falciparum</td><td>26.4 ± 0.1</td><td>29.0 ± 0.3</td><td>30.0 ± 0.3</td><td><strong>33.0 ± 0.7</strong></td><td>34.0 ± 0.7</td><td>36.2 ± 1.2</td><td>37.2 ± NA</td><td>39.0 ± NA</td><td>-</td><td>-</td></tr><tr><td>Positive wells</td><td>3/3</td><td>3/3</td><td>3/3</td><td><strong>5/5</strong></td><td>4/5</td><td>3/5</td><td>1/5</td><td>1/5</td><td>-</td><td>-</td></tr><tr><td>P. vivax</td><td>25.5 ± 0.1</td><td>28.5 ± 0.2</td><td>29.0 ± 0.2</td><td><strong>33.0 ± 0.5</strong></td><td>34.4 ± 0.5</td><td>35.0 ± 0.5</td><td>36.3 ± NA</td><td>36.0 ± 0.1</td><td>36.5 ± NA</td><td>-</td></tr><tr><td>Positive wells</td><td>3/3</td><td>3/3</td><td>5/5</td><td><strong>5/5</strong></td><td>4/5</td><td>3/5</td><td>1/5</td><td>3/5</td><td>1/5</td><td>-</td></tr><tr><td>P. malariae</td><td>32.3 ± 0.3</td><td>34.6 ± 0.4</td><td><strong>34.2 ± 0.6</strong></td><td>35.1 ± 0.8</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Positive wells</td><td>3/3</td><td>3/5</td><td><strong>5/5</strong></td><td>2/5</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>P. ovale</td><td>28.0 ± 0.3</td><td>30.0 ± 0.2</td><td>31.4 ± 0.4</td><td><strong>34.2 ± 1.2</strong></td><td>36.2 ± 0.2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Positive wells</td><td>3/3</td><td>5/5</td><td>3/3</td><td><strong>5/5</strong></td><td>2/4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>P. knowlesi</td><td>25.6 ± 0.1</td><td>27.6 ± 0.1</td><td>30.3 ± 0.1</td><td>33.6 ± 0.5</td><td><strong>34.5 ± 1.0</strong></td><td>35.0 ± NA</td><td>35.5 ± 0.1</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Positive wells</td><td>3/3</td><td>3/3</td><td>2/2</td><td>2/2</td><td><strong>5/5</strong></td><td>1/5</td><td>3/5</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

**Detection of mixed infections.** Cross-reactions between species have been assessed in simplex amplification for each of the five-malaria species. No cross-reaction has been observed.

The impact of a multiplex strategy on PCR efficiency has been done for *P. vivax* amplification (Supplementary Fig. S1). Ct differences observed for each of the dilution points were between – 0.6% and + 1%. These variations were included within the coefficient of variation of the method. Therefore, the results were not impacted by the use of a combination of two or three primer pairs during the multiplex strategy.

The interferences between species have been assessed using artificially prepared mixed-infections from patient infected blood (Table 7) or using plasmids containing *P. falciparum* and *P. vivax* target sequence (Supplementary Table S2). Results from the two types of matrix showed an impact on the Ct values when one of the species was predominant. The identification of the minor one was therefore more difficult. For example, a *P. vivax* infection in minority among *P. falciparum* parasites could not be detected when the major species generated a Ct greater than 29 (Table 7). Based on our experience, in case of a Ct result below 27, it would be judicious to reanalyze each species of the multiplex in simplex to limit the risk of a hidden minor species. This should be defined according to the malaria epidemiology of each area and the probability of having mixed-infection. However, the limit of detection of the method was not impacted when the two species were far from the detection limit of the method. For thirteen samples, we also compared the capacity of real-time amplification or microscopy to detect mixed-infections. As expected, real-time PCR was more sensitive (Supplementary Table S3).
<table><caption>**Table 7.** Detection of *P. falciparum* and *P. vivax* in real mixed infection templates.</caption><thead><tr><th rowspan="3">Artificial blood mixed-infection</th><th rowspan="3">Respective parasitaemia</th><th colspan="3">Cycle threshold ± standard deviation</th></tr><tr><th colspan="2">Duplex-PCR reagents</th><th rowspan="2">Simplex reagents</th></tr><tr><th>P. falciparum</th><th>P. vivax</th></tr></thead><tbody><tr><td>P. falciparum-P. vivax</td><td>100 p/µl-100 p/µl</td><td>31.0 ± 0.3</td><td>31.5 ± 0.3</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>100 p/µl-10 p/µl</td><td>29.5 ± 0.4</td><td>33.0 ± 0.9</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>100 p/µl-1 p/µl</td><td>29.5 ± 0.5</td><td>-</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>10 p/µl-100 p/µl</td><td>32.0 ± 0.3</td><td>30.1 ± 0.2</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>10 p/µl-10 p/µl</td><td>32.3 ± 0.4</td><td>32.0 ± 0.1</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>10 p/µl-1 p/µl</td><td>32.0 ± 0.1</td><td>37.0 ± 0.3</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>1 p/µl-100 p/µl</td><td>41.5 ± 1.1</td><td>29.0 ± 0.3</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>1 p/µl-10 p/µl</td><td>35.0 ± 0.3</td><td>32.0 ± 1.5</td><td></td></tr><tr><td>P. falciparum-P. vivax</td><td>1 p/µl-1 p/µl</td><td>34.5 ± 0.3</td><td>36.3 ± 0.6</td><td></td></tr><tr><td>P. falciparum</td><td>100 p/µl</td><td>30.2 ± 0.5</td><td>-</td><td>30.0 ± 0.4</td></tr><tr><td>P. falciparum</td><td>10 p/µl</td><td>33.0 ± 0.6</td><td>-</td><td>33.0 ± 0.0</td></tr><tr><td>P. falciparum</td><td>1 p/µl</td><td>36.2 ± 0.6</td><td>-</td><td>36.4 ± 1.5</td></tr><tr><td>P. vivax</td><td>100 p/µl</td><td>-</td><td>30.2 ± 1.3</td><td>29.4 ± 0.1</td></tr><tr><td>P. vivax</td><td>10 p/µl</td><td>-</td><td>31.1 ± 0.2</td><td>31.0 ± 0.1</td></tr><tr><td>P. vivax</td><td>1 p/µl</td><td>-</td><td>36.0 ± 0.9</td><td>35.0 ± 0.5</td></tr></tbody></table>

**Robustness of the method.** A depth analysis of the possible reasons of deviation of the results was carried out (freezing/thawing DNA, primers and probe degradation; temperature variation ate the hybridization step ± 1%; pipetting variation during mix preparation ± 10...). The acceptable limit was set at 10%.

A series of 10 freezes/thaws of *P. falciparum* primers and probes has been tested. A drift of 6.3% was observed below the acceptable limit of 10%.

However, the number of freezing/thawing cycles had a major impact on DNA integrity of the samples, including the IQC. A drift of 15% was observed after 14 freezing/thawing cycles even when DNA was stored in aliquots of 30 µl (Fig. 2). Aliquoting IQC in small volume was therefore implemented to avoid several freezing/thawing cycles of each aliquot.

After choosing the best temperature for denaturation (95 °C) and amplification (60 °C), we determined the Maximum Permissible Error (MPE) of the StepOnePlus™ thermocycler in case of a temperature shift through time. The MPE is also required during metrological verification to establish the conformity of the thermocycler. The six parameters have been analyzed at several annealing temperature. By taking the plasmid as a template, we observe variable results with a reduction of Ct values for *P. falciparum*, *P. vivax* and *P. knowlesi* when the annealing temperature increased (Supplementary Fig. S2).

**Discussion**

According to WHO, the microscopic examination of blood smears stained with Giemsa is the gold standard for diagnosis¹. In a context of low transmission and elimination, the skills of microscopists are difficult to maintain and the asymptomatic carriers are difficult to detect. qPCR remains the most sensitive and species-specific methods compared to microscopy, ELISA or LAMP PCR¹⁶.

In this study, we implemented a robust, sensitive and rapid multiplex real-time PCR approach to detect the five human malaria species in venous blood collected on EDTA-vacutainer or capillary blood collected on filter paper in the context of the 15,189-accreditation requirements. The limit of detection (LoD) from DNA extracted on 200 µl of venous blood was 0.25 parasites/µl for *Plasmodium falciparum*, *P. vivax* and *P. knowlesi*, 1 parasite/µL for *P. ovale* and 5 parasites/µl for *P. malariae*. From 4-punches of dried blood spot, we achieve the limits of 5 parasites/µl for *Plasmodium vivax*, 20 parasites/µl for *P. falciparum*, *P. ovale* and *P. knowlesi* and 125 parasites/µl

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Impact of freezing/thawing of *Plasmodium falciparum* sample DNA in 30 µl-aliquots.</figcaption></figure>
for *P. malariae*. For active detection campaigns which could include thousands of samples, we also implement a real time PCR targeting the genus *Plasmodium* spp. to facilitate the first screening for malaria. The detection limits determined for this method were similar to the one observed with the species-specific method for each species. This first round of detection of *Plasmodium* spp. carriers aims to limit the time and cost spending during active case detection campaign. The cost of this method has been evaluated at 20$/sample analyzed in triplicate for *Plasmodium* spp. This cost is competitive compared to the 15–40$/sample analyzed in single well of the different commercial kits or previously published methods¹⁷. The species-specific approach to detect the five human malaria species was around 27$/sample (prices including DNA extraction and qPCR without human resources). However, qPCR will always be an expensive method compare to RDTs but the only performant one to detect the last parasite carriers.

This method has been designed to reduce costs but also sample usage. As such we propose a duplex *P. falciparum*/*P. vivax* and a triplex *P. malariae*/*P. ovale*/*P. knowlesi*, in a final volume of 12 µl, containing 1 µl of DNA, tested in triplicate for each sample. Generally, previously published real time PCR used a final volume chosen of 25 µl containing two to five microliters of DNA¹³,¹⁸,¹⁹. During implementation, 5 µl of DNA has been tested to increase sensitivities but the impact on the detection limits was not high enough comparing to the consumption of the matrix (data not shown). However, we observed a clear impact on the detection limits when the number of targets within each multiplex was reduced. Four species within the same well as Shokopoulos et al. described, generated higher interferences and a drastic drop of the detection limits. Even using a duplex approach, we lost in sensitivity in case of mixed infection. In our transmission settings, French Guiana, the probability of mixed infection *P. falciparum*/*P. vivax* is low, 0.6% ± 0.3 over the past four years (data not shown). We therefore considered to analyze in simplex the second target when the first one presented a Ct value below 27. The other possibility to avoid to reanalyze samples, is to adapt the content of each multiplex. Each laboratory implementing this method would have to define their own diagnosis flowchart according to the epidemiological situation of malaria and the number of mixed-infections in their area. During our validation process, mixed infections with low parasitemia of *P. vivax* (around 1 parasites/µl or below) were not detectable if a *P. falciparum* infection of about 100 parasites/µl was also present. Generally, a 100-fold difference in concentration could impact the detect limit for diagnosis of the minor species.

This validated method also include the amplification of the human β2Megaglobuline gene to control DNA extraction. With a minimum Ct value to validate the results, this allow to confirm the volume of blood analyzed as well as the different steps of the DNA extraction. If the Ct values are not correct, the DNA extraction is repeated. However, in case of leukopenia or hyperleukocytosis, the Ct value will stay outside the range even if the volume of blood is correct. Those cases are rare and additional information such as complete blood count information are required. However, the result for malaria diagnosis is correct.

Malaria elimination is defined by the WHO as the interruption of local transmission of a specific species, in a defined geographical area, following activities specifically carried out for this purpose¹. This multiplex real-time PCR method has been used during a mass screen and treat (MSAT) approach, Palustop, in French Guiana between October and December 2017 in the municipality of St Georges de l'Oyapock¹⁵. Using this real-time PCR, 1501 samples were analyzed from patients living in endemic area, 79% without symptoms and between 2 month and 93-year-old. The sex ratio was 0.88. One hundred positive samples have been identified, 90 with *Plasmodium vivax* and 10 with *Plasmodium falciparum*. No mixed-infection has been observed. The general prevalence of malaria using this method, has been determined at 6.6% [5.3–7.9], 74% of the people were asymptomatic. Among them, only nine samples out of 100 were positive by RDT¹⁵. These results show that RDT is not sensitive enough to detect asymptomatic carriers. In fact, in French Guiana, malaria has decreased for almost 15 years. Since 2005, *P. vivax* predominates but the proportion of *P. falciparum* on the Guiana Shield (region encompassing French Guiana, Guyana, Surinam, part of Brazil, Venezuela and Colombia) is one of the highest of South Americas. *P. malariae* transmission could sporadically occurred in the deep forest related to a wild fauna transmission with a very low prevalence. With a little more than one hundred cases in 2021, the territory targets malaria elimination for 2025. In this context, various operational research projects and active case detection have been initiated since 2017¹⁵,²⁰.

Recently, ultra-sensitive qPCR (us-qPCR) have been developed to detect asymptomatic low-density malaria infections in endemic areas²¹. The term ultra-sensitive refers to the use of genes with higher copy number (i.e. *pfvarATS* for *P. falciparum*, or *pvcox1* for *P. vivax*) and/or increasing the volume of blood to be extracted. These PCR generally target *P. falciparum* and *P. vivax* carriers and aim to detect asymptomatic carriers with low parasitemia acting as parasite reservoir²². However, these methods achieve sensitivities below the threshold of one gametocyte/µl of blood which is considered as the limit to transmit malaria to mosquitoes²³. When an individual is positive by qPCR, if his gametocytemia is below this threshold, he will not contribute to malaria transmission²⁴. However, fluctuation of this parameter is poorly known and a longitudinal follow-up of asymptomatic carriers would be interesting in order to describe gametocytemia through time in asymptomatic carriers. This knowledge would help us to define, the benefit of implementing one diagnosis method or the other according to their sensitivity.

## Conclusion

We established a real time PCR method that can detect and identify the five human malaria parasites (*P. falciparum*, *P. vivax*, *P. malariae*, *P. ovale* and *P. knowlesi*) from whole venous sample and dried blood spot using a small quantity of DNA extract. Composition of each multiplex amplification should be adapted to the local epidemiological situation including the prevalence of mixed-infections.

This PCR requires equipment and well-trained personnel. The DNA extraction step is time consuming and treatments will probably have to be prescribed in the following days. Despite these limitations, this PCR is
accredited in our lab according to the norm ISO 15,189 for medical biology which is strict. This quality process confirms the high performances, robustness and stability of the method. Furthermore, during active malaria detection campaigns the method has demonstrated its feasibility and usefulness in elimination settings.

## Materiel and methods

**Blood samples.** Blood samples were obtained from blood venipuncture and collected either from patients consulting for fever in a health structure in French Guiana, France mainland or during active case detection campaigns. The parasitemia estimation was determined by microscopy and thin smear. We observed at least 50 fields of 250 red blood cells on average and counted the number of parasitized red blood cells using the cell counter. If no parasitized red cells are identified, we observed up to 100 fields then perform a thick film. If parasites are observed, we calculated the percentage of parasitized red blood cells or the parasitaemia following WHO recommendations. The parasitemia was determined by two readers and eventually a third one in case of discrepancies above 20% differences.

**DNA extraction.** DNA was extracted from 200 µl of venous blood or four punches of capillary blood collected on filter paper using the QIAamp® genomic DNA kits (Qiagen) according to the manufacturer's instructions. DNA elution volume was 100 µL in AE buffer when 200 µl of whole blood was extracted and 75 µl of AE buffer when 4-punches of diameter 3 mm of dried blood spots were extracted. The same protocol was used for all samples and external quality samples.

**Primers and probes.** The primers and the probes (TaqMan) have been synthesized by Applied Biosystems and purified by high-performance liquid chromatography when labelled (Table 8). The β2-megaglobulin gene (β2MG) was used as an internal control of DNA extraction.

The targeted sequences of the 18 S rRNA gene of *P. falciparum*, *P. vivax*, *P. malariae*, *P. ovale* and *P. knowlesi* (Accession number M19172.1, U93233.1, M54897.1, KF696363.1 and PKNH_0320900 respectively) were cloned into five different pEX-K4 plasmids (Eurogentec, Belgium) as Internal Quality Control (IQC) to continuously monitor the performances of the methods. IQC were used at a concentration of 1000 copies/µl in order to obtain a Ct around 30. This value has been chosen to observe potential drifts of the real time PCR from one PCR to another. IQC results were analyzed using the standardized WESTGARD rules used in medical biology. Serial dilutions of the plasmids were realized in nuclease free water from 10¹² copies to 1 copy per µl in order to artificially constitute mixed-infections. Negative controls included TE buffer (10 mM Tris-HCl [pH 8], 1 mM EDTA).

<table><caption>**Table 8.** Primers and probes used for identification of *Plasmodium* species targeting 18S genes.</caption><thead><tr><th>Species</th><th>Primer or probe</th><th>Sequence 5'-3'</th><th>Fluorescent label</th><th>Reference</th></tr></thead><tbody><tr><td colspan="5">Duplex</td></tr><tr><td rowspan="3">P. falciparum</td><td>Fal-F</td><td>CCGACTAGGTGTTGGATGAAAGTGTTAA</td><td rowspan="3">VIC-MGBNFQ</td><td rowspan="3">Shokopes et al.<sup>12,18</sup></td></tr><tr><td>Fal-R</td><td>AACCCAAAGACTTTGATTTCTCATAA</td></tr><tr><td>Fal probe</td><td>AGCAATCTAAAAGTCACCTCGAAAGATGACT</td></tr><tr><td rowspan="3">P. vivax</td><td>Viv-F</td><td>CCGACTAGGCTTTGGATGAAAGATTTTA</td><td rowspan="3">NED-MGBNFQ</td><td rowspan="3"></td></tr><tr><td>Viv-R</td><td>AACCCAAAGACTTTGATTTCTCATAA</td></tr><tr><td>Viv probe</td><td>AGCAATCTAAGAATAAACTCCGAAGAGAAAATT</td></tr><tr><td colspan="5">Triplex</td></tr><tr><td rowspan="3">P. malariae</td><td>Mal-F</td><td>CCGACTAGGTGTTGGATGATAGAGTAAA</td><td rowspan="3">FAM-MGBNFQ</td><td rowspan="3">Shokopes et al.<sup>12,18</sup></td></tr><tr><td>Mal-R</td><td>AACCCAAAGACTTTGATTTCTCATAA</td></tr><tr><td>Mal probe</td><td>CTATCTAAAAGAAACACTCAT</td></tr><tr><td rowspan="3">P. ovale</td><td>Ova-F</td><td>CCGACTAGGTTTTGGATGAAAGATTTTT</td><td rowspan="3">VIC-MGBNFQ</td><td rowspan="3"></td></tr><tr><td>Ova-R</td><td>AACCCAAAGACTTTGATTTCTCATAA</td></tr><tr><td>Ova probe</td><td>CGAAAGGAATTTTCTTATT</td></tr><tr><td rowspan="3">P. knowlesi</td><td>Knw-F</td><td>CTAAAATGCGCACAAAGTCGAT</td><td rowspan="3">NED-MGBNFQ</td><td rowspan="3">De Canale et al.<sup>13</sup></td></tr><tr><td>Knw-R</td><td>GCAGTTAAACGCTCGTAGTTGAA</td></tr><tr><td>Knw probe-1</td><td>CGGAGGCATCAGTTAT</td></tr><tr><td colspan="5">Simplex</td></tr><tr><td rowspan="3">Plasmodium spp.</td><td>Spp-F</td><td>AGCTCTTTCTTGATTTCTTGG</td><td rowspan="3">VIC-MGBNFQ</td><td rowspan="3">Hassanpour et al.<sup>14</sup></td></tr><tr><td>Spp-R</td><td>CAGACAAATCATATTCACGAACT</td></tr><tr><td>Spp probe-1</td><td>AAACGGCCATGCATCACCAT</td></tr><tr><td colspan="5">Simplex</td></tr><tr><td rowspan="3">Human β2-megaglobulin</td><td>2βMG-F</td><td>TGAGTATGCCTGCCGTGTGA</td><td rowspan="3">FAM-MGBNFQ</td><td rowspan="3">Shokopes et al.<sup>12,18</sup></td></tr><tr><td>2βMG-R</td><td>ACTCATACACAACTTTCAGCAGCTTAC</td></tr><tr><td>2βMG probe</td><td>CCATGTGACTTTGTCACAGCCCAAGATAGTT</td></tr></tbody></table>
**Real time PCR.** The real-time PCR reaction consisted of 450 nM of each primer, 125 nM of probe, 5X TaqMan Universal Master Mix (Applied Biosystems, USA), and 1 µl of DNA in a 12 µl final volume. Reactions were performed on the StepOnePlus® (Applied Biosystems, USA). The following cycling conditions were applied: 50 °C for 2 min, 95 °C for 10 min, and 45 cycles of 95 °C for 15 s and 60 °C for 1 min. Fluorescence data was collected during the annealing/extension step at 60 °C.

Amplifications were performed in three separate reactions. The detection of *P. falciparum* and *P. vivax* in the first PCR and *P. malariae*, *P. ovale* and *P. knowlesi* in the second one was realized in triplicates to increase sensitivity and robustness of the method.

Cycle Threshold (Ct) values were analyzed by setting the threshold at 0.02 in order to harmonize the data from one experiment to another. The control of DNA extraction was realized using a simplex PCR which amplify the β2Megaglobulin Human gene. The manipulations were carried out in a single well and decision thresholds for venous blood was established at Ct ≤ 24 and for dried blood spot at Ct ≤ 33.

**Repeatability and reproducibility of the real time PCR.** Repeatability (intra-assay variation) have been done on two samples per species (Table 9). For *P. falciparum* and *P. vivax* we used patient's samples with parasitemia expressed with percentage. Due to a valuable and limited matrix, especially for *P. malariae*, *P. ovale* and *P. knowlesi* species, the number of repetitions was limited and External Quality Assessment (EQAs) with parasitemia expressed with parasites per microliters were used.

The reproducibility (inter-assay variation) was determined using the internal quality control, including around 80 values obtained within a range of 7 to 15 consecutive days. All samples were analyzed in triplicates. Coefficients of variation (CVs) were then calculated based on the Mean of the CTs of the triplicate.

**Analytical sensitivities, specificities and detection of mixed infections.** Sensitivity of the method to detect the five species by real-time PCR was evaluated by comparing the expected results with results obtained in EQAs. The specificities were first assayed in silico for the primers and probes by using BLAST software in order to avoid non-specific amplification. We prepared artificial mixed-infections (tenfold dilution factors) using blood of patients infected with *P. falciparum* and *P. vivax* or with plasmid containing the *P. falciparum* and *P. vivax* target sequences. These artificial mixed infections were analyzed in duplicates.

**Real-time PCR limits of detection.** The analytical evaluation was carried out on samples for *P. falciparum* and *P. vivax* or EQA with known parasitemia for *P. malariae*, *P. ovale*, *P. knowlesi*.

For *P. falciparum* and *P. vivax*, blood samples including at least 90% of synchronous and in well-conditions rings were diluted in uninfected whole blood to achieve a concentration range of $1 \times 10^3$ p/µl to 0.005 p/µl before extraction.

External quality controls were provided by WHO through UKNEQAS. The EQAs are received in two types of matrix: lyophilized whole blood and dried blood spot. The lyophilized whole blood was reconstituted in 500 µl of sterile water and extracted with QIAamp® genomic DNA kits (Qiagen). To create the same concentration ranges for *P. malariae*, *P. ovale* and *P. knowlesi*, DNA was directly diluted based on the initial parasitemia of the EQA after considering the concentration factor at the elution step of DNA extraction. Each sample or EQA was tested by real time PCR in 3 wells per parasitemia.

**Robustness and stability of reagents and DNA.** The stability of the reagents has been evaluated using a succession of ten cycles of freezing / thawing for probe/primer set and fifteen cycles for *P. falciparum* DNA sample with a parasitemia of 2.1% and diluted to 1/10,000. A robustness test was also done to determine the Maximum Permissible Error (MPE) of the StepOnePlus™ machine. Several pairs of temperatures were tested for denaturation (94 and 95 °C) and hybridization/elongation (59 to 63 °C). The plasmids of the five spices was used and three negative venous blood sample was used to test the internal control β2-megaglobulin.

**Ethical consideration.** For the present research, authors confirm that all experiments were performed in accordance with French regulations. The parasitology laboratory is receiving samples according to its mandate of National Reference Center since 1989 (Décret no 2016-806 du 16 juin 2016 relatif aux centres nationaux de référence pour la lutte contre les maladies transmissibles, for the 2023-2027 period). According to the French legislation (article L.1211-2 and related of the French Public Health Code), secondary use for scientific purpose of human clinical remaining samples collecting during care or surveillance are possible as long as the corresponding patients are informed and has not given any objection to them. To fulfill this requirement the Hospital brochure informed every patient through a brochure entitled "Information for patients" which detailed the opposition procedure. Solely the samples without immediate and delayed patient opposition, and/or from their legal guardians, have been included in this study.

<table><caption>**Table 9.** Parasitemia of samples used to evaluate the repeatability of the method. *p*/µl parasites/µl of blood.</caption><thead><tr><th></th><th><i>P. falciparum</i></th><th><i>P. vivax</i></th><th><i>P. malariae</i></th><th><i>P. ovale</i></th><th><i>P. knowlesi</i></th></tr></thead><tbody><tr><td>High-parasitemia</td><td>5.6%</td><td>2.2%</td><td>125 p/µl</td><td>200 p/µl</td><td>2000 p/µl</td></tr><tr><td>Low-parasitemia</td><td>0.005%</td><td>0.01%</td><td>20 p/µl</td><td>5 p/µl</td><td>1 p/µl</td></tr></tbody></table>
Authorization associated to the samples collected during the active case detection campaign, Palustop, was approved by the Comité de Protection des Personnes du Sud-Ouest et Outre-Mer 4 N° AM-36/1/CPP15-024 (French ethics committee). The database was anonymized and declared to the French Regulatory Commission (Commission Nationale Informatique et Libertés, CNIL, n°917186). Those samples collected by the National Reference Center were registered by the French Ministry for Research (declaration number DC-2010-1223; collection N°2). Solely the samples without immediate and delayed patient opposition, and/or from their legal guardians, have been included in this study.

## Data availability

The authors declare that the data supporting the findings of this study are available within the paper and its supplementary information.

Received: 22 February 2023; Accepted: 11 July 2023

Published online: 14 July 2023

## References

1. Who Health Organization. World Malaria Report 2022. https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2022 (2022).

2. Cowman, A. F., Healer, J., Marapana, D. & Marsh, K. Malaria: Biology and disease. Cell **167**, 610–624. https://doi.org/10.1016/j.cell.2016.07.055 (2016).

3. Barber, B. E., Grigg, M. J., William, T., Yeo, T. W. & Anstey, N. M. The treatment of *Plasmodium knowlesi* malaria. *Trends Parasitol.* **33**, 242–253. https://doi.org/10.1016/j.pt.2016.09.002 (2017).

4. Baird, J. K. Evidence and implications of mortality associated with acute *Plasmodium vivax* malaria. *Clin. Microbiol. Rev.* **26**, 36–57. https://doi.org/10.1128/CMR.00074-12 (2013).

5. Poespoprodjo, J. R. et al. Vivax malaria: A major cause of morbidity in early infancy. *Clin. Infect. Dis.* **48**, 1704–1712. https://doi.org/10.1086/599041 (2009).

6. Zheng, Z. & Cheng, Z. Advances in molecular diagnosis of malaria. *Adv. Clin. Chem.* **80**, 155–192. https://doi.org/10.1016/bs.acc.2016.11.006 (2017).

7. Snounou, G. et al. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. *Mol. Biochem. Parasitol.* **61**, 315–320. https://doi.org/10.1016/0166-6851(93)90077-B (1993).

8. WHO. 2014: Note d'orientation Sur le Diagnostic du Paludisme.pdf. (2014).

9. Mercereau-Puijalon, O., Barale, J. C. & Bischoff, E. Three multigene families in *Plasmodium parasites*: Facts and questions. *Int. J. Parasitol.* **32**, 1323–1344. https://doi.org/10.1016/s0020-7519(02)00111-x (2002).

10. Bousema, T., Okell, L., Felger, I. & Drakeley, C. Asymptomatic malaria infections: Detectability, transmissibility and public health relevance. *Nat. Rev. Microbiol.* **12**, 833–840. https://doi.org/10.1038/nrmicro3364 (2014).

11. Berry, A., Iriart, X. & Magnaval, J.-F. Nouvelles méthodes de diagnostic du paludisme. *Rev. Franc. Lab.* **2009**, 65–70. https://doi.org/10.1016/S1773-035X(09)70252-7 (2009).

12. Shokopoulos, S. E., Ndao, M., Kowalewska-Grochowska, K. & Yanow, S. K. Multiplexed real-time PCR assay for discrimination of *Plasmodium* species with improved sensitivity for mixed infections. *Methods BM* **47**, 975–980. https://doi.org/10.1128/JCM.01858-08 (2009).

13. De Canale, E. et al. Plasmodium knowlesi malaria in a traveller returning from the Philippines to Italy, 2016. *New Microbiol.* **40**, 291–294 (2017).

14. Hassanpour, G. et al. Simplified pan-species real-time PCR-based detection of *Plasmodium* Spp in blood smear. *Iran J. Parasitol.* **11**, 463–470 (2016).

15. Mosnier, E. et al. Prevalence of *Plasmodium* spp. in the Amazonian border context (French Guiana-Brazil): Associated factors and spatial distribution. *Am. J. Trop. Med. Hyg.* **102**, 130–141. https://doi.org/10.4269/ajtmh.19-0378 (2020).

16. Morris, U. & Aydin-Schmidt, B. Performance and application of commercially available loop-mediated isothermal amplification (LAMP) kits in malaria endemic and non-endemic settings. *Diagnostics* https://doi.org/10.3390/diagnostics11020336 (2021).

17. Zelman, B. W. et al. Costs and cost-effectiveness of malaria reactive case detection using loop-mediated isothermal amplification compared to microscopy in the low transmission setting of Aceh Province, Indonesia. *Malar. J.* **17**, 220. https://doi.org/10.1186/s12936-018-2361-y (2018).

18. Shokopoulos, S. E., Ndao, M., Kowalewska-Grochowska, K. & Yanow, S. K. Multiplexed real-time PCR assay for discrimination of *Plasmodium* species with improved sensitivity for mixed infections. *J. Clin. Microbiol.* **47**, 975–980. https://doi.org/10.1128/JCM.01858-08 (2009).

19. Rougemont, M. et al. Detection of four *Plasmodium* species in blood from humans by 18S rRNA gene subunit-based and species-specific real-time PCR assays. *Methods BM Détect. SNP* **42**, 5636–5643. https://doi.org/10.1128/JCM.42.12.5636-5643.2004 (2004).

20. Douine, M. et al. Malakit: An innovative pilot project to self-diagnose and self-treat malaria among illegal gold miners in the Guiana Shield. *Malar. J.* **17**, 158. https://doi.org/10.1186/s12936-018-2306-5 (2018).

21. Danwang, C., Kirakoya-Samadoulougou, F. & Samadoulougou, S. Assessing field performance of ultrasensitive rapid diagnostic tests for malaria: A systematic review and meta-analysis. *Malar. J.* **20**, 245. https://doi.org/10.1186/s12936-021-03783-2 (2021).

22. Hofmann, N. E. et al. Assessment of Ultra-Sensitive Malaria Diagnosis Versus Standard Molecular Diagnostics for Malaria Elimination: An in-Depth Molecular Community Cross-sectional Study (Elsevier, 2018).

23. Gruenberg, M. et al. qRT-PCR versus IFA-based quantification of male and female gametocytes in low-density plasmodium falciparum infections and their relevance for transmission. *J. Infect. Dis.* **221**, 598–607. https://doi.org/10.1093/infdis/jiz420 (2020).

24. Meibalan, E. et al. *Plasmodium falciparum* gametocyte density and infectivity in peripheral blood and skin tissue of naturally infected parasite carriers in Burkina Faso. *J. Infect. Dis.* **223**, 1822–1830. https://doi.org/10.1093/infdis/jiz680 (2021).

## Acknowledgements

The National Reference Center for Malaria in French Guiana acknowledges its partners involved in diagnosis and care of malaria in French Guiana, who collaborate for several years in malaria surveillance and allow generation of precious data for public health and malaria control. The authors are grateful to Valentin Joste for reviewing the manuscript and fruitful discussion. This work was funded by Santé publique France (French Ministry of Health), European Funds for Regional Development, N° Synergie: GY0012082, and the French Ministry for National Education Higher Education and Research.
**Author contributions**

Y.L. and L.M. contributed to the study design. Y.L., S.P. and L.M. oversaw the conduct of the study. Y.L., B.V. and M.D performed the analysis. S.H. and E.M. contributed to the sample collection. Y.L. and L.M. wrote the original draft. All authors critically reviewed the manuscript during writing. C.F., S.H., S.P., E.M. and L.M. reviewed and edited the manuscript.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-023-38621-9.

**Correspondence** and requests for materials should be addressed to Y.L.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023