JCDR - Register at Journal of Clinical and Diagnostic Research
Journal of Clinical and Diagnostic Research, ISSN - 0973 - 709X
Microbiology Section DOI : 10.7860/JCDR/2018/37316.12307
Year : 2018 | Month : Nov | Volume : 12 | Issue : 11 Full Version Page : DC22 - DC29

An Insight for the Future Development of Diagnostic Tool by Exploiting Novel Leishmania Donovani Recombinant Hypothetical Protein

Fauzia Jamal1, Manish Kumar Singh2, Jagadish Hansa3, Pushpanjali4, Pushkar Shivam5, Sarita Kumari6, Shyam Narayan7, Shubhankar K Singh8

1 Senior Research Fellow, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Patna, Bihar, India.
2 Senior Research Fellow, Department of Biotechnology, Tilka Manjhi Bhagalpur University, Bhagalpur, Bihar, India.
3 Scientist B, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Patna, Bihar, India.
4 Senior Research Fellow, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Patna, Bihar, India.
5 Senior Research Fellow, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Patna, Bihar, India.
6 Senior Research Fellow, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Patna, Bihar, India.
7 Scientist D, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Patna, Bihar, India.
8 Scientist D, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Patna, Bihar, India.


NAME, ADDRESS, E-MAIL ID OF THE CORRESPONDING AUTHOR: Dr. Shubhankar K Singh, Scientist D, Department of Microbiology, Rajendra Memorial Research Institute of Medical Sciences, Agamkuan, Patna-800007, Bihar, India.
E-mail: shubhankark.singh@icmr.gov.in
Abstract

Introduction

Assenting diagnosis of Visceral Leishmaniasis (VL) relies on the detection of Leishmania donovani (L. donovani) in splenic and bone marrow specimens obtained by invasive techniques. Thus, the development of inexpensive, non-invasive serological test encompassing high specificity, sensitivity and diagnostic efficacy is urgently needed.

Aim

To assess the significance of recombinant proteins possessing B-cell epitopes in VL diagnosis.

Materials and Methods

Employing immunoinformatics approach, the B-cell epitope footprint of L. donovani hypothetical proteins (derived from earlier studies) were decrypted holding good antigenic character with numerous epitopes. L. donovani hypothetical proteins XP_003860226.1 and XP_003861271.1 were first time cloned as His-tagged fusion proteins and purified as novel recombinant protein antigens, designated rLdhyb and rLdhyc respectively. Sanger sequencing method was exploited to sequence gene insert (GeneBank accession number MH479406).

Results

B-cell epitopes revealed 100% conservancy with L. infantum. Immunoinformatics data revealed no significant sequence similarity with homo sapien and the causative agent of other diseases like tuberculosis, typhoid, malaria etc., resembling in symptoms to VL. Sequencing chromatogram of cloned gene Ldhyb and Ldhyc revealed 98% and 94% identity with L. donovani. ELISA revealed the absolute specificity with sensitivity of 95.4% for rLdhyb and 91% for rLdhyc. Area under curve for rLdhyb, rLdhyc and SLA were 0.99, 0.99 and 0.961, with standard error 0.002, 0.007 and 0.019 respectively. The in silico data was coherently supported by in vitro result.

Conclusion

Absolute specificity, high sensitivity and diagnostic efficacy (rLdhyb: 98%; rLdhyc: 97%) advocated their excellent biomarker property. The present findings provide some basic insights for the future development of novel hypothetical proteins based non-invasive diagnostic tool for VL detection.

Keywords

B-cell epitopes,Diagnostic efficacy,Immunoinformatic,Sanger sequencing,Sensitivity

Introduction

In addition to clinical symptoms, clinical history, geographic location, and travel history of the patient, the detection and diagnosis of Leishmaniasis still rely on very old and labour-intensive microscopic validation. Leishmaniasis is an economically orphan disease occurring amongst socially neglected people throughout the world. It is caused by a digenetic protozoan parasite transmitted to vertebrate hosts by the bite of blood-sucking female phlebotomine sand flies [1]. The most severe form of this disease is VL, which is associated with high fatality in absence of proper diagnosis and treatment. More than 70% of cases of VL cases from India are reported from Bihar and the current control measures rely solely on chemotherapy [2]. A robust surveillance system, as well as simple and inexpensive diagnostic kit, is utmost required to control this dreadful parasitic disease [3-5].

Though the microscopic diagnosis remains the gold standard for parasite examination [6], the sensitivity varies from bone marrow smears (65%) to splenic aspirates (>95%) [7]. Furthermore, microscopic evaluation of parasite burden is highly sensitive and mostly relying on the expertise of the microscopist. The induction of strong humoral immune response in active VL patients made it an attractive target for serodiagnosis [8-10]. Often, a simple, rapid immunochromatographic test (rK39) is most frequently used for the diagnosis of VL however; 15-32% of the healthy individual living in the endemic region gave false positive results. This might be due to the activation of the B-cells presents in healthy individuals living in the endemic area, which result in the generation of polyclonal antibody against L. donovani antigens. Furthermore, rK39 illustrate false positive results in relapse cases and the sensitivity varies across the ethnic population (85%-97%) [11]. As the patients suffering from VL reveal a significant increase in the parasite-specific immunoglobulin, thus making it an attractive target for the development of serological tests for VL diagnosis [12]. Several attempts have been made to develop recombinant chimeric antigen expressing immuno-dominant B-cell epitopes of Leishmania spp. for serodiagnosis of VL but none of them showed promising results in the clinical trials or made it to the market [13,14].

Previously, authors had reported the recognition of key antigens (XP_003860226.1 and XP_003861271.1) from Circulating Immune Complexes (CIC) of Leishmania donovani (L. donovani) patients using immunoproteomics approach [15]. Concomitantly, the present study was mainly focused to clone, sequence, express and purify novel proteins and explore their immuno-reactivity using immunoinformatics, immunoblotting and Enzyme-Linked Immunosorbent Assay (ELISA). Also, an attempt was made to evaluate their ability to detect the specific antibody in the sera of active VL cases as well as in other disease backgrounds. Hence, the current study was focused on the precise diagnosis of VL.

Materials and Methods

