prevalence, risk factors and genetic diversity of

0 downloads 0 Views 6MB Size Report
Gen 18S rRNA T.equi dan B. caballi yang diasingkan daripada darah dan kutu dianalisa untuk ...... 18S rRNA genotype named genotype C (group C) novel was identified in. Jordanian equids ...... iastate.edu/ Fact sheets/pdfs/equine_piroplasmosis.pdf. Chahan, B. ...... muda forest reserve, Kedah Malaysia. Southeast Asian ...
PREVALENCE, RISK FACTORS AND GENETIC DIVERSITY OF EQUINE PIROPLASMOSIS IN KELANTAN, MALAYSIA

QAES TALB SHUKUR ALSARHAN

DOCTOR OF PHILOSOPHY

2017

8

Prevalence, Risk Factors and Genetic Diversity of Equine Piroplasmosis in Kelantan, Malaysia

by

QAES TALB SHUKUR ALSARHAN

A thesis submitted in fulfillment of the requirements for the degree of Doctor of Philosophy

Faculty of Veterinary Medicine UNIVERSITI MALAYSIA KELANTAN

2017

9

THESIS DECLARATION

I hereby certify that the work embodied in this thesis is the result of the original research and has not been submitted for a higher degree to any other University or Institution. OPEN ACCESS

EMBARGOES

I agree that my thesis is to be made immediately available as hardcopy or on-line open access (full text).

I agree that my thesis is to be made available as hardcopy or on-line (full text) for a period approved by the Post Graduate Committee. Dated from ___________ until ___________

CONFIDENTIAL

(Contains confidential information under the office Official Secret Act 1972)*

RESTRICTED

(Contains restricted information as specified by the organization where research was done) *

I acknowledge that Universiti Malaysia Kelantan reserves the right as follows. 1. The thesis is the property of Universiti Malaysia Kelantan. 2. The library of Universiti Malaysia Kelantan has the right to make copies for the purpose of research only. 3. The library has the right to make copies of the thesis for academic exchange.

_____________________ SIGNATURE _______________________

__________________________ SIGNATURE OF SUPERVISOR _________________________

IC/ PASSPORT NO.

NAME OF SUPERVISOR

Date:

Date:

10

ACKNOWLEDGMENT

First of all, my thanks are to my God. I would like to express my sincere thanks to my supervisor, Associated Professor Dr. Mohd Mokhtar Arshad for his advice, patience and enthusiasm throughout the years of work. I would also like to express my gratitude and thanks to my co-supervisors, Professor Dr. Imad Ibrahim Al-sultan, Professor Dr. Mohd Azam Khan Goriman Khan, Universiti Malaysia Kelantan (UMK) and Dr. Azlinda Abu Bakar, Universiti Sanis Malaysia (USM), for helpful discussions and comments on my thesis. I would also like to extend my gratitude and appreciation to Dr. Maizan Mohammed for helping and advising me on molecular work. I am indebted to the Faculty of Veterinary Medicine, UMK, for making this possible by providing all necessary chemicals and equipment in the laboratory. I am deeply and extremely grateful to all of the UMK laboratory assistants, especially to Mr. Nor Faizull, Mr. Badrul Hisham, Mrs Eizzati, and Miss Nani Izreen for their support in so many ways. I would also like to thank the UMK veterinary clinic staff, especially Dr. Mimi, Mr. Hamid, and Mr. Nizam for their help in the sample collections. Special acknowledgment goest to my friend, Dr. Omer Khazaal, Dr. Ali Saeed, Dr. Ahmad Mahmood, Dr. Maher Mohammed for their direct and indirect assistance during my PhD study.

Qaes

11

TABLE OF CONTENTS

PAGE THESIS DECLARATION

i

ACKNOWLEDGEMENTS

ii

TABLE OF CONTENTS

iii

LIST OF TABLES

xi

LIST OF FIGURES

xv

LIST OF ABBREIVATIONS

xx

LIST OF SYMBOLS

xxiii

LIST OF EQUATIONS

xxiv

ABSTRAK

xxv

ABSTRACT

xxvi

CHAPTER 1 INTRODUCTION 1.1

General introduction

1

1.2

Problem statement

5

1.3

Research questions

6

1.4

Hypothesis

6

1.5

The objectives of study

7

CHAPTER 2 LITERATURE REVIEW 2.1

History of equine piroplasmosis (EP)

8

2.2

Etiological agents and taxonomy

10

12

2.2.1 Theileria equi

12

2.2.1.1 Morphology

12

2.2.1.2 Life cycle

14

2.2.2 Babesia caballi

17

