Genotyping of Coxiella burnetiifrom domestic ruminants and human in Hungary: indication of various genotypes
BMC Veterinary Research volume 10, Article number: 107 (2014)
Information about the genotypic characteristic of Coxiella burnetii from Hungary is lacking. The aim of this study is to describe the genetic diversity of C. burnetii in Hungary and compare genotypes with those found elsewhere. A total of 12 samples: (cattle, n = 6, sheep, n = 5 and human, n = 1) collected from across Hungary were studied by a 10-loci multispacer sequence typing (MST) and 6-loci multiple-locus variable-number of tandem repeat analysis (MLVA). Phylogenetic relationships among MST genotypes show how these Hungarian samples are related to others collected around the world.
Three MST genotypes were identified: sequence type (ST) 20 has also been identified in ruminants from other European countries and the USA, ST28 was previously identified in Kazakhstan, and the proposed ST37 is novel. All MST genotypes yielded different MLVA genotypes and three different MLVA genotypes were identified within ST20 samples alone. Two novel MLVA types 0-9-5-5-6-2 (AG) and 0-8-4-5-6-2 (AF) (Ms23-Ms24-Ms27-Ms28-Ms33-Ms34) were defined in the ovine materials correlated with ST28 and ST37. Samples from different parts of the phylogenetic tree were associated with different hosts, suggesting host-specific adaptations.
Even with the limited number of samples analysed, this study revealed high genetic diversity among C. burnetii in Hungary. Understanding the background genetic diversity will be essential in identifying and controlling outbreaks.
The aetiological agent of Q (Query) fever is Coxiella burnetii, a Gram-negative, obligate intracellular, zoonotic bacterium which is widely distributed throughout the world. Domestic ruminants (cattle, sheep and goats) have the main role in the disease cycle, but C. burnetii has been isolated from other domestic and wild animals too . Domestic ruminants are reservoirs of the agent with usually subclinical manifestation of the disease, but may suffer from reproductive disorders, and abortion and stillbirth can occur. The bacteria are transmitted mainly by inhalation. Additionally, infections may occur through the consumption of raw milk and milk products . In humans, Q fever is typically an acute febrile illness with nonspecific clinical signs, such as atypical pneumonia and hepatitis in roughly 40% of the cases while the 60% of the individuals remain asymptomatic after infection . A small percentage (~5%) of infected people may develop chronic infection with life-threatening valvular endocarditis [4, 5].
Q fever is a notifiable disease in Hungary. Antibodies against C. burnetii were first detected in the sera of abattoir workers in 1950 , and infections were first diagnosed in 1956 in dairy and sheep farms . According to recently published data, seroprevalence among cattle and sheep in Hungary were 38.0% and 6.0% with enzyme-linked immunosorbent assay respectively, which correspond with the European averages [8, 9]. The number of yearly reported acute human infections in Hungary ranged between 36 and 68 during the past five years. Major outbreaks were last registered in the period of 1976–1980 and in 2013 in Hungary (unpublished data).
Genetic characterisation of C. burnetii is required for epidemiological investigations in Q fever outbreaks and for surveillance purposes. Several typing systems exist, including pulsed-field gel electrophoresis , sequence analysis or restriction fragment length polymorphism of single genes (16S ribosomal RNA, icd, com1 or mucZ) [11–15]. The main problems with these methods are that their discriminatory power is poor and their reproducibility and transferability is not always straightforward . Recently, two PCR-based typing methods have come into view: multispacer sequence typing (MST)  and multi-locus variable-number tandem repeat analysis (MLVA) [18, 19]. Both methods have high discriminatory power and are reproducible, allowing interlaboratory comparisons without the need for isolation and cultivation of the organism under biosafety level three conditions.
Multiple sequence types have been identified from the Central European region in ticks and human samples: ST2 from Ukraine, ST16, ST18 from Romania and Slovakia, ST23 from Czech Republic and Slovakia and ST22, ST29 from Slovakia. While ST27 and ST32 were detected in goat samples in Austria [17, 20]. MLVA types L, P and partial MLVA profiles were identified from cow milk samples originating from Slovakia, Austria and Croatia .
The aim of this study was to identify C. burnetii MST and MLVA genotypes occurring in Hungary and to compare them with types from foreign countries.
This study included five ovine and five cattle abortion samples (cotyledons) and a cattle milk sample collected from different parts of Hungary (Table 1, Figure 1). The blood sample of a 72 years-old man suffering from acute Q fever (1:128 IgGII titre with micro-immunofluorescence test, Focus Diagnostics, Cypress, CA) was also characterized. The samples were not related to any outbreak, they were collected through routine diagnostic examinations. All tests were performed in accordance with all applicable institutional and national guidelines and regulations, approved by the ethics committees of the Institute for Veterinary Medical Research (animal) and the National Center for Epidemiology (human) and with the consent of the patient.
Extraction of DNA from the samples was performed with the QIAamp Tissue and Blood kit (Qiagen GmbH, Hilden, Germany) according to the manufacturer’s instructions. To quantify the approximate bacterial load and to determine the cycle threshold (Ct) values, all samples were analysed by TaqMan based real-time polymerase chain reaction (PCR) assay targeting the multicopy IS1111 insertion element of the C. burnetii genome .
MST consisted of amplification followed by sequencing of ten different spacer regions of the C. burnetii genome: Cox2, 5, 18, 20, 22, 37, 51, 56, 57 and 61 . Primer sequences and reaction conditions were described earlier . The compiled sequences of the PCR products were analysed with the BioEdit Sequence Alignment Editor 7.1.11 . Sequence types (STs) were determined using the MST database and previous publications [17, 20, 25]. Unfortunately, the allele codes for ST35 and ST36 are not publicly available therefore we were not able to include them in the phylogenetic analysis . To determine the phylogenetic relationships of our samples to known MST genotypes [17, 20] we used the sequence information of our three STs and constructed a phylogenetic tree with the 112 polymorphisms and methods described in Hornstra et al. . As samples VS16/27/38/42 (sheep abortion) appeared to have a new ST and homoplasies can alter topologies of trees, we also constructed a sub-tree to verify its placement using the same methods and 112 polymorphisms as above but only with STs: 8, 9, 10, 27, 28, 31, and the two sheep abortion STs from this study.
MLVA was performed by single PCRs targeting six variable microsatellite markers as performed in other recent studies [18, 27–30]. Ms27, Ms28 and Ms34 contain repeat units of six base pairs and Ms23, Ms24 and Ms33 contain repeat units of seven base pairs. The 3′-end-labelled forward and reverse primer sequences, and PCR conditions were applied as described before [27–30]. The amplification products were run on an ABI 3100 Genetic Analyser and electropherograms were evaluated with the Peak Scanner Software 2.0 (Applied Biosystems Inc., Foster City, CA). DNA of the Nine Mile strain (RSA 493, Coxevac, Ceva Inc., Budapest, Hungary) was used as a reference. The repeat numbers of each marker were determined by extrapolation using the obtained length of the sample fragments relative to the obtained fragment-length of the reference strain. A modified coding convention was used for Ms33: nine imperfect repeat units were counted in the Nine Mile strain (G. Vergnaud, personal communication). Novel MLVA types were determined if the described allele combinations had not been previously described in any publications applying the same MLVA method [21, 27–29].
All 12 samples were positive with the Coxiella specific real-time PCR with Ct values below 33 and all but one yielded complete MST and MLVA genotypes (Table 1). MST alignment and MLVA profiles were deposited in Dryad (, doi:10.5061/dryad.rc8q4) and in the Coxiella MLVA database (, accession numbers: Hung00-Hung11).
The MST analysis of the six cattle and five ovine samples revealed a novel (proposed ST37) and two previously described (ST20, ST28) STs (Table 1, Figure 2). The allele codes of the newly identified genotype ST37 are 5-4-9-5-8-5-2-3-4-6 for spacers Cox2-Cox5-Cox18-Cox20-Cox22-Cox37-Cox51-Cox56-Cox57-Cox61, respectively. Phylogenetic analysis according to Honstra et al.  resulted in 156 equally parsimonious tree topologies, one of which is shown in Figure 2, with a homoplasy index (excluding uninformative characters) of 0.0769. Sub-tree construction using STs 8, 9, 10, 27, 28, 33, and the two sheep abortion STs from this study resulted in a single tree with no homoplasy whose topology exactly matched the topology for this clade in the tree selected in Figure 2. In clonal groups of organisms with little genetic diversity the homoplasy index has been shown to be the most appropriate measure of phylogenetic accuracy [33, 34], therefore the results of this sub-tree provide confidence in the accuracy of the placement of this novel ST (proposed ST37) in Figure 2.
Five MLVA genotypes were identified among the 11 cattle and ovine samples with complete genotypes, out of which two novel MLVA types were described in this study. Three MLVA types (I, J, M) were detected among the six cattle samples, and the two novel genotypes (AF and AG) were defined from the five ovine samples. In the case of the human sample only a partial genotype was obtained despite the relatively high DNA load (Ct value was 28.31). The identified allele combination in the human sample did not match any previously defined genotypes, implying the discovery of a novel MLVA type as well. The partial MLVA genotype from the human sample showed high similarity to the novel genotype AG defined in sheep, as four out of the five microsatellites were the same.
All samples originating from cattle were identified as ST20, regardless of their MLVA types. Samples from sheep belonged to two MST types: ST28 was detected in a sample with AG MLVA genotype, and the proposed ST37 was described in samples with AF MLVA type. Host specificity of the genotypes was observed as all MLVA genotypes from ST20 were described from cattle only and the two other genotypes were detected in sheep exclusively.
Molecular characterisation of C. burnetii is a useful tool to explore the genotypic diversity in a region and to determine relationships between variants of this bacterium. The main advantage of both MST and MLVA typing methods is their high discriminatory power and that the cultivation of the bacteria is unnecessary for their application.
MST and MLVA genotyping methods revealed genotypic diversity among samples from domestic ruminants and human in Hungary. Three C. burnetii STs (ST20, ST28 and ST37) are present in Hungarian ruminants.
Indicating the broad range of ST20 occurrence, other samples have been reported from The Netherlands, Germany, Spain, in human clinical samples from France and in cow’s milk, soil, and goat placental material from the United States [17, 26, 29, 35, 39]. Microvariants of the ST20 genotype occur as well and have been found elsewhere; MLVA genotypes I, J and M were described before from cow milk and milk products from several countries in Europe (France, the Netherlands, Portugal, Spain and Switzerland) and other parts of the world (e.g. Qatar, Saudi Arabia) . ST20 was also detected in an outbreak among goats recorded from the United Kingdom . In conclusion, although ST20 is most often associated with cattle products, this genotype occasionally infects other species as well.
ST28 samples have been collected from cows, sheep, Hyalomma sp. tick and human blood in Kazakhstan . ST37 genotype is closely related to ST27-28, 31 and 8–10 genotypes that have been collected from humans in Austria, Kazakhstan, Uzbekistan, France and Portugal as well as sheep and goats from Spain, France, Canada and the USA (Figure 2) [17, 20, 28, 29, 39]. The comparison of the novel genotypes (AG and AF) with others described in Europe revealed similarities with genotypes AA and T, differing at 2 and 3 loci respectively and also associated with sheep, goats and human [28, 29]. Again, that these genotypes are most commonly associated with sheep and goats provides further evidence for host preference . We speculate that there may be some association between the Hungarian human and ovine infections (MLVA genotype AG) given genotype similarities and as the geographical distance between the two samples was only 50 km.
No matches were found among the MST and MLVA profiles from Hungary and from the dominant genotype (ST33) associated with the Q fever outbreak of The Netherlands that emerged in 2007 [27, 35]. STs described earlier from other Central European countries differed also from the STs found in Hungary [17, 20].
The present study provides information about the genotypic diversity of C. burnetii occurring in Hungary. From the examined 12 C. burnetii samples originating from cattle, sheep and human three MST types (containing one novel profile) and five MLVA types (including two novel profiles) have been identified. These genotypes cluster with host, suggesting host-specific evolutionary adaptations as shown by Pearson et al.  as well. The additional analyses of samples from a broader range of livestock species are needed to confirm this finding. Our results provide data for the European surveillance and help to increase information about Q fever in this Central European region.
Multispacer sequence typing
Multiple-locus variable number tandem repeat analysis.
Angelakis E, Raoult D: Emergence of Q fever. Iran J Publ Health. 2011, 40: 1-18.
Angelakis E, Raoult D: Q fever. Vet Microbiol. 2010, 140: 297-309. 10.1016/j.vetmic.2009.07.016.
Raoult D, Marrie T, Mege J: Natural history and pathophysiology of Q fever. Lancet Infect Dis. 2005, 5: 219-226. 10.1016/S1473-3099(05)70052-9.
Tilburg JJ, Melchers WJ, Pettersson AM, Rossen JW, Hermans MH, van Hannen EJ, Nabuurs-Franssen MH, de Vries M, Horrevorts AM, Klaassen CHW: Interlaboratory evaluation of different extraction and real-time PCR methods for detection of Coxiella burnetii DNA in serum. J Clin Microbiol. 2010, 48: 3923-3927. 10.1128/JCM.01006-10.
Landais C, Fenollar F, Thuny F, Raoult D: From acute Q fever to endocarditis: serological follow-up strategy. Clin Infect Dis. 2007, 44: 1337-1340. 10.1086/515401.
Farkas E, Gerő S, Takátsy G: A Q láz előfordulása Magyarországon. Orv Hetil. 1950, 91: 717-718.
Romváry J: Incidence of Q-fever in a cattle herd. Magy Állatorv Lapja. 1957, 12: 25-27.
Gyuranecz M, Dénes B, Hornok S, Kovács P, Horváth G, Jurkovich V, Varga T, Hajtós I, Szabó R, Magyar T, Vass N, Hofmann-Lehmann R, Erdélyi K, Bhide M, Dán Á: Prevalence of Coxiella burnetii in Hungary: screening of dairy cows, sheep, commercial milk samples, and ticks. Vector Borne Zoonotic Dis. 2012, 12: 650-653. 10.1089/vbz.2011.0953.
Guatteo R, Seegers H, Taurel AF, Joly A, Beaudeau F: Prevalence of Coxiella burnetii infection in domestic ruminants: a critical review. Vet Microbiol. 2011, 149: 1-16. 10.1016/j.vetmic.2010.10.007.
Jäger C, Willems H, Thiele D, Baljer G: Molecular characterization of Coxiella burnetii isolates. Epidemiol Infect. 1998, 120: 157-164. 10.1017/S0950268897008510.
Hendrix L, Samuel J, Mallavia L: Differentiation of Coxiella burnetii isolates by analysis of restriction-endonuclease-digested DNA separated by SDS-PAGE. J Gen Microbiol. 1991, 137: 269-276. 10.1099/00221287-137-2-269.
Stein A, Saunders NA, Taylor AG, Raoult D: Phylogenic homogeneity of Coxiella burnetii strains as determinated by 16S ribosomal RNA sequencing. FEMS Microbiol Lett. 1993, 113: 339-344. 10.1111/j.1574-6968.1993.tb06537.x.
Nguyen SV, Hirai K: Differentiation of Coxiella burnetii isolates by sequence determination and PCR-restriction fragment length polymorphism analysis of isocitrate dehydrogenase gene. FEMS Microbiol Lett. 1999, 180: 249-254. 10.1111/j.1574-6968.1999.tb08803.x.
Zhang GQ, To H, Yamaguchi T, Fukushi H, Hirai K: Differentiation of Coxiella burnetii by sequence analysis of the gene (com1) encoding a 27-kDa outer membrane protein. Microbiol Immunol. 1997, 41: 871-877. 10.1111/j.1348-0421.1997.tb01943.x.
Sekeyová Z, Roux V, Raoult D: Intraspecies diversity of Coxiella burnetii as revealed by com1 and mucZ sequence comparison. FEMS Microbiol Lett. 1999, 180: 61-67. 10.1016/S0378-1097(99)00462-0.
Massung RF, Cutler SJ, Frangoulidis D: Molecular typing of Coxiella burnetii (Q fever). Coxiella burnetii: Recent advances and new perspectives in research of the Q fever bacterium. Volume 984. Edited by: Toman R, Heinzen RA, Samuel JE, Mege J-L. Dordrecht: Springer Science + Business Media; 2012:381-396.
Glazunova O, Roux V, Freylikman O, Sekeyova Z, Fournous G, Tyczka J, Tokarevich N, Kovacova E, Marrie TJ, Raoult D: Coxiella burnetii genotyping. Emerg Infect Dis. 2005, 11: 1211-1217.
Arricau-Bouvery N, Hauck Y, Bejaoui A, Frangoulidis D, Bodier CC, 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-10.1186/1471-2180-6-38.
Svraka S, Toman R, Skultety L, Slaba K, Homan WL: Estabilishment of a genotyping scheme for Coxiella burnetii. FEMS Microbiol Lett. 2006, 254: 268-274. 10.1111/j.1574-6968.2005.00036.x.
Multi Spacers Typing - Coxiella Database. [http://ifr48.timone.univ-mrs.fr/MST_Coxiella/mst/]
Tilburg JJ, Roest HJ, Nabuurs-Franssen MH, Horrevorts AM, Klaassen CH: Genotyping reveals the presence of a predominant genotype of Coxiella burnetii in consumer milk products. J Clin Microbiol. 2012, 50: 2156-2158. 10.1128/JCM.06831-11.
Loftis AD, Reeves WK, Szumlas DE, Abbassy MM, Helmy IM, Moriarity JR, Dasch GA: Surveillance of Egyptian fleas for agents of public health significance: Anaplasma, Bartonella, Coxiella, Ehrlichia, Rickettsia and Yersinia pestis. Am J Trop Med Hyg. 2006, 75: 41-48.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. 1999, 41: 95-98.
Mediannikov O, Fenollar F, Socolovschi C, Diatta G, Bassene H, Molez J-F, Sokhna C, Trape J-F, Raoult D: Coxiella burnetii in humans and ticks in rural Senegal. PLoS Negl Trop Dis. 2010, 4: e654-10.1371/journal.pntd.0000654.
Hornstra HM, Priestley RA, Georgia SM, Kachur S, Birdsell DN, Hilsabeck R, Gates LT, Samuel JE, Heinzen RA, Kersh GJ, Keim P, Massung RF, Pearson T: Rapid Typing of Coxiella burnetii. PLoS ONE. 2011, 6: e26201-10.1371/journal.pone.0026201.
Tilburg JJ, Rossen JWA, van Hannen EJ, Melchers WJ, Hermans MH, van de Bovenkamp J, Roest HJ, de Bruin A, Nabuurs-Franssen MH, Horrevorts AM, Klaassen CH: Genotypic diversity of Coxiella burnetii in the 2007–2010 Q fever outbreak episodes in The Netherlands. J Clin Microbiol. 2012, 50: 1076-1078. 10.1128/JCM.05497-11.
Santos AS, Tilburg JJ, Botelho A, Barahona MJ, Nuncio MS, Nabuurs-Franssen MH, Klaassen CH: Genotypic diversity of clinical Coxiella burnetii isolates from Portugal based on MST and MLVA typing. Int J Med Microbiol. 2012, 302: 253-256. 10.1016/j.ijmm.2012.08.003.
Astobiza I, Tilburg JJ, Pinero A, Hurtado A, Garcia-Pérez AL, Nabuurs-Franssen MH, Klaassen CH: Genotyping of Coxiella burnetii from domestic ruminants in northern Spain. BMC Vet Res. 2012, 8: 241-10.1186/1746-6148-8-241.
Klaassen CH, Nabuurs-Franssen MH, Tilburg JJ, Hamans MA, Horrevorts AM: Multigenotype Q fever outbreak, the Netherlands. Emerg Infect Dis. 2009, 15: 613-614. 10.3201/eid1504.081612.
Dryad. [https://datadryad.org/, doi:10.5061/dryad.rc8q4]
Coxiella burnetii Cooperative MLVA database. [http://mlva.u-psud.fr/mlvav4/genotyping/]
Pearson T, Hornstra HM, Sahl JW, Schaack S, Schupp JM, Beckstrom-Sternberg SM, O’Neill MW, Priestley RA, Champion MD, Beckstrom-Sternberg JS, Kersh GJ, Samuel JE, Massung RF, Keim P: When outgroups fail; phylogenomics of rooting the emerging pathogen, Coxiella burnetii. Syst Biol. 2013, 62: 752-762. 10.1093/sysbio/syt038.
Pearson T, Okinaka RT, Foster JT, Keim P: Phylogenetic understanding of clonal populations in an era of whole genome sequencing. Infect Genet Evol. 2009, 2009 (9): 1010-1019.
Tilburg JJ, Roest HJ, Buffet S, Nabuurs-Franssen MH, Horrevorts AM, Raoult D, Klaassen CH: Epidemic genotype of Coxiella burnetii among goats, sheep, and humans in the Netherlands. Emerg Infect Dis. 2012, 18: 887-889. 10.3201/eid1805.111907.
Reichel R, Mearns R, Brunton L, Jones R, Horigan M, Vipond R, Vincent G, Evans S: Description of a Coxiella burnetii abortion outbreak in a dairy goat herd, and associated serology, PCR and genotyping results. Res Vet Sci. 2012, 93: 1217-1224. 10.1016/j.rvsc.2012.04.007.
Chmielewski T, Sidi-Boumedine K, Duquesne V, Podsiadly E, Thiery R, Tylewska-Wierzbanowska S: Molecular epidemiology of Q fever in Poland. Pol J Mic. 2009, 58: 9-13.
Mahamat A, Edouard S, Demar M, Abboud P, Patrice JY, La Scola B, Okzande A, Djossou F, Raoult D: Unique clone of Coxiella burnetii causing severe Q fever, French Guiana. Emerg Infect Diseases. 2013, 19: 1102-10.3201/eid1907.130044.
Pearson T, Hornstra HM, Hilsabeck R, Gates LT, Olivas SM, Birdsell D, Hall C, German S, Cook JM, Seymour ML, Priestley R, Kondas AV, Friedman CL, Price EP, Schupp J, Liu C, Price LB, Massung RF, Kersh GJ, Keim P: 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.
This work was supported by the Lendület program (LP2012-22) of the Hungarian Academy of Sciences. HMH, TRP and PSK were supported by the U.S. Department of Homeland Security S&T CB Division Bioforensics R&D Program (HSHQDC-10-C-00139).
The authors declare that they have no competing interests.
KMS genotyped the samples, analysed the data and wrote the manuscript. ZK, AS, ÁD, EB collected the samples, performed the PCR screening and edited the manuscript. HMH, TRP and PSK helped in genotyping and in the writing of the manuscript. MG designed the study and wrote the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Sulyok, K.M., Kreizinger, Z., Hornstra, H.M. et al. Genotyping of Coxiella burnetiifrom domestic ruminants and human in Hungary: indication of various genotypes. BMC Vet Res 10, 107 (2014). https://doi.org/10.1186/1746-6148-10-107
- Coxiella burnetii
- Q fever