The present study was designed to evaluate the diagnostic importance of recombinant hypothetical proteins carrying numerous B-cell epitopes. This study includes initial experimental finding which may serve as a basis for the development of the precise and specific VL diagnosis in future. To achieve this target, proteins carrying B-cell epitopes were validated using immunoinformatics approach, and the, proteins under study were cloned, sequenced, expressed and purified. Finally, diagnostic efficacy of the recombinant proteins was evaluated through ELISA.

Study Population

All clinical investigations were conducted at Rajendra Memorial Research Institute of Medical Sciences (RMRIMS), ICMR, Patna, Bihar, India, in the month of September 2017, on human serum sample after getting written, informed consent as per guideline of the Institutional Ethical Committee, registration number: ECR/480/Inst/BR/2014/RR-17 under rule 122DD of the Drug and cosmetic Rule 1945, RMRIMS, Agamkuan, Patna, in accordance with the Helsinki declaration of 1975 which was revised in 2000. Altogether 99 blood samples from human subjects of both sex and all age groups in between 5 to 45 years were studied. For ELISA, sera samples from 22 (VL-BT) along with 11 samples each from Healthy Endemic (HE), Healthy Non-Endemic (HNE), tuberculosis, viral flu, malaria, asthma and filariasis were procured. The sera were stored at -20°C until further use.

Evaluation of Immunogenicity and Conservancy of the Peptide

