Next Article in Journal
Molecular Pathotyping of Plasmodiophora brassicae—Genomes, Marker Genes, and Obstacles
Previous Article in Journal
First Detection of the West Nile Virus Koutango Lineage in Sandflies in Niger
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Q Fever: Seroprevalence, Risk Factors in Slaughter Livestock and Genotypes of Coxiella burnetii in South Africa

by
Maruping Mangena
1,2,*,
Nomakorinte Gcebe
2,
Rian Pierneef
3,
Peter N. Thompson
1 and
Abiodun A. Adesiyun
1,4
1
Department of Production Animal Studies, Faculty of Veterinary Science, University of Pretoria, Private Bag X04, Onderstepoort, Pretoria 0110, South Africa
2
Agricultural Research Council–Bacteriology and Zoonotic Diseases Diagnostic Laboratory, Onderstepoort Veterinary Research, Private Bag X 05, Onderstepoort, Pretoria 0110, South Africa
3
Agricultural Research Council-Biotechnology Platform, 100 Old Soutpan Road, Onderstepoort, Pretoria 0110, South Africa
4
School of Veterinary Medicine, Faculty of Medical Sciences, University of the West Indies, St. Augustine Campus, St. Augustine, Trinidad and Tobago
*
Author to whom correspondence should be addressed.
Pathogens 2021, 10(3), 258; https://doi.org/10.3390/pathogens10030258
Submission received: 20 January 2021 / Revised: 5 February 2021 / Accepted: 7 February 2021 / Published: 24 February 2021

Abstract

:
Q fever is a neglected zoonosis in South Africa, causing significant losses in livestock and game animals through reproductive disorders. However, there are limited studies on the extent of Coxiella burnetii infections in livestock in South Africa. Further, there is also lack of knowledge about the types of C. burnetii strains that are currently circulating in the country. Therefore, a cross-sectional, abattoir-based study was conducted to determine the seroprevalence of C. burnetii and associated risk factors, and to characterize C. burnetii strains from slaughter livestock at red meat abattoirs in Gauteng, South Africa. Of the 507 animals tested, 6.9% (95% CI: 4.9–9.5%) were positive for antibodies against C. burnetii. The seroprevalence was 9.4% (31/331) in cattle, 4.3% (3/69) in sheep, and 0.9% (1/107) in pigs. Out of the 63 tissue samples from 35 seropositive animals including material from two sheep aborted fetuses from Mangaung district (Free State province), 12.7% (8/63) tested positive by IS1111 PCR. Genotyping of the eight PCR-positive tissues from eight animals by MLVA revealed two novel genotypes, not available in Coxiella MLVA databases. It is concluded that slaughter animals pose a risk of exposing abattoir and farm workers to C. burnetii in South Africa.
Keywords:
LCVs; SCVs; ELISA; tissues; PCR; MLVA

1. Introduction

Q fever is caused by the obligate intercellular bacterium, Coxiella burnetii. The pathogen belongs to gamma subdivision of Proteobacteria, order Legionellales, family Coxiellaceae and C. burnetii is the only species in the genus [1]. The ability to express resistant cell forms enables C. burnetii to survive in harsh environments [2]. Once the bacterium enters the phagolysosome of the eukaryotic cell, it undergoes incompletely uncharacterized life cycle forms [3,4]. These are two different morphological forms, namely large (LCVs) and small cell variants (SCVs) [5]. Large cell variants have similar characteristics as typical gram-positive bacteria during exponential growth. These include distinct outer membrane, periplasmic space, cytoplasmic membrane, and diffuse nucleoid, reaching more than 1 µm in length [5]. Compared to LVCs, SCVs are smaller, 0.2 to 0.5 µm in diameter and have electron-dense, condensed chromatin, and condensed cytoplasm. These SCVs are resistant to osmotic shock, oxidative stress, heat shock, sonication, and pressure, one reason that C. burnetii is able to withstand harsh environmental conditions [6].
Coxiella burnetii infects a wide range of hosts including cattle, sheep, goats, and wild animals such as deer, buffaloes, squirrels, and rabbits [7]. The bacterium is excreted in milk, feces and birth products of infected animals. Animals mainly get infected through inhalation of dried aerosolized particles originating from milk, feces, and birth products such as placenta of infected animals [1].
There is evidence that Q fever infections cause huge economic losses in livestock and wildlife through late abortions, premature births, and low birth weights [8]. For instance, in the Netherlands 50,355, goats and sheep were culled between December 2009 and June 2010 on 89 bulk tank milk positive farms causing huge losses [9]. Furthermore, during the outbreak in the Netherlands, more than 4000 confirmed Coxiella human cases were reported with 20% of the infected individuals hospitalized [10]. Despite this evidence, Q fever is not widely known to be endemic in South Africa and thus the disease is not recognized as a controlled and notifiable disease in terms of the Animal Diseases Act 35 of 1984. Thus, there is little known on the prevalence of Q fever in livestock and wildlife in South Africa. The disease goes unnoticed causing huge economic losses through reproductive disorders such as late abortions as it is not continuously monitored.
Despite a study by [11] recently reporting Q fever seroprevalence of 38.4% in cattle in Mnisi community, Bushbuckridge municipality, South Africa, little is known on the prevalence of Q fever in livestock or wildlife in South Africa, albeit in red meat abattoirs. Furthermore, studies by [11,12] were based only on serological diagnosis of Q fever focusing on detection of antibody production. The study by [13] focused on Coxiella PCR detection from ticks, and in all the above-mentioned studies, no characterization of the Coxiella genotypes was conducted. Thus, there is also no knowledge of genotypes of C. burnetii currently circulating in the country. Due to limited data available on Q fever seroprevalence and lack of knowledge on C. burnetii strains currently circulating in South Africa, we investigated the seroprevalence and associated risk factors, and further characterized C. burnetii strains circulating in slaughter animals at red meat abattoirs in five districts of Gauteng, the most populous province of South Africa, and from a farm in Mangaung district in the Free State province.

2. Materials and Methods

2.1. Study Area

