A portable image-based cytometer for rapid malaria detection and quantification

Abstract

Increasing resistance by malaria parasites to currently used antimalarials across the developing world warrants timely detection and classification so that appropriate drug combinations can be administered before clinical complications arise. However, this is often challenged by low levels of infection (referred to as parasitemia) and presence of predominantly young parasitic forms in the patients’ peripheral blood. Herein, we developed a simple, inexpensive and portable image-based cytometer that detects and numerically counts Plasmodium falciparum infected red blood cells (iRBCs) from Giemsa-stained smears derived from infected blood. Our cytometer is able to classify all parasitic subpopulations by quantifying the area occupied by the parasites within iRBCs, with high specificity, sensitivity and negligible false positives (~ 0.0025%). Moreover, we demonstrate the application of our image-based cytometer in testing anti-malarial efficacy against a commercial flow cytometer and demonstrate comparable results between the two methods. Collectively, these results highlight the possibility to use our image-based cytometer as a cheap, rapid and accurate alternative for antimalarial testing without compromising on efficiency and minimal processing time. With appropriate filters applied into the algorithm, to rule out leukocytes and reticulocytes, our cytometer may also be used for field diagnosis of malaria.

About the Authors:

Dahou Yang

Contributed equally to this work with: Dahou Yang, Gowtham Subramanian

Affiliation: Pillar of Engineering Product Development, Singapore University of Technology and Design, Singapore, Singapore

Gowtham Subramanian

Contributed equally to this work with: Dahou Yang, Gowtham Subramanian

Affiliation: Pillar of Engineering Product Development, Singapore University of Technology and Design, Singapore, Singapore

Jinming Duan

Affiliation: School of Computer Science, University of Nottingham, Nottingham, United Kingdom

Shaobing Gao

Affiliation: School of Life Science and Technology, University of Electronic Science and Technology of China, Chengdu, China

Li Bai

Affiliation: School of Computer Science, University of Nottingham, Nottingham, United Kingdom

Rajesh Chandramohanadas

* E-mail: [email protected] (RC); [email protected] (YA)

Affiliation: Pillar of Engineering Product Development, Singapore University of Technology and Design, Singapore, Singapore

ORCID http://orcid.org/0000-0002-8760-0396

Ye Ai

* E-mail: [email protected] (RC); [email protected] (YA)

Affiliation: Pillar of Engineering Product Development, Singapore University of Technology and Design, Singapore, Singapore

Introduction

Malaria, one of the most devastating infectious diseases around the globe, is caused by protozoan parasites of the genus Plasmodium. There are five major species of plasmodia that infect humans, out of which Plasmodium (P.) falciparum causes the majority of morbidity and mortality in Africa followed by less lethal P. vivax infections across South-East Asia, altogether infecting 200 million people and resulting in over half a million deaths every year [1]. Currently used antimalarials include chloroquine, artemisinin sulfadoxine-pyrimethamine combination, atovaquone and clindamycin. However, parasites have acquired resistance to most of the above-mentioned drugs both in Africa and South-East Asia, rendering them inefficient for future usage [2–4].

Early diagnosis and treatment are required to avoid anemia, organ failure [5] and malaria-associated deaths [6]. Lack of reliable methods and tools in the field settings provide tremendous impedance to early diagnosis in malaria endemic areas such as sub-Saharan Africa. Traditional and widely practised method of malaria diagnosis relies on observable clinical symptoms associated, which are results of general host response to an infection. Undoubtedly, this method foists on several challenges due to the nonspecific nature of the symptoms that could be potentially caused due to other immune challenges, which could ultimately result in inappropriate and unnecessary exploitation of antimalarials [7, 8]. Rapid Diagnostic Tests (RDTs) used in clinics measure the presence of parasitic antigens such as aldolase or lactate dehydrogenase [9, 10]. However, these methods rely on detection of antigens derived from the parasites, rather than detecting the parasites themselves. Though RDTs are often able to differentiate most malarial species by means of their antigenic properties, overall sensitivity of detection is far below the threshold of microscopy-based malaria detection, exhibiting huge variations among the patients [11–15] and failing to guide treatments which could be fatal when mixed infections occur [16]. These are efficient methods; however the results can vary depending on the severity of infection and the high occurrence of young parasitic forms (ring-stage infections) in the peripheral blood, which are only beginning to establish metabolic processes. Due to the recent advancements in technology, several malaria diagnostic techniques such as microarray [17], PCR [18], loop-mediated isothermal amplification (LAMP) [19], flow cytometry [20], hemozoin detection using automated hematology analyser [21] have been developed for efficient malaria diagnosis. Diagnostic methods leveraging PCR, 18s-rRNA detection [22], mitochondrial cytochrome b activity [23, 24], PgMt19 and PfMT869 mitochondrial regions [25], and the Pvr47 and Pfr364 genes [26, 27], have been used for detecting Plasmodium species. However, quantitative amplification of genes demands careful processing of blood to remove the inhibitors of amplification and thermal cycling, making it a cumbersome procedure. Notably, all of these methods demand expensive laboratory facilities and highly trained personnel to conduct complicated analytical procedures and data analyses, which may not be feasible in resource-poor malaria endemic countries.

