Next Article in Journal
Social Relationships, Age and the Use of Preventive Health Services: Findings from the German Ageing Survey
Next Article in Special Issue
Prolongation of Tick-Borne Encephalitis Cycles in Warmer Climatic Conditions
Previous Article in Journal
Evaluation of the Effectiveness of Mobile App-Based Stress-Management Program: A Randomized Controlled Trial
Previous Article in Special Issue
Single Core Genome Sequencing for Detection of both Borrelia burgdorferi Sensu Lato and Relapsing Fever Borrelia Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 16
Issue 21
10.3390/ijerph16214271
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Serological and Molecular Investigation of Coxiella burnetii in Small Ruminants and Ticks in Punjab, Pakistan

by
Qudrat Ullah
1,2,3,
Hosny El-Adawy
3,4,*,
Tariq Jamil
3,
Huma Jamil
2,
Zafar Iqbal Qureshi
2,
Muhammad Saqib
5,
Shakeeb Ullah
1,
Muhammad Kamal Shah
1,
Alam Zeb Khan
1,
Muhammad Zubair
6,
Iahtasham Khan
7,
Katja Mertens-Scholz
3,
Klaus Henning
3 and
Heinrich Neubauer
3
1
Faculty of Veterinary and Animal Sciences, Gomal University, Dera Ismail Khan 29050, Pakistan
2
Department of Theriogenology, University of Agriculture Faisalabad, Faisalabad 38000, Pakistan
3
Friedrich-Loeffler-Institut, Institute of Bacterial Infections and Zoonoses, 07743 Jena, Germany
4
Faculty of Veterinary Medicine, Kafrelsheik University, Kafr El-Sheik 35516, Egypt
5
Department of Clinical Medicine and Surgery, University of Agriculture Faisalabad, Faisalabad 38000, Pakistan
6
Faculty of Veterinary and Animal Sciences, University of Poonch Rawalakot, Azad Kashmir 12350, Pakistan
7
Section of Epidemiology and Public Health, University of Veterinary and Animal Sciences, Lahore, sub-campus Jhang, Jhang 35200, Pakistan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(21), 4271; https://doi.org/10.3390/ijerph16214271
Submission received: 11 October 2019 / Revised: 29 October 2019 / Accepted: 31 October 2019 / Published: 4 November 2019
(This article belongs to the Special Issue Ticks and Tick Vectored Diseases—Biology to Society)
Versions Notes

Abstract

:
Coxiellosis is a zoonotic disease caused by the obligate intracellular bacterium Coxiella burnetii affecting the productive and reproductive capabilities of animals. This study was conducted to gain insight into the seroprevalence of coxiellosis in small ruminants in seven farms of the Punjab, Pakistan. Potential risk factors were assessed. In total, 1000 serum samples (500 from sheep and 500 from goats) and 163 ticks were collected from the ruminants. All these 163 ticks were merged into 55 pools (29 pools for ticks from sheep and 26 pools for ticks from goat). Serum samples were investigated using an indirect ELISA and PCR. Coxiella burnetii DNA was detected in 29 pooled seropositive samples and 11 pooled ticks by real-time qPCR. Serological analysis revealed a prevalence of 15.6% and 15.0% in sheep and goats, respectively. A significant association was found between seropositivity and different variables like district, lactational status, reproductive status, body condition and reproductive disorders. Univariate analysis showed that detection of C. burnetii DNA in tick pools was significantly associated with the presence of ticks on sheep and goats. However, a non-significant association was found for the prevalence of C. burnetii DNA in serum pools. Hence, C. burnetii infection is prevalent in small ruminants and ticks maintained at livestock farms in Punjab, Pakistan.

1. Introduction

Coxiellosis (Q Fever) is a zoonotic bacterial infection affecting various hosts including humans, ruminants, and pets. The pathogenic agent is to be found everywhere with the exception of New Zealand and French Polynesia [1,2]. Coxiellosis is caused by Coxiella (C.) burnetii, an obligate Gram-negative intracellular bacterium [3,4]. Coxiella burnetii can propagate within the acidic phagolysosomal vacuole of mononuclear phagocytes and shows two distinct life cycles and lipopolysaccharide (LPS) phase variation between phase I and phase II [5]. Coxiella burnetii was detected in all animal species and humans but small ruminants are the most important reservoir and common source of human infection [6]. The largest Dutch Q fever outbreak (2007–2010) with more than 4000 reported human cases was primarily linked to dairy goat farms [7]. The bacteria are mainly transmitted by aerosols. Humans usually get Q fever by breathing in contaminated droplets released by birth products of infected animals and drinking raw milk. Coxiellosis may occasionally be transmitted vertically or sexually but these are not common routes. Ticks may play a role in the transmission of disease in animals but this is questioned for human infection [6,8].
Coxiellosis is usually asymptomatic with sub-clinical presentation in animals and usually not considered a problem for animal health except in ruminants where it causes reproductive problems like abortion, premature delivery, stillbirth and weak offspring [9,10,11]. Abortion is the most important clinical manifestation of coxiellosis in sheep and goats with an incidence of 3% to 80% [12]. In cattle, infertility is the main clinical manifestation [13].
Different techniques can be used for the diagnosis of C. burnetii infection in animals but ELISA for serological investigation and PCR for molecular detection of Coxiella DNA are believed to be the most accurate ways to diagnose [8,14].
In Pakistan, coxiellosis (Q fever) is one of the highly neglected diseases in humans and animal species. From 1955 to 2016, there are only six studies reported the prevalence of C. burnetii infection in humans and animals and most of them are based on a complement fixation test (CFT). According to these studies, the prevalence of coxiellosis ranges from 4.6% to 40% in all livestock species and 10.2% to 26.8% in humans [8,15,16,17,18,19]. In view of the steady increase of the prevalence of this disease worldwide, the present study was designed to investigate the prevalence of coxiellosis and its association with potential associated risk factors in small ruminants maintained at livestock farms of Punjab, Pakistan.

2. Materials and Methods

2.1. Description of Study Area