2.2.2.1 Morphology

17

2.2.2.2 Life cycle

18

2.3 The genes commonly targeted in the T. equi and B. caballi

21

2.4 Sequencing and genetic diversity for T. equi and B. caballi

22

2.5

Epidemiology of equine piroplasmosis

25

2.5.1 Geographic distribution

25

2.5.2 Susceptibility to the disease

33

2.6

2.7

2.5.2.1 The susceptibility related to equids factors

33

2.5.2.2 The susceptibility related to environmental and stables factors

36

Transmission of the causative agents

38

2.6.1 Biological transmission

38

2.6.2 Iatrogenic or mechanical transmission

40

2.6.3 Intrauterine or transplacental transmission

41

The tick vectors for equine piroplasms infections

42

2.7.1 Taxonomy of Ixodid ticks

46

2.7.2 Morphology of Ixodid ticks

47

2.7.2.1 Tick family identification

48

2.7.2.2 Tick genus identification

52

13

2.7.2.3 Tick species identification

53

2.7.3 Tick life cycle

53

2.7.4 Overview of the identified Ixodid ticks genera in this study

55

2.7.4.1 Genus Rhipicephalus Koch, 1844 (Soulsby, 1982)

55

2.7.4.2 Genus Haemaphysalis Koch, 1844 (Soulsby, 1982)

58

2.7.4.3 Genus Dermacentor Koch, 1844 (Soulsby, 1982)

59

2.8

Ixodid ticks in Malaysia

60

2.9

Pathogenesis of equine piroplasmosis

62

2.10 Clinical signs of equine piroplasmosis

66

2.11 Clinical pathology of equine piroplasmosis

69

2.11.1 Changes in hematological parameters

69

2.11.2 Changes in serum biochemistry parameters

73

2.11.2.1 Aspartate aminotransferase (AST)

74

2.11.2.2 Alanine aminotransferase (ALT)

74

2.11.2.3 Alkaline phosphatase (ALKP)

75

2.11.2.4 Total bilirubin

75

2.11.2.5 Total protein

76

2.11.2.6 Blood urea nitrogen (BUN)

76

2.11.2.7 Calcium

77

2.11.2.8 Glucose

78

2.11.2.9 Phosphorus

78

2.11.2.10 Creatinine

79

2.11.3 Pathological changes of equine piroplasmosis

14

79

2.11.3.1 Macroscopic findings

80

2.11.3.2 Microscopic findings

80

2.12

Immunity to equine piroplasmosis

81

2.13

Public health of significance of equine piroplasmosis

85

2.14

Diagnosis of equine piroplasmosis

86

2.14.1 Overview of the methods used for diagnosis EP in this study

87

2.14.1.1 Microscopic examination of stained blood smears

87

2.14.1.2 Competitive enzyme linked immunosorbent assay

89

2.14.1.3 Conventional and multiplex polymerase chain reaction

91

2.14.2 Other methods for diagnosis of equine piroplasmosis

94

2.14.2.1 Biological tests

94

2.14.2.2 In vitro culture technique

95

2.14.2.3 Other serological tests

96

2.14.2.4 Other molecular techniques

97

2.15 Differential Diagnosis of equine piroplasmosis

99

2.16 Prognosis of equine piroplasmosis

99

2.17 Treatment of equine piroplasmosis

100

2.18 Control of of equine piroplasmosis

103

2.18.1 Vaccination

103

2.18.2 Control of ticks

104

CHAPTER 3 DIAGNOSIS, PREVALENCE, RISK FACTORS AND VECTOR OF EQUINE PIROPLASMOSIS IN EQUIDS IN KELANTAN 3.1

Introduction

105

15

3.2

Materials and Methods

108

3.2.1 Study area

108

3.2.2 Determining the number of equids

109

3.2.3 Animals and sample collection

109

3.2.4 Epidemiological data collection

110

3.2.5 Climatic data

113

3.2.6 Laboratory analysis

113

3.2.6.1

Microscopic examination of blood smears

113

3.2.6.2

Competitive enzyme linked immunosorbent assay

115

3.2.6.3

Polymerase chain reaction techniques

120

3.2.6.3.1 DNA extraction from equids blood

120

3.2.6.3.2 Determination of DNA concentration and purity

123

3.2.6.3.3 PCR amplification of piroplasms DNA from equids blood samples

124

3.2.7 Evaluation the efficiency of different methods for diagnosis the disease 3.2.8 Statistical analysis

129

130

3.3 Results

131

3.3.1 Morphological and biometerical finding with parasitemia

131

3.3.2 Prevalence of EP by different tests

135