The current gold standard technique for detecting malaria is microscopic examination of Giemsa-stained thin and thick blood smears [28] both in the field as well as in laboratory. This method allows detection of densely stained parasites against a background of lightly stained RBCs and widely accepted owing to cost effectiveness, simplicity and rapidity. Though, this method detects parasitemia levels of up to 1 infected cell in 10⁶ cells [29], the microscopic examination requiring traditional bulky microscope is laborious and often fails when the parasitemia is low, a situation very common in the case of P. vivax infections. In addition, logistic issues and challenges associated with transporting the traditional bulky microscopes to remote and rural malaria endemic regions [30] and fulfilling operational and maintenance requirements could be challenging. Furthermore, manual counting from Giemsa-smears is known to vary depending on the personnel engaged and quality of smears involving a high chance of misinterpreting other microorganisms such as bacteria and fungi as Plasmodium parasites and difficulties associated with identifying different Plasmodium strains [31]. Previously, several laboratories explored image processing based automated cell counting for parasitemia estimation, but most of them suffered from high false positive values [32–34], inferior accuracy [35], inability to differentiate parasitic stages [36] primarily rings [37] and the requirement of fluorescent dyes [38, 39].

To address these problems, we have developed a low-cost, portable image-based cytometer capable of detecting parasitic infections at levels as low as 0.2%. Additionally, our system allows reliable classification of all parasitic life stages at different magnifications. Finally, we demonstrate that the newly developed image-based cytometer can be adopted for antimalarial screening with high accuracy and efficiency and thus can be routinely exploited as a powerful alternative to expensive, technically demanding procedures such as traditional flow cytometry.

Materials and methods

Blood collection and parasite culture

Blood used in this project was purchased from Interstate Blood Bank. Before culturing malaria parasites, blood was transferred to EDTA Tubes (VACUETTE^® EDTA Tubes, Greiner Bio-One), washed three times in RPMI 1640 (Sigma-Aldrich) by centrifuging at 600g for 10 minutes to remove the buffy coat. RBCs were stored at 50% hematocrit in malaria culture medium (MCM), which is RPMI supplemented with bovine serum (Albumax II, Gibco-Singapore), 2.5 μg/mL gentamycin and hypoxanthine. Standard laboratory strains of P. falciparum, 3D7, were used in all experiments. Parasites were cultured in human O⁺ erythrocytes at 2.5% haematocrit in MCM under standard conditions [40]. Parasites were synchronized by magnetic selection at late stages (46–48 h) with the help of a SuperMACS magnet (Miltenyi Biotech, Bergisch Gladbach, Germany) and introduced into fresh RBCs, followed by sorbitol synchronization three hours later to allow a tight window of invasion [41].

Growth inhibition assays and parasitemia determination by flow cytometry and microscopic examination

Parasites at trophozoite stage (24–26 hpi) were treated with a range of concentrations of antimalarial drugs at 2% parasitemia and 2.5% hematocrit. Samples were cultured for 50–52 h until the trophozoites (28–30 hpi) appeared in the next replication cycle and then harvested for further analysis.

Flow cytometry was performed on Accuri C6 flow cytometer (Becton, Dickinson, CA. USA). At least 100,000 RBCs were analysed for each sample to determine parasitemia. To do this, 50 μl aliquots were collected from each sample and fixed with 0.1% glutaraldehyde (Sigma-Aldrich) at 4°C overnight. Fixed cells were collected by centrifuging at 400 g for 3 min, washed in PBS, permeabilized using 0.25% Triton X-100/PBS (Sigma-Aldrich) for 10 min at room temperature and washed again. Further, samples were incubated with 25 μg/ml Hoechst 33342 (Thermo Fisher) for 30 min in dark, and then parasitemia was scored by flow cytometry as reported elsewhere [42]. Data analysis and statistics were performed using GraphPad Prism (GraphPad Software) in accordance with the recommended protocol for nonlinear regression of a log-(inhibitor)-versus-response curve.

Thin smears for microscopic examination were prepared on glass slides and fixed with 100% methanol (Merck). The smears were stained with freshly prepared, filtered 1 in 10 dilution of Giemsa (Sigma-Aldrich) solution in de-ionised water for 10 min. Parasites were observed under an oil immersion lens (100X) using an optical microscope (Leica) for standard experiments.

Design of in-house image-based cytometer and image acquisition

A low-cost, portable image-based cytometer was built for image acquisition from Giemsa-stained blood smears. As illustrated in Fig 1A, it consists of a fine-focus (Z axis) adjustment platform (Dino-Lite, Model RK-10 Rack, Singapore), a 14 MP colour camera (ToupTek, P/N:TP114000A, Hangzhou, China), 20X objective or 100X oil immersion objective, two-dimensional translation stage and a white light source with tunable light intensity. Giemsa-stained smears were prepared carefully to obtain uniform cell distribution and placed on the 3D-printed sample holder for imaging. Images were captured with a commercial software (ToupTek, Hangzhou, China) by manually translating the smear slide on the stage. Only fields containing obvious debris were discarded. A 20X objective was used to capture large fields for imaging at high-throughput. However, the use of 20X objective is not sufficient to identify ring stage, in this case, a 100X oil immersion objective was utilized for ring forms imaging and parasitemia analysis.

[Figure omitted. See PDF.]

Fig 1. The developed image-based cytometer.

(A) Schematic experimental setup of the cytometer. (B) Photograph of the image-based cytometer. (C) An image of the blood smear taken by the cytometer.

https://doi.org/10.1371/journal.pone.0179161.g001

Image analysis algorithm

We used image processing and machine learning algorithms to analyse the smear images captured by the developed image-based cytometer. The image analysis algorithm can be divided into three key steps. First, a pre-processing total generalized variation (TGV) denoising method was used to remove the unavoidable noise in the original images. Second, an efficient and robust local adaptive thresholding approach was used to segment the smoothed image with possible intensity inhomogeneity. Finally, a machine learning method was performed to prevent false detections in the case when the image does not contain any parasites. Details of each step will be discussed in the following.