This study was carried out in small ruminants kept in seven different livestock farms of the Punjab province of Pakistan. Punjab is the largest province, with the highest human and animal population in Pakistan. It is also the second largest province in the country, with an area of 205,344 km2. Geographically, it is located at 31.1704°N and 72.7097°E in the semiarid lowland region. The temperature ranges from −2 to 45 °C but can reach 50 °C (122 °F) in summer and may fall to −10 °C in winter. Mean annual rainfall varies greatly with the highest rainfall in the northern regions [8,18,20]. Agriculture and livestock are the major sources of socio-economic development especially in rural regions of Pakistan. As per the livestock census of Pakistan, Punjab dominates the livestock population with 49%, 65%, 24% and 37% of the cattle, buffaloes, sheep and goats of the country, respectively [21].
Since, no previous studies on C. burnetii have been conducted in these districts so far, it might help us to select a particular region/area to be studied. We choose districts, which reflect the major livestock population in the province with an increased annual incidence of livestock-associated zoonotic diseases (Directorates of Animal and Human Health, Punjab).

2.2. Estimation of Sample Size

This study was conducted during January to June 2016 in seven governmental livestock farms located in the Khanewal, Khushab, Rajanpur, Bhawalpur, Layyah, Okara and Bhakkar districts of Punjab province, Pakistan. Survey Toolbox software (Ausvet, The Australian Biosecurity Cooperative Research Centre for Emerging Infectious Disease, Australia) was used for the random selection of farms and animals [22]. The prevalence of Q fever in small ruminants in the study areas was unknown. For the calculation of sample size, a 95% confidence interval with an expected prevalence of 50% and desired absolute precision of 5% was applied [23]. Three hundred and eighty-four serum samples were required from all farms. In total, 1000 samples (500 from sheep and 500 from goats) were collected. During sampling, data on the individual animal, i.e., species, locality of farm (districts), age, sex, parity, lactational status, reproductive status (pregnant or non-pregnant), body condition, and reproductive disorders, and general farm management were recorded using a structured questionnaire.

2.3. Sample Collection

Approximately 10 mL of blood was collected from the jugular vein of each animal using disposable needles and evacuated blood collection tubes (Improvacuter, Shanghai International Hamburg Holding, GmbH, Germany). Each tube was properly labeled for its identification and information about its origin. Blood collection tubes were placed vertically in a cool box packed with gel freezing pads and shipped to the laboratory for further processing. At the laboratory, centrifugation was performed at 4500 rpm for 10 min for proper separation of the serum from whole blood. The serum was preserved in disposable screw-capped cryovials (Cryo.STM; Greiner Bio-one, GmbH Frickenhausen, Germany) and stored at −20 °C until use for further investigation.

2.4. Collection and Handling of Ticks

During the serum sampling (January to June, 2016), 163 ticks (85 from sheep and 78 from goats) were collected from the same seven governmental livestock farms where serum samples were collected. Out of 163 ticks, 85 ticks were collected from sheep, while 78 were collected from goats. The ticks were placed in properly labeled bottles containing 70% ethyl alcohol as a preservative at the University of Agriculture Faisalabad. Ticks were shipped to Friedrich-Loeffler-Institut (FLI), Jena, Germany, following international regulations for transportation. Ticks were merged into 55 pools (29 pools from sheep and 26 pools from goats). Each pool was composed of 3 ticks (with only pool no. 29 containing 1 tick). These pools were then stored at −20 °C until use for the detection of Coxiella DNA.

2.5. Ethics Statement

Blood samples were collected from small ruminants as per the guidelines of the International Animal Care and Use Committee (IACUC) and after obtaining the consent of the Secretary, Livestock and Dairy Development department, Punjab, Pakistan (wide letter No.SO (I&C)/L7DD/2-6/2016). The samples were processed based on the approval of the Ethical Research Board at the University of Agriculture Faisalabad, Pakistan (wide letter No. ORIC/3253, dated November 16, 2013).

2.6. Diagnostic Tests

2.6.1. Indirect ELISA

All serum samples were analyzed by Indirect ELISA (IDEXX Q-Fever (Coxiella burnetii) Antibody Test Kit, Germany). The plates were read at 450 nm with an ELISA reader (Anthos 2020, Wals, Austria) according to the manufacturer’s recommendation. To calculate the results, the optical densities (ODs) of the sera were considered to be ELISA-positive if they had an OD value of 40% or more, suspected if the value was between 30% and 40%, and were considered negative if the OD value was <30%. Confirmed positive results by indirect ELISA were considered as a final criteria for the calculation of seroprevalence.

2.6.2. Real-Time qPCR

After serological analysis, 29 pools of 153 seropositive and 81 suspected serum samples and 10 pools of seronegative samples were analyzed using TaqMan-based real-time PCR assay targeting the multi copy IS1111 element as well as the single copy icd gene of C. burnetii [24]. Each serum pool was comprised of six serum samples. All molecular procedures were carried out in specially designed cabins to prevent any kind of contamination. DNA was extracted using a high pure PCR template preparation kit (Roche kit, version 20, Roche Diagnostics GmbH, Germany).
The surface of ticks was washed by serial passages in 10% and/or 70% ethyl alcohol and rinsed in sterile water. DNA was extracted from 55 tick pools as described by using the protocol for the purification of genomic DNA from insects using high pure PCR template preparation kit (Roche, Roche Diagnostics GmbH, Germany) and investigated using real-time qPCR targeting multiple copy IS1111 transposase gene as well as the icd gene of C. burnetii. Real-time qPCR performed in duplicate.
In each qPCR run, plasmid standards (for IS1111 from 105 to 101, while for icd 106 to 101) were used as a positive control and DNase and RNase-free water was used as a negative control. Results of qPCR were counted as positive at a cut off ≤35. A standard curve was generated using decimal dilutions of the positive control. The PCR primer sequences and conditions were applied as described by Chakrabartty et al. 2016 [25].

2.7. Statistical Analysis

Minitab17 software (Minitab, LLC, State College, PA, USA) was used for calculation of prevalence percentage and respective 95% confidence intervals. ODD ratios, respective 95% confidence interval and Chi-square values were calculated using IBM SPSS Statistics 13.0 for Windows® (IBM Corporation, Armonk, NY, USA). Chi-square test was applied for the calculation of the significance of association (p < 0.05) between the seroprevalence of coxiellosis and various farm-related variables.

3. Results

3.1. Indirect ELISA Results