3.3.3 Evaluation of cELISA and multiplex PCR for detecting T. equi and B. caballi infections

137

3.3.4 Prevalence of T. equi, B. caballi and both protozoa infections by regions

16

143

3.3.5 Risk factors associated with seroprevalence of EP causative agents

150

3.3.6 Ixodid ticks: identification and infestation rate

163

3.4 Discussion

171

3.4.1 Morphological and biometerical finding with parasitemia

172

3.4.2 Determing the prevalence of EP infections using different methods

173

3.4.3 Evaluating of cELISA and multiplex PCR for detecting T. equi and B. caballi infections

177

3.4.4 Equids factors associated with T. equi, B. caballi and both protozoa

179

3.4.5 Stables factors associated with T. equi, B. caballi and both protozoa

182

3.4.6 Climatic factors associated with T. equi, B. caballi and both protozoa

185

3.4.7 Ixodid ticks: identification and infestation rate

186

3.5 Conclusions

189

CHAPTER 4 EVALUATION OF HEMATOLOGY, BIOCHEMISTRY AND ANTIBODY TITER BETWEEN EQUIDS CLINICALLY AND SUBCLINICALLY INFECTED WITH EQUINE PIROPLASMOSIS 4.1

Introduction

192

4.2

Materials and methods

194

4.2.1 Animals and study area

194

4.2.2 Recording clinical sings

194

4.2.3 Fecal samples collection

194

4.2.4 Determination of antibody titers against T. equi and B. caballi

195

4.2.5 Hematological analysis

195

17

4.2.6 Serum biochemistry analysis

196

4.2.7 Statistical analysis

197

4.3 Results

4.4

198

4.3.1 Seroprevalence of clinical and subclinical forms in equids with the relative risk

198

4.3.2 Level of antibodies titration against piroplasms in equids with clinical and subclinical forms of equine piroplasmosis

198

4.3.3 Hematological and serum biochemistry parameters in equids with clinical and subclinical form of equine piroplasmosis

199

Discussion

206

4.5 Conclusions

212

CHAPTER 5 GENETIC DIVERSITY AND PHYLOGENIC ANALYSES OF THEILERIA EQUI AND BABESIA CABALLI DETECTED IN EQUIDS AND IXODID TICKS IN KELANTAN 5.1

Introduction

214

5.2

Materials and Methods

217

5.2.1 Ticks collection from equids

217

5.2.2 DNA extraction from Ixodid ticks

217

5.2.3 PCR amplification of piroplasms DNA extracted Ixodid ticks

219

5.2.4 DNA sequencing and phylogeny analyses

219

5.2.5 Statistical analyses

221

5.3 Results

221

5.3.1 Detection rate of T. equi, B. caballi and both protozoa in Ixodid ticks using multiplex PCR

221

5.3.2 Similarity rates on detection of T. equi, B. caballi and both protozoa DNAs between equids and Ixodid ticks

222

18

5.3.3 Genotypes of T. equi and B. caballi in equids and Ixodid ticks

226

5.3.4 Similarity within and between the genotypes of T. equi and B. caballi

229

5.3.5 Phylogenic analysis of T. equi and B. caballi sequences

231

5.4 Discussion

240

5.5 Conclusions

244

CHAPTER 6 CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 6.1 General conclusions

246

6.2 Recommrndations

248

6.3 Future work

249

REFERENCES

250

APPENDIX- A

289

APPENDIX- B

293

APPENDIX- C

301

APPENDIX- D

319

LIST OF PUBLICATIONS

334

19

LIST OF TABLES

NO.

PAGE

2.1

Prevalence of T. equi, B. caballi and different countries with the references.

2.2

Hard ticks species with the geographic location, number of the hosts and confirmed or suspected to be transmitted of Theileria equi or Babesia caballi (the causative agents of EP)(Scoles & Ueti, 2015).

45

2.3

Ixodid tick species distribution in different states of Malaysia.

61

2.4

Different methods obtainable for the diagnosis of EP and their purpose (OIE, 2014b).

94

3.1

The Oligonucleotide primers used to amplify the parasites 18S rRNA genes.

125

3.2

PCR program for samples subject to conventional PCR and multiplex PCR.

125

3.3

PCR procedure for DNA extracted samples subject to multiplex PCR and conventional PCR assay (reaction volume 25 μl).

127

3.4

The epidemiological table and formulas were used to copare between tests.

130

3.5

Morphological features, biomedical data and parasitemia of T. equi, B. caballi and both protozoa base on microscopic examination of blood smears.

135

3.6