The potential B- and T-cell epitopes can be determined using immunology focused resources and software, which saves both the time and cost required for laboratory analysis [16,17]. To measure the immunogenicity of the selected peptide, various immunoinformatics algorithms such as BpiPred, ABCpred and EPMLR were used with default threshold value for linear B-cell epitopes prediction [2,18]. BLASTp (http://blast.ncbi.nlm.nih.gov/Blast.cgi) search was performed against the NCBI protein database to retrieve the homologs protein sequences using L. donovani XP_003860226.1 and XP_003861271.1 as query [16]. Epitope conservancy analysis tool from the Immune Epitope Database (IEDB) analysis resource (http://tools.immuneepitope.org/tools/conservancy/iedb_input) was used.

L. Donovani Culture, Isolation of Genomic DNA and Soluble Leishmania Antigen (SLA)

L. donovani culture, genomic DNA isolation and SLA preparation were done according to the protocol mentioned in Singh MK et al., [19].

Amplification and Cloning of L. Donovani rLdhyb and rLdhyc

The rLdhyb and rLdhyc encoding sequences (NCBI Accession no. XP_003860226.1 and NCBI Accession no. XP_003861271.1 respectively) from L. donovani genomic DNA was amplified by PCR. A set of forward (5’TTT GGA TCCATG CTG CGT TTC TGC3’) and reverse (5’AAT AAG CTT TCA GCA CAC CGT CCG3’) primers for were used for rLdhyb and forward (5’-GTC GGA TCC ATG AAG CCG CTA GTG -3’) and a reverse (5’-GTC CTC GAG TCA GCT GTT GTT CTT-3’) primers for rLdhyc. These primers were synthesized commercially (IDT, India). The rLdhyb and rLdhyc gene (XM_003860178.1 and XM_003861223.1) was amplified using the above primers in Thermal cycler (Bio-Rad, USA). Melting temperature (Tm) of forward and reverse primers for rLdhyb was 60.5°C and 59.7°C, and that for rLdhyc were 61.5°C and 58.8°C respectively. The PCR condition consisted of an initial step of denaturation at 95°C for two minutes and followed by 30 cycles of denaturation at 94°C for one minute, with annealing at 58.1°C for one minute, extension at 72°C for one minute and last cycle of final extension at 72°C for 10 minutes. The stock solution was composed with MgCl2, dNTP (Sigma-Aldrich, USA), PCR buffer, forward primer, reverse primer, Taq DNA polymerase (Sigma-Aldrich, USA) and DNA for amplification [Table/Fig-1]. The PCR product was electrophoresed on 1% agarose gel along with 1 Kb ladder (Promega, USA). The amplified product was eluted using the gel extraction kit (Qiagen, Germany). Vector (pET-28a) was isolated from the fresh culture using plasmid isolation kit (Qiagen, Germany). The eluted rLdhyb amplicon of 771 bp was inserted between BamHI and HindIII sites and the rLdhyc amplicon of 1176 bp was inserted in BamHI and XhoI sites of isolated expression vector pET28a (+) (Novagen). Ligation was done using T4 DNA Ligase (Fermentas, Thermo Fisher Scientific, USA). Prior to transfection, chemically competent cells were prepared using 0.1 M CaCl2, followed by heat shock at 42°C. Freshly prepared competent Escherichia coli DH5α cells were transfected with the recombinant plasmid. The transformation was confirmed by colony PCR and restriction double digestion using respective restriction enzymes.

Components of PCR.

PCR ComponentsStock SolutionWorking SolutionTotal reaction Volume (25 μL)
MgCl225 mM1 mM2.5 μL
dNTP10 mM200 μM2.0 μL
PCR buffer10X1X2.5 μL
Forward Primer10 μM1 μM1 μL
Reverse Primer10 μM1 μM1 μL
Taq Polymerase2.5 u/μL2.5 u/100 μL0.5 μL
DNAVariable500 ng to 1 μg2 μL
H2O--13.5 μL

DNA Sequence Analysis and Alignment

Gene inserts of rLdhyb and rLdhyc clone were sequenced by Sanger sequencing methods. The cloned genes insert were isolated from vector pET28a with the help of respective restriction enzymes. To verify certainty of the gene cloned in vector, specific forward and reverse primers for rLdhyb and for rLdhyc were used. DNA sequencing was done commercially by Eurofins Genomics. Chromatogram analysis was performed by means of Finch TV 1.4.0 software. Homology searches and sequence alignment were performed using the nucleotide BLAST program http://www.ncbi.nlm.nih.gov/BLAST.

Expression and Purification of rLdhyb and rLdhyc

pET-28a-rLdhyb and pET-28a-rLdhyc construct were isolated from DH5-α and retransformed in competent Escherichia coli BL-21 (DE3, an expression host for recombinant construct). Recombinant proteins were produced by Isopropyl β-D-thiogalactoside (IPTG) induction at 0.8 mM (rLdhyb) and 1 mM (rLdhyc) for 4 hours at 25°C and at 1 mM for overnight at 37°C for rLdhyb and rLdhyc respectively. Expression was confirmed by SDS-PAGE following the protocol of Jamal F et al., [15].

A 600 μL Ni-NTA bead (Qiagen, Germany) in a column was used to purify proteins. The purity of eluted pET-28a-rLdhyb and pET-28a-rLdhyc was analysed by 10% SDS-PAGE (Laemmli, 1970). The Concentration of purified recombinant protein was estimated by the Bradford method using Bovine Serum Albumin (BSA) as standard. SDS-PAGE and western blotting was performed according to the protocol mentioned by Jamal F et al., [15].

ELISA for Detection of rLdhyb and rLdhyc Antibody in VL Patients

For determining most appropriate antigen concentration and antibody dilution, titration was performed at different antigen concentration and antibody serial dilution. Antibody concentration was determined from 1: 50 to 1: 51,000 dilutions. Antigen in form of rLdhyb and rLdhyc was taken in the concentration of 0.5 μg, 1 μg and 1.5 μg per well. ELISA was performed following the protocol mentioned by Jamal F et al., [18].

Statistical Analysis

Statistical analysis was performed using Graph Pad Prism (version 6.0). The lower limits of positivity (cut-off) for the diagnostic antigens were established for optimal sensitivity and specificity. ROC curves were used to determine the ELISA cut-off, sensitivity, specificity and AUC. The cut-off values (dotted line) for negative and positive sample discrimination were calculated using the mean±three times the standard deviation of all negative samples.

Result

In Silico B-Cell Epitopes Prediction

A large number of linear B-cell epitopes were found to be present in recombinant proteins rLdhyb and rLdhyc (XP_003860226.1 and XP_003861271.1 respectively), analysed from different web-servers. [Table/Fig-2a,b], demonstrates numerous B-cell epitopes present in the proteins under study utilising IEDB analysis with the default threshold value set at -0.5. Common best eight epitopes from different servers were selected [Table/Fig-3,4]. The proteins carrying epitopes with good blast score were carried forward for evaluation as a diagnostic candidate. In-depth conservancy analysis of epitopes revealed ~absolute value within species causing the visceralized form of Leishmaniasis [Table/Fig-5,6]. Comparatively B-cell epitopes exhibited meager similarity with Homo sapien [Table/Fig-5,6]. On the other pointer, in silico analysis of whole Ldhyb protein of L. donovani depicted 99% identity with L. infantum, whereas L. braziliensis and L. panamensis shows 68% and 67% identity with 100% query cover. L. major and L. maxicana depicted 92% and 88% identity with 67% query cover. In silico data showed that the target peptide might produce antibodies not only against L. donovani but also against other Leishmania species. Likewise, Ldhyb of L. donovani revealed 99% identity with L. infantum.Ldhyc also revealed 79%-91% identity to other Leishmania spp. with 100% query cover. Ldhyb displayed only 50% identity to Leptomonas spp. with 86-98% query cover. It shows only 32% identity with 40% query cover with Trypanosoma vivax. There was no similarity between Homo sapiens and causative agents of other diseases. Hence, these sequences were conserved among different species of Leishmania.

Linear B-cell epitopes predicted in Ldhyb (XP_003860226.1): (a) and Ldhyc (XP_003861271.1); (b) proteins at default threshold value set at -0.5 using IEDB analysis.

Linear B-cell epitopes prediction of rLdhyb (XP_003860226.1) from different servers.

RankPositionABCPred (Sequence)BepiPredEPMLR
137-52TGARSTARERRYYTQPTGARSTARERRY----TGARSTARERRYYTQP
260-75TSTIVGRQKSDADAAATSTIVGRQKSDADAAATSTIVGRQKSDADAAA
3100-115GHPRSTTRCRTAHQSPGHPRSTTRCRTAHQSPGHPRSTTRCRT-----
753-68LRTPQYGTSTIVGRQKLRTPQYGTSTIVGRQKLRTPQYGTSTIVGRQK
972-87DAAATWSGEASSPVAEDAAATWSGEASSPVAEDAAATWSGEASSPVAE
945-60ERRYYTQPLRTPQYGTERRY-TQPLRTPQYGTERRYYTQPLRTPQYGT
13227-242REIALHTLCNVFAYVMREIALH---REIALH-LCNVFAYVM
16151-166QLRLCTAILAAFYITKQL…….AFYITKQLRLCTAILAA----K

Linear B-cell epitopes prediction of rLdhyc (XP_003861271.1) from different servers.

PositionABCPred (Sequence)BepiPredEPMLR
311-326TVTKTRPVEVGEPLYSTVTKTRPVEVGEPLYSTVTKTRPVEVGEP
123-138QRSSDSASGCDARASSQRSSDSASGCDARASSQRSSD---------SS
147-162LEERHTECFPDPRNDALEER-TECFPDPRNDALEERHT-PDPRNDA
170-185PSSSSSTSASDSVDTAPSSSSSTSASDSVDTAP--SSTSASDSVDTA
294-309GPASVDPPAAGAITDEGPASVDPPAAGAITDEGPASVDPPAAGAITDE
284-299TRKRTVPLHFGPASVDGPASVDTRKRTVPLHFGPASVD
190-205TAAPPSSSPPAAVTPPTAAPPSSSPPAAVTPPAAPPSSSPPAAVTPP
129-144ASGCDARASSVTPLSPASGCDARASSVTPLSP------------SSVTPLSP

Conservancy analysis of predicted B-cell epitopes of rLdhyb.

PositionPeptideL. donovaniL. infantumHomo sapien
37-52TGARSTARERRYYTQP100 (100)100 (100)73 (68)
60-75TSTIVGRQKSDADAAA100 (100)100 (100)71 (75)
100-115GHPRSTTRCRTAHQSP100 (100)100 (100)75 (50)
53-68LRTPQYGTSTIVGRQK100 (100)94 (100)83 (75)
72-87DAAATWSGEASSPVAE100 (100)100 (100)58 (75)
45-60ERRYYTQPLRTPQYGT100 (100)100 (100)78 (62)
227-242REIALHTLCNVFAYVM100 (100)100 (100)67 (93)
151-166QLRLCTAILAAFYITK100 (100)100 (100)70 (68)

*Identity (Query cover)


Conservancy analysis of predicted B-cell epitopes of rLdhyc.

S No.PositionPeptideL. donovaniL. infantumHomo sapien
1311-326TVTKTRPVEVGEPLYS100 (100)100 (100)88 (50)
2123-138QRSSDSASGCDARASS100 (100)100 (100)64 (68)
3147-162LEERHTECFPDPRNDA100 (100)100 (100)59 (87)
4170-185PSSSSSTSASDSVDTA100 (100)100 (100)80 (56)
5294-309GPASVDPPAAGAITDE100 (100)100 (100)89 (56)
6284-299TRKRTVPLHFGPASVD100 (100)100 (100)89 (58)
7190-205TAAPPSSSPPAAVTPP100 (100)100 (100)60 (75)
8129-144ASGCDARASSVTPLSP100 (100)100 (100)64 (87)

*Identity (Query cover)


Cloning, Expression and Purification of Hypothetical Proteins rLdhyb and rLdhyc

rLdhyb and rLdhyc genes were successfully amplified from genomic DNA isolated from the reference strain of L. donovani (MHOM/IN/83/AG83). Amplification of gene rLdhyb and rLdhyc appeared at ~771bp and ~1.2 Kb region [Table/Fig-7a,8a: Lane 2], which was evaluated in comparison to 1 Kb ladder [Table/Fig-7a,8a: Lane 1]. Double restriction digestion was done with their respective enzymes to confirm cloning [Table/Fig-7b and 8b: Lane 1 and 2]. Western blotting was performed with an anti-His monoclonal antibody, which confirmed the expression and purification of recombinant proteins as a single over-expressed band. Confirmation of purification was shown in [Table/Fig-9] whereas [Table/Fig-7b (Lane 3), 8b (Lane 3)] being broad range (10-250 kDa) molecular weight marker (PureGene, Genetix Brand) respectively.

PCR and confirmation of transformation: a) PCR showing Ldhyb gene amplification. Gene was amplified using a forward 5’TTT GGA TCCATG CTG CGT TTC TGC3’ and reverse 5’AAT AAG CTT TCA GCA CAC CGT CCG3’ primers for rLdhyb using L. donovani genomic DNA, lane 1: Molecular weight marker (1 Kb DNA ladder); lane 2: Amplified PCR gene product; 7b) Double restriction digestion with restriction enzymes to confirm transformation.