The results of iELISA showed that 153 serum samples were seropositive, and 81 serum samples were considered as suspected. Univariate analysis was conducted to investigate the relationship of seven livestock farms and associated variables with the prevalence of C. burnetii infection (Table 1 and Table 2). The overall prevalence of coxiellosis in seven districts of the Punjab province was found to be 15.3%. The highest prevalence (26.2%) was recorded in the district Layyah, while the lowest prevalence (5.8%) was reported in animals of the Rajanpur and Khanewal district. Statistical analysis showed a significant (p = 0.000) difference for the prevalence in various districts. There is no significant difference between seroprevalence in sheep (15.6%) and goats (15.0%). Male animals showed higher seroprevalence (18.7%) than females (14.9%) (p = 0.302). However, seropositivity for C. burnetii infection was non-significantly associated with age (p = 0.085) and parity (p = 0.838). A higher significant prevalence (20.0%) was recorded in non-lactating animals than in lactating animals (8.5%) (p = 0.000). Regarding the reproductive status, the anti-Coxiella antibodies were detected in 17.3% and 10.6% of non-pregnant and pregnant animals, respectively (p = 0.002).
The prevalence of C. burnetii infection in animals with a history of various reproductive problems was significantly different (p = 0.000). The seroprevalence in animals with a previous history of abortion, premature delivery, stillbirth and repeated breeding was 51.6%, 32.4%, 31.8% and 25%, respectively.
As far as body condition was concerned, higher prevalence, (44.9%) was recorded in small ruminants with weak body condition, while lower prevalence (37.7%) was found in ewes with the delivery of a weak lamb. Statistical analysis revealed that the prevalence of the disease was significantly different (p = 0.000) between weak animals and animals with a weak lamb or kid delivery.

3.2. Real-Time qPCR Results

3.2.1. Serum Pools

Out of 29 serum pools tested with qPCR, C. burnetii DNA was not detected in seropositive or seronegative pools. Coxiella burnetii DNA was amplified in five (17.2%) serologically suspected serum pools. Statistical analysis revealed that the prevalence of C. burnetii DNA in serum pools was not significantly (p = 0.333) different between sheep (13.3%) and goats (21.4%) (Table 3).

3.2.2. Ticks Pools

Coxiella burnetii DNA was detected in 20% of tick pools. Univariate analysis of data revealed a significant difference in the prevalence of C. burnetii DNA in tick pools collected from sheep (31.0%) and goats (7.7%) (p = 0.031, df = 1, χ2 = 4.668) (Table 4).

4. Discussion