A serological study was conducted in Gauteng province which is situated in the Highveld region of South Africa with a total area of 18,176 km2. The province consisted of six districts, namely City of Tshwane located in the north-eastern part of the province covering an area of 6298 km2; City of Johannesburg, 1645 km2; Sedibeng district covering most of south-eastern and western Johannesburg with an area of 4173 km2; Ekurhuleni district in eastern and north-eastern part of Johannesburg with an area of 1975 km2; West Rand, which covers the south-western part of the province with an area of 4087 km2; and Metsweding, which was later merged with Tshwane, located in the north-east and western part of the province with total area of 1643 km2 (Figure 1). For IS1111 PCR and MLVA genotyping, we also included pooled organs from two sheep fetuses (liver, spleen, and lungs) which were obtained from a farm in Mangaung district in the Free State province. The Free State province has with an area of 6284 km2. A questionnaire was used to obtain demographic information on the animals since they originated from different provinces in South Africa. Other data such as sex, origin of animals (auctions/feed lots), and species were obtained from respective abattoir managers using questionnaires.

2.2. Sample Size and Study Population

The sample size for serological tests was estimated using the formula of [15]: n = 1.962 × p × (1 − p)/d2 where p = estimated prevalence, d = precision, n = estimated sample size. Considering that there are no current data on the prevalence of Q fever in livestock sampled in red meat abattoirs in South Africa, we used an estimated prevalence of 50% i.e., 0.5 and a precision of 4.5% i.e., 0.045, giving a required minimum sample size of 475. However, for the study, a total of 507 serum samples and 1018 corresponding reproductive tissues were collected at during slaughter at 19 randomly selected red meat abattoirs in Gauteng province (Table 1). These abattoirs consisted of 16 high-throughput (HT) (daily output ≥ 25 carcasses) and three low-throughput (LT) red meat abattoirs in five districts of Gauteng (excluding City of Johannesburg). The 1018 reproductive tissues collected from 507 slaughter animals during slaughter in Gauteng red meat abattoirs comprised 355 penises and 355 testes. A further 80 ovaries, 4 oviducts, 79 mammary glands, and 145 uterus tissue samples were collected from slaughter livestock in Gauteng red meat abattoirs during slaughter (Table 1). In addition, there were also two pooled aborted sheep tissue samples (liver, spleen, and lungs) from diagnostic samples collected from a farm in the Manguang district, Free State province (Table 1). Each abattoir was visited once and based on the study design 30 and 20 serum samples were collected from HT and LT abattoirs respectively, whenever the number of slaughtered animals was available. Blood samples were collected using Bencton, Dickinson and Company (BD)-Vacutainer® SSTTM II Advance 10 mL serum collection tubes (Franklin Lakes, NJ, USA) and transported to laboratory at room temperature. Serum was harvested same day by centrifuging the clotted blood in collection tubes at 1000× g for 10 min and stored at −20 °C until analyzed. Tissues were collected in zip lock bags and transported to the laboratory the same day at 4 °C to be stored at −20 °C until analysis.

2.3. Serological Testing

For detection of IgG antibodies against C. burnetii, the IDEXX Q fever 2/strip antibody test kit was used according to manufacturer’s instructions (IDDEX Laboratories, Liebelfld-Bern, Switzerland). All reagents from the IDEXX Q Fever 2/strip antibody test kit together with the frozen serum samples were brought to room temperature and concentrated wash buffer diluted 10 times with distilled water. The diluted buffer was used to dilute the negative, positive controls, as well as the serum samples 400-fold. Exactly 100 µL of the diluted control and collected serum samples were transferred into the C. burnetii antigen-coated plate wells and incubated at 37 °C for one hour (h). The plate was then washed three times with 300 µL of the 10-fold diluted wash buffer using a BioTek ELx50 automated microplate washer (BioTEk, Winooski, VT, USA), 100 µL of conjugate was added to each well, and the plates incubated at 37 °C for 1 h. The plate was washed again three times with 300 µL of the 10-fold diluted wash buffer, 100 µL of 3,3′,5,5′-Tetramethylbenzidine (TMB) substrate was then added to each well and the plates incubated at room temperature for 15 min, away from direct light. After 15 min, 100 µL of stop solution was added to each well and the plates immediately read at 450 nm using Thermo Labsystems Multiskan MS Original microplate reader (Thermo Fischer Scientific, Waltham, MA, USA). The assay was validated as follows: For the assay to be valid, the average optical density value of the two negative controls (NCx) at 450 nanometers (A450) should be less than or equals to (≤) 0.500. The average value of the two positive controls (PCx) at 450 nanometers (A450) should be less than or equals to (≤) 2.500. Then PCx-NCx (A450) should be greater than or equals to (≥) 0.300. Sample to positive (S/P) ratio was calculated according to the formula:
S / P   % = 100 × Sample ( A 450 ) NCx ( A 450 ) PCx ( A 450 ) NCx ( A 450 )
For interpretation of results, S/P% < 30% represented a negative result, 30% ≤ S/P% < 40% a suspect, while S/P% ≥ 40% represented a positive result. Seven samples gave suspect results. These samples were repeated once using three different aliquots of the same sample and all of them were negative. The identified 63 reproductive tissues from Gauteng red meat abattoirs corresponding to positive serum samples, including two pooled fetus tissue samples, were subjected to PCR confirmation and MLVA typing.

3. Molecular Characterization

3.1. DNA Extraction and PCR for Detection of C. burnetii

PCR confirmation was conducted on reproductive tissues from all seropositive animals from this study as well as for diagnostic tissue samples. The tissue samples consisted of 61 tissues from 35 animals (12 mammary glands, 7 uterus, 18 penis, 21 testes, and 3 ovaries) and two pooled aborted sheep fetus tissues from two animals (spleen, liver, and lungs) (Table 1). Tissue samples were cut into small pieces and 10 g from each sample in placed 10 mL ice cold buffered phosphate saline (PBS) pH 7.4 in 50 mL bead ruptor homogenizing tubes containing 2.8 mm ceramic beads. The tissue samples were then homogenized using the automated BEAD RUPTOR ELITE Bead Mill homogenizer (Omni International, Kennesaw, GA, USA). The tissues were homogenized at a speed of three meters per second (3 m/s) for 90 s (s). DNA extraction from the homogenates was conducted using the Qiagen DNeasy® blood and tissue kit as previously described by [16]. The homogenates were centrifuged for 15 min at 4000 rpm and 200 µL of the supernatant transferred to 2 mL centrifuge tubes. To the tubes, 180 µL of tissue lysis (ATL) buffer, 20 µL of proteinase K were added, suspension vortexed, and incubated at 56 °C overnight. After overnight incubation, 200 µL of lysis (AL) buffer was added and the suspension vortexed for 15 sand incubated at 70 °C for 10 min. Absolute ethanol (200 µL) was added to the mixture, vortexed, and transferred to DNeasy® spin columns. The columns were then washed twice with wash buffers; AW1 and AW2, respectively. DNA was eluted from the columns with 200 µL of elution buffer (AE).
PCR for detection of C. burnetti in tissues of the seropositive animals was conducted in a 50 µL reaction targeting the multi-copy transposase gene in insertion element; IS1111 [17] using primers described in Table 2 The Coxiella gene fragment (gblock) from Integrated DNA Technologies (Coralville, IA, USA) and water were used as positive and negative controls in the reaction, respectively. The reaction mixture contained 400 nM of each primer as listed in Table 2 (IS1111F and IS1111R), 25 µL of the Ampliqon 2× Taq DNA polymerase Master Mix Red (Ampliqon A/S, Odense, Denmark) and 10 µL of the extracted DNA. PCR amplification was conducted using BIO RAD T100™ thermal cycler (BIO RAD, Hercules, CA, USA). Cycling conditions consisted of initial denaturation at 95 °C for 15 min, 35 cycles of denaturation 95 °C for 30 s, annealing at 60 °C for 30 s, and extension at 72 °C for 60 s for 35 cycles. Final extension was carried out at 72 °C for 10 min and amplicons visualized on a 1.5% w/v ethidium bromide-stained agarose gel with an expected size of 146 bp [17] estimated using Quick-Load® 100 bp DNA Ladder (New England Biolabs, Ipswich, MA, USA).