Overall prevalence of equine piroplasmosis (singly T. equi, singly B. caballi and both protozoa) in equids in Kelantan by microscopic examination, cELISA and multiplex PCR.

136

3.7

Prevalence of T. equi, B. caballi and both protozoa in equids in Kelantan by microscopic examination, cELISA and multiplex PCR.

137

20

both protozoa in the

29

3.8

The sensitivity, specificity, accuracy, positive predictive value and negative predictive value of cELISA for ditecting T. equi antibodies in comparison to microscopic examination (Gold standard).

139

3.9

The sensitivity, specificity, accuracy, positive predictive value and negative predictive value of multiplex PCR for detecting T. equi DNAs in comparison to microscopic examination (Gold standard).

139

3.10

The sensitivity, specificity, accuracy, positive predictive value and negative predictive value of cELISA for detecting B. caballi antibodies in comparison to microscopic examination (Gold standard).

140

3.11

The sensitivity, specificity, accuracy, positive predictive value and negative predictive value of multiplex PCR for detecting B. caballi DNAs in comparison to microscopic examination (Gold standard).

140

3.12

Comparison between conventional PCR and multiplex PCR techniques for detection piroplasms DNA in equids blood (n= 306).

141

3.13

Prevalence of T. equi, B. caballi and both protozoa in sampling stables in different regions and sub-regions in Kelantan state.

145

3.14

Relative risk of EP associated with the type of protozoa in equids.

154

3.15

Relative risk of equids factors associated with seropositivity of T. equi, B. caballi and both protozoa.

155

3.16

Relative risk of regional factors associated with seropositivity of T. equi, B. caballi and both protozoa.

157

3.17

Relative risk of monthly factors associated with seropositivity of T. equi, B. caballi and both protozoa.

158

3.18

Relative risk of management and ticks factors associated with seropositivity of T. equi, B. caballi and both protozoa.

160

3.19

Relative risk of climatic factors associated with seropositivity of T. equi, B. caballi and both protozoa.

161

21

3.20

Percentage of sables infested with Ixodid ticks by Kelantan regions.

166

3.21

Details and number of identified Ixodid tick species infesting equids and nearby animals related to the prevalence of EP in each field.

167

3.22

Distribution of identified Ixodid tick species and their predilection sites on equids and nearby animals.

169

3.23

The abundance and percentage of Ixodid ticks on equids and nearby animals.

170

4.1

Health status factor equids associated with seroprevalence of T. equi, B. caballi and both protozoa.

201

4.2

Relative risk of of equids health status factors associated with the T. equi, B. caballi and both protozoa.

201

4.3

The clinical symptom parameters in the equids with clinical form compared to the equids with subclinical form and healthy group.

204

4.4

Haematological changes in the equids with clinical form compared to the equids with subclinical form and healthy group.

205

4.5

Serum biochemistry changes in the equids with clinical form compared to the equids with subclinical infectedform and healthy group.

206

5.1

Detection of T. equi and B. caballi in Ixodid ticks (n= 31) using multiplex PCR.

222

5.2

Similarity on detection of T. equi and B. caballi in the equids and Ixodid ticks using multiplex PCR.

224

5.3

Detection rate of T. equi and B. caballi genotypes in equids base on individual BLASTn analysis of positive samples.

227

5.4

Detection rate of T. equi and B. caballi genotypes in Kelantan regions.

228

5.5

Detection rate of T. equi and B. caballi genotypes in Ixodid ticks (n=31) base on individual BLASTn analysis of sequences.

229

22

5.6

Similarity within and between T. equi genotypes using multiple sequence alignment- CLUSTALW (GenomeNet).

230

5.7

Similarity within and between B. caballi genotypes using multiple sequence alignment- CLUSTALW (GenomeNet).

230

5.8

GeneBank accession numbers of Kelantan T. equi genotypes in the equids and ticks.

234

5.9

GeneBank accession numbers of Kelantan B. caballi genotypes in the equids and ticks. Homology between obtained sequences (GeneBank accession numbers) of T. equi and GeneBank database using online sequence BLASTn.

235

Homology between obtained sequences (GeneBank accession numbers) of B. caballi and GeneBank database using online sequence BLASTn.

237

5.10

5.11

23

236

LIST OF FIGURES

NO.

PAGE

2.1

A) Different forms of T. equi inside the erythrocytes and the characteristic Maltese cross form in a Giemsa‘s stained blood smear 100X (S. Kumar & R. Kumar, 2007); B) Microschizontes and macroschizontes (Koch‘s blue bodies) stages of Theileria spp inside the lymphocyte in a Giemsa‘s stained lymph smear 100X (Oryan et al., 2013).