Coxiellosis is a zoonotic disease caused by C. burnetii. The diagnosis of C. burnetii infection in animals is of great importance not only to identify the infected flocks but also to determine the risk of disease transmission to humans [8,26]. The epidemiological information about the geographic distribution of C. burnetii infection is very limited in humans and livestock populations in Pakistan. To the best of our knowledge, this is the first epidemiological survey investigating the prevalence of coxiellosis in small ruminants and ticks maintained at different governmental livestock farms in Punjab, Pakistan. In the current study, the prevalence of C. burnetii infection in small ruminants has been investigated using ELISA and real-time qPCR. While in ticks, real-time PCR assay was used for the detection of C. burnetii DNA in pooled samples.
In this study, the prevalence of coxiellosis varied greatly in the seven districts of the Punjab province. The prevalence ranged from 26.2% (95% CI: 21.1–31.7) in district Layyah to 5.8% (95% CI: 3.2–9.5) in districts Rajanpur and Khanewal. Univariate analysis proved a significant (p = 0.000) association of prevalence and district. Ezatkhah et al. (2015) recorded a seroprevalence of 26.4% in small ruminants of five counties of the southeast region of Iran, which ranged from 5% in Sarbaz to 39.2% in Iranshahr [27]. The higher prevalence of the disease might be due to the prevailing climatic and weather conditions in the region studied [28]. This variation in the prevalence of infection within flock and different geographical areas might be associated with the farm hygienic measures, routine management practices and environmental factors such as vegetation, soil moisture and the presence of infected animals in the surroundings [7,29]. These management and environmental factors might be responsible for the higher seroprevalence of coxiellosis in different districts. A significant high prevalence of 69.4% and 75.0% was reported previously by Zahid et al. (2016) in the Layyah and Muzaffargarh districts of Punjab province, respectively [19]. On the other hand, a very low seroprevalence (5.8%) was recorded in the districts Rajanpur and Khanewal, which might be attributed to better hygienic practices.
Results of this study revealed an overall seroprevalence in small ruminants of 15.3% (95% CI: 12.3–18.7), which varied from 15.6% (95% CI: 12.07–18.33) to 15.0% (95% CI: 12.62–18.98) between sheep and goats (p = 0.792). Rizzo et al., 2016, reported a similar flock level seroprevalence (14.5%) in sheep and goats [7]. However, a few previous studies reported higher seropositivity in sheep [30,31]. Controversially, other studies found higher prevalence in goats [7,32]. Variations in intrinsic susceptibility to this pathogen among small ruminant species have not been described in detail [32]. Thus, the species-wise variations in the prevalence of disease may be more likely attributed to different farm management practices and variation in flock density in sheep and goat farms [19,33]. This survey revealed that small ruminants are important reservoirs of C. burnetii infection in the study area.
Male animals showed a higher seroprevalence (18.7%) than females (14.9%) (p = 0.302). However, seropositivity for C. burnetii infection was non-significantly associated with age (p = 0.085) and parity (p = 0.838).
In the current study, the prevalence of C. burnetii infection varied greatly between non-lactating (20.0%) and lactating animals (8.5%) (p = 0.000). While in Kenya, the seroprevalence of coxiellosis was higher in lactating animals as compared to non-lactating animals, although the difference was not statistically significant [28]. Paul et al. (2012) reported a significantly (p < 0.005) higher prevalence of coxiellosis in dairy breeds as compared to the beef breeds of cattle [29]. The target sites for proliferation of C. burnetii are the placenta and mammary glands of animals [34,35]. In ruminants, immediately after inoculation, the organism reaches the predilection site through the blood stream and resides in the supramammary lymph nodes, mammary glands and placenta of pregnant animals, where it replicates and may contribute to higher prevalence in lactating animals [14].
In this study, the non-pregnant small ruminants showed higher seroprevalence (17.3%) than pregnant (10.6%) (p = 0.002), which is in agreement with the study conducted by Abushahba et al. (2017), which observed a non-significant difference regarding reproductive status in small ruminants and the prevalence of Q fever in Egypt [36]. The results of this study were in contrast to a study conducted previously in the Netherlands that reported a higher seroprevalence of Coxiella in pregnant and periparturient small ruminants when compared with non-pregnant ruminants and those in the early gestation period [37].
In this study, the prevalence of coxiellosis was higher (11.1%) in young non-pregnant animals aged from 1 to 2.5 years than in pregnant animals (6.8%). In middle-aged small ruminants (2.5 to 4 years), the prevalence was 19.7% and 16% in non-pregnant and pregnant, respectively. The older non-pregnant animals (≥4 years) showed a higher prevalence (17.5%) than pregnant animals of the same age (9.6%).
Interestingly, the prevalence of anti-Coxiella antibodies was 20% in young females (≤1 year).
Results of the current study revealed that seroprevalence of coxiellosis was significantly (p = 0.000) different in animals with a previous history of reproductive problems including abortion (51.6%), premature delivery (32.4%), stillbirth (31.8%) and repeated breeding (25%). These findings are in agreement with the results of García-Pérez et al. (2009), who observed a significantly higher prevalence of coxiellosis in animals with a history of reproductive problems [38]. Controversially, Ruiz-Fons et al. (2010) and Asadi et al. (2013) did not record any association between the prevalence of coxiellosis and a history of reproductive disorders [39,40]. Antibodies against C. burnetii found in small ruminants with a history of reproductive problems do not necessarily mean that C. burnetii is the only cause of reproductive loss in these animals. Infection with other pathogens such as Brucella melitensis and Toxoplasma gondii may be possible confounders. Furthermore, malnutrition in pregnant females might be a reason for abortion. During the periparturient period, the massive multiplication of these pathogens occurs within trophoblast cells, causing necrotic suppurative placentitis, which ultimately leads to pregnancy failure in the form of abortion, stillbirth, premature delivery and birth of weak offspring [41,42].
A significant (p = 0.000, df = 3, χ2 = 124.868) difference in prevalence was found between C. burnetii seropositivity and animals with poor body conditions. Highest seropositivity was recorded in weak animals (45.7%) followed by the does that gave birth to weak kids (44.9%) and ewes that delivered weak lambs (37.7%). Agerholm (2013) also reported that C. burnetii-infected ewes gave birth to noticeably small and underweight lambs [43]. Similarly, Saegerman et al. (2013) found an association between C. burnetii seropositivity and delivery of weak calves (OR = 2.14 with 95% CI: 1.05–4.39) [44].
Infected goats with Q fever after abortion showed signs of endometritis. Full-term kids were weak, with low body weight and high mortality. Many of the apparently healthy kids infected with C. burnetii were suffering from digestive and respiratory tract problems [45].
The higher seropositivity in females with weak lamb and kid deliveries might be due to the massive bacterial multiplication that occurs during the periparturient period in pregnant animals [42].
Real-time PCR assay targeting a single copy of the icd gene was used for the detection of C. burnetii DNA [24]. Results of this study revealed that C. burnetii DNA was not detected in the seropositive samples. However, the DNA was detected in five (17.2%) pools of suspected serum samples.
The presence of C. burnetii DNA in suspected serum pools indicated that infection was probably in the incubation period or in the early acute phase, as C. burnetii DNA cannot be detected in serum after two weeks of the onset of clinical infection [46,47,48,49].
Schneeberger et al. (2010) also recorded PCR positive results in 98% seronegative samples collected from acute Q fever patients [48]. Vincent et al. (2015) found that early serum samples from acute Q fever patients are infectious and an important source of viable C. burnetii [50]. Apart from serological investigation, C. burnetii-specific PCR can be effectively used for the genomic detection of coxiellosis in the early acute phase of Q fever, but contradictory sensitivities have also been reported [48,51]. Hence, qPCR can also be an additional tool to achieve a better understanding of the disease prevalence in the veterinary medicine practice.
Coxiella burnetii DNA was detected in 31.0% (95% CI: 15.3–50.8) and 7.7% (95% CI: 0.9–25.1) of ticks collected from sheep and goats, respectively. Univariate analysis revealed a significant (p = 0.031, df = 1, χ2 = 4.668) difference in the prevalence of C. burnetii DNA in tick pools collected from sheep and goats.
Knobel et al. (2013) showed that C. burnetii DNA can be detected in 2.5% and 20.0% of ticks collected from cattle and dogs, respectively [52]. In another study conducted in Ethiopia, a significant number of C. burnetii DNA (25.0%) was found in samples from questing ticks collected from domestic animals [53]. In contrast, Sprong et al. (2012) reported only 0.2% C. burnetii DNA in questing Ixodes ricinus ticks and ticks in wildlife, pets and domestic animals by using a multiplex q-PCR [54]. Although, some pioneering studies have focused on C. burnetii in ticks, the role of ticks in the epidemiology of Q fever remains uncertain. Even the highly virulent Nine Mile reference strain of C. burnetii was isolated from a guinea pig upon Dermacentor andersoni ticks had fed [55]. The occasional reports of an unexpectedly high prevalence of C. burnetii DNA in ticks may reflect their role as a vector for the transmission of Q fever [55].
The marker IS1111 transposable element was routinely targeted during the epidemiological investigation of C. burnetii prevalence in ticks [55], while C. burnetii detection assays based only on IS1111 may lead to misidentification with Coxiella-like endosymbionts [56]. Reeves et al. 2006 were able to amplify a 612 bp icd fragment, displaying 93.0% homology with C. burnetii, from a Coxiella-like bacterium that infects ticks [57]. In contrast, Elsa et al., 2015, proved that icd was not amplified from the endosymbiont of Ornithodoros capensis ticks and highlights the fact that the amplification of a specific genetic marker strongly depends on the PCR method and the primer sequences used [56].
In this study, we used a combination of primers targeting different markers, IS1111 and icd, to guarantee the specificity of C. burnetii detection.
Standardizing a methodology for the detection of C. burnetii DNA across laboratories is essential to make sure that no cross-reaction exists not only with other abortive agents but also with Coxiella-like organisms.

5. Conclusions