The denoising of the original image was implemented by the TGV method [43], which is a powerful image pre-processing tool that has been extensively used in image processing community [44–46]. The TGV regularisation has the capability of representing image characteristics up to an arbitrary order of differentiation (piecewise constant, piecewise affine, piecewise quadratic etc.). The TGV model is of the following form(1)where E(u,p) means the energy functional with respect to two variables u and p, u is the denoised/smoothed image, p = (p1,p2) is a symmetrised gradient vector which is closely related to , f is the input original image with intrinsic noise, and the operator ε(p) is the symmetrised derivative which is defined as 0.5(∇p + ∇p^T). The first energy term on the right hand side of Eq 1 is the data fidelity term which constrains the smoothed image u to be similar to the original image f. The second term preserves edges of objects in u and removes the noise from f in the meantime. The third term imposes smoothness on u and also eliminates the staircase artefact produced by the second term. The minimum of Eq 1 is taken over all the symmetrised gradient of the deformation field p = (p1,p2) on the image space Ω. The positive λ coefficients and α in Eq 1 balance the data fidelity term (the first energy term), the first order derivative (the second energy term) and the second order derivative (the third energy term). The value of α is normally set to 2 and this setting is suitable for most applications [44] and does not need to be tuned. In our experiment, we therefore set α = 2 and only vary λ to achieve different smoothness scales.

As TGV takes both first and second derivatives into consideration, it can highlight the edges of objects and smooth image without creating additional artefacts. This property also enables TGV to overcome the problem of intensity inhomogeneity that has been found to widely exist in our captured images. Moreover, since the TGV functional of Eq 1 is convex, it guarantees a global optimal solution as well as allows the use of powerful modern optimization techniques. Finally, TGV is translationally invariant and rotationally invariant, meaning that the denoising results are not affected by the viewpoint of the images taken from different angles. In spite of the outstanding performance of TGV for image denoising, it is difficult to minimise such functional due to its non-linear and non-smooth nature. In this study, we minimise TGV with a fast numerical algorithm based on split Bregman.[47] For a complete implementation on this algorithm, we refer the reader to a previous study [48].

The smoothed image can be then segmented via a simple local adaptive thresholding approach [49], which is given as(2)where T is the binary image, I = u in this case, and ls(I,ws) means that I is convolved with a suitable operator, i.e. the mean or median filter. ws is the window size of the filter and C is a user-defined threshold value. The adaptive thresholding produces binary segmentation with fewer isolated points, giving a better result than a simple high-pass threshold. It is worth mentioning that this approach can segment large sized images (e.g. 2000 × 2000 pixels) with real-time computational speed, making it an ideal segmentation tool for the images used in the experiments. In addition, previous studies have shown that it is robust against inhomogeneity in medical images and can also obtain higher accuracy [50].

To prevent false detections in images that do not contain any parasites, we performed the machine-learning algorithm to automatically classify the parasite and non-parasite images. Specifically, 50 images were first selected manually as a training dataset, and divided into two groups, one containing 25 parasite images and the other group with 25 non-parasite images. The selected representative images were then segmented separately. For each RGB channel of the segmented image, the average intensity of all the pixels that fall in the segmented cells was calculated, thus forming a three-dimensional feature that can be effectively used for training a precise classification model. In this study, we choose Support Vector Machine (SVM) classifier with the linear kernel [51] to train these extracted features. The SVM classifies two different datasets by finding an optimal hyperplane that has the largest margin distance between them. It is more robust and accurate than other machine learning techniques. We therefore used the trained model by the SVM to automatically classify the rest of parasite and non-parasite images.

Total RBC number and infected RBC number quantification process

As shown in Fig 2, the flowchart of the algorithm is composed of three sections: 1) to estimate the total number of all RBCs; 2) to estimate the number of infected RBCs (iRBCs); and 3) to classify distinct developmental stages of the parasite. For images acquired under objective 20X, we used the mean filter with the window size ws = 120 and set C = 0.03, λ = 3 for the segmentation of all RBCs, which produces binary segmentation with isolated RBCs and clustered RBCs. The number of RBCs in the clustered region was estimated by dividing the whole clustered area by the average area of a single RBC. For iRBCs segmentation, a high threshold value C = 0.2 was used with ws = 100, λ = 3. Similarly, under 100X oil immersion objective, ws = 300, C = 0.01 and λ = 5 were used for RBCs segmentation and ws = 300, C = 0.15 and λ = 10 were set for iRBCs segmentation. It should be noted that the C value can be affected by the light intensity. After standardizations, we kept a constant, optimum light intensity and fixed exposure time of camera for our further assays. To address the possible risk of counting multiply infected parasites in one single iRBC as multiple iRBCs, we performed the following calculation: if the distance between two segmented parasites is less than 150 pixels, these parasites will be considered as single iRBCs. It might still cause the false classification when two overlapping iRBCs are observed, but this was an extremely rare case due to the careful preparation of thin and uniform smear slides.

[Figure omitted. See PDF.]

Fig 2. Flowchart of the developed image analysis algorithm.

Image A was taken by objective 20X and image B was taken by oil immersion objective 100X.

https://doi.org/10.1371/journal.pone.0179161.g002

Classification of parasites’ developmental stage

