Description: Cross section study estimating the current herd-level prevalence of C. burnetii in the Irish dairy herd and to identify associated herd-level risk factors using bulk tank milk ELISA results.
Authors
First Author (and Corresponding Author)
Katie Corridan
School of Public Health, Physiotherapy and Sports Science
University College Dublin
Bellfield
katie.corridan@ucd.ie
Second Author
Patrick Wall
Emeritus Professor
School of Public Health, Physiotherapy and Sports Science
University College Dublin
Bellfield
Dublin 4
patrick.wall@ucd.ie
Third Author
Guy McGrath
Centre for Veterinary Epidemiology and Risk Analysis
University College Dublin
Bellfield
Dublin 4
guy.mcgrath@ucd.ie
Fourth Author
Conor McAloon
Centre for Veterinary Epidemiology and Risk Analysis
University College Dublin
Bellfield
Dublin 4
conor.mcaloon@ucd.ie
Key words:
Coxiella burnetii
Q fever
Coxiellosis, herd-level prevalence
bulk tank milk
ELISA
zoonosis
List of abbreviations
BMT bulk milk tank
EBI Economic Breeding Index
ELISA enzyme-linked immunosorbent assay
IBR infectious bovine rhinotracheitis
IFA immunofluorescent assay
Title: Q Fever in the Irish dairy herd
Abstract
Introduction: Coxiella burnetii, the causative agent of Q fever, is a notifiable zoonotic pathogen in Ireland. While typically subclinical in ruminants, infection is associated with reproductive losses. In humans, disease can range from asymptomatic to more serious complications. Ruminants have been identified as the main reservoir for human infection. Ireland’s dairy industry has expanded substantially in recent years, yet current data on the national prevalence of Coxiella burnetii in dairy herds are limited. Understanding the herd-level prevalence and associated risk factors is essential for informing disease management and control strategies.
Methods: Bulk milk tank testing results from an Irish dairy cooperative herd health programme were analysed to determine the apparent and true prevalence of Coxiella burnetii antibodies, using 2022 data. Further analysis was conducted to determine the
relationship between Coxiella burnetii prevalence and co-morbid disease and herd
characteristics.
Results: 2,691 dairy herds were included in the sample. The true prevalence of Coxiella burnetii antibodies was 62.9%. Coxiella burnetii prevalence was associated with increasing herd size and replacement rate.
Conclusion: This study provides updated data, revealing the highest herd-level prevalence of Coxiella burnetii antibodies reported to date. The association with larger herd size is particularly relevant with the substantial growth in the national dairy herd over the past decade. These findings reinforce the need for further research into transmission dynamics, impact on production and zoonotic risk.
A
Background:
Coxiellosis or Q fever is a notifiable zoonotic disease caused by the bacterium Coxiella burnetii. Ruminants have been identified as the main reservoir for human infection (1). In animals, coxiellosis can pose a significant risk to herd health and loss of production through decreased fertility, abortions, and anorexia (2). In humans, Q fever most commonly presents as a self-limiting flu-like illness but may cause pneumonia and hepatitis in rare cases and has been associated with spontaneous abortion in pregnant women (3). Infection in humans is frequently subclinical, leading to a consensus in both medical and veterinary literature that Q fever is likely underdiagnosed and underreported (4, 5).
In Q-fever outbreaks, human transmission occurs through inhalation of aerosol dust particles (6). Q fever can also be an occupationally acquired disease transmitted to farmers, their families, staff, and veterinarians, through inhalation of contaminated particles, ingestion or direct contact with birth fluids or placenta. The organism is also shed in milk, urine and faeces. Factors such as the proximity and intensity of exposure to infected animals or their products, personal hygiene practices, and individual susceptibility further influence transmission risk (7).
To accurately assess the potential human health risk, it is essential to understand the prevalence of Coxiella burnetii in the cattle population. Previous studies have indicated that Q fever is endemic in the Irish dairy herd, with a reported herd-level prevalence of 19.5% from 1489 randomly selected bulk milk tank samples — slightly above the estimated European average of 15.1% (8). However, this data is now a decade old. Since then, average herd size in Ireland has increased significantly, and changes in herd management, biosecurity practices, and environmental factors may have influenced disease dynamics. Updated prevalence estimates are therefore needed to better assess current animal and public health risks. Bulk tank milk (BTM) sampling combined with enzyme-linked immunosorbent assay (ELISA) has been shown to be a cost-effective method for estimating national herd-level prevalence (9).
The objective of this study was to estimate the current herd-level prevalence of C. burnetii in the Irish dairy herd and to identify associated herd-level risk factors using bulk tank milk ELISA results. These findings will support veterinary and public health decision-making, particularly in the context of occupational and environmental exposure.
Methods
Programme/dataset description.
This study utilised BMT samples collected as part of an Irish dairy cooperative herd health programme. The study sample comprises all participants in an Irish dairy cooperative Herd Health Programme in 2022. By their nature, dairy cooperatives are regionalised. Herd owners subscribe to the Herd Health Programme annually for a subscription fee. The programme is open to all Irish milk producers, regardless of their co-operative affiliation.
Sample testing.
A
The samples used for BMT testing are taken concurrently with routine milk constituent payment testing. BMT samples were collected aseptically from each herd and stored at 4°C until testing, which occurred within 72 hours. BMT samples were tested using the ID Screen Q Fever Indirect Multi-species (Innovative Diagnostics, Grabels), an indirect ELISA designed to detect antibodies against
Coxiella burnetii in milk samples. Testing was performed according to the manufacturer’s protocol.
Data processing
Test results were dichotomised to positive or negative test status based on their signal-to-positive ratio (S/P ratio) according to the manufacturer’s guidelines. Herd level variables (e.g. herd size, herd average economic breeding index (EBI), percentage of herd at first lactation) were extracted from the Irish Cattle Breeding Federation for the year 2022 (10).
All data were anonymised and linked to unique herd identification numbers before transfer to the researchers.
Data analysis
A herd was deemed positive if any one of the BTM samples in that calendar year tested positive. Apparent herd-level prevalence was calculated as the proportion of positive herds among all tested herds in that year. Confidence intervals were estimated using the Clopper-Pearson method (11). Herd-level true prevalence was estimated using the Rogan-Gladen estimator, with confidence intervals estimated using using the delta method (12). For these calculations, test sensitivity and specificity values of 0.86 and 0.99, respectively, were used based on previous evaluation of the ID Screen ELISA Q Fever test by Paul et al. (13), who estimated these parameters using bulk tank milk samples.
Herd-level risk factors for Coxiella burnetii antibody positivity were identified using logistic regression with the herd level status as the dependent variable. Independent variables were herd size, average herd Economic Breeding Index (EBI), percentage of herd at first lactation, exposure to infectious bovine rhinotracheitis (IBR), salmonellosis (Salmonella enterica), leptospirosis (Leptospira Hardjo), Neospora caninum, Ostertagia ostertagi (stomach worms), and Fasciola hepatica (liver fluke). Each variable was first trialled in a univariate analyses, and variables with P < 0.20 were subsequently offered to a multivariable model. A backward stepwise selection approach was used to remove variables in order of P value, with variables with the largest p-value removed first. Variables with a P < 0.05 were retained in the final model. All statistical analysis was conducted in R studio version 4.3.0. (R Core Team (14).
Ethics exemption was granted from UCD Taught Masters Research Ethics Committee reference TMREC-22-108-CORRIDAN-WALL.
Results
Descriptive statistics
The final dataset for analysis consisted of 2,691 herds comprising a total of 227,706 individual cows. Descriptive statistics for herd characteristics included in the sample and the national average as reported by the Irish Cattle Breeding Federation can be seen in Table 1. The median herd size for this sample was 116 lactating females. The median herd percentage at first lactation was 20%, with a range of 0% to 100%. The median herd average EBI was €151, with a range of -€82 to €239.
Table 1
Descriptive statistics for herd characteristics
| | Study Median | Study Range | National Average |
|---|
Herd Size | 116 | 1-1175 | 107 |
% Herd at first lactation | 20% | 0–20% | 19% |
EBI | €151 | -€82 - +€239 | €151 |
None of the herds included in the sample had been vaccinated for Q fever, but there was a high level of vaccination across the population sample for other diseases. 1,867 herds (69.4%) were vaccinated for IBR, 1877 herds (69.6%) were vaccinated for of leptospirosis, and 1540 herds (57.2%) were vaccinated for Salmonella. The majority of herds also demonstrated a low parasitic burden. 85.3% of the sample was negative for Liver Fluke and 78% of herds returned either a negative or low positive result for stomach worms.
Prevalence estimates
Of the 2,691 herds tested, 1,465 (54.4%; 95% CI: 52.5%–56.3%) tested positive on bulk tank milk. After adjusting for test sensitivity (86%) and specificity (99%) using the Rogan-Gladen estimator, the true herd-level prevalence was estimated at 62.9% (95% CI: 60.7%–65.1%).
A
There was no significant association found between testing positive for
Coxiella burnetii on BMT sampling and stomach worm burden (P = 0.58). As such, this variable was not included in the multivariate analysis.
The variables which were deemed significant in the above univariate analysis were then carried through to a multivariable analysis for backward stepwise regression. Only one variable, neospora, was removed during this process due to lack of statistical significance (P = 0.88). The results of the multivariable analysis are included in Table 2.
Table 2: Multivariable Logistic Regression Analysis |
|---|
Variable | Odds Ratio | 95% confidence interval | P value |
IBR Vaccinated | Referent | | |
IBR Unvaccinated Negative | 1.024 | (0.810–1.295) | 0.844 |
IBR Unvaccinated Positive | 0.722 | (0.57–0.914) | 0.01 |
Leptospirosis Vaccinated | Referent | | |
Leptospirosis Unvaccinated Negative | 0.718 | (0.546–0.944) | 0.018 |
Leptospirosis Low Positive | 1.438 | (0.845–2.447) | 0.177 |
Leptospirosis Medium Positive | 0.813 | (0.582–1.136) | 0.2222 |
Leptospirosis High Positive | 1.213 | (0.891–1.652) | 0.2130 |
Salmonella Vaccinated | Referent | | |
Salmonella Unvaccinated Negative | 0.817 | (0.665–1.004) | 0.0481 |
Salmonella Unvaccinated Positive | 1.276 | (0.981–1.660) | 0.0688 |
Liver Fluke | 0.690 | (0.548–0.868) | 0.0015 |
Herd size | | (1.005–1.008) | p < 0.001 |
EBI | 0.995 | (0.993–0.997) | p < 0.001 |
First Lactation | 0.259 | (0.110–0.612) | 0.0020 |
In the multivariable logistic regression model, several significant associations were identified. Herds which were vaccinated for IBR had significantly (P = 0.01) higher odds of testing positive for Coxiella burnetii on bulk tank sampling than herds which were unvaccinated testing positive for IBR antibodies. There was no significant difference between herds which were vaccinated for IBR and those which were unvaccinated negative. Herds which were vaccinated for leptospirosis had significantly (P = 0.018) higher odds of testing positive for Coxiella burnetii antibodies in bulk tank samples than herds who were unvaccinated negative for leptospirosis. There was no significant difference between herds who were vaccinated for leptospirosis and those who were unvaccinated testing positive. Herds which were vaccinated for salmonella had significantly (P = 0.048) lower odds of testing positive for Coxiella burnetii antibodies in bulk milk sampling than those herds who were unvaccinated negative. There was no significant difference noted between herds who were vaccinated for salmonella and those who were unvaccinated testing positive for salmonella antibodies.
Herds which were positive for liver fluke had significantly (P = 0.0015) lower odds of testing positive for Coxiella burnetii antibodies in bulk tank samples.
Increasing herd size had significantly (P < 0.001) increased the odds of testing positive for Coxiella burnetii on bulk tank sampling. Increasing herd average EBI (P < 0.001) and proportion of the herd in first lactation (P = 0.0020) significantly decreased the odds of testing positive for Coxiella burnetii on bulk tank sampling.
Discussion
The large study sample population of 2,691 dairy herds included in this study represent nearly 15% of the total number of dairy herds in Ireland as per the Irish Cattle Breeding Federation. The sample population was relatively well matched with the national average in terms of, herd size, replacement rate and genetic merit (10).
In this study the apparent prevalence for Coxiella burnetii was 54.4% rising to 62.9% after adjusting for ELISA sensitivity and specificity (13). Earlier Irish bulk-tank ELISA studies reported a true prevalence of 19.5% in 1,489 randomly selected dairy herds (8) and an apparent prevalence of 37.9% in 290 self-selected dairy herds (9).
Coxiella burnetii is enzootic worldwide, yet estimations of prevalence among both human and bovine populations vary widely in the literature. Although Coxiella burnetii is a notifiable disease in 29 EU countries, surveillance is passive meaning most data is case based (6).
The EU average for Coxiella burnetii antibody prevalence among all cattle is reported to be 15.1% (6). Ryan et al. (9) had reported a significantly higher risk of Coxiella burnetii among dairy cattle than beef cattle which could account, at least in part, for the higher prevalence noted in this study. Pandit et al. (15) reported a prevalence of 69.3% using ELISA BMT testing in Western France. Wood et al. (16) used immunofluorescent assay testing (IFA) to estimate a seroprevalence of 5.2% among extensively managed beef cattle in Australia. Evidently, prevalence estimations differ widely depending on geographical location and farm management system.
Increasing herd size was found to be associated with significantly (p < 0.0001) higher odds of BMT Coxiella burnetii antibody positivity, echoing prior Irish studies (8, 9). This is particularly pertinent for the Irish context as the Irish dairy industry has undergone a period of significant growth and expansion since the abolition of EU milk quotas in 2015. During this period, the Irish national dairy herd has increased by 25%, demonstrating an increase each year. Furthermore, in the same period the Irish suckler herd (breeding beef cows) population has declined by over 18%. With previous research indicating that dairy cows are significantly more likely to test positive for Coxiella burnetii, this highlights the potential increasing prevalence of Q fever among Irish dairy cattle due to demographic shifts.
Herds that were unvaccinated and tested positive for IBR had significantly lower odds of Q fever positivity compared to vaccinated herds. Conversely, herds which were unvaccinated and tested negative for leptospirosis or salmonella, had significantly lower odds of Q fever positivity compared to vaccinated herds. Given the observational nature of the data, it is not possible to determine causality between vaccination status and Q fever risk. These opposing patterns likely reflect differences in herd management—such as greater animal movement for replacement stock or contract rearing—rather than direct vaccine effects. The cross-sectional design of this study cannot infer causality and that apparent associations may be due to confounding factors such as broader management strategies, biosecurity or co-infections.
Many of the large human outbreaks of Q fever have been traced to small ruminant farms such as the Netherlands 2007–2010 outbreak (17). Ireland has a thriving dairy industry, but minimal dairy goat or dairy sheep farms. The differing ruminant population patterns between countries may confer differing Q fever risks. Whole Genome Sequencing is an interesting development in the field of Q fever dynamics and zoonotic disease overall. Genotyping research from the United States identified a single genotype, ST20, in bovine milk samples and different genotype, ST8, in caprine milk (18), bringing the possibility of subspecies adaptations into question. Given the most recent large human outbreaks have been linked to a small ruminant source, perhaps the genotype associated with bovine Q fever has more limited transmissibility to humans but to explore this further genotyping studies are required in human cases. Evidently this issue is not limited to Q fever, a One Health pathogen genotype database would enable rapid and accurate source detection enabling faster implementation of outbreak control measures as well as prevention measures.
The high prevalence of Q fever in dairy herds warrants further investigation. While clinical cases in individual cattle are usually treated by vets, many animals remain asymptomatic, making BMT testing valuable for detecting hidden infections. Shedding of bacteria can occur without obvious symptoms, posing ongoing zoonotic risks. For example, Ladbury et al. (19) linked a large human outbreak in the Netherlands to an infected goat farm that showed no animal health issues. A Teagasc survey found 92% of Irish farmers were unaware of the zoonotic risk from seemingly healthy animals, highlighting the need for farmer education (20).
The occupational risk of Q fever is well documented with transmission via birthing fluids, urine and faeces. Given the high dairy herd prevalence, human Q fever may be underdiagnosed and underreported, particularly among farmers and vets. Critically, vulnerable groups, like pregnant women or those with vascular conditions, should avoid animals during calving season (21). Hygiene measures and personal protective equipment during calving are critical.
The large sample population included in this study is a significant strength. The high prevalence detected in this study may prompt physicians to include Q fever as a differential diagnosis when patients are presenting with non-specific symptoms, potentially revealing a greater burden of disease than currently recognised. The high prevalence reported in this study also needs to be highlighted to farmers, farm workers and their families so that hygiene and mitigation practices can protect those at risk.
A possible limitation of this study is the regional scope of the testing data, because of this, the results may not fully reflect national-level prevalence or risk factors. Bulk milk tank testing detects herd exposure rather than active shedding so it may overestimate infection risk. In addition, sampling bias may exist as herds included in this study are paid participants in a herd health programme whose management practices may differ from those who chose not to enrol in such a programme. Furthermore, herds that have enrolled will have received specialist veterinary advice to maximise herd health and productivity so that perhaps vaccination or dosing regimens may differ for herds enrolled in the programme.
Future research should include a spatiotemporal analysis of human Q fever notifications in Ireland and compare with the bovine prevalence and distribution. Inclusion of small ruminant data may further quantify the risk posed by Q fever in Ireland. The contrasting vaccination patterns between herds warrant further investigation to elucidate their impact on C. burnetii infection dynamics. Finally, genotypic analyses of isolates from animal and human sources would enhance understanding of transmission pathways and inform One Health strategies.
Conclusion
A
In conclusion, this study revealed a high herd-level prevalence of
Coxiella burnetii, with larger herd size and certain vaccination practices emerging as significant risk factors. Bulk milk tank testing proved an efficient, cost-effective method. Future work should integrate human Q fever case notifications, include small ruminant populations, and apply genotypic analyses to clarify transmission pathways. Strengthened collaboration across veterinary and human health sectors, supported by ongoing research, is essential to guide targeted interventions and reduce zoonotic Q fever risk.
Declarations
Ethics
Ethics exemption was granted from UCD Taught Masters Research Ethics Committee reference TMREC-22-108-CORRIDAN-WALL.
A
Data Availability
The datasets generated and/or analysed during the current study are not publicly available but are available from the corresponding author on reasonable request.
A
Author Contribution
All authors contributed to the conception and design of the study. All authors contributed to the interpretation of data and the preparation of the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors would like to thank an industry partner for providing access to data used in this study, in accordance with their request to remain anonymous. This research was conducted as part of a Master of Public Health degree at University College Dublin.
References
1.Woldehiwet Z. Q fever (coxiellosis): epidemiology and pathogenesis. Res Vet Sci. 2004;77(2):93–100.
2.Rodolakis A. Q fever in dairy animals. Ann N Y Acad Sci. 2009;1166(1):90–3.
3.Quijada SG, Terán BM, Murias PS, Anitua AA, Cermeño JLB, Frías AB. Q fever and spontaneous abortion. Clin Microbiol Infect. 2012;18(6):533–8.
4.Kaufman HW, Chen Z, Radcliff J, Batterman HJ, Leake J. Q fever: an under-reported reportable communicable disease. Epidemiol Infect. 2018;146(10):1240–4. 10.1017/S0950268818001395.
5.van Loenhout JAJ, Wielders CCH, Morroy G, Cox MJ, van der Hoek W, Hautvast JLA, et al. Severely impaired health status of non-notified Q fever patients leads to underestimation of true burden of disease. Epidemiol Infect. 2015;143(12):2580–7.
6.European Centre for Disease Prevention and Control. Q fever annual epidemiological report for 2019 [Internet]. Stockholm: ECDC. 2021 [cited 2023 Feb 27]. Available from: https://www.ecdc.europa.eu/sites/default/files/documents/Q-fever-annual-epidemiological-report-for-2019.pdf
7.Guattéo R, Czaplicki G, Mainil J, Saegerman C, Porter SR. Q Fever: Current State of Knowledge and Perspectives of Research of a Neglected Zoonosis. Int J Microbiol. 2011;2011:92–113.
8.Ryan ED, Wrigley K, Hallinan A, McGrath G, Clegg TA. Antibodies to Coxiella burnetii in Irish bulk tank milk samples. Vet Rec. 2018;182(19):550.
9.Ryan ED, Kirby M, Collins DM, Sayers R, Mee JF, Clegg T. Prevalence of Coxiella burnetii (Q fever) antibodies in bovine serum and bulk-milk samples. Epidemiol Infect. 2011;139:1413–7.
10.Irish Cattle Breeding Federation. Dairy statistics [Internet]. ICBF: 2022 [cited 2023 Feb 27]. Available from: https://www.icbf.com/dairy-statistics/
11.Clopper-Pearson method. In. A Dictionary of Statistics. Oxford: Oxford University Press; 2014.
12.Rogan WJ, Gladen B. Estimating prevalence from the results of a screening test. Am J Epidemiol. 1978;107(1):71–6. 10.1093/oxfordjournals.aje.a112510.
13.Paul S, Toft N, Agerholm JS, Christoffersen A, Agger JF. Bayesian estimation of sensitivity and specificity of Coxiella burnetii antibody ELISA tests in bovine blood and milk. Prev Vet Med. 2013;109(3–4):258–63.
14.R Core Team. R: A language and environment for statistical computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2023 [cited 2025 Oct 13]. Available from: https://www.R-project.org/
15.Pandit P, Hoch T, Ezanno P, Beaudeau F, Vergu E. Spread of Coxiella burnetii between dairy cattle herds in an enzootic region: modelling contributions of airborne transmission and trade. Vet Res. 2016;47(1):1–16.
16.Wood CM, Perkins NR, Tozer SJ, Johnson W, Barnes TS, McGowan M, et al. Prevalence and spatial distribution of Coxiella burnetii seropositivity in northern Australian beef cattle adjusted for diagnostic test uncertainty. Prev Vet Med. 2021;189:105282. 10.1016/j.prevetmed.2021.105282.
17.van Leuken JP, Havelaar AH, van der Hoek W, Ladbury GA, Hackert VH, Swart AN, et al. A model for the early identification of sources of airborne pathogens in an outdoor environment. PLoS ONE. 2013;8(12):e80412.
18.Pearson T, Hornstra HM, Hilsabeck R, Gates LT, Olivas SM, Birdsell DM, et al. High prevalence and two dominant host-specific genotypes of Coxiella burnetii in U.S. milk. BMC Microbiol. 2014;14:41. 10.1186/1471-2180-14-41.
19.Ladbury GA, van Leuken JP, Swart A, Vellema P, Schimmer B, ter Schegget R, et al. Integrating interdisciplinary methodologies for One Health: goat farm re-implicated as the probable source of an urban Q fever outbreak, the Netherlands, 2009. BMC Infect Dis. 2015;15:1–11.
20.Mahon MM, Sheehan MC, Kelleher PF, Johnson AJ, Doyle SM. An assessment of Irish farmers’ knowledge of the risk of spread of infection from animals to humans and their transmission prevention practices. Epidemiol Infect. 2017;145(12):2424–35.
21.Parker NR, Barralet JH, Bell AM. Q fever Lancet. 2006;367(9511):679–88.