The findings of the present study for the detection of C. burnetii infection in small ruminants in different geographical regions of Pakistan indicates that coxiellosis is endemic in this country. Since there is very limited awareness of the disease among veterinarian and farm workers of the studied livestock farms, an awareness-raising campaign should be launched to educate livestock farmers and professionals about proper preventive and control measures for coxiellosis.
This survey highlights the importance of clinical herd health management to avoid economical losses in small ruminants.

Author Contributions

Data curation, Q.U. and H.J.; Formal analysis, Q.U.; Investigation, Q.U., H.J., Z.I.Q., M.S., S.U., M.K.S., A.Z.K., M.Z., T.J., K.M.-S. and K.H.; Methodology, Q.U. and H.J.; Supervision, H.E.-A. and M.S.; Writing—original draft, Q.U., H.E.-A., H.J., Z.I.Q., M.S., S.U., M.K.S., A.Z.K., M.Z. and I.K.; Writing—review and editing, H.E.-A., T.J., K.M.-S. and H.N.

Funding

This research received no external funding.

Acknowledgments

The authors thank Muhammad Hammad Hussain at Animal Research Center, Ministry of Agriculture and Fisheries, Muscat, Oman, for helping in sample collection and data analysis. The authors are thankful to the Higher Education Commission of Pakistan for supporting this study (Scholar PIN No. 112-34911-2AV1-119). The authors thank Federal Foreign Office, Germany, “German Biosecurity Programme”, in collaboration with Friedrich-Loeffler-Institut, Institute of Bacterial Infections and Zoonoses, Jena, Germany, for their support in this work. The authors are thankful to the Livestock and Dairy Development Department of Punjab and the local veterinarians and farmers for facilitating in serum sample and tick collection from small ruminants maintained at livestock farms.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Million, M.; Raoult, D. Recent advances in the study of Q fever epidemiology, diagnosis and management. J. Infect. 2015, 71, S2–S9. [Google Scholar] [CrossRef] [PubMed]
  2. Musso, D.; Broult, J.; Parola, P.; Raoult, D.; Fournier, P.E. Absence of antibodies to Rickettsia spp., Bartonella spp., Ehrlichia spp. and Coxiella burnetii in Tahiti, French Polynesia. BMC Infect. Dis. 2014, 14, 255–258. [Google Scholar] [CrossRef] [PubMed]
  3. Bielawska-Drozd, A.; Cieslik, P.; Mirski, T.; Gawel, J.; Michalski, A.; Niemcewicz, M.; Bartoszcze, M.; Zakowska, D.; Lasocki, K.; Knap, J.; et al. Prevalence of Coxiella burnetii in environmental samples collected from cattle farms in Eastern and Central Poland (2011–2012). Vet. Microbiol. 2014, 174, 600–606. [Google Scholar] [CrossRef] [PubMed]
  4. Van Leuken, J.P.G.; Swart, A.N.; Brandsma, J.; Terink, W.; Van de Kassteele, J.; Droogers, P.; Sauter, F.; Havelaar, A.H.; Van der Hoek, W. Human Q fever incidence is associated to spatiotemporal environmental conditions. One Health 2016, 2, 77–87. [Google Scholar] [CrossRef] [Green Version]
  5. Van Schaik, E.J.; Chen, C.; Mertens, K.; Weber, M.M.; Samuel, J.E. Molecular pathogenesis of the obligate intracellular bacterium Coxiella burnetii. Nat. Rev. Microbiol. 2013, 11, 561–573. [Google Scholar] [CrossRef]
  6. OIE. Manual of the Diagnostic Tests and Vaccines for Terrestrial Animals; Office International des Epizootics: Paris, France, 2016; pp. 1–15. [Google Scholar]
  7. Rizzo, F.; Vitale, N.; Ballardini, M.; Borromeo, V.; Luzzago, C.; Chiavacci, L.; Mandola, M.L. Q fever seroprevalence and risk factors in sheep and goats in northwest Italy. Prev. Vet. Med. 2016, 130, 10–17. [Google Scholar] [CrossRef]
  8. Ullah, Q.; Jamil, H.; Iqbal Qureshi, Z.; Saqib, M.; Neubauer, H. Sero-epidemiology of Q Fever (Coxiellosis) in small ruminants kept at government livestock farms of Punjab, Pakistan. Pak. J. Zool. 2018, 51, 135–140. [Google Scholar] [CrossRef]
  9. EFSA. EFSA panel on animal health and welfare (AHAW); Scientific opinion on Q Fever. EFSA J. 2010, 8, 1–114. [Google Scholar] [CrossRef]
  10. Lang, G.H. Q fever: An emerging public health concern in Canada. Can. J. Vet. Res. 1989, 53, 1–6. [Google Scholar]
  11. Carbonero, A.; Guzman, L.T.; Montano, K.; Torralbo, A.; Arenas-Montes, A.; Saa, L.R. Coxiella burnetii seroprevalence and associated risk factors in dairy and mixed cattle farms from Ecuador. Prev. Vet. Med. 2015, 118, 427–435. [Google Scholar] [CrossRef]
  12. Arricau-Bouvery, N.; Rodolakis, A. Is Q fever an emerging or re-emerging zoonosis? Vet. Res. 2005, 36, 327–349. [Google Scholar] [CrossRef] [PubMed]
  13. Roest, H.J. Coxiella Burnetii in Pregnant Goats; Central Veterinary Institute of Wageningen UR: Lelystad, The Netherlands, 2013. [Google Scholar]
  14. Hadush, A.; Kandi, V.; Pal, M. Epidemiology and public health implications of Q fever. Perspect. Med. Res. 2016, 4, 42–46. [Google Scholar]
  15. Ahmed, I.P. A serological investigation of Q fever in Pakistan. J. Pak. Med. Assoc. 1987, 37, 126–129. [Google Scholar] [PubMed]
  16. Ayaz, M.; Bari, A.; Humayun, A. Coxiellosis in man and animals in northern parts of Pakistan. Proc. Pak. Congr. Zool. 1993, 13, 425–431. [Google Scholar]
  17. Kaplan, M.M.; Bertagna, P. The geographical distribution of Q fever. Bull. World Health Organ. 1955, 13, 829–860. [Google Scholar] [PubMed]
  18. Shabbir, M.Z.; Akram, S.; Hassan, Z.U.; Hanif, K.; Rabbani, M.; Muhammad, J.; Chaudhary, M.H.; Abbas, T.; Ghori, M.T.; Rashid, H.; et al. Evidence of Coxiella burnetii in Punjab province, Pakistan. Acta Trop. 2016, 163, 61–69. [Google Scholar] [CrossRef]
  19. Zahid, M.U.; Hussain, M.H.; Saqib, M.; Neubauer, H.; Abbas, G.; Khan, I.; Mansoor, M.K.; Asi, M.N.; Ahmad, T.; Muhammad, G. Seroprevalence of Q Fever (Coxiellosis) in small ruminants of two districts in Punjab, Pakistan. Vector Borne Zoonotic Dis. 2016, 16, 449–454. [Google Scholar] [CrossRef]
  20. Khattak, M.S.; Ali, S. Assessment of temperature and rainfall trends in Punjab province of Pakistan for the period 1961–2014. J. Himal. Earth Sci. 2015, 48, 42–61. [Google Scholar]
  21. Pakistan Livestock Census 2006. Available online: http://www.pbs.gov.pk/sites/default/files/aco/publications/pakistan-livestock-cencus2006/All%20Pakistan%20Report%20%28Introduction%29.pdf (accessed on 1 November 2019).
  22. Cameron, A. Survey Toolbox for Livestock Diseases; Australian Center for Internationl Agricultural Research: Bruce, Australia, 1999. [Google Scholar]
  23. Thrusfield, M. Veterinary Epidemiology, 3rd ed.; Blackwell Sci Ltd.: Hoboken, NJ, USA, 2007. [Google Scholar]
  24. Klee, S.R.; Tyczka, J.; Ellerbrok, H.; Franz, T.; Linke, S.; Baljer, G.; Appel, B. Highly sensitive real-time PCR for specific detection and quantification of Coxiella burnetii. BMC Microbiol. 2006, 6, 2. [Google Scholar] [CrossRef]
  25. Chakrabartty, A.; Bhattacharjee, P.K.; Sarker, R.R.; Rahman, A.K.M.A.; Henning, K.; Neubauer, H.; Rahman, M.S. Prevalence of Coxiella burnetii infection in cattle, black bengal goats and ticks in Bangladesh. Bangladesh J. Vet. Med. 2016, 14, 65–68. [Google Scholar] [CrossRef]
  26. OIE. Q Fever OIE Terrestrial Manual 2015; Office International des Epizootics: Paris, France, 2015; pp. 1–15. [Google Scholar]
  27. Ezatkhah, M.; Alimolaei, M.; Khalili, M.; Sharifi, H. Seroepidemiological study of Q fever in small ruminants from Southeast Iran. J. Infect. Public Health 2015, 8, 170–176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Muema, J.; Thumbi, S.M.; Obonyo, M.; Wanyoike, S.; Nanyingi, M.; Osoro, E.; Bitek, A.; Karanja, S. Seroprevalence and factors associated with Coxiella burnetii infection in small ruminants in Baringo County, Kenya. Zoonoses Public Health 2017, 64, e31–e43. [Google Scholar] [CrossRef] [PubMed]
  29. Paul, S.; Agger, J.F.; Markussen, B.; Christoffersen, A.B.; Agerholm, J.S. Factors associated with Coxiella burnetii antibody positivity in Danish dairy cows. Prev. Vet. Med. 2012, 107, 57–64. [Google Scholar] [CrossRef] [PubMed]
  30. Anastacio, S.; Tavares, N.; Carolino, N.; Sidi-Boumedine, K.; da Silva, G.J. Serological evidence of exposure to Coxiella burnetii in sheep and goats in central Portugal. Vet. Microbiol. 2013, 167, 500–505. [Google Scholar] [CrossRef] [PubMed]
  31. Kshash, Q.H. Prevalence of Q-fever in small ruminant in Al-Qassim city. Basrah J. Vet. Res. 2012, 11, 342–348. [Google Scholar] [CrossRef]
  32. Klaasen, M.; Roest, H.; van der Hoek, W.; Goossens, B.; Secka, A.; Stegeman, A. Coxiella burnetii seroprevalence in small ruminants in the Gambia. PLoS ONE 2014, 9, e85424. [Google Scholar] [CrossRef]
  33. Khaled, H.; Sidi-Boumedine, K.; Merdja, S.; Dufour, P.; Dahmani, A.; Thiery, R.; Rousset, E.; Bouyoucef, A. Serological and molecular evidence of Q fever among small ruminant flocks in Algeria. Comp. Immunol. Microbiol. Infect. Dis. 2016, 47, 19–25. [Google Scholar] [CrossRef]
  34. Jung, B.Y.; Seo, M.G.; Lee, S.H.; Byun, J.W.; Oem, J.K.; Kwak, D. Molecular and serologic detection of Coxiella burnetii in native Korean goats (Capra hircus coreanae). Vet. Microbiol. 2014, 173, 152–155. [Google Scholar] [CrossRef]
  35. Shapiro, A.J.; Bosward, K.L.; Heller, J.; Norris, J.M. Seroprevalence of Coxiella burnetii in domesticated and feral cats in Eastern Australia. Vet. Microbiol. 2015, 177, 154–161. [Google Scholar] [CrossRef]
  36. Abushahba, M.F.N.; Abdelbaset, A.E.; Rawy, M.S.; Ahmed, S.O. Cross-sectional study for determining the prevalence of Q fever in small ruminants and humans at El Minya Governorate, Egypt. BMC Res. Notes 2017, 10, 538. [Google Scholar] [CrossRef]
  37. Van den Brom, R.; Moll, L.; van Schaik, G.; Vellema, P. Demography of Q fever seroprevalence in sheep and goats in The Netherlands in 2008. Prev. Vet. Med. 2013, 109, 76–82. [Google Scholar] [CrossRef] [PubMed]
  38. Garcia-Perez, A.L.; Astobiza, I.; Barandika, J.F.; Atxaerandio, R.; Hurtado, A.; Juste, R.A. Short communication: Investigation of Coxiella burnetii occurrence in dairy sheep flocks by bulk-tank milk analysis and antibody level determination. J. Dairy Sci. 2009, 92, 1581–1584. [Google Scholar] [CrossRef] [PubMed]
  39. Asadi, J.; Kafi, M.; Khalili, M. Seroprevalence of Q fever in sheep and goat flocks with a history of abortion in Iran between 2011 and 2012. Vet. Ital. 2013, 49, 163–168. [Google Scholar] [PubMed]
  40. Ruiz-Fons, F.; Astobiza, I.; Barandika, J.F.; Hurtado, A.; Atxaerandio, R.; Juste, R.A.; García-Pérez, A.L. Seroepidemiological study of Q fever in domestic ruminants in semi-extensive grazing systems. BMC Vet. Res. 2010, 6, 1–3. [Google Scholar] [CrossRef]
  41. Arserim, N.B.; Yesilmen, S.; Tel, O.Y.; Ozekinci, T.; Keskin, O.; Pulat, H.; Vural, A. Seroprevalance of coxiellosis in cows, sheep, goats and humans in Diyarbakir region of Turkey. Afr. J. Microbiol. Res. 2011, 5, 2041–2043. [Google Scholar]
  42. Sanchez, J.; Souriau, A.; Buendia, A.J.; Arricau-Bouvery, N.; Martinez, C.M.; Salinas, J.; Rodolakis, A.; Navarro, J.A. Experimental Coxiella burnetii infection in pregnant goats: A histopathological and immunohistochemical study. J. Comp. Pathol. 2006, 135, 108–115. [Google Scholar] [CrossRef]
  43. Agerholm, J.S. Coxiella burnetii associated reproductive disorders in domestic animals—A critical review. Acta Vet. Scand. 2013, 55, 13–25. [Google Scholar] [CrossRef]
  44. Saegerman, C.; Speybroeck, N.; Dal Pozzo, F.; Czaplicki, G. Clinical indicators of exposure to Coxiella burnetii in dairy herds. Transbound. Emerg. Dis. 2015, 62, 46–54. [Google Scholar] [CrossRef]
  45. Ganter, M. Zoonotic risks from small ruminants. Vet. Microbiol. 2015, 181, 53–65. [Google Scholar] [CrossRef]
  46. Howe, G.B.; Loveless, B.M.; Norwood, D.; Craw, P.; Waag, D.; England, M.; Lowe, J.R.; Courtney, B.C.; Pitt, M.L.; Kulesh, D.A. Real-time PCR for the early detection and quantification of Coxiella burnetii as an alternative to the murine bioassay. Mol. Cell. Probes 2009, 23, 127–131. [Google Scholar] [CrossRef]
  47. Niemczuk, K.; Szymanska-Czerwinska, M.; Smietanka, K.; Bocian, L. Comparison of diagnostic potential of serological, molecular and cell culture methods for detection of Q fever in ruminants. Vet. Microbiol. 2014, 171, 147–152. [Google Scholar] [CrossRef] [PubMed]
  48. Schneeberger, P.M.; Hermans, M.H.; van Hannen, E.J.; Schellekens, J.J.; Leenders, A.C.; Wever, P.C. Real-time PCR with serum samples is indispensable for early diagnosis of acute Q fever. Clin. Vaccine Immunol. 2010, 17, 286–290. [Google Scholar] [CrossRef] [PubMed]
  49. Wielders, C.C.H.; Wijnbergen, P.C.A.; Renders, N.H.M.; Schellekens, J.J.A.; Schneeberger, P.M.; Wever, P.C.; Hermans, M.H.A. High Coxiella burnetii DNA load in serum during acute Q fever is associated with progression to a serologic profile indicative of chronic Q fever. J. Clin. Microbiol. 2013, 51, 3192–3198. [Google Scholar] [CrossRef] [PubMed]
  50. Vincent, G.A.; Graves, S.R.; Robson, J.M.; Nguyen, C.; Hussain-Yusuf, H.; Islam, A.; Fenwick, S.G.; Stenos, J. Isolation of Coxiella burnetii from serum of patients with acute Q fever. J. Microbiol. Methods 2015, 119, 74–78. [Google Scholar] [CrossRef] [PubMed]
  51. Fournier, P.E.; Raoult, D. Comparison of PCR and serology assays for early diagnosis of acute Q fever. J. Clin. Microbiol. 2003, 41, 5094–5098. [Google Scholar] [CrossRef]
  52. Knobel, D.L.; Maina, A.N.; Cutler, S.J.; Ogola, E.; Feikin, D.R.; Junghae, M.; Halliday, J.E.B.; Richards, A.L.; Breiman, R.F.; Cleaveland, S.; et al. Coxiella burnetii in humans, domestic ruminants, and ticks in rural Western Kenya. Am. J. Trop. Med. Hyg. 2013, 88, 513–518. [Google Scholar] [CrossRef]
  53. Kumsa, B.; Socolovschi, C.; Almeras, L.; Raoult, D.; Parola, P. Occurrence and genotyping of Coxiella burnetii in ixodid ticks in Oromia, Ethiopia. Am. J. Trop. Med. Hyg. 2015, 93, 1074–1081. [Google Scholar] [CrossRef]
  54. Sprong, H.; Tijsse-Klasen, E.; Langelaar, M.; De Bruin, A.; Fonville, M.; Gassner, F.; Takken, W.; Van Wieren, S.; Nijhof, A.; Jongejan, F.; et al. Prevalence of Coxiella burnetii in ticks after a large outbreak of Q fever. Zoonoses Public Health 2012, 59, 69–75. [Google Scholar] [CrossRef]
  55. Duron, O.; Sidi-Boumedine, K.; Rousset, E.; Moutailler, S.; Jourdain, E. The importance of ticks in Q Fever transmission: What has (and has not) been demonstrated? Trends Parasitol. 2015, 31, 536–552. [Google Scholar] [CrossRef]
  56. Elsa, J.; Duron, O.; Séverine, B.; González-Acuña, D.; Sidi-Boumedine, K. Molecular methods routinely used to detect Coxiella burnetii in ticks cross-react with Coxiella-like bacteria. Infect. Ecol. Epidemiol. 2015, 5, 29230. [Google Scholar] [CrossRef]
  57. Reeves, W.K.; Loftis, A.D.; Sanders, F.; Spinks, M.D.; Wills, W.; Denison, A.M.; Dasch, G.A. Borrelia, Coxiella, and Rickettsia in Carios capensis (Acari: Argasidae) from a brown pelican (Pelecanus occidentalis) rookery in South Carolina, USA. Exp. Appl. Acarol. 2006, 39, 321–329. [Google Scholar] [CrossRef] [PubMed]
