Eight new genomes and synthetic controls increase the accessibility of rapid melt-MAMA SNP typing of Coxiella burnetii

Phylogeography of Human and Animal Coxiella burnetii Strains: Genetic Fingerprinting of Q Fever in Belgium

Q fever is a zoonotic disease caused by the bacteria Coxiella burnetii. Domestic ruminants are the primary source for human infection, and the identification of likely contamination routes from the reservoir animals the critical point to implement control programs. This study shows that Q fever is detected in Belgium in abortion of cattle, goat and sheep at a different degree of apparent prevalence (1.93%, 9.19%, and 5.50%, respectively). In addition, and for the first time, it is detected in abortion of alpaca (Vicugna pacos), raising questions on the role of these animals as reservoirs. To determine the relationship between animal and human strains, Multiple Locus Variable-number Tandem Repeat Analysis (MLVA) (n=146), Single-Nucleotide Polymorphism (SNP) (n=92) and Whole Genome Sequencing (WGS) (n=4) methods were used to characterize samples/strains during 2009-2019. Three MLVA clusters (A, B, C) subdivided in 23 subclusters (A1-A12, B1-B8, C1-C3) and 3 SNP types (SNP1, SNP2, SNP6) were identified. The SNP2 type/MLVA cluster A was the most abundant and dispersed genotype over the entire territory, but it seemed not responsible for human cases, as it was only present in animal samples. The SNP1/MLVA B and SNP6/MLVA C clusters were mostly found in small ruminant and human samples, with the rare possibility of spillovers in cattle. SNP1/MLVA B cluster was present in all Belgian areas, while the SNP6/MLVA C cluster appeared more concentrated in the Western provinces. A broad analysis of European MLVA profiles confirmed the host-species distribution described for Belgian samples. In silico genotyping (WGS) further identified the spacer types and the genomic groups of C. burnetii Belgian strains: cattle and goat SNP2/MLVA A isolates belonged to ST61 and genomic group III, while the goat SNP1/MLVA B strain was classified as ST33 and genomic group II. In conclusion, Q fever is widespread in all Belgian domestic ruminants and in alpaca. We determined that the public health risk in Belgium is likely linked to specific genomic groups (SNP1/MLVA B and SNP6/MLVA C) mostly found in small ruminant strains. Considering the concordance between Belgian and European results, these considerations could be extended to other European countries.

Variant site strain typer (VaST): efficient strain typing using a minimal number of variant genomic sites

BACKGROUND

Targeted PCR amplicon sequencing (TAS) techniques provide a sensitive, scalable, and cost-effective way to query and identify closely related bacterial species and strains. Typically, this is accomplished by targeting housekeeping genes that provide resolution down to the family, genera, and sometimes species level. Unfortunately, this level of resolution is not sufficient in many applications where strain-level identification of bacteria is required (biodefense, forensics, clinical diagnostics, and outbreak investigations). Adding more genomic targets will increase the resolution, but the challenge is identifying the appropriate targets. VaST was developed to address this challenge by finding the minimum number of targets that, in combination, achieve maximum strain-level resolution for any strain complex. The final combination of target regions identified by the algorithm produce a unique haplotype for each strain which can be used as a fingerprint for identifying unknown samples in a TAS assay. VaST ensures that the targets have conserved primer regions so that the targets can be amplified in all of the known strains and it also favors the inclusion of targets with basal variants which makes the set more robust when identifying previously unseen strains.

RESULTS

We analyzed VaST's performance using a number of different pathogenic species that are relevant to human disease outbreaks and biodefense. The number of targets required to achieve full resolution ranged from 20 to 88% fewer sites than what would be required in the worst case and most of the resolution is achieved within the first 20 targets. We computationally and experimentally validated one of the VaST panels and found that the targets led to accurate phylogenetic placement of strains, even when the strains were not a part of the original panel design.

CONCLUSIONS

VaST is an open source software that, when provided a set of variant sites, can find the minimum number of sites that will provide maximum resolution of a strain complex, and it has many different run-time options that can accommodate a wide range of applications. VaST can be an effective tool in the design of strain identification panels that, when combined with TAS technologies, offer an efficient and inexpensive strain typing protocol.

