2.1 Ethical approval
This study was revised and approved by the Ethical Research Committee, Department of Biological Science, Shaqra University, according to the ethical principles of human and animal research (Approval no. SH 06-2018).
2.2 Study areas
The investigation was conducted from November 2018 to August 2019 in two provinces of Saudi Arabia.The Riyadh province of Saudi Arabia has an area of 404,240 km² and is located in the central part of Saudi Arabia between 24˚.38˚N and 46˚.43˚E (Figure 1). This province is characterised by very hot summers with an average high temperature of 45˚C in July. Winters are cold, the overall climate is arid, receiving very little annual rainfall (21.4 mm), with the relative humidity ranging from 10% to 47% throughout the year. Riyadh province is also known to have many dust storms (http://www.pme.gov.sa). Asir province, has an area of 76,690 km² and is located in the southwestern part of Saudi Arabia between 19°0′N and 43°0′E (Figure 1). Asir province is situated on a high plateau that receives more rainfall than the rest of the country and contains the country's highest peaks, which rise to almost 3,000 m. Asir has a tropical and subtropical climate and the average annual rainfall in the highlands is expected torange from 300 to 500 ml across two rainy seasons. As a result, there is much more natural vegetation and forests (http://www.pme.gov.sa).
2.3.Sampling of dogs and cats and blood collections
Stray dogs and cats from the two provinces were trapped by a live bait traps (Havahart®) and selected randomly.A total of 74 dogs from Riyadh province (48 males and 26 females) and a total of 70 dogs and 44 cats from Asir province (dogs; 43 males and 27 females, cats; 30, male and 14 females). Dogs and cats varied in ages from ≤6 months to> 6 months and were examined forhaemoparasites. Most of the dogs and cats appeared healthy at the time of blood collection. Blood samples were collected from each animal (0.5–3 ml) from the cephalic vein into vacutainer tubes (BD Vacutainer® Tube, Gribbles Pathology, VIC, Australia) and transported to the parasitology laboratory, Department of Biological Sciences, Faculty of Science and Humanities, Shaqra University, for DNA extraction.
2.4 DNA extraction
Total genomic DNA (gDNA) was isolated from the blood samples using the Wizard® Genomic DNA Purification Kit (Promega, Madison,WI, USA) and eluted into 50 μl or 100 μl of elution buffer as per the manufacturer’s instruction. An aliquot between 50 μl and 100 µl of gDNA from each of the samples was stored at −80°C prior to being sent to Veterinary Pathology Diagnostic Services (VPDS), Sydney School of Veterinary Science, The University of Sydney for PCR analysis. Upon arrival atVPDS gDNA was stored at -20 ˚C for up to 1 month prior to molecular diagnostics.
2.5. Commercial Multiplexed Tandem PCR (MT-PCR) for Small Animal Anaemia
A commercial diagnostic MT-PCR panel for small animal anaemia was performed using the mini-plex 12 system (R910738, AusDiagnostics Pty. Ltd., Australia) as per the manufacturers’ instruction. The MT-PCR assay is a two-step nested PCR assay simultaneously targeting B. gibsoni, B. vogeli, M. haematoparvum, M. Haemocanis and A. platys and was run on the Easy-PlexTM platform (AusDiagnostics Pty. Ltd., Australia). The assay was run in duplicate using 10 µlundiluted samples (n=188). Each run included controls to detect PCR inhibition (SPIKE) and sample adequacy control (ANONO) as per the manufacturers’ protocol (AusDiagnostics Pty. Ltd., Australia).
The positive samples were tested using 16S rRNA primers (S0697/S0698) for Anaplasma platys [24], 18S rRNA primers (S0701/S0702) for the universal amplification of Mycoplasma spp. [25] and 18S rRNA primers for Babesia spp.[26]. All PCR reactions were run using MyTaqTMRedMix(Bioline, Australia) in a Veriti Thermal Cycler (Life Sciences, Australia). Primers were included at a final concentration of 400nM/µl and 2 µl of template DNA was used per reaction. All non-nested PCR were run using the following cycling conditions: 95 °C for 1 min and 35 cycles of 95 °C for 15 s, 50 °C for 15 s and 72 °C for 10 s followed by 72 °C for 5 min. PCR reactions with Babgen-F/Babgen-R primers were run with an annealing temperature of 55 °C, and all other conditions as previously described. The first round of the nested PCR was run using the following cycling conditions: 95 °C for 1 min and 35 cycles of 95 °C for 15 s, 55 °C for 15 s and 72 °C for 10 s followed by 72 °C for 5 min. The second round of the nested PCR was run using a 60 °C annealing temperature with all other conditions as described for the first round. A negative and positive control was included in all assays. All PCR products were separated by electrophoresis in 2% agarose gel stained with GelRedTM (Biotium, USA) and visualised using UV light. Discrete bands of expected size were submitted for bidirectional sequencing using amplification primers (Macrogen, South Korea). Sequences were assembled and compared to closely related sequences using CLC Main Workbench 6.8.1 (CLC bio, Denmark).
2.6 Real-time PCR: Bartonella spp.
A multiplex TaqMan probe real time PCR assay targeting the gltAgene of Rickettsia spp. was multiplexed with an assay targeting the ssrA gene of Bartonella spp. [27, 28, 29]. PCR reactions were run using SensiFASTTM Probe No-ROX Kit (Bioline, Australia) on a CFX95 TouchTM Real-Time PCR detection system (BioRad Laboratories Inc., Australia).
2.7.Statistical analysis
Statistical analyses were performed with the statistics package SPSS (v.17.0; IBM, New York, New York). Positive PCR test was set as an outcome variable and the independent variables were age, gender and, health status. The effect of independent variables on the outcome variables were evaluated by chi-square and Fisher's exact test and Odds Ratio (OR) calculation. Differences were considered significant if the P value was <0.05.