Next Article in Journal
Automatic Segmentation of Pelvic Cancers Using Deep Learning: State-of-the-Art Approaches and Challenges
Next Article in Special Issue
Challenges Encountered When Evaluating an Antibody-Detecting Point-of-Care Test for Taeniosis in an Endemic Community in Zambia: A Prospective Diagnostic Accuracy Study
Previous Article in Journal
Automated COVID-19 and Heart Failure Detection Using DNA Pattern Technique with Cough Sounds
Previous Article in Special Issue
Assessment of the Sedia HIV Self-Test Device: Usability and Performance in the Hands of Untrained Users in Johannesburg, South Africa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of Quantitative Rapid Isothermal Amplification Assay for Leishmania donovani

1
Institute of Animal Hygiene and Veterinary Public Health, University of Leipzig, An den Tierkliniken 1, 04103 Leipzig, Germany
2
Nutrition and Clinical Services Division, International Centre for Diarrheal Disease Research, Dhaka 1212, Bangladesh
3
Institute of Microbiology & Virology, Brandenburg Medical School, 01968 Senftenberg, Germany
4
Laboratory Sciences and Services Division, International Centre for Diarrheal Disease Research, Dhaka 1212, Bangladesh
*
Authors to whom correspondence should be addressed.
Submission received: 28 September 2021 / Revised: 12 October 2021 / Accepted: 18 October 2021 / Published: 22 October 2021
(This article belongs to the Special Issue Point-of-Care Diagnostics for Low-Resource Settings)

Abstract

:
Quantification of pathogen load, although challenging, is of paramount importance for accurate diagnosis and clinical management of a range of infectious diseases in a point-of-need testing (PONT) scenario such as in resource-limited settings. We formulated a quantification approach to test the standard-curve based absolute quantification ability of isothermal recombinase polymerase amplification (RPA) assay. As a test of principle, a 10-fold dilution series of Leishmania donovani (LD) genomic DNA prepared in nuclease-free-water (NFW), and from culture-spiked-blood (CSB) were tested, and a 15 min assay was performed. A modified algorithm was formulated to derive the detection outcome. The threshold-record times (Tr) in seconds thus obtained were plotted against the initial load of parasite genomes for log-linear regression analysis. The quantitative RPA (Q-RPA) assay was further evaluated against a LD quantitative (q)-PCR assay with DNA extracted from visceral and post-Kala-azar dermal leishmaniasis case specimens and stratified into different ranges of threshold cycle (Ct). The best-fitted regression models were found linear with mean r2/root mean square error (RMSE) values of residual points (in seconds) estimated as 0.996/8.063 and 0.992/7.46 for replicated series of NFW and CSB, respectively. In both series, the lower limit of detection reached less than 0.1 parasite genome equivalent DNA. Absolute agreement between Q-RPA and LD-qPCR was found for test positivity, and strong positive correlations were observed between the Tr and Ct values (r = 0.89; p < 0.0001) as well as between the absolute parasite loads (r = 0.87; p < 0.0001) quantified by respective assays. The findings in this very first Q-RPA assay for leishmaniasis are suggestive of its potential in monitoring LD load in clinical specimens, and the development of rapid Q-RPA assays for other infectious diseases.

1. Introduction

Modern nucleic acid testing (or NAT) approaches such as polymerase chain reaction (PCR) have led to a considerable increase in diagnostic sensitivity, which often outperforms culture-based approaches [1]. However, the fundamental yet unresolved issue of whether mere detection of pathogens can be interpreted for the clinical outcome, either in a well-established diagnostic microbiological laboratory or point-of-need testing (PONT) scenario, remains a challenge. Furthermore, monitoring of pathogen load post-infection can be relevant for clinical management [2], deployment of treatment initiatives [3], drug efficacy and pathogen resistance [4,5], pathogen proliferation kinetics for early detection of disease [6], and therefore relevant for clinical outcome [7,8,9]. Thus, while qualitative diagnostics are the mainstay of diagnosis of most infectious diseases, the development and implementation of quantitative molecular diagnostic assays is of prime importance for a growing number of treatable infectious diseases.
Real-time quantitative (q)-PCR, which employs a time-to-signal based relationship of a specific threshold concentration of amplified products with the initial copy number of the template (i.e., number of cycles or Ct), is the most widely studied approach in the quantification of the pathogen load in biological specimen. Estimation of the amount of initial template concentration by real-time qPCR can be performed either by relative or absolute quantification [10]. For quantitative nucleic acid analysis methods to be used at the point of need, such as in resource-limited settings, among desired characteristics are simplicity in operation, low-cost and feasibility. However, due to specific requirements like cyclic heating and cooling systems and thermally stringent materials, PCR reactions within an enclosed microfluidic system- be it in a conventional or real-time apparatus- render an imperfect solution for incorporation into PONT platforms. Although, efforts have gone into necessary optimization and miniaturization of the assay as well as detection systems onboard [11,12,13], it would be advantageous to focus on detection methods that are less expensive and inherently simple, miniaturizable, robust, and portable for PONT molecular diagnosis. Thus, research has turned towards isothermal nucleic acid amplification techniques [14] that mainly differ from a typical PCR reaction in the strand separation mechanism. Importantly, determination of pathogen load is a desirable characteristic of isothermal technologies for PONT integration [15] and this in response, led to progresses in quantitative strategies for isothermal assay formats such as nucleic acid sequence based amplification (NASBA), loop-mediated amplification (LAMP) and recombinase polymerase amplification (RPA) [16,17,18,19,20].
Leishmania is an intracellular parasite that causes a range of pathologies involving cutaneous and visceral organs of human and animals. Among several clinically important species, L. donovani (LD) typically causes visceral leishmaniasis (VL/Kala-azar) which is often fatal if left untreated, and post-Kala-azar dermal leishmaniasis (PKDL)- a skin complication that may lead to deformity. Although these neglected tropical disease forms are mainly endemic in the Indian subcontinent, however increasing evidences suggest the presence of co-endemic and imported Leishmania species [21,22,23]. Over the last decade a variety of molecular tools including conventional and real-time qPCR, and isothermal based methods such as LAMP, RPA, and NASBA assays have been developed in the search of appropriate tools for the detection of different species of Leishmania [20,24,25], but quantification has been confined to real-time qPCR. The quantification of pathogen load in leishmaniasis has been suggested to be helpful to monitor drug response, parasite kinetics and infection dynamics. For instance, a prospective study to determine the efficacy of two formulations of amphotericin B suggested significant correlation between parasite load and clinical outcome in VL patients, while monitoring of parasite kinetics signified definitive cure with a parasite load as low as below 2 parasites/mL of blood [26]. On the other hand, a threshold load of 1 parasite/mL of blood was observed to distinguish asymptomatic carriage from disease development [27], whereas, a cut-off of 10 parasites/mL blood or lower can indicate an increased chance of disease relapse post-treatment [26,27,28,29]. In an earlier study, we showed that a real-time qPCR assay for LD (LD-qPCR) could successfully differentiate cured VL and PKDL from active and/or relapsed disease forms- rendering the potential of quantification of parasite load as a prognostic marker of VL and PKDL [30]. Nevertheless, RPA offers several advantages over other molecular tools in terms of feasibility, time to result, simplicity in primer design, capability to target longer sequences, tolerance to inhibitors, cost and dispensability of a heating source. It facilitated development of detection assays for a number of pathogens of human and animal importance [31]. For the detection of LD in clinical specimens, with the broad scope of its deployment in the endemic settings where diagnostic establishment of PCR based applications can be challenging, we have previously developed a qualitative LD-RPA assay [24], which has recently completed its phase II diagnostic trial [32]. Monitoring of LD load in patients during diagnosis and in treatment follow-ups can be clinically significant (e.g., drug dosage and efficacy testing, relapse and cure assessment etc.) and may facilitate assessing the prevalence of the potentially transmitting population in resource-poor endemic settings, however implementation of qPCR for this purpose remains farfetched.
Quantitative LD-RPA as a plausible alternative for PONT diagnosis and prognosis in field conditions could be a very powerful tool. Here we describe an absolute RPA quantification assay for LD. Overall, our observations are in support of the potential of the quantitative RPA assay to detect and monitor LD infection in leishmaniasis patients.