After estimating the number of iRBCs, the algorithm performed the stage classification based on the area occupied by individual parasites within an iRBC. Fig 2 depicts the criteria of parasites’ stage classification under 20X and 100X magnifications. For 100X oil immersion objective, if the area occupied by the parasite is larger than 2830 pixels and less than 4500 pixels, then it was identified as a ring-stage parasite. It was classified as trophozoite stage when the area is greater than 4500 pixels and less than 16000 pixels. If the parasite-occupied area is greater than 16000 pixels and less than 48000 pixels, it was considered as a schizont. For 20X objective, the magnification was insufficient to image the rings, thus we only identify the two late stages under this low magnification. All these threshold pixel numbers for stage classification were verified independently based on conventional cytological examination. The developed system was used to determine parasitemia from seven independent samples. Each sample set included three smears and a minimum number of 30 fields were captured continuously resulting in counting at least 4000 RBCs for each smear. Imaging 30 fields for each smear took about 10 minutes, meaning it took 20 seconds to image each field.

Results and discussion

Image analysis results under 20X objective

The 20X objective was used for high-throughput imaging, and the analyses of images captured by our image-based cytometer were performed automatically using our image analysis algorithm, as shown in Fig 3. Table 1 shows the results from six randomly selected images from the same smear by using the image analysis algorithm and naked eye manual counting. Errors of total RBCs, iRBCs and parasitic classification were less than 5%. These negligible errors of total RBC number were primarily in regions where RBCs were clustered, minimally compromising the counting accuracy. Estimation of the number of iRBCs and life cycle stages were mildly influenced by debris. Nevertheless, error values less than 5% showed the excellent performance and reliability of our image-based cytometer compared to manual counting.

[Figure omitted. See PDF.]

Table 1. Comparison of cell counting and developmental stage classification (Objective 20X).

https://doi.org/10.1371/journal.pone.0179161.t001

[Figure omitted. See PDF.]

Fig 3. Images automatically analyzed by the developed image-based cytometer.

(A) Original image taken by objective 20X. (B) Representation of all the extracted RBCs. (C) Extracted schizont stage iRBCs. (D) Extracted trophozoite stage iRBCs.

https://doi.org/10.1371/journal.pone.0179161.g003

Image analysis results under 100X oil immersion objective

Since 20X objective is insufficient to image the ring stages with high accuracy, a 100X oil immersion objective was used for parasitemia detection and classification. As illustrated in Fig 4, image-based cytometer equipped with 100X allowed us to reliably classify all the three parasitic stages. Table 2 shows the inspection results of RBCs, iRBCs segmentation and stage-specific classification. The ability of our proposed cytometer to classify and document different developmental stages may lead to potential pharmacological applications in identifying the potential blockers of egress, invasion and parasitic development inhibitory molecules [38]. In a single image, errors of a specific parasite stage may be large (due to limited number of parasites in the field), but the total error of randomly selected 6 images was well below 5%, which shows the robustness of our system. Giemsa-stained images taken at 20X and 100X objective lens in the aforementioned studies are shown in S1 and S2 Figs of supplementary information, respectively.

[Figure omitted. See PDF.]

Table 2. Comparison of cell counting and developmental stage classification (Immersion objective 100X).

https://doi.org/10.1371/journal.pone.0179161.t002

[Figure omitted. See PDF.]

Fig 4. Images automatically analyzed by the developed image-based cytometer.

(A) Original image taken by immersion objective 100X. (B) All the extracted RBCs. (C) All the extracted iRBCs. (D) Extracted ring stage. (E) Extracted trophozoite stage. (F) Extracted schizont stage.

https://doi.org/10.1371/journal.pone.0179161.g004

Estimation of false positives

There are several factors that could account for the sensitivity of the microscopic examination of thin blood smears such as quality of the smears (uniform dispersion of single cells vs clump of cells), duration and dilution of Giemsa staining, and debris and dirt arising from Giemsa solution [52]. Previous reports have shown a huge variation in parasitemia levels between 6.3% and 45.8% based on 48 microscopic slides circulated between four laboratories in Mpumalanga province of South Africa [53]. We therefore set up trial experiments to standardize the parameters such as the dilution and duration of Giemsa staining, uniform distribution of RBCs to reduce and prevent the false positive values. From our trial experiments, we noticed that uniform distribution of RBCs stained with 1 in 10 dilution of Giemsa in distilled water for 10 minutes and washed with distilled water showed the best results. To determine the false positive value of image-based cytometer, healthy RBCs (hRBCs) were diluted in MCM and same number of RBCs (ranging from 2000 to 10000) was scored using both flow cytometer and our image-based cytometer. As expected, there was an absolute zero detection of false positive value in the flow cytometer and the corresponding flow plots are shown in Fig 5B and the image-based cytometer showed 0.0025% false positive detection for 8000 and 10000 cells. As shown in Fig 5A, there were no false positives detected using image-based cytometer when calculating cells less than 5000. In order to calculate the false positives, we performed two independent experiments in duplicates. In one of the replicates of the first independent experiment, we noticed 1 false positive while counting 8000 and 10000 cells, which accounts for 0.0125% and 0.01% false positives. We rounded off our results to decimal places, getting 0.01% false positive value and 0.0025% average false positive value in total for two independent experiments in duplicates and the differences between the values of image-based cytometer and flow cytometer were found to be nonsignificant using two tailed Wilcoxon T test. The result obtained from this experiment was evaluated using Bland-Altman method of comparison and it is shown in Fig 6A. It should be noted that the values are within the 95% limits of agreement (lower 95% limit of agreement = -0.004079 and upper 95% limit of agreement = 0.001579).

[Figure omitted. See PDF.]

Fig 5. Comparison of false positives in hRBCs using image-based cytometer and flow cytometer.

