Single nucleotide polymorphisms (SNPs) are ideal signatures for subtyping monomorphic pathogens such as Bacillus anthracis. Here we report the use of next-generation sequencing technology to investigate the historical, geographic and genetic diversity of Bacillus anthracis in France. 122 strains isolated over a 60-years period throughout the country were whole-genome sequenced and comparative analyses were carried out with a focus on SNPs discovery to discriminate regional sub-groups of strains.

Brucella sp. are the causative agents of brucellosis. One of the main characteristics of the Brucella genus concerns its very high genetic homogeneity. To date, classical bacteriology typing is still considered as the gold standard assay for direct diagnosis of Brucella. Molecular approaches are routinely used for the identification of Brucella at the genus level. However, genotyping is more complex, and to date, no method exists to quickly assign a strain into species and biovar levels, and new approaches are required. Next generation sequencing (NGS) opened a new era into the diagnosis of bacterial diseases. In this study, we designed a high-resolution melting (HRM) method for the rapid screening of DNA and direct assignment into one of the 12 species of the Brucella genus. This method is based on 17 relevant single nucleotide polymorphisms (SNPs), identified and selected from a whole genome SNP (wgSNP) analysis based on 988 genomes (complete and drafts). These markers were tested against the collection of the European Reference Laboratory (EU-RL) for brucellosis (1440 DNAs extracted from Brucella strains). The results confirmed the reliability of the panel of 17 SNP markers, allowing the differentiation of each species of Brucella together with biovars 1, 2, and 3 of B. suis and vaccine strain Rev1 (B. melitensis) within 3 h, which is a considerable gain of time for brucellosis diagnosis. Therefore, this genotyping tool provides a new and quick alternative for Brucella identification based on SNPs with the HRM-PCR assay.

In France, tularemia is caused by Francisella tularensis subsp. holarctica and is a sporadic disease affecting mainly wildlife animals and humans. F. tularensis species presents low genetic diversity that remains poorly described in France, as only a few genomes of isolates from the country are available so far. The objective of this study was to characterize the genetic diversity of F. tularensis in France and describe the phylogenetic distribution of isolates through whole-genome sequencing and molecular typing. Whole genomes of 350 strains of human or animal origin, collected from 1947 to 2018 in France and neighboring countries, were sequenced. A preliminary classification using the established canonical single nucleotide polymorphism (canSNP) nomenclature was performed. All isolates from France (except four) belonged to clade B.44, previously described in Western Europe. To increase the resolution power, a whole-genome SNP analysis was carried out. We were able to accurately reconstruct the population structure according to the global phylogenetic framework, and highlight numerous novel subclades. Whole-genome SNP analysis identified 87 new canSNPs specific to these subclades, among which 82 belonged to clade B.44. Identifying genomic features that are specific to sublineages is highly relevant in epidemiology and public health. We highlighted a large number of clusters among a single clade (B.44), which shows for the first time some genetic diversity among F. tularensis isolates from France, and the star phylogeny observed in clade B.44-subclades revealed that F. tularensis biodiversity in the country is relatively recent and resulted from clonal expansion of a single population. No association between clades and hosts or clinical forms of the disease was detected, but spatiotemporal clusters were identified for the first time in France. This is consistent with the hypothesis of persistence of F. tularensis strains found in Western Europe in the environment, associated with slow replication rates. Moreover, the presence of identical genotypes across long periods of time, and across long distances, supports this hypothesis but also suggests long-distance dispersal of the bacterium.

Bacillus anthracis, the causative agent of anthrax, has been widely described as a clonal species. Here we report the use of both canonical SNP analysis and whole-genome sequencing to characterize the phylogenetic lineages of B. anthracis from the Netherlands. Eleven strains isolated over a 25-years period (1968-1993) were paired-end sequenced using parallel sequencing technology. Five canSNP groups or lineages, i.e. A.Br.001/002 (n=6), A.Br.Aust94 (n=2), A.Br.008/011 (n=1), A.Br.011/009 (n=1) and A.Br.Vollum (n=1) were identified. Comparative analyses, with a focus on SNPs discovery, were carried out using a total of 52 B. anthracis genomes. A phylogeographic "Dutch" cluster within the dominant A.Br.001/002 group was discovered, involving isolates from a single outbreak. Diagnostic SNPs specific to the newly identified sub-groups were developed into high-resolution melting SNP discriminative assays for the purpose of rapid molecular epidemiology. Phylogenetic relationships with strains from other parts of the world are discussed.

Despite Brucella suis biovar 2’s (BSB2) active circulation in wildlife, no canine infections have been reported. The present paper is the first to describe two cases of BSB2 infections in French dogs. The first case occurred in 2020 and concerned a 13-year-old male neutered Border Collie with clinical signs of prostatitis. The urine culture revealed the excretion of significant levels of Brucella in the sample. The second case concerned a German Shepherd with bilateral orchitis, in which it was possible to detect Brucella colonies following neutering. HRM-PCR and classical biotyping methods classified both isolated strains as BSB2, in contrast to expected B. canis, which is usually the etiological agent of canine brucellosis in Europe. The wgSNP and MLVA analyses highlighted the genetic proximity of two isolates to BSB2 strains originating from wildlife. No pig farms were present in the proximity of either dog’s residence, ruling out potential spill over from infected pigs. Nevertheless, the dogs used to take walks in the surrounding forests, where contact with wildlife (i.e., wild boars or hares, or their excrements) was possible. These cases highlight the importance of adopting a One Health approach to control the presence of zoonotic bacteria in wild animals and avoid spillovers into domestic animals and, potentially, humans.

Anaplasma phagocytophilum is a zoonotic tick-borne intracellular alpha-proteobacterium causing tick-borne fever, which leads to significant economic losses in domestic ruminants in Europe. Its epidemiological cycles are complex and reservoir host species of bovine strains have not yet been identified. Given that little genetic information is available on strains circulating within a defined bovine environment, our objective was to assess the genetic diversity of A. phagocytophilum obtained from the same farms over time. Blood samplings were performed several times in two European herds. In the French herd, 169 EDTA-blood samples were obtained from 115 cows (32 were sampled two to four times). In the German herd, 20 cows were sampled six times (120 EDTA-blood samples). The presence of A. phagocytophilum DNA was assessed using a qPCR targeting msp2. The positive DNA samples underwent MLST at nine genetic markers (typA, ctrA, msp4, pleD, recG, polA, groEL, gyrA, and ankA). For each locus, sequences were aligned with available bacterial sequences derived from cattle, horse, dog, and roe deer hosts, and concatenated neighbor joining trees were constructed using three to six loci. Around 20% (57/289) of samples were positive. Forty positive samples from 23 French and six German cows (11 of them being positive at two time points) were sequenced. Six loci (typA, ctrA, msp4, pleD, recG, and polA) allowed to build concatenated phylogenetic trees, which led to two distinct groups of bovine variants in the French herd (hereafter called A and B), whereas only group A was detected in the German herd. In 42% of French samples, double chromatogram peaks were encountered in up to four loci. Eleven cows were found infected three weeks to 17 months after first sampling and harboured a new variant belonging to one or the other group. Our results demonstrate the occurrence of two major bovine strain groups and the simultaneous infection of single cows by more than one A. phagocytophilum strain. This challenges the role of cattle as reservoirs for A. phagocytophilum. This role may be facilitated via long-term bacterial persistence in individual cows and active circulation at the herd scale.