14

2.2

The life cycle of T. equi Illustration by Massaro Ueti (Wise et al., 2013).

16

2.3

Different forms of B. caballi inside the erythrocytes and the characteristic pair of joint at their posterior ends organisms forming in a Giemsa‘s stained blood smear (S. Kumar & R. Kumar, 2007).

18

2.4

The life cycle of Babesia caballi. Illustration by Massaro Uet (Wise et al., 2013).

21

2.5

The internal structure of the hard tick (40X) by Edward et al. (2009).

48

2.6

Parts of a generalized hard tick; A) Dorsal view of capitulum (mouthpart); B) Ventral view of capitulum; C) Dorsal view of female with body parts; D) Ventral view of male body parts; E) Segment of the leg; F) Ventral view of coxae (Wall & Shearer, 2001).

50

2.7

Diagrammatic dorsal view of the capitulum of seven genera of hard ticks: (A) Ixodes, (B) Hyalomma, (C) Dermacentor, (D) Amblyomma, (E) Boophilus, (F) Rhipicephalus and (G) Haemaphysalis (Wall & Shearer, 2001).

52

2.8

The life cycle of Ixodid ticks (adapted from www.life cycle of ticks family Ixodidae png.)

54

2.9

The type of ixodid ticks depending on the number of hosts attached (P. Jain & A. Jain, 2006).

55

24

3.1

Geographical map of Kelantan state showing the locations of equids sampling stables using Map Window arc GIS 10 program.

108

3.2

The clinical examination card.

112

3.3

A-F) Demonstration of thin and thick blood smears preparation (http://www.slideshare.net/drAjayAgale/05-peripheral-bloodsmear-examination).

115

3.4

Various shapes and stages with measurement of T. equi in stained blood smears with 5% Giemsa, examined under oil immersion lens (100X); A) Pyriform (pair of joint) and single pyriform ; B) Maltese cross and anaplasmoid shape ; C) Rod shape; D) The schizonte stages: microschizontes and macroschizontes (Koch‘s blue bodies) of T. equi within lymphocytes.

133

3.5

Various shapes with measurement of B.caballi in stained blood smears with 5% Giemsa, examined under oil immersion lens (100X); A) Double pear acute and obtuse angle, single pear and round shape; B) Signet ring shape; C) Amoeboid shape.

134

3.6

Gel electrophoresis image showing: lanes M) Exact Mark 1001500bp DNA ladder; Lane 1-9) Conventional PCR technique detected Theileria spp. and Babesia spp. using ‗catch-all‘ primers in approximately band size 496 bp; Lane N) DNA extracted from piroplasms-free horse used as negative control.

141

3.7

Gel electrophoresis image showing: lane M) Exact Mark 1001500bp DNA ladder; Lane P) DNA extracted from clinically infected case used as positive controls for T. equi and B. caballi; Lane 1-6) Multiplex PCR technique detected only T. equi in approximately band size 360 bp; Lane N) DNA extracted from piroplasms-free horse used as negative control.

142

3.8

Gel electrophoresis image showing: lane M) Exact Mark 1001500bp DNA ladder; Lane P) DNA extracted from clinically infected case used as positive controls for T. equi and B. caballi; Lane 1-6) Multiplex PCR technique detected only B. caballi in approximately band size 650 bp; Lane N) DNA extracted from piroplasms-free horse used as negative control.

142

25

3.9

Gel electrophoresis image showing: lane M) Exact Mark 1001500bp DNA ladder; Lane P) DNA extracted from clinically infected case used as positive controls for T. equi and B. caballi; Lane 1-6) Multiplex PCR technique detected both protozoa using specific primers in approximately band size 360 bp and 650 bp respectively; Lane N) DNA extracted from piroplasms-free horse used as negative control.

143

3.10

Geographical map of Kelantan state in Malaysia showing the distribution of T. equi infection. The different marks show infection rate in each stable at different sub-districts using Map Window arc GIS 10 program.

147

3.11

Geographical map of Kelantan state in Malaysia showing the distribution of B. caballi infection. The different marks show infection rate in each stable at different sub-districts using Map Window arc GIS 10 program.

148

3.12

Geographical map of Kelantan state in Malaysia showing the distribution of both protozoa infections. The different marks show infection rate in each stable at different sub-districts using Map Window arc GIS 10 program.

149

3.13

Seroprevalence of T. equi, B. caballi and both protozoa by months; *) Values significantly different (P < 0.05) compared to November month.

159

3.14

Seroprevalence of T. equi, B. caballi and both protozoa by means of monthly temperature (ºC); *) Values significantly different (P