(A) Briefly, freshly drawn RBCs were counted and scored using image-based cytometer and flow cytometer, respectively. There were no false parasites obtained until 5000 RBCs were scored using image-based cytometer. Our developed cytometer exhibited 0.0025% false positive values while calculating more than 5000 cells. The figure represents the average of two independent experiments performed in duplicates and the error bars represent the standard deviation of the average values. (B) The figure represents the histogram plots for false positives estimated using flow cytometer while counting 2000 (a), 5000 (b), 8000 (c) and 10000 (d) events.

https://doi.org/10.1371/journal.pone.0179161.g005

[Figure omitted. See PDF.]

Fig 6. Evaluation of the performance and reliability of the image-based cytometer using Bland-Altman method of comparison.

The differences between the false positives of flow cytometer and image-based cytometer versus their average is plotted in (A). Ratio of parasitemia detected through flow cytometer and image-based cytometer and ratio of IC₅₀ values identified through flow cytometer and image-based cytometer are shown in (B) and (C) respectively. All values are within the 95% limits of agreement, which are shown as two dotted lines in the figures.

https://doi.org/10.1371/journal.pone.0179161.g006

Pharmacological testing using the developed image-based cytometer

Reliable, easy and affordable methods of inexpensive and high throughput screening techniques are needed for quantitative evaluation of new classes of antimalarials. Quantification of radioactive substance uptake such as [³H]hypoxanthine [54], [³H]isoleucine [55] and [³H]ethanolamine [56] by malaria parasites in the presence of the drugs have been measured to assess the drug efficacy. Several tools including flow cytometry and fluorescence microscopy are commonly employed for antimalarial tests in research laboratories. Here, we investigated the application of our image-based cytometer in pharmacological assays using three well-established antimalarials: chloroquine, artemisinin and cycloheximide in comparison to a more laborious flow cytometry method. These three antimalarials were selected since they act through distinct mechanisms and halt the parasitic growth and development under different stages [57].

Prior to proceeding with drug testing, we evaluated the potency of our image-based cytometer over a range of parasitemia and the results were compared against flow cytometry. Mixed stage parasites were diluted to obtain parasitemia ranging between 0.2% and 9%. Giemsa-stained smears were taken for assessing the parasitemia (in a blind manner to not reveal sample identity) using image-based cytometer in parallel to flow cytometry-based assay. As shown in Fig 7, parasitemia scored using our image-based cytometer were totally comparable to the values obtained from flow cytometric experiments as confirmed through two-tailed paired T test. Corresponding flow cytometer plots are shown in S3 Fig of the supplementary information. The results obtained through this parasitemia determination experiment is evaluated using Bland-Altman method of comparison and it is shown in Fig 6B. It should be noted that the values are within the 95% limits of agreement (lower 95% limit of agreement = 0.6539 and upper 95% limit of agreement = 1.493).

[Figure omitted. See PDF.]

Fig 7. Comparison of parasitemia values using image-based cytometer and flow cytometer.

Mixed stage malaria parasites were diluted in hRBCs to obtain parasitemia ranging between 0.2% and 9% and scored using image-based cytometer. The results were comparable with flow cytometer values with no significant difference. The figure represents the average of three independent measurements performed and the error bars represent the standard deviation of the average values.

https://doi.org/10.1371/journal.pone.0179161.g007

Next, dose-response curves for chloroquine, artemisinin and cycloheximide estimated through image-based cytometer and standard flow cytometer are shown in Fig 8. Despite minor differences, IC₅₀ values defined as half maximal inhibitory concentration calculated through these methods were similar (and mostly in agreement with reported values elsewhere [57–59]) as listed in Table 3. These results demonstrate the feasibility of our proposed system for inexpensive, highly sensitive and reliable, high throughput screening of antimalarials. The IC₅₀ values obtained by the image-based cytometer was evaluated using Bland-Altman method of comparison and it is shown in Fig 6C. It should be noted that the values are within the 95% limits of agreement (lower 95% limit of agreement = 0.9970 and upper 95% limit of agreement = 1.170).

[Figure omitted. See PDF.]

Table 3. Comparison of IC₅₀ determined by two different methods.

https://doi.org/10.1371/journal.pone.0179161.t003

[Figure omitted. See PDF.]

Fig 8. Comparison of dose response curves for the three antimalarial drugs generated using image-based cytometer and flow cytometer.

Well-known anti-malarials; (A) Chloroquine, (B) Artemisinin and (C) Cycloheximide were screened against the laboratory strain, 3D7. Standard growth inhibition assay was performed and IC₅₀ values determined using our image-based cytometer and flow cytometry in parallel [Table 3], to obtain comparable results. The error bars represent standard errors.

https://doi.org/10.1371/journal.pone.0179161.g008

Conclusions

In this study, we have developed an image-based cytometer for scoring and classifying the stages of malaria parasites. This portable system, which is a combination of commercial parts and in-house 3D-printed parts, has a compact dimension of 22 cm (L) x 15 cm (W) x 23 cm (H) and weights less than 2000 grams and can be used in routine laboratory drug screening and potential on field-diagnosis. By comparing the results obtained by our developed method to manual counting and flow cytometry, we have validated its robustness and accuracy to differentially estimate all parasitic stages. For the parasitemia test and IC₅₀ determination assays, results were comparable with flow cytometry method. Several advantages of our portable image-based cytometer are summarized as follows: 1) it provides comparable results to flow cytometry and allows stage-specific scoring of parasitemia; 2) it is a non-fluorescent method with simple preparation and operation, exhibiting high computational speed less than 10 seconds per image (Dimensions: 4096 x 3286) using commonly used laptops; 3) this image-based cytometer is low-cost (less than 600 USD) and compact, therefore making it a reliable and affordable tool to (i) detect and score level of infection and (ii) facilitate screening of antimalarials from natural or synthetic products in resource-poor settings.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Six randomly selected images from the same smear taken and analysed using the developed image-based cytometer under 20X objective lens.