3.2. Sequence Confirmation of C. burnetii

PCR confirmation of tissues from two animals (n = 2) STCN17 and 40241 from Gauteng province red meat abattoirs and a farm in Mangaung, Free State province, respectively, was conducted using Sanger sequencing. The IS1111 PCR products of the two tissues were sent to Inqaba Biotechnical industries (Pty) Ltd (Pretoria, South Africa). For sequencing and sequences manually edited using the BioEdit Sequence alignment editor (version 7.2.5). PCR positive DNA from eight seropositive tissue samples from the eight animals (n = 8) was then subjected to MLVA typing.

3.3. MLVA Typing

Eight PCR positive DNA samples from eight animals were genotyped using the Dutch six-locus MLVA panel as previously described by [18,19,20]. This panel utilizes two sets of microsatellite markers. Panel one consists of hexanucleotide repeats (MS27, MS28, and MS34) while panel two is made up of heptanucleotide repeats (MS23, MS24, and MS33) as in Table 2. PCR amplification of the microsatellite markers was carried out in total reaction volume of 20 µL reaction containing 5 µL of the extracted DNA, 10 µL of the Ampliqon 2× Tag DNA polymerase Master Mix Red (Ampliqon A/S, Odense, Denmark. The microsatellite primer pairs are as described in Table 2. The PCR reactions were conducted using BIO RAD T100™ thermal cycler (BIO RAD, Hercules, CA, USA). PCR conditions consisted of initial denaturation at 95 °C for 10 min, 40 cycles of denaturation 95 °C for 30 s, annealing at 60 °C for 30 s, and extension at 72 °C for one min for 40 cycles. Final extension for at 72 °C for 10 min. Microsatellite panel one amplicons visualized on a 2% w/v ethidium bromide-stained agarose gel while panel two amplicons were visualized on a 3% w/v agarose gel [21,22]. Amplicon sizes for each marker were estimated using O’RangeRuler 20 bp DNA ladder (Thermo Fisher Scientific, Waltham, MA, USA).
Table 2. List of PCR primers used in Coxiella IS11111 PCR and MLVA typing.
Table 2. List of PCR primers used in Coxiella IS11111 PCR and MLVA typing.
NamePrimer SequenceReference(s)
IS1111F5′ CGCAGCACGTCAAACCG3′[17]
IS1111R5′TATCTTTAACAGCGCTTGAACGTC3′[17]
MS23F5′CGCMTAGCGACACAACCAC3′[18,20]
MS23R5′GACGGGCTAAATTACACCTGCT3′[18,20]
MS24F5′TGGAGGGACTCCGATTAAAA3′[18,20]
MS24R5′GCCACACAACTCTGTTTTCAG3′[18,20]
MS27F5′TCTTTATTTCAGGCCGGAGT3′[18,20]
MS27R5′GAACGACTCATTGAACACACG3′[18,20]
MS28F5′AGCAAAGAAATGTGAGGATCG3′[18,20]
MS28R5′GCCAAAGGGATATTTTTGTCCTTC3′[18,20]
MS33F5′TCGCGTAGCGACACAACC3′[18,20]
MS33R5′GTAGCCCGTATGACGCGAAC3′[18,20]
MS34F5′TTCTTCGGTGAGTTGCTGTG3′[18,20]
MS34R5′GCAATGACTATCAGCGACTCGAA3′[18,20]

4. Data Analysis

The data were coded using Microsoft Excel and analyzed in Stata 15 (StataCorp, College Station, TX, USA). Exact binomial confidence intervals were calculated for proportions. Univariate associations of sex, species, breed, abattoir throughput, district, and origin of animals with C. burnetii seropositivity were assessed using a two-tailed Fisher’s exact test. All variables were then entered into a mixed-effects multiple logistic regression model with abattoir as a random effect; however, due to collinearity with species, breed was not considered in the multivariable model. The model was developed by backward elimination until all remaining variables were significant (Wald p < 0.05). Model fit was assessed using the Hosmer–Lemeshow goodness-of-fit test.
Coxiella burnetii Sanger sequencing data were analyzed using basic local alignment tool; BLAST; NCBI [23]. MLVA data were analyzed using Bionumerics 7.6 software (Applied Maths, Sint-Martens-Latem, Belgium). Copy numbers were for each marker (MS23, MS24, MS27, MS28, MS33, and MS34) were obtained by extrapolating fragments sizes for each marker with those of the Nine Mile RSA493. From the comparisons of the obtained fragments with the Nine Mile RSA493, a sequence number representing number of repeats for each locus and sample was determined, resulting into the defining of the MLVA profile of the particular sample [24]. The relationship between the genotypes on the study and previously on described genotypes was determined by comparing them with C. burnetii genotypes found in the C. burnetii 2014 Nijmegen database created by J.J.H.C. Tillburg [24]. Unweighted Pair Group Method with Arithmetic mean (UPGMA) was used to calculate genetic distances between the genotypes and used to generating dendrograms and minimum spanning trees. First, distance matrices were calculated using the “daisy” function with the “gower” parameter specified to determine Gower distances with the R package “cluster” [25]. Thereafter, UPGMA trees were constructed and visualized with ggtree [26]. Minimum spanning trees were calculated using the “ape” package [27] with the “mst” function and visualized using “igraph” [28] and “ggnetwork” [29]. R v.4.0.2 [30] was used for the construction of dendrograms and minimum spanning trees.

5. Results

5.1. Serology

Of the 507 animals tested, 6.9% (95% CI: 4.9–9.5%) were positive for antibodies against C. burnetii. The seroprevalence in livestock by species was 9.4% (31/331) in cattle, 4.3% (3/69) in sheep, and 0.9% (1/107) in pigs. The difference in seroprevalence between sexes was most pronounced in cattle, where it was 17/81 (21%) in females vs. 14/250 (6%) in males (Table 3). Significant univariate associations with C. burnetii seropositivity were seen for species, sex, breed, district, and animal origin (Table 3). However, sex and district were no longer significant in the final multivariable logistic regression model (Table 4). The Hosmer–Lemeshow test showed adequate model fit (p = 0.670). Animals from auctions were more likely to have been exposed to C. burnetii than animals from farms and feedlots (OR = 5.7; 95%CI: 2.6–12.4; p < 0.001). Although not significant in the univariate analysis, the multivariable model showed that the odds of Q fever seropositivity in LT abattoirs was significantly higher than in HT abattoirs (OR 4.1; 95%CI: 1.2–14.0; p = 0.023) (Table 4).

5.2. Molecular Detection and MLVA Typing of C. burnetii

Molecular detection of C. burnetii in 63 tissues from 35 seropositive animals, positive by IS1111 PCR showed that 12.7% (8/63) were positive with amplicons approximately 146 bp in size (Figure 2). The eight positive tissues comprised 9.5% (6/63) that originated from cattle in red meat abattoirs of Gauteng province while 3.2% (2/63) were pooled organs from aborted sheep fetuses originating from a farm in Mangaung district, Free State province (Table 1).
The edited sequences from two isolates using BLAST revealed that the two amplicons, STFCNC17 and 40241, had a 99.1% similarity with the partial coding sequence of C. burnetii strain 54T1 transposase gene (MT268532.1). There was uniform amplification for all the six cattle tissue samples and two sheep samples for markers MS23, MS24, MS27, MS28, and MS33 but only the two sheep samples from Mangaung district amplified for MS34. MLVA profiles of six cattle and two sheep isolates were determined and compared with C. burnetii genotypes found in the C. burnetii 2014 Nijmegen database created by J.J.H.C. Tillburg found on the database: (http://mlva.u-psud.fr/mlvav4/genotyping/view.php (accessed on 14 January 2021)). The comparisons revealed that the six cattle tissues isolates (STFC17, RTSC8, RNBRC16, KMLDC8, MGC11, and MCM21) from Gauteng red meat abattoirs were novel genotypes, not available on public databases. These genotypes shared a close relationship with a human blood and valve isolate 20090317Frankrijk004 from Marseille, France (Figure 3) with a distance of three map units (m.u). The 2009031Frankrijk004 isolate from Marseille, France, belongs to MLVA genotype 92 and MST 1 and does not originate from an abattoir (Figure 4). The two sheep isolates from Mangaung district (40241 and 40242) were also previously uncharacterized genotypes not available in the public databases. They also shared a close relationship with a human blood and valve isolate 20090317Frankrijk044 from Marseille, France with a distance of three map units (3 m.u; Figure 5). The Frankrijk044 isolate belong to MLVA genotype 92 and MST 1 and does not originate from an abattoir. Minimum spanning trees (Figure 4 and Figure 6) show that by geographic location, all the eight isolates were closely related to human isolates from France, although they were isolated from slaughter livestock.

6. Discussion and Conclusions

This is the first abattoir-based study to show serological evidence of antibodies to C. burnetii in livestock using indirect ELISA in South Africa. We observed a seroprevalence of 6.9% with the highest found in cattle and the lowest in pigs. The serology results are consistent with a study conducted by [31] in Northern Ireland. This study reported 64.5% Q fever seroprevalence in dairy cattle as compared to 21% in sheep using indirect immmunofluorescence assay. Similarly, a considerably higher seroprevalence of Q fever (39%) was detected in cattle in Zimbabwe [32] also using the indirect immmunofluorescence assay compared to 9.4% in cattle observed in our study using indirect ELISA. Other studies have found little differences in Q fever seroprevalence between cattle (21.7%) and sheep (28.4%) in the Volta region of Ghana [33], which is consistent with our study as there was no significant difference in seroprevalence between cattle and sheep in the multivariable model (p = 0.369).
Another study by [34] reported highest seroprevalence of Q fever in pigs (11.3%) than in both cattle (4.3%) and sheep (0.0) in Trinidad using capillary agglutination test. This is in contrast with findings in the present study where we observed the lowest Q fever seroprevalence in pigs (0.9%) followed by sheep (4.3%) with cattle having the highest prevalence (9.4%).
In this study, breed differences were noted, with highest seroprevalence (27%) detected in Nguni cattle compared to other breeds. There is therefore a possibility that Nguni breed of cattle may have increased exposure potential for C. burnetii in the country and possibly may expose abattoir workers and veterinarians to Q fever. This is because one of the Q fever transmission routes from animals to humans is direct contact with infected carcasses and birth products [1].
We observed a higher seroprevalence of C. burnetii in females (11.8%) than in males (4.8%), although the differences were no longer significant in the multivariable model (p = 0.258). This result is consistent with a study done on Danish cattle by [35], which found a higher true seroprevalence of Q fever in females (9.4%) than in males (2.6%). This the higher Q fever seroprevalence in females than males can be explained in part that an infected male could infect several females during reproductive cycles. Coxiella burnetii have been detected in bulls’ semen, indicating sexual transmission of the disease [36].
We observed in the study that livestock originating from auctions were more likely to have been exposed to C. burnetii compared to livestock from farms and feedlots. This could be partly explained by the fact that animals purchased from auctions may have originated from a wide range of farms with possible differences in exposure to C. burnetii, together with possible close contact of animals in pens at auctions, thus facilitating the spread of the pathogen among animals.
The current data showed statistically significant (p = 0.003) differences in prevalence of antibodies against C. burnetii between districts where slaughter livestock were sampled. This might not represent the effect of geographical locations of farms and abattoirs in Gauteng province on the prevalence of antibodies to C. burnetii. This is because animals slaughtered in abattoirs may originate from other provinces across the country. However, based on the locations of the abattoirs, our findings display differences in seroprevalence of Q fever ranging from 0.0% in the West Rand to 12.2% in Tshwane districts of Gauteng province. The significantly high seroprevalence of Q fever in the then South-Eastern Transvaal (now Mpumalanga) observed by [12] corresponds with a report by [11] who observed a higher Q fever seroprevalence in the Bushbuckridge area of Mpumalanga province as compared to the findings in the study. Similarly, Ref. [12] reported a significantly higher seroprevalence of Q fever in the then South-Eastern Transvaal (now Mpumalanga) than in the then Western Transvaal (now North West Province), which was attributed to the differences in the distribution of the blue tick, Boophilus decoloratus. This tick was a known important vector of the pathogen in 1985, which was a dry year [12]. During that time in 1985, the abundance of this tick was low in the Western Transvaal corresponding with low Q fever seroprevalence in that area. This may suggest that B. decoloratus is involved in the transmission of Q fever. Another tick species, Haemaphysalis leachi, a known Q fever carrier in South Africa had a distribution closely like that of Q fever in the same year, 1985 [12]. This suggests that this tick species might also be involved in Q fever transmission among cattle in Gauteng province.
We observed Q fever prevalence of 12.7% by PCR in the study. This is relatively low compared to other studies such as [13] who reported 41% prevalence. Other similar studies have also reported higher Q fever prevalence elsewhere using PCR. For instance, Ref. [22] reported C. burnetii prevalence of 16.6% from ruminants and wildlife in Portugal using IS1111 PCR. In the study by [37], goats had the highest prevalence with 23.5%, cattle with 20.8%, and sheep with 10%. These observations by [37] differ with the findings in our study where cattle reported highest Q fever prevalence with 9.5% (6/63) while sheep reported 3.2% (2/63) using IS1111 PCR.
Molecular characterization of the C. burnetii isolates using MLVA showed that the eight isolates are novel genotypes. However, they are closely related to genotypes from France, Europe, albeit human isolates, with genetic distance of 3 m.u. This could suggest that two genotypes were discovered, one that could be specific to cattle and the other to sheep. This observation also suggests that the isolates from this study in South Africa and the ones from France might have originated from different locations but shared a common ancestor.
Based on the seroprevalence and PCR detection of Q fever, cattle had the highest level of exposure to the pathogen compared to sheep and pigs. There is also the possibility that cattle pose the highest risk of abattoir workers being exposed to the pathogen at slaughter and processing. Although the prevalence was not high, the detection of widespread exposure to C. burnetii in slaughter livestock cannot be ignored. Hence, the pathogen may be of economic importance to the livestock industry and of zoonotic significance to personnel in the farming sector, including livestock farmers, animal attendants, and veterinarians, and to consumers of under-cooked products such as milk. It is therefore imperative to conduct more studies on Q fever in livestock in other provinces of South Africa. These studies should specifically focus on the isolation and molecular, genomic, and proteomic analyses of circulating strains of C. burnetii in the country.
In conclusion, we have documented the seroprevalence and risk factors associated with C. burnetii in red meat abattoirs of Gauteng province, South Africa. In addition, we have detected C. burnetii by PCR in slaughter animals as well as from aborted fetuses. Using MLVA we detected two novel genotypes in sheep and cattle, respectively.

Author Contributions

Conceptualization, A.A.A. and N.G.; methodology, M.M., A.A.A., N.G., and P.N.T.; software, P.N.T., R.P., and N.G.; Validation, M.M., A.A.A., N.G., and P.N.T.; formal analysis, M.M., A.A.A., N.G., P.N.T., and R.P.; investigation. M.M., A.A.A., N.G., and P.N.T.; resources, A.A.A., N.G., and P.N.T.; data curation, M.M., A.A.A., N.G., P.N.T., and R.P.; writing—original draft preparation, M.M.; writing—review and editing A.A.A., N.G., P.N.T., R.P.; visualization, M.M., A.A.A., N.G., P.N.T., and R.P.; supervision, A.A.A., N.G., and P.N.T.; project administration, A.A.A. and N.G.; funding acquisition, A.A.A. and N.G. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for the study was obtained from National Research Foundation (Grant # SFH150702122812), the Red Meat Research and Development in South Africa (RMRD-SA), Gauteng Department of Agriculture and Rural Development (GDARD), CSIR IBS-HCD Bursary and the Department of Trade and Industry (DTI)-THRIP.

Institutional Review Board Statement

Permission to conduct the study was obtained from the Department of Agriculture, Forestry and Fisheries (DAFF) in terms of the Animal Diseases Act 35 of 1984, Reference number 12/11/1/1/6. The study was approved by the University of Pretoria Animal Ethics Committee (V085-17). Ethical clearance to perform the study was also obtained from the Agricultural Research Council-Onderstepoort Veterinary Research, Animal Use Ethics Committee (AEC 12.16).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Sequence confirmation of Coxiella amplicon sequences in the study was conducted by comparing the sequences with those available on BLAST. This information can be found on https://blast.ncbi.nlm.nih.gov/Blast.cgi? (accessed on 14 January 2021) The Coxiella MLVA Publicly available datasets were used for comparisons with genotypes discovered in the study. This data can be found on http://mlva.u-psud.fr/mlvav4/genotyping/view.php (accessed on 14 January 2021).

Acknowledgments

We are grateful to Charles Katsande, Veterinary Public Health officials and red meat abattoirs from Gauteng Province for assistance in obtaining permission from the abattoir owners to conduct the study at their facilities, and for scheduling the visits to the abattoirs to collect samples. We also acknowledge Olivier Mathee and Brian Peba of the Agricultural Research Council-Epidemiology, Parasites and Vectors (ARC-EPV) Laboratory, for technical assistance with ELISA and equipment. We also thank Kudakwashe Jambwa and the whole Bacterial PCR Laboratory team at the Agricultural Research Council-Onderstepoort Veterinary Research (ARC-OVR) campus for technical assistance with PCR assays. We also acknowledge Henriette Van Heerden for assistance in facilitating scheduling of abattoir sampling visits. Our gratitude also go to Solomon Odemuyiwa from University of Missouri for assistance with Coxiella gblocks and PCR methodology. We also thank Enerst Ngoepe for revision of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Angelakis, E.; Raoult, D. Q Fever. Vet. Microbiol. 2010, 140, 297–309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Seshadri, R.; Paulsen, I.T.; Eisen, J.A.; Read, T.D.; Nelson, K.E.; Nelson, W.C. Complete genome sequence of the Q fever pathogen Coxiella burnetii. Proc. Natl. Acad. Sci. USA 2003, 100, 5455–5460. [Google Scholar] [CrossRef] [Green Version]
  3. Davis, G.E.; Cox, H.R. A filter-passing infectious agent isolated from ticks. IN. Isolation from Dermacentor andersonii, reactions in animals and filtration. Public Health Rep. 1938, 53, 2259. [Google Scholar] [CrossRef]
  4. Heinzen, R.A.; Hackstadt, T.; Samuel, J.E. Developmental biology of Coxiella burnetii. TIMS 1999, 7, 149–154. [Google Scholar] [CrossRef]
  5. McCaul, T.F.; Williams, J.C. Developmental cycle of Coxiella burnetii: Structure and morphogenesis of vegetative and sporogenic differentiations. J. Bacteriol. 1981, 147, 1063–1076. [Google Scholar] [CrossRef] [Green Version]
  6. Amano, K.J.; Williams, C.; McCaul, T.E.; Peacock, M.G. Biochemical and immunological properties of Coxiella burnetii cell wall andpeptidoglycan-protein complex fractions. J. Bacteriol. 1984, 160, 982–988. [Google Scholar] [CrossRef] [Green Version]
  7. Marrie, T.J.; Schlech, W.F., 3rd; Williams, J.C.; Yates, L. Q fever pneumonia associated with exposure to wild rabbits. Lancet 1986, 1, 427–429. [Google Scholar] [CrossRef]
  8. Honarmand, H. Q fever: An old but still a poorly understood disease. Interdiscip. Perspect Infect Dis. 2012, 131932. [Google Scholar] [CrossRef] [Green Version]
  9. Whelan, J.; Schimmer, B.; Schneeberger, P.; Meekelenkamp, J.; Ijff, A.; Van der Hoek, W.; Van Beest, M.R. Q fever among culling workers, the Netherlands, 2009–2010. Emerg. Infect Dis. 2011, 17, 1719–1723. [Google Scholar] [CrossRef] [PubMed]
  10. Van der Hoek, W.; Dijkstra, F.; Schimmer, B.; Schneeberger, P.M.; Vellema, P.; Wijkmans, C.; Ter Schegget, R.; Hackert, V.; Van Duynhoven, Y. Q fever in the Netherlands: An update on the epidemiology and control measures. Euro Surveill 2010, 15. Available online: http://www.njmonline.nl/getpdf.php?t=a&id=10000657 (accessed on 14 January 2021).
  11. Adesiyun, A.A.; Knobel, D.L.; Thompson, P.N.; Wentzel, J.; Kolo, F.B.; Kolo, A.O.; Conan, A.; Simpson, G.J.G. Sero-Epidemiological Study of Selected Zoonotic and Abortifacient Pathogens in Cattle at a Wildlife-Livestock Interface in South Africa. Vector Borne Zoonotic Dis. 2020, 4, 258–267. [Google Scholar] [CrossRef]
  12. Gummow, B.; Poerstamper, N.; Herr, S. The incidence of Coxiella burnetii antibodies in cattle in the Transvaal. Onderstepoort J. Vet. Res. 1987, 54, 569–571. [Google Scholar] [PubMed]
  13. Mtshali, K.; Khumalo, Z.; Nakao, R.; Grab, D.J.; Sugimoto, C.; Thekisoe, O. Molecular detection of zoonotic tick-borne pathogens from ticks collected from ruminants in four South African provinces. J. Vet. Med. Sci. 2016, 77, 1573–1579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Google Maps. Available online: https://maps.google.com/maps (accessed on 1 February 2021).
  15. Thrusfield, M. Veterinary Epidemiology; John Wiley & Sons: Hoboken, NJ, USA, 2013; p. 1289. [Google Scholar]
  16. Ozkaraca, M.; Ceribasi, S.; Ceribasi, A.O.; Kilic, A.; Altun, S.; Comakli, S.; Ongor, H. Determination of Coxiella burnetii in bovine foetuses using PCR and immunohistochemistry. Vet. Med. Czech 2017, 61, 421–427. [Google Scholar] [CrossRef] [Green Version]
  17. De Bruin, A.; Van Alphen, P.T.; Van der Plaats, R.Q.; De Heer, L.N.; Reusken, C.B.; Van Rotterdam, B.J.; Janse, I. Molecular typing of Coxiella burnetii from animal and environmental matrices during Q fever epidemics in the Netherlands. BMC Vet. Res. 2012, 8, 165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Arricau-Bouvery, N.; Hauck, Y.; Bejaoui, A.; Frangoulidis, D.; Bodier, C.; Souriau, A.; Meyer, H.; Neubauer, H.; Rodolakis, A.; Vergnaud, G. Molecular characterization of Coxiella burnetii isolates by infrequent restriction site-PCR and MLVA typing. BMC Microbiol. 2006, 6, 38. [Google Scholar] [CrossRef] [Green Version]
  19. Santos, A.S.; Tilburg, J.J.; Botelho, A.; Barahona, M.J.; Nuncio, M.S.; Nabuurs-Franssen, M.H.; Klaassen, C.H. Genotypic diversity of clinical Coxiella burnetii isolates from Portugal based on MST and MLVA typing. Int. J. Med. Microbiol. 2012, 302, 253–256. [Google Scholar] [CrossRef]
  20. Svraka, S.; Toman, R.; Skultety, L.; Slaba, K.; Homan, W.L. Establishment of a genotyping scheme for Coxiella burnetii. FEMS Microbiol. Lett. 2006, 254, 268–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Ceglie, L.; Guerrini, E.; Rampazzo, E.; Barberio, A.; Tilburg, J.J.; Hagen, F.; Lucchese, L.; Zuliani, F.; Marangon, S.; Natale, A. Molecular characterization by MLVA of Coxiella burnetii strains infecting dairy cows and goats of North-eastern Italy. Microbes. Infect. 2015, 17, 776–781. [Google Scholar] [CrossRef] [PubMed]
  22. Glazunova, O.; Roux, V.; Freylikman, O.; Sekeyova, Z.; Fournous, G.; Tyczka, J.; Tokarevich, N.; Kovacova, E.; Marrie, T.J.; Raoult, D. Coxiella burnetii genotyping. Emerg. Infect Dis. 2005, 11, 1211–1217. [Google Scholar] [PubMed]
  23. Boratyn, G.M.; Camacho, C.; Cooper, P.S.; Coulouris, G.; Fong, A.; Ma, N.; Madden, T.L.; Matten, W.T.; McGinnis, S.D.; Merezhuk, Y.; et al. BLAST: A more efficient report with usability improvements. Nucleic Acids Res. 2013, 41, W29–W33. [Google Scholar] [CrossRef] [Green Version]
  24. Tilburg, J.J.; Rossen, J.W.; Van Hannen, E.J.; Melchers, W.J.; Hermans, M.H.; Van de Bovenkamp, J.; Roest, H.J.; De Bruin, A.; Nabuurs-Franssen, M.H.; Horrevorts, A.M.; et al. H Genotypic diversity of Coxiella burnetii in the 2007–2010 Q fever outbreak episodes in The Netherlands. J. Clin. Microbiol. 2012, 50, 1076–1078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Maechler, M. Finding Groups in Data: Cluster Analysis Extended Rousseeuw et. R package Version 2; Wiley-Interscience: Hoboken, NJ, USA, 2019. [Google Scholar]
  26. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T.Y. gtree: An R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  27. Paradis, E.; Schliep, K. Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 2019, 35, 52628. [Google Scholar] [CrossRef] [PubMed]
  28. Csardi, G.; Nepusz, T. The igraph software package for complex network research. InterJ. Complex Syst. 2006, 1695, 1–9. [Google Scholar]
  29. Briatte, F. Ggnetwork: Geometries to Plot Networks with ’ggplot2’. 2016. R package version 0.5.1. pp. 27–29. Available online: https://github.com/briatte/ggnetwork (accessed on 14 January 2021).
  30. Team R Core. R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria. 2013. Available online: http://www.R-project.org (accessed on 29 October 2020).
  31. McCaughey, C.; Murray, L.J.; McKenna, J.P.; Menzies, F.D.; McCullough, S.J.; O′Neill, H.J.; Wyatt, D.E.; Cardwell, C.R.; Coyle, P.V. Coxiella burnetii (Q fever) seroprevalence in cattle. Epidemiol. Infect. 2010, 138, 21–27. [Google Scholar] [CrossRef] [Green Version]
  32. Kelly, P.J.; Matthewman, L.A.; Mason, P.R.; Raoult, D. Q fever in Zimbabwe. A review of the disease and the results of a serosurvey of humans, cattle, goats and dogs. S. Afr. Med. J. 1993, 83, 21–25. [Google Scholar] [PubMed]
  33. Johnson, S.A.; Kaneene, J.B.; Asare-Dompreh, K.; Tasiame, W.; Mensah, I.G.; Afakye, K.; Simpson, S.V.; Addo, K. Seroprevalence of Q fever in cattle, sheep and goats in the Volta region of Ghana. Vet. Med. Sci. 2019, 5, 402–411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Adesiyun, A.A.; Cazabon, E.P. Seroprevalences of Brucellosis, Q fever and toxoplasmosis in slaughter livestock in Trinidad. Rev. Elev. Med. Vet. Pays Trop 1996, 49, 28–30. [Google Scholar]
  35. Paul, S.; Agger, J.F.; Agerholm, J.S.; Markussen, B. Prevalence and risk factors of Coxiella burnetii seropositivity in Danish beef and dairy cattle at slaughter adjusted for test uncertainty. Prev. Vet. Med. 2014, 113, 504–511. [Google Scholar] [CrossRef]
  36. Kruszewska, D.; Tylewska-Wierzbanowska, S. Isolation of Coxiella burnetii from bull semen. Res. Vet Sci. 1997, 62, 299–300. [Google Scholar] [CrossRef] [PubMed]
  37. Cumbassá, A.; Barahona, M.J.; Cunha, M.V.; Azórin, B.; Fonseca, C.; Rosalino, L.M.; Tilburg, J.; Hagen, F.; Santos, A.S.; Botelho, A. Coxiella burnetii DNA detected in domestic ruminants and wildlife from Portugal. Vet. Microbiol. 2015, 180, 136–141. [Google Scholar]
Figure 1. Map of Gauteng province districts where blood and reproductive tissues were collected from red meat abattoirs [14].
Figure 1. Map of Gauteng province districts where blood and reproductive tissues were collected from red meat abattoirs [14].
Pathogens 10 00258 g001
Figure 2. Gel electrophoresis of seropositive tissue DNA Using IS1111 PCR. The first lane is Quick-Load® 100 bp DNA Ladder (New England Biolabs, Ipswich, MA, USA). The positive control is Coxiella gene fragment (gblock) from IDT (Coralville, IA, USA), negative control is distilled water and blank is an empty lane. RTSC11 Mammary gland, STFNC17 ovary, RNBRC16 penis, KMFC1 penis, MGC11 testes, and MCMC21 penis are cattle samples from Gauteng red meat abattoirs while F40241 and 40242 are pooled sheep fetuses samples from a farm in Mangaung district. The expected PCR product is approximately 146 bp as shown in the figure.
Figure 2. Gel electrophoresis of seropositive tissue DNA Using IS1111 PCR. The first lane is Quick-Load® 100 bp DNA Ladder (New England Biolabs, Ipswich, MA, USA). The positive control is Coxiella gene fragment (gblock) from IDT (Coralville, IA, USA), negative control is distilled water and blank is an empty lane. RTSC11 Mammary gland, STFNC17 ovary, RNBRC16 penis, KMFC1 penis, MGC11 testes, and MCMC21 penis are cattle samples from Gauteng red meat abattoirs while F40241 and 40242 are pooled sheep fetuses samples from a farm in Mangaung district. The expected PCR product is approximately 146 bp as shown in the figure.
Pathogens 10 00258 g002
Figure 3. Phylogenic tree of C. burnetii genotypes of six tissues from six animals (STFC17, RTSC8, RNBRC16, KMLDC8, MGC11, and MCM21) originating from Gauteng red meat abattoirs based on MLVA-6 database (http://mlva.u-psud.fr/mlvav4/genotyping/view.php (accessed on 14 January 2021)) created by Dr J.J.H.C. Tillburg. The genotype did not amplify for all tissues for MS34 (data not shown) and Nine Mile RSA493 with copy numbers 9-27-4-6-9-5 was used reference (data not shown).
Figure 3. Phylogenic tree of C. burnetii genotypes of six tissues from six animals (STFC17, RTSC8, RNBRC16, KMLDC8, MGC11, and MCM21) originating from Gauteng red meat abattoirs based on MLVA-6 database (http://mlva.u-psud.fr/mlvav4/genotyping/view.php (accessed on 14 January 2021)) created by Dr J.J.H.C. Tillburg. The genotype did not amplify for all tissues for MS34 (data not shown) and Nine Mile RSA493 with copy numbers 9-27-4-6-9-5 was used reference (data not shown).
Pathogens 10 00258 g003
Figure 4. Minimum spanning tree showing the relationship between the eight tissue samples from Gauteng tissue samples from Gauteng province red meat abattoirs and genotypes from the C. burnetii 2014 Nijmegen database by geographical location show that the samples are closely related with human isolate 20090317Frankrijk004 from Marseille, France. The eight tissues from Gauteng province red meat abattoirs share the same MLVA-6 copy numbers as the isolate 20090317Frankrijk004.
Figure 4. Minimum spanning tree showing the relationship between the eight tissue samples from Gauteng tissue samples from Gauteng province red meat abattoirs and genotypes from the C. burnetii 2014 Nijmegen database by geographical location show that the samples are closely related with human isolate 20090317Frankrijk004 from Marseille, France. The eight tissues from Gauteng province red meat abattoirs share the same MLVA-6 copy numbers as the isolate 20090317Frankrijk004.
Pathogens 10 00258 g004
Figure 5. Phylogenic tree of C. burnetii genotypes of two sheep tissues (pooled organs) from two animals (40241 and 40242) originating a farm in the Mangaung district based on MLVA-6 database (http://mlva.u-psud.fr/mlvav4/genotyping/view.php (accessed on 14 January 2021)) created by Dr J.J.H.C. Tillburg. This genotype amplified for all MLVA-6 markers Nine Mile RSA493 with copy numbers 9-27-4-6-9-5 was used reference (data not shown).
Figure 5. Phylogenic tree of C. burnetii genotypes of two sheep tissues (pooled organs) from two animals (40241 and 40242) originating a farm in the Mangaung district based on MLVA-6 database (http://mlva.u-psud.fr/mlvav4/genotyping/view.php (accessed on 14 January 2021)) created by Dr J.J.H.C. Tillburg. This genotype amplified for all MLVA-6 markers Nine Mile RSA493 with copy numbers 9-27-4-6-9-5 was used reference (data not shown).
Pathogens 10 00258 g005
Figure 6. Minimum spanning tree showing the relationship between the two sheep tissue samples (40241 and 40242) from a farm in Mangaung district and genotypes from the C. burnetii 2014 Nijmegen database by geographical location show the isolates samples are closely related with human isolate 20090317Frankrijk044 from Marseille, France.
Figure 6. Minimum spanning tree showing the relationship between the two sheep tissue samples (40241 and 40242) from a farm in Mangaung district and genotypes from the C. burnetii 2014 Nijmegen database by geographical location show the isolates samples are closely related with human isolate 20090317Frankrijk044 from Marseille, France.
Pathogens 10 00258 g006
Table 1. Types of samples collected at red meat abattoirs in Gauteng province and tissue samples from a farm in Free State province used in the study.
Table 1. Types of samples collected at red meat abattoirs in Gauteng province and tissue samples from a farm in Free State province used in the study.
ProvinceSample TypeSample NameNumber (n)
GautengSerumSerum507
GautengTissuePenis355
GautengTissueTestes355
GautengTissueUterus145
GautengTissueMammary gland79
GautengTissueOvary80
GautengTissueOviduct4
Free statePooled tissuesSpleen, liver, lung2
Table 3. Seroprevalence of Q fever in livestock slaughtered at red meat abattoirs in Gauteng, South Africa, and univariate analysis of associated factors.
Table 3. Seroprevalence of Q fever in livestock slaughtered at red meat abattoirs in Gauteng, South Africa, and univariate analysis of associated factors.
VariableLevelPrevalence (%)95%CI*p-Value
SpeciesBovine9.46.5–13.00.003
Ovine4.30.9–12.0
Porcine0.90.002–5.0
SexMale4.82.8–7.60.007
Female11.87.2–18.1
BreedBonsmara6.63.8–10.5<0.001
Jersey6.91.4–19.0
Nguni26.714.7–42.0
Dorper4.30.9–12.0
Large white0.90.2–5.0
DistrictTshwane12.27.9–17.80.003
Ekurhuleni3.40.7–9.7
Metsweding101.2–32.0
Sedibeng4.31.7–8.7
West Rand00.0–7.1
Abattoir throughputHigh7.14.8–9.91
Low5.81.6–143.9
Origin of animalsFarm/feedlot4.82.9–7.3<0.001
Auction16.79.6–26.0
Total 6.94.9–9.5
CI* Confidence interval.
Table 4. Multivariable analysis of risk factors associated with seropositivity to C. burnetii in livestock at red meat abattoirs in Gauteng, South Africa.
Table 4. Multivariable analysis of risk factors associated with seropositivity to C. burnetii in livestock at red meat abattoirs in Gauteng, South Africa.
VariableLevelOdds Ratio95% CI*p-Value
SpeciesBovine1 *0.2–2.00.369
Ovine0.60.006–0.40.003
Porcine0.04
Abattoir throughputHigh1 *1.2–14.00.023
Low4.1
Origin of animalsFarm/feedlot1 *2.6–12.4<0.001
Auction5.7
CI* Confidence interval. * Reference level.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mangena, M.; Gcebe, N.; Pierneef, R.; Thompson, P.N.; Adesiyun, A.A. Q Fever: Seroprevalence, Risk Factors in Slaughter Livestock and Genotypes of Coxiella burnetii in South Africa. Pathogens 2021, 10, 258. https://doi.org/10.3390/pathogens10030258

AMA Style

Mangena M, Gcebe N, Pierneef R, Thompson PN, Adesiyun AA. Q Fever: Seroprevalence, Risk Factors in Slaughter Livestock and Genotypes of Coxiella burnetii in South Africa. Pathogens. 2021; 10(3):258. https://doi.org/10.3390/pathogens10030258

Chicago/Turabian Style

Mangena, Maruping, Nomakorinte Gcebe, Rian Pierneef, Peter N. Thompson, and Abiodun A. Adesiyun. 2021. "Q Fever: Seroprevalence, Risk Factors in Slaughter Livestock and Genotypes of Coxiella burnetii in South Africa" Pathogens 10, no. 3: 258. https://doi.org/10.3390/pathogens10030258

APA Style

Mangena, M., Gcebe, N., Pierneef, R., Thompson, P. N., & Adesiyun, A. A. (2021). Q Fever: Seroprevalence, Risk Factors in Slaughter Livestock and Genotypes of Coxiella burnetii in South Africa. Pathogens, 10(3), 258. https://doi.org/10.3390/pathogens10030258

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