2. Materials and Methods

2.1. Theoretical Considerations

The RPA process initiates in a phase separated system formulated by a viscous crowding agent which prompts the reactions to occur within colloidal globules that localize RPA enzymes, nucleotides and template molecules [33]. We assume that this compartmentalized reaction principle of RPA can be harnessed to explain template-time relationship from unperturbed amplification data unlike the default (classical) RPA method that incorporates a mixing step [34]. We further considered that the template-time relation is more likely to be synchronized regardless of template load at the early growth phase of RPA reaction where the limited supply of DNA determines the initial rate of amplification and initiates a roughly exponential trend. Due possibly to the diffusion mediated effects [35,36] on reaction space and time, the template-time relationship becomes progressively nonsynchronous in the continuity of growth phase (Figure S1). Several methodological variations have been formulated to date in order to increase RPA amplification efficiency. The comparative features of these approaches as well as the theoretical considerations behind our methodological approach are highlighted in the Supplementary Materials (Supplementary S1 and Figure S2).

2.2. Sample Source

The samples that were tested in this study comprise of: (1) LD-qPCR positive archived DNA extracted from specimens collected from VL (blood; n = 30) and PKDL (skin biopsy; n = 17) cases; (2) LD-qPCR negative archived DNA samples extracted from specimens collected from asymptomatic (blood; n = 30) subjects, who were positive in rk39 rapid immunochromatographic test (ICT; InBios International Inc., Seattle, WA, USA) or direct agglutination test (DAT) for LD bodies; (3) standard DNA of a WHO reference strain of LD (L. donovani MHOM/IN/80/DD8); and (4) a blood sample spiked with VL patient-derived and in-house cultured promastigotes in a serial dilution series. Archived samples were obtained from the Ethical Review Committee, icddr,b, approved study (PR-18023), in which VL and PKDL patients were enrolled and provided informed consent regarding the use of the anonymous samples in further research.

2.3. DNA Standard Preparation

Serial dilutions of standard DNA of L. donovani MHOM/IN/80/DD8 strain were prepared in nuclease free water (NFW) from 1 ng to 1 fg of parasite DNA corresponding to 104 to 10−2 parasites per reaction of RPA (100 fg equals 1 parasite). For preparation of promastigote culture spiked in blood (CSB), approximately 106 LD parasites were taken and centrifuged at 6000 rpm for 5 min. One ml of PBS (0.02 M, pH 7.2) was added to wash the pellet followed by centrifugation (6000 rpm for 5 min). The pellet was dissolved in 200 µL of PBS (i.e., 106 parasites/200 µL). Twenty µL was added to 180 µL of LD-free human blood to make a final volume containing 105 parasites, which was serially diluted (1:10) in blood to obtain 104 to 100 parasites in 200 µL of blood. This corresponded to 2.5 × 102 to 2.5 × 10−2 parasite equivalent genomic DNA per reaction of RPA (described below). DNA from spiked samples were extracted with DNeasy Blood and Tissue kit (Qiagen, Hilden, Germany).

2.4. RPA and Real-Time qPCR Assays

The classical RPA (hereby referred to as or C-RPA) was performed using the TwistAmp exo kit (TwistDx Ltd., Cambridge, UK) in a 50 μL volume (5 μL template DNA and 45 μL mastermix) as per manufacturer’s recommendation, and according to our previous protocol for LD-RPA development [24]. In brief, after preparing the mastermix, all components were added in the lid of the reaction tube leaving the pellet undissolved. The tubes were then vortexed and centrifuged at once and placed into the Axxin T8 fluorescence reader (Axxin, Fairfield, VIC, Australia) at 42 °C/15 min with 10 or 20 s read-time intervals. After 4 min, the tubes were removed for mixing and brief centrifugation, and then placed back into the reader. Each run included one reaction with molecular grade water as a negative control. Fluorescence was detected in the FAM channel (470 nm excitation and 520 nm emission). In addition to the classical protocol, a modified protocol was performed (hereby termed as quantitative RPA or Q-RPA). In brief, the temperature was set to 40 °C for an uninterrupted assay run. Dissolving of the lyophilized enzyme and the nucleotide mix in each tube was ensured before placing it into the Axxin device. The primer and probe sequences used in either assay targeted the kinetoplast minicircle DNA (GenBank accession no. Y11401.1) which has around 10,000 copies in a single LD cell.
LD-qPCR was performed as described [30]. To quantify the parasite load in each sample, every run included a standard curve generated by assaying serially diluted standard LD DNA that correspond to 104–10−1 parasites per reaction.