Table
Table 1. Seroprevalence of coxiellosis in small ruminants in seven districts of Punjab, Pakistan.
Table 1. Seroprevalence of coxiellosis in small ruminants in seven districts of Punjab, Pakistan.
LocationPositive/TestedPrevalence, %95% CIOR #95% CI *Statistics
Layyah74/28326.221.1–31.75.791.75–19.11χ2 = 49.689
df = 6
p = 0.000
Bhakkar37/20817.812.8–23.73.541.04–11.9
Khushab6/4513.35.1–26.82.520.59–10.69
Okara15/13211.46.5–18.02.10.58–7.56
Bhawalpur4/3710.83.0–25.41.980.42–9.43
Khanewal14/2435.83.2–9.510.28–3.61
Rajanpur3/525.81.2–15.91-
Overall153/100015.312.3–18.7--
Individual prevalence was significantly different in the farms investigated, χ2 = 49.689, df = 6, p = 0.000. # OR (odds ratio); * CI (confidence interval).
Table
Table 2. Analysis for the serological results of coxiellosis in small ruminants in seven farms in the Punjab, Pakistan, in relation to potential risk factors.
Table 2. Analysis for the serological results of coxiellosis in small ruminants in seven farms in the Punjab, Pakistan, in relation to potential risk factors.
VariableCategoryPos./TestedPrev. %95% CIOR95% CIp-Value
SpeciesSheep78/50015.612.5–19.11.0470.74–1.47χ2 = 0.069
df = 1
p = 0.792
Goat75/50015.012.0–18.41-
AgeUp to 122/13316.510.0–23.11.0820.63–1.85χ2 = 6.603
df = 3
p = 0.085
>1–2.520/18610.86.1–15.61-
>2.5–447/24119.514.9–25.11.4060.93–2.12
>464/44014.511.6–18.40.6530.38–1.14
SexMale20/10718.711.8–27.41.3140.78–2.21χ2 = 1.064
df = 1
p = 0.302
Female133/89314.912.6–17.41-
ParityNulliparous29/19814.611.7–20.21.2030.63–2.29χ2 = 0.353
df = 2
p = 0.838
Multiparous82/58214.112.7–18.71.1910.65–2.18
Primiparous22/11319.47.5–21.41-
LactationLactating 34/3998.516.8–24.02.8001.84–4.25χ2 = 23.1064
df = 1
p = 0.000
Non-lactating 99/49420.05.8–11.41-
Reproductive StatusPregnant 43/32010.614.2–20.61.7471.15–2.65χ2 = 10.483
df = 1
p = 0.002
Non-pregnant 99/57317.37.5–14.61-
Rep. DisorderAbortion48/9351.641.0–62.110.1746.37–16.26χ2 = 133.984
df = 4
p = 0.000
Premature delivery12/3732.418.0–49.84.5782.21–9.47
Stillbirth14/4431.818.6–47.64.4512.26–8.75
Repeat breeding1/425.00.6–80.63.1790.33–30.93
No78/8229.57.6–11.71-
Body ConditionWeak32/7045.733.7–58.17.9844.73–13.49χ2 = 124.868
df = 3
p = 0.000
Does with weak kid delivery22/4944.930.7–59.87.7254.20–14.20
Ewes with weak lamb delivery20/5337.724.8–52.15.7463.15–10.49
Apparently good79/8289.57.6–11.71-
Table
Table 3. Detection of Coxiella burnetii DNA in 29 pooled serum samples.
Table 3. Detection of Coxiella burnetii DNA in 29 pooled serum samples.
Serum PoolsPositive/TestedPrevalence, %95% CIOR95% CIStatistics
Sheep 2/1513.31.7–40.51-χ2 = 0.333
df = 1
p = 0.564
Goat3/1421.44.7–50.81.610.10–4.05
Overall 5/2917.26.60–34.2--
Table
Table 4. Prevalence of C. burnetii DNA in 55 pools of ticks collected from small ruminants.
Table 4. Prevalence of C. burnetii DNA in 55 pools of ticks collected from small ruminants.
Tick PoolsPositive/TestedPrevalence %95% CIOR95% CIStatistics
Sheep9/2931.015.3–50.84.030.05–1.22χ2 = 4.668
df = 1
p = 0.031
Goat 2/267.70.9–25.11-
Overall 11/5520.010.99–32.10--