PCR and confirmation of transformation: a) PCR showing Ldhyc gene amplification. Gene was amplified using a forward (5’-GTC GGA TCC ATG AAG CCG CTA GTG -3’) and a reverse (5’-GTC CTC GAG TCA GCT GTT GTT CTT-3’) primers for rLdhyc using L. donovani genomic DNA; b) Double restriction digestion with restriction enzymes to confirm insert into vector (lane 2) and only vector (lane 3).

Western blotting showing confirmation of purification with anti-His antibody for rLdhyc, rLdhyb antigenic proteins and broad range MW marker (10-250 kDa) in lane 1, 3 and 2 respectively.

Sequence Analysis of Cloned Genes Inserts using Sanger DNA Sequencing Method

The gene insert in rLdhyb and rLdhyc were sequenced and analysed. The nucleotide sequences of rLdhyb have been deposited in GeneBank under accession number “MH479406” and sequence for rLdhyc was being deposited in GenBank. Both clones were double digested with their corresponding restriction enzyme. Hypothetical protein Ldhyb is an exon region encoded on chromosome number 19; NC_018246.1, encrypted by mRNA (LDBPK_190480) complete cds, comprising 771 bp linear mRNA. Chromatogram acquired by forward primer sequencing revealed 98% identity with 100% query cover (97AGAAGCT-AGCAGC526) with L. donovani strain MHOM/IN/1983/AG83 isolate late passage chromosome 19 accession number CP019526.1 [Table/Fig-10a]. Similarly, reverse primer chromatogram also exhibited 98% identity with 100% query cover (70GGCGCC-TATAGT600) from L. donovani strain CP019526.1 [Table/Fig-10b]. Only 1% gaps were revealed in both ways sequencing owing to nucleotide deletion. Likewise, hypothetical protein Ldhyc is an exon region encoded on chromosome number 24 and encrypted by mRNA (LDBPK_241810) complete cds, comprising 1176 bp linear mRNA. Chromatogram peaks revealed 94% identity with 99% (161CTTGTA-CGAAGG709) and 100% (146ACGCCG-TCATCT887) query cover respectively from forward and reverse primer sequencing respectively, with L. donovani strain Pasteur chromosome number 24 complete sequence; sequence ID: CP022639.1, with 5% Gaps representing nucleotide deletion and insertion [Table/Fig-11a,b]. One and 22 nucleotide insertion site has been identified in cloned gene rLdhyc forward and reverse sequencing respectively. Percentage nucleotide sequence identity between different species has been represented in [Table/Fig-12].

Alignment analysis of cloned gene of rLdhyb showing sequence similarity with L. donovani strain MHOM/IN/1983/AG83 isolate late passage chromosome 19 sequence, Sequence ID: CP019526.1: 10a) DNA sequencing of the cloned gene using the forward primer; 10b) DNA sequencing of the cloned gene using the reverse primer.

Alignment analysis of cloned gene of rLdhyb showing sequence similarity with L. donovani strain pasteur chromosome 24, complete sequence, Sequence ID: CP022639.1. 11a) DNA sequencing of cloned gene using forward primer. 11b) DNA sequencing of cloned gene using reverse primer.

Similarity of cloned sequence with different Leishmania species.

Leishmania speciesCloned gene (Percentage Identity)
rLdhybrLdhyc
Forward Primer SequencingReverse Primer SequencingForward Primer SequencingReverse Primer Sequencing
Leishmania donovani98989494
Leishmania infantum97989494
Leishmania major92938990
Leishmania maxicana9089--

ELISA