(Quantitative analyses are shown in the Table 1).

https://doi.org/10.1371/journal.pone.0179161.s001

(TIFF)

S2 Fig. Six randomly selected images from the same smear taken and analysed using the developed image-based cytometer under 100X oil immersion objective lens.

(Quantitative analyses are shown in the Table 2).

https://doi.org/10.1371/journal.pone.0179161.s002

(TIFF)

S3 Fig. Comparison of true positives (parasitemia) measured using image-based cytometer and flow cytometer.

The figure represents the histogram plots of parasite values calculated using flow cytometer. Figures (a) to (g) represents the different samples from 1 to 7 showed in Fig 7.

https://doi.org/10.1371/journal.pone.0179161.s003

(TIFF)

Acknowledgments

DY and GS gratefully acknowledge the President’s Graduate Fellowship awarded by the Ministry of Education (MOE), Singapore. YA and RC gratefully acknowledge the support provided by SUTD-MIT International Design Centre (IDG11300101). RC also acknowledges Astar-India Collaboration Grant (RGAST1503) for funding support. YA also acknowledges the support from the Singapore Ministry of Education Academic Research Fund Tier 2 (T2MOE1603).

Author Contributions

1. Conceptualization: DY GS RC YA.

2. Data curation: DY GS JD.

3. Formal analysis: DY GS JD.

4. Funding acquisition: YA RC.

5. Investigation: DY GS.

6. Methodology: DY GS JD.

7. Project administration: YA RC.

8. Resources: YA RC.

9. Software: JD DY SG LB.

10. Supervision: YA RC.

11. Validation: YA RC.

12. Visualization: DY GS YA RC.

13. Writing – original draft: DY GS JD.

14. Writing – review & editing: YA RC.

Citation: Yang D, Subramanian G, Duan J, Gao S, Bai L, Chandramohanadas R, et al. (2017) A portable image-based cytometer for rapid malaria detection and quantification. PLoS ONE 12(6): e0179161. https://doi.org/10.1371/journal.pone.0179161

References

1. Organization WHO. World Malaria Report 2015. 2015.

2. Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin Resistance in Plasmodium falciparum Malaria. New England Journal of Medicine. 2009;361(5):455–67. pmid:19641202

3. Noedl H, Socheat D, Satimai W. Artemisinin-Resistant Malaria in Asia. New England Journal of Medicine. 2009;361(5):540–1. pmid:19641219

4. Peters W. Drug resistance in Plasmodium berghei. Vincke and Lips, 1948. I. Chloroquine resistance. Experimental parasitology. 1965;17(1):80–9. pmid:5843289

5. Dondorp AM, Angus BJ, Chotivanich K, Silamut K, Ruangveerayuth R, Hardeman MR, et al. Red blood cell deformability as a predictor of anemia in severe falciparum malaria. American Journal of Tropical Medicine and Hygiene. 1999;60(5):733–7. pmid:10344643

6. Tangpukdee N, Duangdee C, Wilairatana P, Krudsood S. Malaria Diagnosis: A Brief Review. Korean Journal of Parasitology. 2009;47(2):93–102. pmid:19488414

7. Mwangi TW, Mohammed M, Dayo H, Snow RW, Marsh K. Clinical algorithms for malaria diagnosis lack utility among people of different age groups. Tropical Medicine & International Health. 2005;10(6):530–6. pmid:15941415

8. Reyburn H, Mbatia R, Drakeley C, Carneiro I, Mwakasungula E, Mwerinde O, et al. Overdiagnosis of malaria in patients with severe febrile illness in Tanzania: a prospective study. British Medical Journal. 2004;329(7476):1212–5. pmid:15542534

9. Doderer C, Heschung A, Guntz P, Cazenave J-P, Hansmann Y, Senegas A, et al. A new ELISA kit which uses a combination of Plasmodium falciparum extract and recombinant Plasmodium vivax antigens as an alternative to IFAT for detection of malaria antibodies. Malaria journal. 2007;6(1):19.

10. Makler MT, Ries JM, Williams JA, Bancroft JE, Piper RC, Gibbins BL, et al. Parasite lactate-dehydrogenase as an assay for Plasmodium falciparum drug-sensitivity. American Journal of Tropical Medicine and Hygiene. 1993;48(6):739–41. pmid:8333566

11. John SM, Sudarsanam A, Sitaram U, Moody AH. Evaluation of OptiMAL (TM), a dipstick test for the diagnosis of malaria. Annals of Tropical Medicine and Parasitology. 1998;92(5):621–2. pmid:9797836

12. Lema OE, Carter JY, Nagelkerke N, Wangai MW, Kitenge P, Gikunda SM, et al. Comparison of five methods of malaria detection in the outpatient setting. American Journal of Tropical Medicine and Hygiene. 1999;60(2):177–82. pmid:10072132

13. Cho D, Kim KH, Park SC, Kim YK, Lee KN, Lim CS. Evaluation of rapid immunocapture assays for diagnosis of Plasmodium vivax in Korea. Parasitology Research. 2001;87(6):445–8. pmid:11411942