Critical Aspects for Detection of Coxiella burnetii

Coxiella burnetii is a globally distributed zoonotic γ-proteobacterium with an obligatory intracellular lifestyle. It is the causative agent of Q fever in humans and of coxiellosis among ruminants, although the agent is also detected in ticks, birds, and various other mammalian species. Requirements for intracellular multiplication together with the necessity for biosafety level 3 facilities restrict the cultivation of C. burnetii to specialized laboratories. Development of a novel medium formulation enabling axenic growth of C. burnetii has facilitated fundamental genetic studies. This review provides critical insights into direct diagnostic methods currently available for C. burnetii. It encompasses molecular detection methods, isolation, and propagation of the bacteria and its genetic characterization. Differentiation of C. burnetii from Coxiella-like organisms is an essential diagnostic prerequisite, particularly when handling and analyzing ticks.

From Q Fever to Coxiella burnetii Infection: a Paradigm Change

Coxiella burnetii is the agent of Q fever, or "query fever," a zoonosis first described in Australia in 1937. Since this first description, knowledge about this pathogen and its associated infections has increased dramatically. We review here all the progress made over the last 20 years on this topic. C. burnetii is classically a strict intracellular, Gram-negative bacterium. However, a major step in the characterization of this pathogen was achieved by the establishment of its axenic culture. C. burnetii infects a wide range of animals, from arthropods to humans. The genetic determinants of virulence are now better known, thanks to the achievement of determining the genome sequences of several strains of this species and comparative genomic analyses. Q fever can be found worldwide, but the epidemiological features of this disease vary according to the geographic area considered, including situations where it is endemic or hyperendemic, and the occurrence of large epidemic outbreaks. In recent years, a major breakthrough in the understanding of the natural history of human infection with C. burnetii was the breaking of the old dichotomy between "acute" and "chronic" Q fever. The clinical presentation of C. burnetii infection depends on both the virulence of the infecting C. burnetii strain and specific risks factors in the infected patient. Moreover, no persistent infection can exist without a focus of infection. This paradigm change should allow better diagnosis and management of primary infection and long-term complications in patients with C. burnetii infection.

Molecular epidemiology of Coxiella burnetii in French livestock reveals the existence of three main genotype clusters and suggests species-specific associations as well as regional stability

Q fever is a worldwide zoonosis caused by the bacterium Coxiella burnetii. In domestic ruminants, Q fever main clinical manifestations are abortions. Although the clinical signs may differ between ruminant species, C. burnetii's genetic diversity remains understudied in enzootic areas. Here, we focused on France, where Q fever is enzootic, with the aims to (a) identify potential associations between C. burnetii genotypes and ruminant host species; (b) assess the distribution of C. burnetii genotypes both within French farms and across France's major livestock-farming regions; and (c) suggest a subset of markers for future genotypic studies. We used DNA samples collected between 2006 and 2015 from 301 females (160 cows, 76 ewes, 65 goats) aborted of Q fever within 7 different farming regions. C. burnetii diversity was determined using a multiple-locus variable-number of tandem repeat analysis (MLVA) considering 17 markers. Using a phylogenetic approach, we identified 3 main genotypic clusters divided into 12 sub-clusters. These clusters were significantly associated with ruminant species: almost all the cattle genotypes were found in a "cattle-specific" cluster whereas small ruminants genotypes essentially grouped into the two other clusters. The clusters also proved stable over space and time, some genotypes being more specifically observed in certain farming regions. We also observed some within-farm diversity but this diversity was restricted to a same genotypic cluster. Finally, we identified 6 MLVA markers that maximized the representativeness of the diversity described. Overall, we highlighted that molecular epidemiology is a relevant approach to assess C. burnetii's genetic diversity and to reveal the existence of species-specific associations and regional stability. These results will be valuable in the field to trace genotype circulation among ruminants and from ruminants to humans. Ultimately, the potential links between genotypes and virulence traits need to be investigated to adapt control measures in livestock farms.

Introduction and persistence of tularemia in Bulgaria

Introduction