Indirect ELISA was used to determine the reactivity of the recombinant proteins, rLdhyb and rLdhyc against the serological panel. The peptide reactivity was measured with sera from the VL-BT and control groups. Titration curves were performed with ELISA to determine the most appropriate antigen concentration and antibody dilution to be used. One microgram of rLdhyb and rLdhyc were found to be appropriate. Frequency distribution data revealed that ODs for the VL subjects in case of rLdhyb range from 0.39 to as high as 1.019. ODs value for healthy endemic and healthy non-endemic subjects range from 0.11-0.48 and 0.04-0.45, respectively. ODs value for viral flu, leprosy, tuberculosis, malaria, and asthma was 0.21-0.51, 0.07-0.50, 0.047-0.53, 0.15-0.48, 0.17-0.44, respectively, [Table/Fig-13a]. Frequency distribution data revealed that ODs for the VL subjects in case of rLdhyc range from 0.6 to as high as 1.53. ODs value for healthy endemic and healthy non-endemic subjects range from 0.17-0.69 and 0.16-0.67, respectively. ODs value for viral flu, leprosy, tuberculosis, malaria, and asthma were 0.31-0.78, 0.22-0.84, 0.23-0.77, 0.15-0.67, 0.23-0.88 respectively [Table/Fig-13b]. In case of SLA these parameters vary from 0.96-2.00, 0.11-0.48, 0.07-0.4, 0.45-1.16, 0.26-0.97, 0.31-0.84, 0.44-1.12 and 0.38-1.08 [Table/Fig-13c]. The lower limits of positivity (cut-off) for the diagnostic antigens rLdhyb, rLdhyc and SLA were 0.61, 0.95 and 1.31 [Table/Fig-14a-c]. The p-value for each test was <0.0001. Statistical analysis (unpaired t-test) considering VL-BT subject as a test group, whereas healthy and other diseases subjects as control group. The mean, standard deviation, standard error of mean for rLdhyb was calculated to be (0.78, 0.30), (0.15, 0.11), and (0.33, 0.01) respectively [Table/Fig-14a]; for rLdhyc as (1.23, 0.5), (0.24, 0.15) and (0.05, 0.02) respectively [Table/Fig-14b], and for SLA as (1.44, 0.54), (0.38, 0.27) and (0.08, 0.03) respectively [Table/Fig-14c]. Unpaired test was significant with p-value <0.0001. Area under curve for rLdhyb, rLdhyc and SLA were 0.99 and 0.99 and 0.961, with standard error 0.002, 0.007 and 0.019 respectively [Table/Fig-15a-c]. The ELISA revealed the specificity 100%; SpCl (84.56 to 100%); sensitivity 95.4%; Se Cl 95% (87.23-98.57%) for rLdhyb using software-based statistical analysis by GraphPad Prism (version 6.0). For rLdhyc, ELISA revealed the specificity 100%; SpCl (84.56-100%); sensitivity 91%; Se Cl95% (82.16-96.27%), whereas for SLA these data are 100%; SpCl (84.56 to 100%) and 78%; Se Cl95% (67.02-86.58%) respectively. A good diagnostic efficacy of rLdhyb and rLdhyc were 98% and 97%.

ELISA results showing the presence of antibody against rLdhyb, rLdhyc and SLA. Fig 13a, 13b and 13c depicts the frequency distribution graph of antibody for different study groups. The wells were coated with a concentration of 1 μg/well in 100 μL portion of rLdhyb, rLdhyc and SLA. Primary antibody was added in the concentration of 1:100. HRP tagged anti-human secondary antibody was administered in 1:2000 dilution. The ODs was recorded at 450 nm. The ODs value for the VL subjects was two-fold higher than healthy and other disease subjects in case of rLdhyb and rLdhyc.

ELISA results representing the lower limit of positivity. Cut-off value for rLdhyb, rLdhyc and SLA were 0.61, 0.95 and 1.31, respectively. The cut-off values (dotted line) for negative and positive sample discrimination were calculated using the mean±three times the standard deviation of all negative samples.

ROC curve. The ELISA revealed the specificity 100% and sensitivity of 95.4% for rLdhyb using software-based statistical analysis by GraphPad Prism (version 6.0). For rLdhyc, ELISA revealed the specificity 100% and a sensitivity of 91%, whereas for SLA these data were 100% and 78% respectively, with the area under curve 0.99, 0.99 and 0.96 respectively, represented by ROC curve (Fig. 15a-c).

Discussion

Labour-intensive practices such as microscopy still endure the mainstay of parasite detection in laboratories. However, over the past few decades, these diagnostic methods have proven to be highly unreliable. Therefore, there is a need for research emphasising the development of more certain tests that do not sacrifice sensitivity and specificity and can be utilised in both clinical settings as well as in poor resource field settings. The arenas of diagnostic parasitology, treatment, and vaccines are undergoing vivid change. In recent years, there has been a tremendous effort to replace the microscopic examination of invasive splenic aspiration, with newer diagnostic methods focusing on serological, molecular, and proteomic approaches.

In the current scenario, the overall sensitivity of available VL diagnostic methods appears highly unsatisfactory. The procedure like splenic bone marrow and the lymph node aspiration are highly invasive and the sensitivity of the parasitological method on aspirates of bone marrow and the lymph node is much lower, and up to 50% of cases of VL are expected to remain untreated if the treatment is conditional on diagnosis confirmation [20]. Upon activation of B-cells specific antibodies are produced which persists comparatively for a long time in patient’s sera and may give false positive results in serodiagnosis [21,22]. Concerning serodiagnosis, it has also been reported that when soluble antigens are used, serological tests are less specific due to cross-reactivity with other parasitic diseases [23]. In the last decades, an increasing number of recombinant protein candidates have been proposed to replace the crude Leishmania antigen (SLA) for the serodiagnosis of leishmaniasis [23]. Major technological advances in recombinant antigens as reagents for the serological diagnosis of VL have led to high sensitivity and specificity of these serological tests. In the existing system, rK39 is the widely used recombinant antigen for VL diagnosis. However, it detects anti-rk39 antibodies in 20-32% of endemic healthy individuals [24]. Further, a report from Sudan has shown the 67% sensitivity of the rK-39 test in immuno-compromised patients [25]. Thus, many efforts have been done to develop more sensitive and specific tests for the detection of VL helping in epidemiology and to control the disease. Hence, there have been the waves of gusto preferring one technique over others.