14. Richardson DC, Ciach M, Zhong KJY, Crandall I, Kain KC. Evaluation of the Makromed dipstick assay versus PCR for diagnosis of Plasmodium falciparum malaria in returned travelers. Journal of Clinical Microbiology. 2002;40(12):4528–30. pmid:12454146

15. Iqbal J, Khalid N, Hira PR. Comparison of two commercial assays with expert microscopy for confirmation of symptomatically diagnosed malaria. Journal of Clinical Microbiology. 2002;40(12):4675–8. pmid:12454171

16. Genton B, D'Acremont V, Rare L, Baea K, Reeder JC, Alpers MP, et al. Plasmodium vivax and mixed infections are associated with severe malaria in children: A prospective cohort study from Papua New Guinea. Plos Medicine. 2008;5(6):881–9. pmid:18563961

17. Palacios G, Quan PL, Jabado OJ, Conlan S, Hirschberg DL, Liu Y, et al. Panmicrobial oligonucleotide array for diagnosis of infectious diseases. Emerging Infectious Diseases. 2007;13(1):73–81. pmid:17370518

18. Morassin B, Fabre R, Berry A, Magnaval JF. One year's experience with the polymerase chain reaction as a routine method for the diagnosis of imported malaria. American Journal of Tropical Medicine and Hygiene. 2002;66(5):503–8. pmid:12201583

19. Poon LLM, Wong BWY, Ma EHT, Chan KH, Chow LMC, Abeyewickreme W, et al. Sensitive and inexpensive molecular test for falciparum malaria: Detecting Plasmodium falciparum DNA directly from heat-treated blood by loop-mediated isothermal amplification. Clinical Chemistry. 2006;52(2):303–6. pmid:16339303

20. Shapiro HA, Mandy F. Cytometry in malaria: Moving beyond Giemsa. Cytometry Part A. 2007;71A(9):643–5. pmid:17712779

21. de Langen AJ, van Dillen J, de Witte P, Mucheto S, Nagelkerke N, Kager P. Automated detection of malaria pigment: feasibility for malaria diagnosing in an area with seasonal malaria in northern Namibia. Tropical Medicine & International Health. 2006;11(6):809–16.

22. Singh B, Bobogare A, Cox-Singh J, Snounou G, Abdullah MS, Rahman HA. A genus- and species-specific nested polymerase chain reaction malaria detection assay for epidemiologic studies. American Journal of Tropical Medicine and Hygiene. 1999;60(4):687–92. pmid:10348249

23. Steenkeste N, Incardona S, Chy S, Duval L, Ekala MT, Lim P, et al. Towards high-throughput molecular detection of Plasmodium: new approaches and molecular markers. Malaria Journal. 2009;8. pmid:19402894

24. Steenkeste N, Rogers WO, Okell L, Jeanne I, Incardona S, Duval L, et al. Sub-microscopic malaria cases and mixed malaria infection in a remote area of high malaria endemicity in Rattanakiri province, Cambodia: implication for malaria elimination. Malaria Journal. 2010;9. pmid:20409349

25. Polley SD, Mori Y, Watson J, Perkins MD, Gonzalez IJ, Notomi T, et al. Mitochondrial DNA Targets Increase Sensitivity of Malaria Detection Using Loop-Mediated Isothermal Amplification. Journal of Clinical Microbiology. 2010;48(8):2866–71. pmid:20554824

26. Demas A, Oberstaller J, DeBarry J, Lucchi NW, Srinivasamoorthy G, Sumari D, et al. Applied Genomics: Data Mining Reveals Species-Specific Malaria Diagnostic Targets More Sensitive than 18S rRNA. Journal of Clinical Microbiology. 2011;49(7):2411–8. pmid:21525225

27. Fontecha GA, Mendoza M, Banegas E, Poorak M, De Oliveira AM, Mancero T, et al. Comparison of molecular tests for the diagnosis of malaria in Honduras. Malaria Journal. 2012;11. pmid:22513192

28. Payne D. Use and limitations of light-microscopy for diagnosing malaria at the primary health-care level. Bulletin of the World Health Organization. 1988;66(5):621–6. pmid:2463112

29. Murray CK, Gasser RA, Magill AJ, Miller RS. Update on rapid diagnostic testing for malaria. Clinical Microbiology Reviews. 2008;21(1):97-+. pmid:18202438

30. Gascoyne P, Satayavivad J, Ruchirawat M. Microfluidic approaches to malaria detection. Acta Tropica. 2004;89(3):357–69. pmid:14744562

31. Erdman LK, Kain KC. Molecular diagnostic and surveillance tools for global malaria control. Travel medicine and infectious disease. 2008;6(1):82–99.

32. Tek FB, Dempster AG, Kale I. Parasite detection and identification for automated thin blood film malaria diagnosis. Computer Vision and Image Understanding. 2010;114(1):21–32.

33. Memeu DM, Kaduki KA, Mjomba A, Muriuki NS, Gitonga L. Detection of plasmodium parasites from images of thin blood smears. Open Journal of Clinical Diagnostics. 2013;3:183–94.

34. Tsai M-H, Yu S-S, Chan Y-K, Jen C-C. Blood Smear Image Based Malaria Parasite and Infected-Erythrocyte Detection and Segmentation. Journal of medical systems. 2015;39(10):1–14.

35. Ross NE, Pritchard CJ, Rubin DM, Duse AG. Automated image processing method for the diagnosis and classification of malaria on thin blood smears. Medical & Biological Engineering & Computing. 2006;44(5):427–36. pmid:16937184

36. Somasekar J, Reddy BE. Segmentation of erythrocytes infected with malaria parasites for the diagnosis using microscopy imaging. Computers & Electrical Engineering. 2015;45:336–51.