Outbreaks of the zoonotic disease tularemia occurred in north-east Bulgaria in the 1960s. Then came 30 years of epidemiological silence until new outbreaks occurred in west Bulgaria in the 1990s. To investigate how bacterial strains of Francisella tularensis causing tularemia in wildlife and humans in the 1960s and the 1990s were related, we explored their genetic diversity.

Material and methods

Ten F. tularensis genomes from the 1960s (n=3) and the 1990s (n=7) were sequenced, assigned to canonical single-nucleotide polymorphism (canSNP) clades, and compared to reference genomes. We developed four new canSNP polymerase chain reaction (PCR) assays based on the genome sequence information.

Results and discussion

The genetic analysis showed that the outbreaks in the 1960s as well as in the 1990s involved multiple clones and new genetic diversity. The smallest genetic difference found between any of the Bulgarian strains was five SNPs between the strains L2 and 81 isolated 43 years apart, indicating that F. tularensis may persist locally over long time periods without causing outbreaks. The existence of genetically highly similar strain-pairs isolated the same year in the same area from different hosts supports a hypothesis of local expansion of clones during outbreaks. Close relationship (two SNPs) was found between one strain isolated 1961 in northeast Bulgaria and one strain isolated 5 years before in USSR. Historical data coinciding with the actual time point describe the introduction of water rats from USSR into the Bulgarian outbreak area, which may explain the close genetic relationship and the origin of the outbreak.

Conclusion

Genome analysis of strains from two outbreaks in the 1960s and the 1990s provided valuable information on the genetic diversity and persistence of F. tularensis in Bulgaria.

HuMiChip2 for strain level identification and functional profiling of human microbiomes

With the massive data generated by the Human Microbiome Project, how to transform such data into useful information and knowledge remains challenging. Here, with currently available sequencing information (reference genomes and metagenomes), we have developed a comprehensive microarray, HuMiChip2, for strain-level identification and functional characterization of human microbiomes. HuMiChip2 was composed of 29,467 strain-specific probes targeting 2063 microbial strains/species and 133,924 sequence- and group-specific probes targeting 157 key functional gene families involved in various metabolic pathways and host-microbiome interaction processes. Computational evaluation of strain-specific probes suggested that they were not only specific to mock communities of sequenced microorganisms and metagenomes from different human body sites but also to non-sequenced microbial strains. Experimental evaluation of strain-specific probes using single strains/species and mock communities suggested a high specificity of these probes with their corresponding targets. Application of HuMiChip2 to human gut microbiome samples showed the patient microbiomes of alcoholic liver cirrhosis significantly (p < 0.05) shifted their functional structure from the healthy individuals, and the relative abundance of 21 gene families significantly (p < 0.1) differed between the liver cirrhosis patients and healthy individuals. At the strain level, five Bacteroides strains were significantly (p < 0.1) and more frequently detected in liver cirrhosis patients. These results suggest that the developed HuMiChip2 is a useful microbial ecological microarray for both strain-level identification and functional profiling of human microbiomes.

[Not Available]

Q fever is a zoonosis of worldwide distribution with the exception of New Zealand. It is caused by an intracellular bacterium, Coxiella burnetii. The disease often goes underdiagnosed because the main manifestation of its acute form is a general self-limiting flu-like syndrome. The Dutch epidemics renewed attention to this disease, which was less considered before. This review summarizes the description of C. burnetii (taxonomy, intracellular cycle, and genome) and Q fever disease (description, diagnosis, epidemiology, and pathogenesis). Finally, vaccination in humans and animals is also considered.

ConStrains identifies microbial strains in metagenomic datasets

An important fraction of microbial diversity is harbored in strain individuality, so identification of conspecific bacterial strains is imperative for improved understanding of microbial community functions. Limitations in bioinformatics and sequencing technologies have to date precluded strain identification owing to difficulties in phasing short reads to faithfully recover the original strain-level genotypes, which have highly similar sequences. We present ConStrains, an open-source algorithm that identifies conspecific strains from metagenomic sequence data and reconstructs the phylogeny of these strains in microbial communities. The algorithm uses single-nucleotide polymorphism (SNP) patterns in a set of universal genes to infer within-species structures that represent strains. Applying ConStrains to simulated and host-derived datasets provides insights into microbial community dynamics.