2.5. Curation of Q-RPA Amplification Data for Detection and Quantification

The observable template amplification phase denotes positive change in the acquisition of total fluorescence number over time for a consecutive period of read-points, denoted hereby as Ap. We assume that the Ap separates the background/noise floor. Noise-unbiased sample fluorescence (Fn) at any read-point was thus derived and it was used to calculate the fold change ratio in fluorescence acquisition rate (Rfold) for any given read-point above the noise floor. For the reaction to assert exponential behavior, Rfold must be valued >1.0 for a period of consecutive read-points, denoted hereafter as E-phase. The threshold record-time (Tr) would thus correspond to a time point that belongs to E-phase. The optimal cut-off duration (in seconds) of E-phase was then determined by using LD-qPCR positive (n = 30; Ct range: 30.8 to 38.9) and LD-qPCR negative (n = 30) DNA to differentiate between a true template amplification from template negative outcome. To generate standard curve from serial dilution series of parasite genome, Tr value of each dilution was estimated from an arbitrary Fn value (Fn(arbitrary), which is an integer) that is reached within the E-phase by all the respective dilutions. Standard curves were generated from the log-linear regression analysis after plotting Tr values against the parasite loads that represent corresponding dilution. A total of 37 DNA samples extracted from VL (n = 20; blood) and PKDL (n = 17; skin biopsy) clinical specimens were blindly evaluated by C-RPA, Q-RPA and LD-qPCR for comparison. In the Q-RPA and LD-qPCR assay, threshold times and absolute template number were extrapolated from the standard curve.
Furthermore, for generalization of the exponential amplification that is assumed to initiate about the Tr under a hypothetically constant reaction condition, the normalized fluorescence reads starting about the Fn(arbitrary) value were transformed as a function of Tr to fit to a deviated model of what has been previously described mathematically [35]. The details of the detection and quantification algorithm are given in the Supplementary Materials (Supplementary S2).

2.6. Statistical Analysis