37. Lee SA, Leitao R, Zheng G, Yang S, Rodriguez A, Yang C. Color capable sub-pixel resolving optofluidic microscope and its application to blood cell imaging for malaria diagnosis. PloS one. 2011;6(10):e26127. pmid:22022535

38. Moon S, Lee S, Kim H, Freitas-Junior LH, Kang M, Ayong L, et al. An image analysis algorithm for malaria parasite stage classification and viability quantification. PloS one. 2013;8(4):e61812. pmid:23626733

39. Duffy S, Avery VM. Development and Optimization of a Novel 384-Well Anti-Malarial Imaging Assay Validated for High-Throughput Screening. American Journal of Tropical Medicine and Hygiene. 2012;86(1):84–92. pmid:22232455

40. Trager W, Jensen JB. Human malaria parasites in continuous culture. Science. 1976;193(4254):673–5. pmid:781840

41. Lambros C, Vanderberg JP. Synchronization of Plasmodium falciparum erythrocytic stages in culture. Journal of Parasitology. 1979;65(3):418–20. pmid:383936

42. Subramanian G, Rajeev CPB, Mohan CD, Sinha A, Chu TTT, Anusha S, et al. Synthesis and in vitro evaluation of hydrazinyl phthalazines against malaria parasite, Plasmodium falciparum. Bioorganic & Medicinal Chemistry Letters. 2016;26(14):3300–6. pmid:27261180

43. Bredies K, Kunisch K, Pock T. Total generalized variation. SIAM Journal on Imaging Sciences. 2010;3(3):492–526.

44. Duan J, Lu W, Tench C, Gottlob I, Proudlock F, Samani NN, et al. Denoising optical coherence tomography using second order total generalized variation decomposition. Biomedical Signal Processing and Control. 2016;24:120–7.

45. Liu RW, Shi L, Huang W, Xu J, Yu SCH, Wang D. Generalized total variation-based MRI Rician denoising model with spatially adaptive regularization parameters. Magnetic resonance imaging. 2014;32(6):702–20. pmid:24746774

46. Knoll F, Bredies K, Pock T, Stollberger R. Second order total generalized variation (TGV) for MRI. Magnetic resonance in medicine. 2011;65(2):480–91. pmid:21264937

47. Goldstein T, Osher S. The split Bregman method for L1-regularized problems. SIAM journal on imaging sciences. 2009;2(2):323–43.

48. Lu W, Duan J, Qiu Z, Pan Z, Liu RW, Bai L. Implementation of high-order variational models made easy for image processing. Mathematical Methods in the Applied Sciences. 2016.

49. Ding Y, Pardon MC, Agostini A, Faas H, Duan J, Ward WO, et al. Novel Methods for microglia segmentation, feature extraction and classification. IEEE/ACM Transactions on Computational Biology and Bioinformatics. 2016.

50. Duan J, Tench C, Gottlob I, Proudlock F, Bai L. Automated Segmentation of Retinal Layers from Optical Coherent Tomography Images Using Geodesic Distance. arXiv preprint arXiv:160902214. 2016.

51. Scholkopf B, Sung K-K, Burges CJ, Girosi F, Niyogi P, Poggio T, et al. Comparing support vector machines with Gaussian kernels to radial basis function classifiers. IEEE transactions on Signal Processing. 1997;45(11):2758–65.

52. Frean J, Poonsamy B, Shandukani B, Moonasar D, Raman J. Case management of malaria: Diagnosis. Samj South African Medical Journal. 2013;103(10):789–93. pmid:24079635

53. Durrheim DN, Becker PJ, Billinghurst K, Brink A. Diagnostic disagreement—The lessons learnt from malaria diagnosis in Mpumalanga. South African Medical Journal. 1997;87(5):609–11. pmid:9254820

54. Desjardins RE, Canfield C, Haynes J, Chulay J. Quantitative assessment of antimalarial activity in vitro by a semiautomated microdilution technique. Antimicrobial agents and chemotherapy. 1979;16(6):710–8. pmid:394674

55. Iber P, Pavanand K, Wilks N, Colwell E. Evaluation of in vitro drug sensitivity of human Plasmodium falciparum by incorporation of radioactive isoleucine. Journal of the Medical Association of Thailand. 1975;58(11):559–66. pmid:1107459

56. Elabbadi N, Ancelin ML, Vial HJ. Use of radioactive ethanolamine incorporation into phospholipids to assess invitro antimalarial activity by the semiautomated microdilution technique. Antimicrobial Agents and Chemotherapy. 1992;36(1):50–5. pmid:1590699

57. Wilson DW, Langer C, Goodman CD, McFadden GI, Beeson JG. Defining the Timing of Action of Antimalarial Drugs against Plasmodium falciparum. Antimicrobial Agents and Chemotherapy. 2013;57(3):1455–67. pmid:23318799

58. Lee YQ, Goh ASP, Ch'ng JH, Nosten FH, Preiser PR, Pervaiz S, et al. A High-Content Phenotypic Screen Reveals the Disruptive Potency of Quinacrine and 3 ', 4 '-Dichlorobenzamil on the Digestive Vacuole of Plasmodium falciparum. Antimicrobial Agents and Chemotherapy. 2014;58(1):550–8. pmid:24217693

59. Delves M, Plouffe D, Scheurer C, Meister S, Wittlin S, Winzeler EA, et al. The Activities of Current Antimalarial Drugs on the Life Cycle Stages of Plasmodium: A Comparative Study with Human and Rodent Parasites. Plos Medicine. 2012;9(2). pmid:22363211