Ethical approval
This study was reviewed and approved by the Ethics Committee of the Department of Biological Sciences at Shaqra University, according to the ethical principles of animal research (protocol SH 2-2017).
Study areas
The investigation was conducted from January 2018 to May 2019 in Al-Qaseem province and Riyadh city, Saudi Arabia. Al-Qassim province is located in the central part of Saudi Arabia (latitude 25°–23° N and longitude 42°–24° E). It has an area of about 58,046 km2 and in 2017, it had a population of 1,423,000 people [19]. Al-Qassim province is known as an agricultural region and it has a typical desert climate, with an average temperature of 13˚C and a hot summer (an average temperature of 35˚C). It has low annual rainfall (214 mm) and low humidity that ranges from 25% to 76% (http://www.pme.gov.sa). Conversely, Riyadh city is the capital of Saudi Arabia (latitude 24°–08°N and longitude 47°–18° E). It covers an area of about 1,798 km2 and in 2017, was inhabited by approximately seven million people [19]. Riyadh is characterised by very hot summers with an average temperature of 45˚C in July, whereas winters are cold. The overall climate is arid, with scarce annual rainfall (21.4 mm) and a relative humidity that ranges from 10% to 47% throughout the year. Riyadh also has many dust storms (http://www.pme.gov.sa) (Fig. 1).
Patient biopsy tissue collection and DNA extraction
A total of 27 patients who were suspected to be infected by Leishmania species were seen in either King Saud Medical City in Riyadh city (n =16) or Buraydah Central Hospital (n =11) in Al-Qaseem province. The presence of Leishmania was investigated in all samples, which were collected after clinical and microscopy examination [20]. Briefly, skin biopsies (i.e., 5-10 mm in diameter) were taken under sterile conditions from the border of the ulcerous and cutaneous lesions and DNA samples were extracted from all biopsies by MagNaA pure DNA extraction through use of a Pure LC DNA Isolation Kit (Roche Applied Science, Germany) according to the manufacturer’s instructions. The extracted DNA was quantified by Nanodrop spectrophotometer (Thermo, USA). The DNA concentration differed from sample to sample, but it ranged from 18ng/ul to 33ng/ul. An aliquot (100µl of DNA from each sample) was stored at -20°C prior to nPCR amplification and analysis.
Sampling of stray dogs
From January 2018 to May 2019, 311 stray dogs were trapped in Al-Qaseem province by bait traps (Havahart®) (Fig. 1) and were examined physically for canine leishmaniasis skin lesions in the field. Seven dogs were suspected of infection with canine leishmaniasis due to the presence of cutaneous nodules or ulcerated lesions on the skin (Table 1). Skin biopsies (5mm in diameter) were collected under sterile conditions from the borders of the ulcers and were inoculated into M199 medium (Gibco, Life Technologies, Germany), which was supplemented with 25 mmol/l of 4-(2-hydroxyethyl)-1-piperazineethanesulphonic acid (HEPES) (pH:7.5) and 20% fetal bovine serum (Gibco, Life Technologies, Germany). These samples were then incubated at 24 °C. Ten days after sample incubation, parasites were harvested and washed with ice-cold phosphate-buffered saline (10X PBS, pH: 7.4) and stored at -20°C before DNA isolation. DNA from parasite cultures was isolated by use of the ReliaPrep™ gDNA Tissue Miniprep System Kit (Promega, Madison, United States), following the manufacturer’s instructions.
Leishmania nested PCR
The specific external CSB2XF primers (5ˊ-ATTTTTCGCGATTTTCGCAGAAACG-3ˊ) and CSB1XR (5ˊ-CGAGTAGCAGAAACTCCCGTTCA-3ˊ) were used initially. In the second step, specific internal 13Z primers (5ˊ-ACTGGGGGTTGGTGTAAAATAG-3ˊ) and LiR (5ˊ-TCGCAGAACGCCCCT-3ˊ) were applied [21]. The specificity and sensitivity of these primers are reported to be 92% and 100%, respectively [21]. These primers were able to track and multiply the variable part of all forms of the Leishmania kDNA. Amplified fragments of L. infantum were 680bp in length and fragments of L. tropica and L. major were 750bp and 560bp in length, respectively [21].
The first step of the addition of the PCR master mix, which included CSB2XF and CSB1XR, was performed using the AccuPower® PCR PreMix kit (Bioneer, Daejeon, Korea). The prepared PCR pre-mix volumes contained potassium chloride (KCl) at a concentration of 30mM, magnesium chloride (MgCl2) at 1.5mM, tris (hydroxymethyl) aminomethane hydrochloride (Tris-HCL at pH 9.0) at 10mM, Taq DNA polymerase, and deoxynucleoside triphosphate (dNTP) were adjusted to 2µl. In addition, 1μl of each initial CSB1XR and CSB2XF primer at concentrations of 10pmol/ul (Bioneer, Daejeon, Korea) and 3μl of DNA were added to the reaction mixture. Finally, 13μl of deionised water (ddH2O) were added up to a total volume of 20μl for reaction. Negative control was included in the final nPCR. The reaction was performed in a thermal cycler (Techne TC-3000, US) according to the following conditions: initial denaturation temperature of 94°C for 5min; 30 cycles of denaturation at 94°C for 30s, annealing at 55°C for 60s and extension at 72°C for 60s; final extension at 72°C for 7min; and then the reaction was held at 4°C. The second step of PCR involved 13Z and LiR primers and the same PCR master mix except that 3µl of template PCR product were used. The second round, PCR products that were obtained were electrophoresed on a 1.5% agarose gel that contained 1μl Syber safe (Thermo Scientific™, Nalgene, UK) in tris-acetate–ethylenediaminetetraacetic acid (EDTA) buffer (50X) at 100V for 45min and visualised under a UV imaging system (ImageQuant Laz4000, GE Healthcare Life Science, Hammersmith, UK). The size of each sample was estimated by comparison with a 100bp DNA Ladder Marker (Solis BioDyne OU, Estonia).
Leishmania kDNA sequencing and BLAST analysis
Positive amplified products of Leishmania species were sent to Macrogen (South Korea) for sequencing, and the results were compared with the sequences that were available in the GenBank database through use of BLAST (http://blast.ncbi.nlm. nih.gov/). The obtained sequences were aligned with a set of reference sequences that were available in GenBank using CLUSTALW in MEGA software version 7.0 [22]. The phylogenetic tree was constructed using the maximum-likelihood method and with the Hasegawa-Kishino-Yano (HKY) model with 2000 bootstrap replicates in MEGA 7.0 software, using Trypanosoma cruzi (AJ748063) as outgroup [22, 23].