An endeavour of this study was to ascertain an L. donovani specific antigen which would be highly sensitive and specific. In our earlier study, we identified the B-cell epitopes of L. donovani circulating antigens [15], which specified the presence of B-cell epitopes of two L. donovani hypothetical proteins XP_003860226.1 and XP_003861271.1. This led to the discovery of rLdhyb and rLdhyc proteins as a diagnostic antigen. In the present study employing antigen drawn from blood circulation, hypothetical proteins proved to be very promising with the high degree of specificity and sensitivity. Epitomisation of hypothetical proteins in term of their B-cell epitopes possession exposed a large number of potential B-cell epitopes. Although, majority of B-cell epitopes are conformational or discontinuous; but linear B-cell epitopes are widely accepted for experimental purpose. B-cell epitopes are potentially recognisable by the immune system having the capacity to be bound by antibodies. Towards target, the first phase of the study revealed the epitope nature of targeted proteins. Eight best-matched peptides shortlisted from different servers revealed an excellent B-cell epitopes nature for both rLdhyb and rLdhyc protein [Table/Fig-3,4]. A similar assumption was advocated by Immune Epitope Database (IEDB) analysis Resource, Bpipred, BCpred etc.

Sequence homology analysis of rLdhyb and rLdhyc B-cell epitopes displayed 100% identity with L. infantum [Table/Fig-5,6]. On the other hand, whole protein homology analysis also revealed 99% similarity with L. infantum, mentioned in the result section thus, exhibiting its potential for VL diagnosis. Data also advocated the conservancy of these proteins among different species of Leishmania. Moreover, no any significant similarity could be derived between humans and the causative agent of other diseases resembling in symptoms to VL like tuberculosis, typhoid, malaria etc., henceforth no cross-reactivity and false positive result leading to its absolute specificity.

For the applied investigation hypothetical proteins were successfully cloned, expressed and purified. Further, the logical step was to characterise and prepare recombinant antigens. Western blot result [Table/Fig-9] with anti-His tag monoclonal antibody further authenticated the purified proteins. DNA sequencing data of cloned gene Ldhyb and Ldhyc revealed 98% and 94% identity with L. donovani Gene sequence ID CP019526.1 and CP022639.1 respectively. Cloned gene Ldhyb revealed high similarity with Leishmania spp. hence the likelihood of novel strain of same species. Ldhyc revealed less similarity due to a high percentage of nucleotide deletion and insertion. The results indicated that clones are homologous to the theoretical gene deposited in the Genbank. Also, sequencing result from both the primers revealed the same gene, which data further validated the sequencing results. Sanger sequencing data analysis for different Leishmanial spp. also revealed sequence similarity within species [Table/Fig-12]. To date, several Leishmania antigens have been genetically and antigenically characterised and recombinant technology has been used for the development of novel immunoassays for serological diagnosis of infections [14,26-29]. This study is the novel work where these hypothetical proteins of L. donovani have been characterised and sequenced and reported for the first time.

The study is among the first few works that verify biomarkers from L. donovani hypothetical recombinant proteins for VL patient’s diagnosis. In blood circulation of VL patient, more than 50% of the L. donovani antigens were represented as hypothetical protein [30]. Therefore, a complete hypothetical protein containing a large number of B-cell epitopes might be a better candidate for diagnosis. In silico data was validated by using the recombinant form of proteins. A pertinent and operative approach for the production of large quantities of purified antigens to be used in serological tests is the prokaryotic or eukaryotic systems based on recombinant DNA technology. Thus, providing the best alternative for the serodiagnosis of VL in place of crude antigens used in DAT.

Diagnostic value of immunoinformatics data was authenticated with indirect ELISA, even though it is well known that the short peptides may induce a specific response while large peptide may evoke a non-specific response due to the specificity of antigen. Both the targeted protein under study demonstrated absolute specificity (100%: SpCl (84.56-100)%) on the basis of determined cut-off, due to the presence of enormous epitopes. This excellent specificity is more than any other diagnostics for Kala-azar including rK39 (97%), DAT-FD (Direct Agglutination Test based on freeze-dried-94%), DAT-FDF (Direct Agglutination Test based on freeze-dried using filter paper-94%), Katex (99%), except rK26 which demonstrated 100% value [31]. Although rK39 is associated with high sensitivity, its specificity is a problem in the Indian subcontinent. A large proportion of the healthy population, living in the endemic areas may mimic the signs and symptoms of VL, in case of other diseases like malaria, typhoid fever, tuberculosis, etc., and give rK39 positive result [31]. Thus, there is a crucial need of a diagnostic that could precisely detect active VL.

The ELISA revealed the sensitivity of 95.4%; SeCl95% (87.23-98.57%) for rLdhyb which is better than kinesin based rK26 (22.6%) and antigen detection based KAtex (73%). The sensitivity of rBHUP1 was mentioned as 96.5% however, its low specificity restricted its utility [24]. Earlier, rA2 and BHUP2 proteins were also recognised as a potential target for diagnosis but failed to establish as diagnostics [23,32]. For rLdhyc, ELISA shows sensitivity 91%; SeCl95% (82.16-96.27%). Whereas, for SLA the data is 100%; SpCl (84.56-100%) and 78%; SeCl95% (67.02-86.58%) respectively. Present targeted antigens show better potential than previously documented antigens. The lower sensitivity might be due to non-accessible epitopes region of the cloned protein. Recombinant proteins under study reveal the diagnostic efficacy of 98% and 97% respectively. Hence, further strengthen these antigens’ biomarker property. The recombinant proteins presented a possible applicable approach for VL serological diagnosis since these proteins gave a positive result in the immune assay using many different sera from VL patients and negative results with control groups. Thus, these recombinant proteins with high specificity and sensitivity provide the initial findings for the development of VL specific diagnostics.

Limitation