Genome sequence of Coxiella burnetii strain Namibia

We present the whole genome sequence and annotation of the Coxiella burnetii strain Namibia. This strain was isolated from an aborting goat in 1991 in Windhoek, Namibia. The plasmid type QpRS was confirmed in our work. Further genomic typing placed the strain into a unique genomic group. The genome sequence is 2,101,438 bp long and contains 1,979 protein-coding and 51 RNA genes, including one rRNA operon. To overcome the poor yield from cell culture systems, an additional DNA enrichment with whole genome amplification (WGA) methods was applied. We describe a bioinformatics pipeline for improved genome assembly including several filters with a special focus on WGA characteristics.

Genome Sequence of Coxiella burnetii Strain AuQ01 (Arandale) from an Australian Patient with Acute Q Fever

Coxiella burnetii strain AuQ01 was isolated from the serum of an Australian acute Q fever patient and represents the first whole genome from this historical Q fever country. This new genome shows distinct differences from existing genomic data and will enhance the understanding of this query pathogen.

Genotyping of French Bacillus anthracis strains based on 31-loci multi locus VNTR analysis: epidemiology, marker evaluation, and update of the internet genotype database

Background

Bacillus anthracis is known to have low genetic variability. In spite of this lack of diversity, multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA) and single nucleotide polymorphisms (SNPs) including the canonical SNPs assay (canSNP) have proved to be highly effective to differentiate strains. Five different MLVA schemes based on a collection of 31 VNTR loci (MLVA8, MLVA15, MLVA20, MLVA25 and MLVA31) with increased resolving power have been described.

Results

MLVA31 was applied to characterize the French National Reference Laboratory collection. The total collection of 130 strains is resolved in 35 genotypes. The 119 veterinary and environmental strains collection in France were resolved into 26 genotypes belonging to three canSNP lineages and four MLVA clonal complexes (CCs) with particular geographical clustering. A subset of seven loci (MLVA7) is proposed to constitute a first line assay. The loci are compatible with moderate resolution equipment such as agarose gel electrophoresis and show a good congruence value with MLVA31. The associated MLVA and SNP data was imported together with published genotyping data by taking advantage of major enhancements to the MLVAbank software and web site.

Conclusions

The present report provides a wide coverage of the genetic diversity of naturally occurring B. anthracis strains in France as can be revealed by MLVA. The data obtained suggests that once such coverage is achieved, it becomes possible to devise optimized first-line MLVA assays comprising a sufficiently low number of loci to be typed either in one multiplex PCR or on agarose gels. Such a selection of seven loci is proposed here, and future similar investigations in additional countries will indicate to which extend the same selection can be used worldwide as a common minimum set. It is hoped that this approach will contribute to an efficient and low-cost routine surveillance of important pathogens for biosecurity such as B. anthracis.

Draft Genome Sequences of Six Ruminant Coxiella burnetii Isolates of European Origin

Coxiella burnetii is responsible for Q fever, a worldwide zoonosis attributed to the inhalation of aerosols contaminated by livestock birth products. Six draft genome sequences of European C. burnetii isolates from ruminants are presented here. The availability of these genomes will help in understanding the potential host specificity and pathogenicity and in identifying pertinent markers for surveillance and tracing.

CanSNPer: a hierarchical genotype classifier of clonal pathogens

SUMMARY

Advances in typing methodologies have recently reformed the field of molecular epidemiology of pathogens. The falling cost of sequencing technologies is creating a deluge of whole genome sequencing data that burdens bioinformatics resources and tool development. In particular, single nucleotide polymorphisms in core genomes of pathogens are recognized as the most important markers for inferring genetic relationships because they are evolutionarily stable and amenable to high-throughput detection methods. Sequence data will provide an excellent opportunity to extend our understanding of infectious disease when the challenge of extracting knowledge from available sequence resources is met. Here, we present an efficient and user-friendly genotype classification pipeline, CanSNPer, based on an easily expandable database of predefined canonical single nucleotide polymorphisms.

AVAILABILITY AND IMPLEMENTATION

All documentation and Python-based source code for the CanSNPer are freely available at http://github.com/adrlar/CanSNPer.