Share and Cite

MDPI and ACS Style

Ullah, Q.; El-Adawy, H.; Jamil, T.; Jamil, H.; Qureshi, Z.I.; Saqib, M.; Ullah, S.; Shah, M.K.; Khan, A.Z.; Zubair, M.; et al. Serological and Molecular Investigation of Coxiella burnetii in Small Ruminants and Ticks in Punjab, Pakistan. Int. J. Environ. Res. Public Health 2019, 16, 4271. https://doi.org/10.3390/ijerph16214271

AMA Style

Ullah Q, El-Adawy H, Jamil T, Jamil H, Qureshi ZI, Saqib M, Ullah S, Shah MK, Khan AZ, Zubair M, et al. Serological and Molecular Investigation of Coxiella burnetii in Small Ruminants and Ticks in Punjab, Pakistan. International Journal of Environmental Research and Public Health. 2019; 16(21):4271. https://doi.org/10.3390/ijerph16214271

Chicago/Turabian Style

Ullah, Qudrat, Hosny El-Adawy, Tariq Jamil, Huma Jamil, Zafar Iqbal Qureshi, Muhammad Saqib, Shakeeb Ullah, Muhammad Kamal Shah, Alam Zeb Khan, Muhammad Zubair, and et al. 2019. "Serological and Molecular Investigation of Coxiella burnetii in Small Ruminants and Ticks in Punjab, Pakistan" International Journal of Environmental Research and Public Health 16, no. 21: 4271. https://doi.org/10.3390/ijerph16214271

APA Style

Ullah, Q., El-Adawy, H., Jamil, T., Jamil, H., Qureshi, Z. I., Saqib, M., Ullah, S., Shah, M. K., Khan, A. Z., Zubair, M., Khan, I., Mertens-Scholz, K., Henning, K., & Neubauer, H. (2019). Serological and Molecular Investigation of Coxiella burnetii in Small Ruminants and Ticks in Punjab, Pakistan. International Journal of Environmental Research and Public Health, 16(21), 4271. https://doi.org/10.3390/ijerph16214271

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