The studied hypothetical proteins of Leishmania donovani were successfully purified and evidenced their diagnostic potency but in vitro functional characterization of the novel proteins is incumbent. Further large scale validation and authentication of same proteins on different ethnic population is obligatory.

Conclusion

The uniqueness of these L. donovani antigens might be a reason for its high specificity and diagnostic efficacy. As these proteins displayed conservancy of protein sequence within species, these might be used for diagnosing Leishmanial disease other than VL. Moreover, these proteins both in silico and in vitro data established no significant similarity between humans and causative agents of any other diseases. In the present investigation, rLdhyb and rLdhyc antigens revealed an increased ability to discriminate diseases, in comparison to SLA. Though the absolute specificity of rLdhyb and rLdhyc antigens for VL would have made them good diagnostic biomarkers, nevertheless, rLdhyc antigen sensitivity limit its utility. Further evaluation of the different ethnic population may justify their sensitivity. The present findings provide some basic insights for the future development of hypothetical proteins based non-invasive diagnostic tool.

Financial or Other Competing Interests

None.

References

[1]WHO. World Health Organization. The disease. What is Leishmaniaisis. 2017  [Google Scholar]

[2]Dikhit MR, Purkait B, Singh R, Sahoo BR, Kumar A, Kar RK, Activity of a novel sulfonamide compound 2-nitro-N-(pyridin-2-ylmethyl) benzenesulfonamide against Leishmania donovani Drug Design, Development and Therapy 2016 10:175310.2147/DDDT.S9665027307706  [Google Scholar]  [CrossRef]  [PubMed]

[3]Huda MM, Chowdhury R, Ghosh D, Dash AP, Bhattacharya SK, Mondal D, Visceral leishmaniasis-associated mortality in Bangladesh: a retrospective cross-sectional study BMJ Open 2014 4(7):e00540810.1136/bmjopen-2014-00540825052172  [Google Scholar]  [CrossRef]  [PubMed]

[4]Das VN, Pandey RN, Siddiqui NA, Chapman LA, Kumar V, Pandey K, Longitudinal study of transmission in households with visceral leishmaniasis, asymptomatic infections and PKDL in highly endemic villages in Bihar, India PLoS Neglected Tropical Diseases 2016 10(12):e000519610.1371/journal.pntd.000519627974858  [Google Scholar]  [CrossRef]  [PubMed]

[5]NVBDCP. National Road Map for Kala-Azar Elimination. Directorate of National Vector Borne Disease Control Programme, 2014. Available from: http://nvbdcp.gov.in/Doc/Road-map-KA_2014.pdf  [Google Scholar]

[6]Boelaert M, Bhattacharya S, Chappuis F, El Safi SH, Hailu A, Mondal D, Evaluation of rapid diagnostic tests: visceral leishmaniasis Nature Reviews Microbiology 2007 5:S30-39.10.1038/nrmicro1766  [Google Scholar]  [CrossRef]

[7]Sundar S, Rai M, Laboratory diagnosis of visceral leishmaniasis Clinical and Diagnostic Laboratory Immunology 2002 9(5):951-58.10.1128/CDLI.9.5.951-958.200212204943  [Google Scholar]  [CrossRef]  [PubMed]

[8]Terán-Angel G, Schallig H, Zerpa O, Rodríguez V, Ulrich M, Cabrera M, The direct agglutination test as an alternative method for the diagnosis of canine and human visceral leishmaniasis Biomédica 2007 27(3):447-53.10.7705/biomedica.v27i3.20718320110  [Google Scholar]  [CrossRef]  [PubMed]

[9]Jensen AT, Gasim S, Moller T, Ismail A, Gaafar A, Kemp M, Serodiagnosis of Leishmania donovani infections: assessment of enzyme-linked immunosorbent assays using recombinant L. donovani gene B protein (GBP) and a peptide sequence of L. donovani GBP Transactions of the Royal Society of Tropical Medicine and Hygiene 1999 93(2):157-60.10.1016/S0035-9203(99)90291-2  [Google Scholar]  [CrossRef]

[10]Mandal J, Khurana S, Dubey ML, Bhatia P, Varma N, Malla N, Evaluation of direct agglutination test, rk39 Test, and ELISA for the diagnosis of visceral leishmaniasis The American Journal of Tropical Medicine and Hygiene 2008 79(1):76-78.10.4269/ajtmh.2008.79.7618606766  [Google Scholar]  [CrossRef]  [PubMed]

[11]Chappuis F, Rijal S, Jha UK, Desjeux P, Karki BM, Koirala S, Field validity, reproducibility and feasibility of diagnostic tests for visceral leishmaniasis in rural Nepal Tropical Medicine & International Health 2006 11(1):31-40.10.1111/j.1365-3156.2005.01533.x16398753  [Google Scholar]  [CrossRef]  [PubMed]

[12]Kumar A, Misra P, Sisodia B, Shasany AK, Sundar S, Dube A, Mass spectrometry-based proteomic analysis of Leishmania donovani soluble proteins in Indian clinical isolate Pathogens and Disease 2014 70(1):84-87.10.1111/2049-632X.1210324115687  [Google Scholar]  [CrossRef]  [PubMed]

[13]De Souza RF, dos Santos YL, de Souza Vasconcellos R, Borges-Pereira L, Caldas IS, de Almeida MR, Recombinant Leishmania (Leishmania) infantum Ecto-Nucleoside Triphosphate Diphosphohydrolase NTPDase-2 as a new antigen in canine visceral leishmaniasis diagnosis Actatropica 2013 125(1):60-66.10.1016/j.actatropica.2012.09.01123022017  [Google Scholar]  [CrossRef]  [PubMed]

[14]Boarino A, Scalone A, Gradoni L, Ferroglio E, Vitale F, Zanatta R, Development of recombinant chimeric antigen expressing immunodominant B epitopes of Leishmania infantum for serodiagnosis of visceral leishmaniasis Clinical and Diagnostic Laboratory Immunology 2005 12(5):647-53.10.1128/CDLI.12.5.647-653.200515879027  [Google Scholar]  [CrossRef]  [PubMed]

[15]Jamal F, Shivam P, Kumari S, Singh MK, Sardar AH, Murugesan S, Identification of Leishmania donovani antigen in circulating immune complexes of visceral leishmaniasis subjects for diagnosis PloS one 2017 12(8):e018247410.1371/journal.pone.018247428820895  [Google Scholar]  [CrossRef]  [PubMed]

[16]Saha S, Raghava GP, Prediction of continuous B-cell epitopes in an antigen using recurrent neural network Proteins: Structure, Function, and Bioinformatics 2006 65(1):40-48.10.1002/prot.2107816894596  [Google Scholar]  [CrossRef]  [PubMed]

[17]Saha S, Raghava GP, BcePred: prediction of continuous B-cell epitopes in antigenic sequences using physico-chemical properties International Conference on Artificial Immune System 2004 :197-204.10.1007/978-3-540-30220-9_16  [Google Scholar]  [CrossRef]

[18]Jamal F, Dikhit MR, Singh MK, Shivam P, Kumari S, Pushpanjali S, Identification of B-cell Epitope of Leishmania donovani and its application in diagnosis of visceral leishmaniasis Journal of Biomolecular Structure and Dynamics 2017 35(16):3569-80.10.1080/07391102.2016.126324027892844  [Google Scholar]  [CrossRef]  [PubMed]

[19]Singh MK, Jamal F, Dubey AK, Shivam P, Kumari S, Ahmed G, Co-factor-independent phosphoglyceratemutase of Leishmania donovani modulates macrophage signalling and promotes T-cell repertoires bearing epitopes for both MHC-I and MHC-II Parasitology 2018 145(3):292-306.10.1017/S003118201700149429140228  [Google Scholar]  [CrossRef]  [PubMed]

[20]Rafati S, Hassani N, Taslimi Y, Movassagh H, Rochette A, Papadopoulou B, Amastin peptide-binding antibodies as biomarkers of active human visceral leishmaniasis Clinical and Vaccine Immunology 2006 13(10):1104-10.10.1128/CVI.00188-0617028214  [Google Scholar]  [CrossRef]  [PubMed]

[21]Londner MV, Feinsod FM, Faris R, Rosen G, El Said S, Saar AJ, Persistence of human leishmanial antibodies in an endemic area of visceral leishmaniasis in El Agamy, Egypt European Journal of Epidemiology 1988 4(4):473-76.10.1007/BF001464013203729  [Google Scholar]  [CrossRef]  [PubMed]

[22]Gidwani K, Picado A, Ostyn B, Singh SP, Kumar R, Khanal B, Persistence of Leishmania donovani antibodies in past visceral leishmaniasis cases in India Clinical and Vaccine Immunology 2011 18(2):346-48.10.1128/CVI.00473-1021159922  [Google Scholar]  [CrossRef]  [PubMed]

[23]Farahmand M, Nahrevanian H, Application of recombinant proteins for Serodiagnosis of visceral Leishmaniasis in humans and dogs Iran Biomed J 2016 20(3):128  [Google Scholar]

[24]Kumar D, Tiwary P, Dube A, Chakravarty J, Rai M, Sundar S, Cloning, Expression and Purification of L. Donovani Specific Antigen for Serodiagnosis of Visceral Leishmaniasis Journal of Molecular Biomarkers & Diagnosis 2013 4(1):100014110.4172/2155-9929.1000141  [Google Scholar]  [CrossRef]

[25]Zijlstra EE, Nur Y, Desjeux P, Khalil EA, El-Hassan AM, Groen J, Diagnosing visceral leishmaniasis with the recombinant K39 strip test: experience from the Sudan Tropical Medicine & International Health 2001 6(2):108-13.10.1046/j.1365-3156.2001.00680.x  [Google Scholar]  [CrossRef]

[26]Maalej IA, Chenik M, Louzir H, Salah AB, Bahloul C, Amri F, Comparative evaluation of ELISAs based on ten recombinant or purified Leishmania antigens for the serodiagnosis of Mediterranean visceral leishmaniasis The American Journal of Tropical Medicine and Hygiene 2003 68(3):312-20.10.4269/ajtmh.2003.68.31212685637  [Google Scholar]  [CrossRef]  [PubMed]

[27]Martin SK, Thuita-Harun L, Adoyo-Adoyo M, Wasunna KM, A diagnostic ELISA for visceral leishmaniasis, based on antigen from media conditioned by Leishmania donovani promastigotes Annals of Tropical Medicine & Parasitology 1998 92(5):571-77.10.1080/00034983.1998.11813314  [Google Scholar]  [CrossRef]

[28]Raj VS, Ghosh A, Dole VS, Madhubala RE, Myler PJ, Stuart KD, Serodiagnosis of leishmaniasis with recombinant ORFF antigen The American Journal of Tropical Medicine and Hygiene 1999 61(3):482-87.10.4269/ajtmh.1999.61.48210497995  [Google Scholar]  [CrossRef]  [PubMed]

[29]Ali D, Abbady AQ, Kweider M, Soukkarieh C, Cloning, expression, purification and characterization of Leishmania tropica PDI-2 protein Open Life Sciences 2016 11(1):166-76.10.1515/biol-2016-0022  [Google Scholar]  [CrossRef]

[30]Abass E, Bollig N, Reinhard K, Camara B, Mansour D, Visekruna A, rKLO8, a novel Leishmania donovani–derived recombinant immunodominant protein for sensitive detection of visceral leishmaniasis in Sudan PLoS Neglected Tropical Diseases 2013 7(7):e232210.1371/journal.pntd.000232223875052  [Google Scholar]  [CrossRef]  [PubMed]

[31]Sundar S, Singh RK, Bimal SK, Gidwani K, Mishra A, Maurya R, Comparative evaluation of parasitology and serological tests in the diagnosis of visceral leishmaniasis in India: a phase III diagnostic accuracy study Tropical Medicine & International Health 2007 12(2):284-89.  [Google Scholar]

[32]Kumar S, Kumar D, Chakravarty J, Sundar S, Identification and characterization of a novel, 37-kilodalton Leishmania donovani antigen for diagnosis of Indian visceral leishmaniasis Clinical and Vaccine Immunology 2011 18(5):772-75.10.1128/CVI.00559-1021411602  [Google Scholar]  [CrossRef]  [PubMed]