The sensitivity and specificity of the Q-RPA assay against C-RPA and LD-qPCR assays were calculated by standard formula [37]. Depending on the distribution of data, pairwise comparisons were performed between threshold points of C-RPA (Tt) and Q-RPA (Tr), and between the absolute parasite loads estimated by Q-RPA and LD-qPCR. Further, correlation analysis was performed among Tr (Q-RPA), Ct (LD-qPCR) and Tt (C-RPA) and respective correlation coefficients (r-value) were determined. Variability in the assay dilution series was reported as the coefficient of variation (CV; shown as the ratio of mean to standard deviation (SD) in Tr values for each of the dilution points of the replicated dilution series. A p-value < 0.05 was considered as an indication of statistically significant result. Statistical analysis and regeneration of plots and curves were performed using SPSS v20 and Graphpad Prism v8.0.

3. Results

3.1. Generalization, Linearity and Reproducibility of Q-RPA Amplification

The reader generated fluorescence values were used to estimate the acquisition period (Ap), based on which normalization was performed to synchronize noise level. This was followed by estimation of an exponential amplification phase (E-phase) and threshold record-time (Tr) that lies within the early E-phase (i.e., Tr ∈ E-phase ∈ Ap). In determining linearity of the template-time relationship, several unique Fn(arbitrary) values within the E-phase of individual dilutions were tested to estimate Tr, which were then plotted against the initial parasite genome number. The log-linear model was chosen based on the goodness-of-fit (r2) of the Tr curve with the least root mean square error (RMSE) values of the residual Tr points.
A simplified model of the RPA amplification process was thus generated by transformation of normalized data under the assumption that Tr is the starting point of exponential growth and, under ideal assay conditions, each of the reactions regardless of dilution is consistent in amplicon generation rate until reaching the fluorescence maxima. In this case, the relative fluorescence accumulated at the terminal end of the reaction by the lowest dilution was the maximal fluorescence number. The RPA amplicon generation rate (k) was derived by taking the average of Rfold values that belong to the E-phase of all the dilution series samples, and thus a simplistic plot of the RPA amplification was obtained (Figure 1).
The standard curve generated by a 10-fold serial dilution series of reference strain DNA in NFW had a linear dynamic range over 7 log-steps (104 to 10−2 parasite genomes) with mean r2/RMSE values estimated as 0.996/8.063. On the other hand, for the CSB dilution series, linearity was found for 5 log-steps (2.5 × 102 to 2.5 × 10−2 parasite genome) with mean r2/RMSE values of 0.992/7.46 (Figure 2). There was no significant difference (p-value = 0.16) in the slope between NFW and CSB dilution series and neither among the replicates of individual series. However, significant difference (p-value < 0.0001) was observed in the elevation between the NFW and CSB series. Reproducibility of the linearity was assessed as inter-assay variation of Tr values for the same dilution series in independent runs. Inter-assay mean coefficients of variations (CV) were found to be 1.19% (95% CI: 0.07–2.32%) and 1.84% (95% CI: 0.46–3.21%) for NFW and CSB series, respectively, indicating high reproducibility (Table 1).

3.2. Threshold E-Phase and Limit of Detection (LOD)

The LD-qPCR positive case samples (n = 30) distributed around high to low Ct values and the LD-qPCR negative asymptomatic samples (n = 30) were used to test whether the E-phase (i.e., the longest consecutive period where Rfold > 1.0) can be a classifier for positive RPA detection in the sample pool (Table S1). It was found that LD-qPCR positive samples had a significantly higher duration of the E-phase at any period of reaction than LD-qPCR negative samples (median values: 552.5 vs. 75, p < 0.0001). The estimated threshold of 330 s predicted both sensitivity and specificity to reach 100% (95% CI: 88.65–100%) (Figure 3). If this threshold is applied to the standard dilution series, the lower LOD for the NFW dilution series would be 0.1 parasite equivalent genomic DNA, while it would remain 0.025 parasite equivalent genomic DNA for CSB series (Table 1).

3.3. Q-RPA Quantification and Comparison with C-RPA and LD-qPCR

Reproducibility of the Q-RPA assay and its quantitative ability were assessed blindly. For this, 17 out of 30 LD-qPCR positive samples that were assayed in the E-phase cut-off determination step were taken randomly along with 20 additional LD-qPCR positive samples ranging from high to low Ct values (Table S2). A 10-fold dilution series of standard genomic DNA (L. donovani MHOM/IN/80/DD8) was contemporarily assayed. The dilution series thus prepared for both Q-RPA and LD-qPCR assays with adjustment of reaction volumes resulted in a detection limit of 0.1 parasite equivalent genomic DNA per reaction in either assay (not shown). The resultant sample Tr each corresponded to the unique Fn(arbitrary) value set to generate the standard linear curve of the dilution series.
The Q-RPA assay successfully detected all the 37 test samples, rendering 100% (95% CI: 88.43–100.00%) sensitivity against LD-qPCR. When C-RPA was performed on the given sample pool, it was found as well to have perfect concordance (Cohen’s kappa: 1.00; 95% CI: 0.93 to 1.00) with Q-RPA for the detection of LD. No significant difference (Wilcoxon matched-pairs signed rank test, p = 0.70) between the threshold times of C-RPA (Tt) and Q-RPA (Tr) was found, however, poor association (Spearman correlation, r = 0.36; 95% CI: 0.02–0.62) was observed between them. On the other hand, a significant pairwise association (p < 0.0001) was observed between Tr and Ct with an r-value of 0.89 (Spearman correlation; 95% CI: 0.80–0.95). When absolute parasite load (as genome copy number) per reaction was extrapolated from the standard curves of respective Q-RPA and LD-qPCR assays, there was a significant difference (Wilcoxon matched-pairs signed rank test, p < 0.0001) between the parasite loads reported by the two assays, however, strong positive association (Spearman correlation, r = 0.87; p < 0.0001) was observed between respective parasite loads, which indicates that the absolute quantitative ability of Q-RPA assay based on standard curve is comparable to that of quantitative PCR (Figure 4).

4. Discussion

Monitoring of parasite load in clinical samples can be advantageous for both clinical management and epidemiological assessment of leishmaniasis, a neglected tropical disease. In this regard, however, establishment of a field-deployable quantitative molecular tool in the resource-limited endemic settings remains a challenge. In this study, a quantitative RPA assay was developed and its performance assessed with clinical samples. Quantification was assessed using genomic DNA of LD and archived DNA extracted from specimens derived from VL and PKDL patients. In order to perform the Q-RPA assay, several modifications were introduced in the assay procedure as well as in the analytical method to derive results. The Q-RPA assay was allowed to record fluorescence without interruption for the entire assay duration omitting one additional step of mixing and centrifugation. As another deviation from the C-RPA approach [24], the exponential trend in amplification was used to define the output (positive or negative) rather than arbitrarily assigning a threshold in the derivative plot. Furthermore, continuous fluorescence recording allowed to derive a standard Q-RPA curve from the noise-normalized data to extrapolate initial template load of test samples. Quantification features (threshold time and parasite load/reaction) were evaluated on archived clinical samples against our previously developed LD-qPCR assay [30]. Strong correlations in threshold time and parasite load between Q-RPA and LD-qPCR assays as observed in this study indicate the potential of the Q-RPA assay as an alternative approach to quantify LD in clinical specimens.
A likely challenge with regards to our approach is defining background/noise and exponential phase in the reaction- two key elements in the polymerase based quantification principle. In the recombinase filament-induced delay time period preceding the growth phase, the signal is actually small but non-zero. This noise floor is generated by random noise processes possibly arising from instrumentation-specific stochastic factors, inefficient fluorophore quenching, primer-dependent artefacts, and due to sensitivity of the signal to matrix effects [35,38]. Our earlier observation suggests that noisy fluorescence is prone to variability even between replicates of the same assays (results not shown). Therefore, before estimating the exponential reporting window, sample-to-sample (or reaction-to-reaction) variation in noise level was essentially minimized by normalization of individual noise to true background noise derived from the non-template control. Next, the Q-RPA reporting window should have a measurable trend of amplification (i.e., fluorescence level) for template-contained samples that is quantitatively different from that of non-template samples. In this study, this trend was defined for the given reporting interval period (i.e., 10/20 s) as the rate of amplicon generation (i.e., fluorescence acquisition) at a given read-point surpassing the rate obtained at its former read-point in a continuous manner for a cut-off reporting window (i.e., exponential phase or E-phase). Although one of the challenges of the mixing free reaction has been attributed to the delay in producing a result [39], it was shown in the present study that by defining the exponential trend within the context of a reporting window of the E-phase, template load of around 100 copies of LD kDNA in the NFW dilution series could be detected within 11 min.
Since the goal of this study was to evaluate RPA quantification in clinical samples of leishmaniasis, the cut-off reporting window of the exponential phase (E-phase) was determined empirically using LD-qPCR positive (high to low Ct values) and LD-qPCR negative clinical samples. The estimated E-phase cut-off suggests that exponential amplification of at least 330 s can differentiate between a positive and negative outcome with good confidence. However, low number of samples used in the estimation of the E-phase cut-off has possibly downgraded the detection limit of NFW series to 0.1 parasite equivalent genomic DNA per reaction in this study, as the standard curves were considerably linear up to 0.01 parasite equivalent genomic DNA. This limit is nevertheless in congruence with our previously developed and clinically evaluated qPCR assay for LD [30]. Furthermore, the LOD of the CSB dilution series was investigated to evaluate the effect of non-template background DNA and sample matrix on the linearity of the RPA amplification trend, particularly for low template reactions. In the CSB series it was observed that RPA amplification trends follow the template-time relationship in a similar manner as in the NFW series (no significant difference in standard curve slopes). However, the LOD of the CSB series was a bit higher (0.025 parasite equivalent genomic DNA per reaction) than that of the NFW series. Nonetheless, the CSB series was restrained by the use of promastigotes spiked in blood, and therefore further dilution up to a fraction of a promastigote was not possible. Interestingly, the CSB series with a similar dilution range tended to reach the Tr much faster resulting in significantly higher elevation of the standard curve than that of the NFW series. This hints at possible strain-specific differences in kDNA copy numbers (i.e., reference LD strain in NFW series vs. in-house clinical strain in CSB series). At the same time, the RPA template-time relationship does not seem to be perturbed in the presence of sample matrix/background DNA. Moreover, there were no significant differences among the respective slopes for the replicate runs of the NFW and CSB dilution series.
A number of real-time qPCR assays that target different genetic markers of Leishmania species and sample types have been developed to date for use in veterinary or human medicine, as discussed elsewhere [40]. Owing to the conserved sequence of the 18S rDNA, and the heterogeneity as well as the high copy number of minicircle kDNA, these target sequences allow detection at the genus or subgenus level with resolution as low as fractions of parasite equivalent genomic DNA [41]. This is particularly beneficial for quantification of parasite load in cases where enumeration of parasitemia is necessary to differentiate clinical manifestation. Minicircle kDNA can exist at around 10,000 copies per LD cell and this may vary among Leishmania species or between isolates of a single species [29,42,43]. But importantly, the copy number of minicircle kDNA is relatively stable in amastigotes infecting a single patient [27,44]. Furthermore, no significant stage specific differences in copy number between amastigotes and promastigotes is evident, which justifies the use of standard curves generated from promastigote-derived DNA for the quantification of amastigotes of that species or strain [43]. In this study, the same standard genomic DNA of LD was used in both LD-qPCR and Q-RPA assays for comparative evaluation of quantitative features. Since all the test samples could not be assayed against a single dilution series due to limitation of slots in both PCR and RPA readers, Ct and Tr as well as template load per reaction was extrapolated from the respective standard series curves. This produced small differences in correlation coefficient values between threshold points (r = 0.89) and respective template loads (r = 0.87). DNA extracted from clinical specimens that involved two different sample matrices (blood for VL, skin biopsy for PKDL) yielded no differences in sensitivity. The absolute sensitivity of the Q-RPA method was determined against the real time LD-qPCR method, whereas absolute agreement of the method with C-RPA was observed for the outcome of the clinical test samples. Notably, we observed a significantly higher parasite load quantified by the Q-RPA assay than that by the LD-qPCR assay (Figure 4C). This is plausible since we considered the equivalent weighted measure of 1 parasite as 100 fg standard genomic DNA as per our established LD-qPCR assay that targets a repetitive region on the Leishmania chromosome [30], whereas there could also be difference in relative copy number of the extrachromosomal kDNA sequences between standard and clinical samples. Nevertheless, the correlation matrices between the two assays suggest that as an alternative to the LD-qPCR, the quantitative RPA assay has the potential to be used for monitoring of parasitemia levels in individual patients with different clinical manifestations. In the next phase of the study, evaluation with expanded number of samples from individuals with different clinical manifestations and transition (such as relapsed or cured) will be performed to validate the quantitative Q-RPA assay.

5. Conclusions

Quantitative metrics may be used to estimate infectious disease burden and/or offer a longitudinal assessment of response to treatment regimens used in disease mitigation efforts. Our target in this study was to evaluate the quantification capability of RPA to convey clinical significance in PONT field conditions and in particular resource-poor settings avoiding the use of any supplementary machinery such as continuous stirrer, vortex, and centrifuge. In this study, all the calculations in derivation of outcome/quantification from raw fluorescence data were performed using excel sheets and simple statistical formulas only (Figure S3). By defining the amplification trend and heuristically evaluating it with relative genomic and clinical samples, we showed that the Q-RPA assay has quantification potential for monitoring parasitemia. Further studies with a larger sample size and potentially with RPA assays for other pathogens of interest can be undertaken for fine tuning and validation of the approach described here. Moreover, the mathematical considerations to attain an outcome at each read-point potentiates the Q-RPA data analysis in a simple and automatically programmable interface for real time monitoring and end point quantification.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/diagnostics11111963/s1, Supplementary S1: Theoretical considerations, Supplementary S2: Curation of Q-RPA amplification data for detection and quantification, Figure S1: Schematic of onset of an RPA reaction step in a globule, Figure S2: Schematic representation of amplification trends in various RPA approaches, where reaction globules are either template-positive (yellow) or template-negative (green), Figure S3: Q-RPA variable generation from machine read data in MS-excel sheet, Table S1: Information of test samples used in E-phase threshold determination, Table S2: Test sample information and results obtained in LD-qPCR, Q-RPA and C-RPA assays.

Author Contributions

Conceptualization and methodology, M.A.A.K.; investigation, M.A.A.K., K.F. and R.C.; formal analysis, M.A.A.K., K.F., R.C. and P.G.; data curation, M.A.A.K., K.F., R.C., P.G. and F.H.; writing—original draft preparation, M.A.A.K.; writing—review and editing, M.W., A.A.E.W. and D.M.; visualization and validation, M.A.A.K., A.A.E.W. and M.W.; supervision, A.A.E.W. and D.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

Md Anik Ashfaq Khan is a Ph.D. student of the Faculty of Medicine at the University of Leipzig and his position is funded by the German Academic Exchange Service (DAAD). We thank the participants of the original study and project staffs of the International Centre for Diarrheal Disease Research, Bangladesh (icddr,b). We also thank icddr,b’s core donors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carman, B. Molecular techniques should now replace cell culture in diagnostic virology laboratories. Rev. Med. Virol. 2001, 11, 347–349. [Google Scholar] [CrossRef]
  2. Limaye, A.P.; Huang, M.-L.; Leisenring, W.; Stensland, L.; Corey, L.; Boeckh, M. Cytomegalovirus (CMV) DNA Load in Plasma for the Diagnosis of CMV Disease before Engraftment in Hematopoietic Stem-Cell Transplant Recipients. J. Infect. Dis. 2001, 183, 377–382. [Google Scholar] [CrossRef]
  3. Cesaro, S.; Murrone, A.; Mengoli, C.; Pillon, M.; Biasolo, M.A.; Calore, E.; Tridello, G.; Varotto, S.; Alaggio, R.; Zanesco, L.; et al. The real-time polymerase chain reaction-guided modulation of immunosuppression enables the pre-emptive management of Epstein–Barr virus reactivation after allogeneic haematopoietic stem cell transplantation. Br. J. Haematol. 2005, 128, 224–233. [Google Scholar] [CrossRef]
  4. Beshir, K.B.; Hallett, R.L.; Eziefula, A.C.; Bailey, R.; Watson, J.; Wright, S.G.; Chiodini, P.L.; Polley, S.D.; Sutherland, C.J. Measuring the efficacy of anti-malarial drugs in vivo: Quantitative PCR measurement of parasite clearance. Malar. J. 2010, 9, 312. [Google Scholar] [CrossRef] [Green Version]
  5. Pas, S.D.; Noppornpanth, S.; Eijk, A.A.v.d.; Man, R.A.d.; Niesters, H.G.M. Quantification of the newly detected lamivudine resistant YSDD variants of Hepatitis B virus using molecular beacons. J. Clin. Virol. 2005, 32, 166–172. [Google Scholar] [CrossRef] [PubMed]
  6. Lion, T.; Baumgartinger, R.; Watzinger, F.; Matthes-Martin, S.; Suda, M.; Preuner, S.; Futterknecht, B.; Lawitschka, A.; Peters, C.; Pötschger, U.; et al. Molecular monitoring of adenovirus in peripheral blood after allogeneic bone marrow transplantation permits early diagnosis of disseminated disease. Blood 2003, 102, 1114–1120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Biedermann, K.; Dandachi, N.; Trattner, M.; Vogl, G.; Doppelmayr, H.; Moré, E.; Staudach, A.; Dietze, O.; Hauser-Kronberger, C. Comparison of Real-Time PCR Signal-Amplified In Situ Hybridization and Conventional PCR for Detection and Quantification of Human Papillomavirus in Archival Cervical Cancer Tissue. J. Clin. Microbiol. 2004, 42, 3758–3765. [Google Scholar] [CrossRef] [Green Version]
  8. Snijders, P.J.; van den Brule, A.J.; Meijer, C.J. The clinical relevance of human papillomavirus testing: Relationship between analytical and clinical sensitivity. J. Pathol. 2003, 201, 1–6. [Google Scholar] [CrossRef]
  9. Watzinger, F.; Suda, M.; Preuner, S.; Baumgartinger, R.; Ebner, K.; Baskova, L.; Niesters, H.G.M.; Lawitschka, A.; Lion, T. Real-Time Quantitative PCR Assays for Detection and Monitoring of Pathogenic Human Viruses in Immunosuppressed Pediatric Patients. J. Clin. Microbiol. 2004, 42, 5189–5198. [Google Scholar] [CrossRef] [Green Version]
  10. Freeman, W.M.; Walker, S.J.; Vrana, K.E. Quantitative RT-PCR: Pitfalls and Potential. BioTechniques 1999, 26, 112–125. [Google Scholar] [CrossRef]
  11. Neuzil, P.; Zhang, C.; Pipper, J.; Oh, S.; Zhuo, L. Ultra fast miniaturized real-time PCR: 40 cycles in less than six minutes. Nucleic Acids Res. 2006, 34, e77. [Google Scholar] [CrossRef]
  12. Quan, P.-L.; Sauzade, M.; Brouzes, E. dPCR: A Technology Review. Sensors 2018, 18, 1271. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Xiang, Q.; Xu, B.; Li, D. Miniature real time PCR on chip with multi-channel fiber optical fluorescence detection module. Biomed. Microdevices 2007, 9, 443–449. [Google Scholar] [CrossRef]
  14. Zhao, Y.; Chen, F.; Li, Q.; Wang, L.; Fan, C. Isothermal Amplification of Nucleic Acids. Chem. Rev. 2015, 115, 12491–12545. [Google Scholar] [CrossRef]
  15. Craw, P.; Balachandran, W. Isothermal nucleic acid amplification technologies for point-of-care diagnostics: A critical review. Lab Chip 2012, 12, 2469–2486. [Google Scholar] [CrossRef]
  16. Crannell, Z.A.; Rohrman, B.; Richards-Kortum, R. Quantification of HIV-1 DNA Using Real-Time Recombinase Polymerase Amplification. Anal. Chem. 2014, 86, 5615–5619. [Google Scholar] [CrossRef] [Green Version]
  17. Hsieh, K.; Patterson, A.S.; Ferguson, B.S.; Plaxco, K.W.; Soh, H.T. Rapid, Sensitive, and Quantitative Detection of Pathogenic DNA at the Point of Care through Microfluidic Electrochemical Quantitative Loop-Mediated Isothermal Amplification. Angew. Chem. Int. Ed. 2012, 51, 4896–4900. [Google Scholar] [CrossRef] [Green Version]
  18. Nixon, G.; Garson, J.A.; Grant, P.; Nastouli, E.; Foy, C.A.; Huggett, J.F. Comparative Study of Sensitivity, Linearity, and Resistance to Inhibition of Digital and Nondigital Polymerase Chain Reaction and Loop Mediated Isothermal Amplification Assays for Quantification of Human Cytomegalovirus. Anal. Chem. 2014, 86, 4387–4394. [Google Scholar] [CrossRef]
  19. Shen, F.; Davydova, E.K.; Du, W.; Kreutz, J.E.; Piepenburg, O.; Ismagilov, R.F. Digital Isothermal Quantification of Nucleic Acids via Simultaneous Chemical Initiation of Recombinase Polymerase Amplification Reactions on SlipChip. Anal. Chem. 2011, 83, 3533–3540. [Google Scholar] [CrossRef] [Green Version]
  20. van der Meide, W.F.; Schoone, G.J.; Faber, W.R.; Zeegelaar, J.E.; de Vries, H.J.C.; Özbel, Y.; Lai A Fat, R.F.M.; Coelho, L.I.A.R.C.; Kassi, M.; Schallig, H.D.F.H. Quantitative Nucleic Acid Sequence-Based Assay as a New Molecular Tool for Detection and Quantification of Leishmania Parasites in Skin Biopsy Samples. J. Clin. Microbiol. 2005, 43, 5560–5566. [Google Scholar] [CrossRef] [Green Version]
  21. Khan, M.A.A.; Chowdhury, R.; Nath, R.; Hansen, S.; Nath, P.; Maruf, S.; Abd El Wahed, A.; Mondal, D. Imported cutaneous leishmaniasis: Molecular investigation unveils Leishmania major in Bangladesh. Parasites Vectors 2019, 12, 527. [Google Scholar] [CrossRef] [Green Version]
  22. Maurício, I.L.; Stothard, J.R.; Miles, M.A. The Strange Case of Leishmania chagasi. Parasitol. Today 2000, 16, 188–189. [Google Scholar] [CrossRef]
  23. Siriwardana, H.V.Y.D.; Karunanayake, P.; Goonerathne, L.; Karunaweera, N.D. Emergence of visceral leishmaniasis in Sri Lanka: A newly established health threat. Pathog. Glob. Health 2017, 111, 317–326. [Google Scholar] [CrossRef]
  24. Mondal, D.; Ghosh, P.; Khan, M.A.A.; Hossain, F.; Böhlken-Fascher, S.; Matlashewski, G.; Kroeger, A.; Olliaro, P.; Abd El Wahed, A. Mobile suitcase laboratory for rapid detection of Leishmania donovani using recombinase polymerase amplification assay. Parasites Vectors 2016, 9, 281. [Google Scholar] [CrossRef] [Green Version]
  25. Nzelu, C.O.; Kato, H.; Peters, N.C. Loop-mediated isothermal amplification (LAMP): An advanced molecular point-of-care technique for the detection of Leishmania infection. PLoS Negl. Trop. Dis. 2019, 13, e0007698. [Google Scholar] [CrossRef] [Green Version]
  26. Sudarshan, M.; Weirather, J.L.; Wilson, M.E.; Sundar, S. Study of parasite kinetics with antileishmanial drugs using real-time quantitative PCR in Indian visceral leishmaniasis. J. Antimicrob. Chemother. 2011, 66, 1751–1755. [Google Scholar] [CrossRef]
  27. Mary, C.; Faraut, F.; Drogoul, M.P.; Xeridat, B.; Schleinitz, N.; Cuisenier, B.; Dumon, H. Reference values for Leishmania infantum parasitemia in different clinical presentations: Quantitative polymerase chain reaction for therapeutic monitoring and patient follow-up. Am. J. Trop. Med. Hyg. 2006, 75, 858–863. [Google Scholar] [CrossRef] [Green Version]
  28. Bossolasco, S.; Gaiera, G.; Olchini, D.; Gulletta, M.; Martello, L.; Bestetti, A.; Bossi, L.; Germagnoli, L.; Lazzarin, A.; Uberti-Foppa, C.; et al. Real-Time PCR Assay for Clinical Management of Human Immunodeficiency Virus-Infected Patients with Visceral Leishmaniasis. J. Clin. Microbiol. 2003, 41, 5080–5084. [Google Scholar] [CrossRef] [Green Version]
  29. Mary, C.; Faraut, F.; Lascombe, L.; Dumon, H. Quantification of Leishmania infantum DNA by a Real-Time PCR Assay with High Sensitivity. J. Clin. Microbiol. 2004, 42, 5249–5255. [Google Scholar] [CrossRef] [Green Version]
  30. Hossain, F.; Ghosh, P.; Khan, M.A.A.; Duthie, M.S.; Vallur, A.C.; Picone, A.; Howard, R.F.; Reed, S.G.; Mondal, D. Real-time PCR in detection and quantitation of Leishmania donovani for the diagnosis of Visceral Leishmaniasis patients and the monitoring of their response to treatment. PLoS ONE 2017, 12, e0185606. [Google Scholar] [CrossRef]
  31. Li, J.; Macdonald, J.; von Stetten, F. Review: A comprehensive summary of a decade development of the recombinase polymerase amplification. Analyst 2019, 144, 31–67. [Google Scholar] [CrossRef] [Green Version]
  32. Ghosh, P.; Sharma, A.; Bhattarai, N.R.; Abhishek, K.; Nisansala, T.; Kumar, A.; Böhlken-Fascher, S.; Chowdhury, R.; Khan, M.A.A.; Faisal, K.; et al. A Multi-Country, Single-Blinded, Phase 2 Study to Evaluate a Point-of-Need System for Rapid Detection of Leishmaniasis and Its Implementation in Endemic Settings. Microorganisms 2021, 9, 588. [Google Scholar] [CrossRef]
  33. Powell, M.L.; Bowler, F.R.; Martinez, A.J.; Greenwood, C.J.; Armes, N.; Piepenburg, O. New Fpg probe chemistry for direct detection of recombinase polymerase amplification on lateral flow strips. Anal. Biochem. 2018, 543, 108–115. [Google Scholar] [CrossRef]
  34. TwistDx. TwistAmp® Combined Instruction Manual; TwistDx: Cambridge, UK, 2014. [Google Scholar]
  35. Gordon, M.I.; Klemer, D.P.; Fuller, S.L.; Chang, J.H.; Klemer, D.R.; Putnam, M.L. Mathematical modeling of a real-time isothermal amplification assay for Erwinia amylovora. Eng. Rep. 2019, 1, e12047. [Google Scholar] [CrossRef] [Green Version]
  36. Moody, C.; Newell, H.; Viljoen, H. A mathematical model of recombinase polymerase amplification under continuously stirred conditions. Biochem. Eng. J. 2016, 112, 193–201. [Google Scholar] [CrossRef]
  37. Parikh, R.; Mathai, A.; Parikh, S.; Chandra Sekhar, G.; Thomas, R. Understanding and using sensitivity, specificity and predictive values. Indian J. Ophthalmol. 2008, 56, 45–50. [Google Scholar] [CrossRef]
  38. Piepenburg, O.; Williams, C.H.; Stemple, D.L.; Armes, N.A. DNA Detection Using Recombination Proteins. PLoS Biol. 2006, 4, e204. [Google Scholar] [CrossRef]
  39. Kalsi, S.; Valiadi, M.; Tsaloglou, M.-N.; Parry-Jones, L.; Jacobs, A.; Watson, R.; Turner, C.; Amos, R.; Hadwen, B.; Buse, J.; et al. Rapid and sensitive detection of antibiotic resistance on a programmable digital microfluidic platform. Lab Chip 2015, 15, 3065–3075. [Google Scholar] [CrossRef] [PubMed]
  40. Galluzzi, L.; Ceccarelli, M.; Diotallevi, A.; Menotta, M.; Magnani, M. Real-time PCR applications for diagnosis of leishmaniasis. Parasites Vectors 2018, 11, 273. [Google Scholar] [CrossRef]
  41. Kocher, A.; ValiÈRe, S.; BaÑUls, A.-L.; Murienne, J. High-throughput sequencing of kDNA amplicons for the analysis of Leishmania minicircles and identification of Neotropical species. Parasitology 2017, 145, 585–594. [Google Scholar] [CrossRef]
  42. Alonso, D.P.; Costa, D.L.; de Mendonça, I.L.; Costa, C.H.N.; Ribolla, P.E.M. Heterogeneity of Leishmania infantum chagasi Kinetoplast DNA in Teresina (Brazil). Am. Soc. Trop. Med. Hyg. 2010, 82, 819–821. [Google Scholar] [CrossRef] [PubMed]
  43. Weirather, J.L.; Jeronimo, S.M.B.; Gautam, S.; Sundar, S.; Kang, M.; Kurtz, M.A.; Haque, R.; Schriefer, A.; Talhari, S.; Carvalho, E.M.; et al. Serial Quantitative PCR Assay for Detection, Species Discrimination, and Quantification of Leishmania spp. in Human Samples. J. Clin. Microbiol. 2011, 49, 3892–3904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Jara, M.; Adaui, V.; Valencia, B.M.; Martinez, D.; Alba, M.; Castrillon, C.; Cruz, M.; Cruz, I.; Van der Auwera, G.; Llanos-Cuentas, A.; et al. Real-Time PCR Assay for Detection and Quantification of Leishmania (viannia) Organisms in Skin and Mucosal Lesions: Exploratory Study of Parasite Load and Clinical Parameters. J. Clin. Microbiol. 2013, 51, 1826–1833. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. RPA amplification curves of NFW dilution series (104 to 10−2 genome/reaction) in terms of relative fluorescence acquisition in a 15 min (900 s) run protocol. Plots were generated by using (a) machine read data; (b) normalized data; (c) transformed data to generalize the amplification.
Figure 1. RPA amplification curves of NFW dilution series (104 to 10−2 genome/reaction) in terms of relative fluorescence acquisition in a 15 min (900 s) run protocol. Plots were generated by using (a) machine read data; (b) normalized data; (c) transformed data to generalize the amplification.
Diagnostics 11 01963 g001
Figure 2. Performance of Q-RPA for obtaining standard curve. Tr values are plotted from replicated assay runs against parasite load equivalent genomic DNA. The plot represents Tr values and parasite load fitting a linear function. (a) DNA extracted from reference strain and diluted in NFW in the range of 104 to 10−2 parasite per reaction and assayed in duplicate (b) DNA extracted from parasite spiked blood to give a range of 2.5 × 102 to 2.5 × 10−2 parasite per reaction and assayed in triplicate.
Figure 2. Performance of Q-RPA for obtaining standard curve. Tr values are plotted from replicated assay runs against parasite load equivalent genomic DNA. The plot represents Tr values and parasite load fitting a linear function. (a) DNA extracted from reference strain and diluted in NFW in the range of 104 to 10−2 parasite per reaction and assayed in duplicate (b) DNA extracted from parasite spiked blood to give a range of 2.5 × 102 to 2.5 × 10−2 parasite per reaction and assayed in triplicate.
Diagnostics 11 01963 g002
Figure 3. E-phase threshold determination. (a) Comparison between LD-qPCR positive and LD-qPCR negative samples in E-phase duration (in seconds) of Q-RPA. (b) Receiver-operating-characteristic curve (AUC = 1.00, p < 0.0001) for E-phase as the predictor of positive results by Q-RPA. Coordinates of the curve (green) are indicated by black dots.
Figure 3. E-phase threshold determination. (a) Comparison between LD-qPCR positive and LD-qPCR negative samples in E-phase duration (in seconds) of Q-RPA. (b) Receiver-operating-characteristic curve (AUC = 1.00, p < 0.0001) for E-phase as the predictor of positive results by Q-RPA. Coordinates of the curve (green) are indicated by black dots.
Diagnostics 11 01963 g003
Figure 4. Quantitative performance of Q-RPA assay. (a) Threshold time difference between C-RPA (Tt) and Q-RPA (Tr); (b) Correlation plots for Tr vs. Tt and Tt vs. Ct; (c) Absolute parasite load per reaction extrapolated from LD-qPCR and Q-RPA standard dilution curves, **** p < 0.0001; (d) Correlation plot for Tr vs. Ct and relative parasite loads between LD-qPCR and Q-RPA assays.
Figure 4. Quantitative performance of Q-RPA assay. (a) Threshold time difference between C-RPA (Tt) and Q-RPA (Tr); (b) Correlation plots for Tr vs. Tt and Tt vs. Ct; (c) Absolute parasite load per reaction extrapolated from LD-qPCR and Q-RPA standard dilution curves, **** p < 0.0001; (d) Correlation plot for Tr vs. Ct and relative parasite loads between LD-qPCR and Q-RPA assays.
Diagnostics 11 01963 g004
Table 1. Reproducibility of the Q-RPA assay.
Table 1. Reproducibility of the Q-RPA assay.
Inter-Assay Variation of Tr in Seconds (NFW)Inter-Assay Variation of Tr in Seconds (CSB)
Parasite loadAssay 1Assay 2MeanSDCV (%)Parasite loadAssay 1Assay 2Assay 3MeanSDCV (%)
1 × 104301.32302.02301.670.500.172.5 × 102262.67255.33273.10263.708.933.39
1 × 103360.45357.02358.732.430.682.5 × 101304.77297.56311.60304.647.022.30
1 × 102413.19393.14403.1614.183.522.5 × 100341.96354.85350.90349.246.601.89
1 × 101452.70462.87457.797.191.572.5 × 10−1395.85403.17395.08398.034.461.12
1 × 100496.73500.12498.422.400.482.5 × 10−2459.90455.29457.57457.582.300.50
1 × 10−1555.79556.92556.360.800.14
1 × 10−2610.96626.92618.9411.291.82
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khan, M.A.A.; Faisal, K.; Chowdhury, R.; Ghosh, P.; Hossain, F.; Weidmann, M.; Mondal, D.; Abd El Wahed, A. Development of Quantitative Rapid Isothermal Amplification Assay for Leishmania donovani. Diagnostics 2021, 11, 1963. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics11111963

AMA Style

Khan MAA, Faisal K, Chowdhury R, Ghosh P, Hossain F, Weidmann M, Mondal D, Abd El Wahed A. Development of Quantitative Rapid Isothermal Amplification Assay for Leishmania donovani. Diagnostics. 2021; 11(11):1963. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics11111963

Chicago/Turabian Style

Khan, Md Anik Ashfaq, Khaledul Faisal, Rajashree Chowdhury, Prakash Ghosh, Faria Hossain, Manfred Weidmann, Dinesh Mondal, and Ahmed Abd El Wahed. 2021. "Development of Quantitative Rapid Isothermal Amplification Assay for Leishmania donovani" Diagnostics 11, no. 11: 1963. https://0-doi-org.brum.beds.ac.uk/10.3390/diagnostics11111963

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop