Core Genome Multilocus Sequence Typing and Single Nucleotide Polymorphism Analysis in the Epidemiology of Brucella melitensis Infections

Molecular epidemiology of brucellosis in Asia: insights from genotyping analyses

Brucellosis infects humans and animals worldwide but is particularly prevalent in Asia. In many Asian countries, molecular diagnostic tools for accurate molecular diagnostics and molecular epidemiology are lacking. Nonetheless, some countries have conducted in-depth molecular epidemiological studies. The objective of this study was to reveal the genetic relationships, geographic origins, and distributions of Brucella strains across Asia for two primary species, B. abortus and B. melitensis. For this, we systematically searched genotyping data from published studies on the molecular epidemiology of Brucella species for both humans and livestock in Asia. We used data from multilocus sequence typing (MLST), multiple-locus variable-number tandem repeat analysis (MLVA), and whole genome sequencing analysis of Brucella strains. We also analyzed the MLVA genotypes of 129 B. abortus isolates and 242 B. melitensis isolates with known origins in Asia from an online MLVA database using MLVA-11 data in minimum spanning trees and MLVA-16 data in neighbor-joining trees. We found that the B. melitensis East Mediterranean lineage is predominant across the continent, with only a small number of samples from the Africa and Americas lineages, and none from the West Mediterranean lineage. The "abortus C" genotype was the most common group of B. abortus in Asia, with limited genetic variation for this species. Several studies also reported that Near Eastern countries frequently encounter human brucellosis cases of B. abortus from genotypes 42 and 43. Our study highlights the inconsistent collection of genetic data for Brucella species across Asia and a need for more extensive sampling in most countries. Finally, a consistent nomenclature is necessary to define various groupings of strains within a lineage (i.e., clade) so uniform terminology should denote particular genetic groups that are understood by all researchers.

An Evaluation of the Lineage of Brucella Isolates in Turkey by a Whole-Genome Single-Nucleotide Polymorphism Analysis

Brucellosis is a disease caused by the Brucella (B.) species. It is a zoonotic disease that affects farm animals and causes economic losses in many countries worldwide. Brucella has the ability to persist in the environment and infect the host at low doses. Thus, it is more important to trace brucellosis outbreaks, identify their sources of infection, and interrupt their transmission. Some countries already have initial data, but most of these data are based on a Multiple-Locus Variable-Number Tandem-Repeat Analysis (MLVA), which is completely unsuitable for studying the Brucella genome. Since brucellosis is an endemic disease in Turkey, this study aimed to examine the genome of Turkish Brucella isolates collected between 2018 and 2020, except for one isolate, which was from 2012. A total of 28 strains of B. melitensis (n = 15) and B. abortus (n = 13) were analyzed using a core-genome single-nucleotide polymorphism (cgSNP) analysis. A potential connection between the Turkish isolates and entries from Sweden, Israel, Syria, Austria, and India for B. melitensis was detected. For B. abortus, there may be potential associations with entries from China. This explains the tight ties found between Brucella strains from neighboring countries and isolates from Turkey. Therefore, it is recommended that strict measures be taken and the possible effects of uncontrolled animal introduction are emphasized.

Delay in the diagnosis of Brucella abortus bacteremia in a nonendemic country: a case report

BACKGROUND

It is challenging to diagnose brucellosis in nonendemic regions because it is a nonspecific febrile disease. The accurate identification of Brucella spp. in clinical microbiology laboratories (CMLs) continues to pose difficulties. Most reports of misidentification are for B. melitensis, and we report a rare case of misidentified B. abortus.

CASE PRESENTATION

A 67-year-old man visited an outpatient clinic complaining of fatigue, fever, and weight loss. The patient had a history of slaughtering cows with brucellosis one year prior, and his Brucella antibody tests were negative twice. After blood culture, the administration of doxycycline and rifampin was initiated. The patient was hospitalized due to a positive blood culture. Gram-negative coccobacilli were detected in aerobic blood culture bottles, but the CML's lack of experience with Brucella prevented appropriate further testing. Inaccurate identification results were obtained for a GN ID card of VITEK 2 (bioMérieux, USA) and matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) using a MALDI Biotyper (Bruker, Germany). The strain showed 100.0% identity with Brucella spp. according to 16S rRNA sequencing. MALDI-TOF MS peaks were reanalyzed using the CDC MicrobeNet database to determine Brucella spp. (score value: 2.023). The patient was discharged after nine days of hospitalization and improved after maintaining only doxycycline for six weeks. The isolate was also identified as Brucella abortus by genomic evidence.

CONCLUSION

Automated identification instruments and MALDI-TOF MS are widely used to identify bacteria in CMLs, but there are limitations in accurately identifying Brucella spp. It is important for CMLs to be aware of the possibility of brucellosis through communication with clinicians. Performing an analysis with an additional well-curated MALDI-TOF MS database such as Bruker security-relevant (SR) database or CDC MicrobeNet database is helpful for quickly identifying the genus Brucella.

Molecular Epidemiological Study of a Human Brucellosis Outbreak - Weihai City, Shandong Province, China, 2022

What is already known about this topic?

Brucellosis, mainly caused by Brucella melitensis (B. melitensis), is regarded as a significant zoonotic disease in China. In Weihai, located at the eastern end of the Shandong Peninsula, brucellosis has been in a low epidemic phase for the past five years.

What is added by this report?

This was the initial report of a brucellosis outbreak in the last five years. Strains of B. melitensis bv. 3 from Weihai and other cities showed a close genetic relationship, suggesting a potential common ancestry.

What are the implications for public health practice?

Epidemiological investigations depend on standardized and effective molecular typing methods and analysis tools for public health laboratories to identify and trace outbreaks. Understanding the circulation patterns of livestock in free-range households in heavily affected areas is essential for controlling the spread of brucellosis.

New insights into the pathogenesis and transmission of Brucella pinnipedialis: systemic infection in two bottlenose dolphins (Tursiops truncatus)

IMPORTANCE

Brucella spp. are zoonotic pathogens that can affect both terrestrial and marine mammals. Brucella ceti has been identified in various cetacean species, but only one sequence type (ST27) has been reported in humans. However, it is important to conduct surveillance studies to better understand the impact of marine Brucella species on marine mammals, a typically understudied host group. Here, we describe a systemic infection by two related strains of Brucella pinnipedialis (ST25) in a couple of live-stranded bottlenose dolphins, with more severe lesions in the younger animal. Furthermore, B. pinnipedialis was first detected in milk from a female cetacean that stranded with its offspring. Our study reveals novel insights into the epidemiology and pathological consequences of B. pinnipedialis infections in cetaceans, emphasizing the crucial importance of ongoing surveillance and accurate diagnosis to understand the impact of this pathogen on marine mammal populations.

Genotypic peculiarities of a human brucellosis case caused by Brucella suis biovar 5

Human brucellosis cases are rare in non-endemic countries, such as Germany, where infections are predominantly caused by Brucella melitensis. The German National Reference Laboratory for Bovine, Porcine, Ovine and Caprine Brucellosis received a suspected Brucella sp. isolate from a patient for identification. Bacteriological tests and PCR-based diagnostics showed the isolate to be B. suis, but did not yield cohesive results regarding the biovar. Whole genome sequencing and subsequent genotyping was employed for a detailed characterization of the isolate and elucidating the reason for failure of the diagnostic PCR to correctly identify the biovar. The isolate was found to be B. suis bv. 5, a rare biovar with limited geographical distribution primarily found in the Northern Caucasus. Due to a deletion in one of the target regions of the diagnostic PCR, the isolate could not be correctly typed. Based on in silico genotyping it could be excluded that the isolate was identical to one of the B. suis bv. 5 reference strains. Here, we report a rare case of a B. suis bv. 5 field isolate. Furthermore, by reporting this finding, we want to make practitioners aware of possible misinterpretation of PCR results, as it cannot be excluded that the detected deletion is common among the B. suis bv. 5 community, as there is currently a lack of field isolates.

Native circulating Brucella melitensis lineages causing a brucellosis epidemic in Qinghai, China

Since 2010, the cases and incidences of human brucellosis have been increasing annually in Qinghai (QH) Province. Molecular epidemiology and phylogenetic analyses of strains from this region are crucial to better understand the transmission of the disease and the evolutionary patterns of Brucella strains. In this study, classical bio-typing assay, multilocus variable-number tandem repeat analysis, and the whole-genome sequencing-single-nucleotide polymorphism approach were used to illustrate the epidemiological and evolutionary patterns of Brucella melitensis. A total of 54 B. melitensis bv. 3 strains were isolated and molecularly characterized, with all strains belonging to the East Mediterranean lineages. Cross-regional transmission events (i.e., between counties) were caused by common sources of infection, suggesting that predominant circulating genotypes are endemic in different regions. Strengthening surveillance in animal brucellosis and controlling infected animals' cross-border movement are necessary. Two strains isolated from humans and marmots were clustered in the same sub-clade, implying the possible existence of direct and/or indirect contact between sheep (and goats) and wildlife (marmots), but this needs to be verified by further investigations. The global-scale phylogenetic analysis indicated that 54 strains sorted into six subclades, four of which formed independent lineages, suggesting that the increase in the incidence rate of human brucellosis may be caused by local circulating lineages. Further strengthening the serology and pathogen surveillance of animals (wildlife) and humans will contribute to an in-depth understanding of the transmission chain of human brucellosis in this region.

Genotype diversity of brucellosis agents isolated from humans and animals in Greece based on whole-genome sequencing

BACKGROUND

Brucellosis is a zoonotic disease whose causative agent, Brucella spp., is endemic in many countries of the Mediterranean basin, including Greece. Although the occurrence of brucellosis must be reported to the authorities, it is believed that the disease is under-reported in Greece, and knowledge about the genomic diversity of brucellae is lacking.

METHODS

Thus, 44 Brucella isolates, primarily B. melitensis, collected between 1999 and 2009 from humans and small ruminants in Greece were subjected to whole genome sequencing using short-read technology. The raw reads and assembled genomes were used for in silico genotyping based on single nucleotide substitutions and alleles. Further, specific genomic regions encoding putative virulence genes were screened for characteristic nucleotide changes, which arose in different genotype lineages.

RESULTS

In silico genotyping revealed that the isolates belonged to three of the known sublineages of the East Mediterranean genotype. In addition, a novel subgenotype was identified that was basal to the other East Mediterranean sublineages, comprising two Greek strains. The majority of the isolates can be assumed to be of endemic origin, as they were clustered with strains from the Western Balkans or Turkey, whereas one strain of human origin could be associated with travel to another endemic region, e.g. Portugal. Further, nucleotide substitutions in the housekeeping gene rpoB and virulence-associated genes were detected, which were characteristic of the different subgenotypes. One of the isolates originating from an aborted bovine foetus was identified as B. abortus vaccine strain RB51.

CONCLUSION

The results demonstrate the existence of several distinct persistent Brucella sp. foci in Greece. To detect these and for tracing infection chains, extensive sampling initiatives are required.

Isolation of Brucella inopinata from a White's tree frog (Litoria caerulea): pose exotic frogs a potential risk to human health?

Introduction

Cold-blooded hosts, particularly exotic frogs, have become a newly recognized reservoir for atypical Brucella species and strains worldwide, but their pathogenicity to humans remains largely unknown. Here we report the isolation and molecular characterization of a B. inopinata strain (FO700662) cultured from clinical samples taken from a captive diseased White's Tree Frog (Litoria caerulea) in Switzerland. The isolation of B. inopinata from a frog along with other reports of human infection by atypical Brucella raises the question of whether atypical Brucella could pose a risk to human health and deserves further attention.

Methods

The investigations included histopathological analysis of the frog, bacterial culture and in-depth molecular characterization of strain FO700662 based on genome sequencing data.

Results and Discussion

Originally identified as Ochrobactrum based on its rapid growth and biochemical profile, strain FO700622 was positive for the Brucella- specific markers bcsp31 and IS711. It showed the specific banding pattern of B. inopinata in conventional Bruce-ladder multiplex PCR and also had identical 16S rRNA and recA gene sequences as B. inopinata. Subsequent genome sequencing followed by core genome-based MLST (cgMLST) analysis using 2704 targets (74% of the total chromosome) revealed only 173 allelic differences compared to the type strain of B. inopinata BO1^T, while previously considered the closest related strain BO2 differed in 2046 alleles. The overall average nucleotide identity (ANI) between the type strain BO1^T and FO700622 was 99,89%, confirming that both strains were almost identical. In silico MLST-21 and MLVA-16 also identified strain FO700662 as B. inopinata. The nucleotide and amino acid-based phylogenetic reconstruction and comparative genome analysis again placed the isolate together with B. inopinata with 100% support. In conclusion, our data unequivocally classified strain FO700622, isolated from an exotic frog, as belonging to B. inopinata.

Whole-genome sequencing for genetic diversity analysis of Iranian Brucella spp. isolated from humans and livestock

Brucellosis is one of the most common zoonoses in the Middle East. It is causing economic losses to the livestock industry and has a great public health concern. Little is known about the genetic diversity and distribution of brucellae in Iran. Therefore, forty Brucella spp. strains (B. abortus and B. melitensis) isolated from animals and humans were analyzed by whole genome sequencing (WGS) technology using single nucleotide polymorphism (SNP) analysis and core genome multilocus sequence typing (cgMLST). Brucella isolates were obtained from lymph nodes (cows and camels), milk (cows, camels and sheep), and aborted foetus samples (sheep and goats), as well as cerebrospinal fluid and blood of humans. The isolates were originating from thirteen provinces of Iran and isolated between 2015 and 2020. According to in-silico MLST, ST8 and ST2 were the most frequent sequence types in B. melitensis and B. abortus, respectively. Based on phylogeographic reconstruction using cgSNP analysis, the investigated Iranian B. melitensis strains belonged to the American and Mediterranean lineages of the B. melitensis phylogeny. Furthermore, cgSNP analysis revealed a similarity between Iranian B. abortus isolates and strains from Iraq and Egypt. Therefore, the origin of the Iranian strains can be suggested to be strains from neighboring and Middle East countries. Moreover, cgMLST analysis showed that the Iranian B. melitensis strains were closely relative to strains recovered from sheep and humans in Iraq, Afghanistan, Syria, Turkmenistan, and Pakistan. In the current panel of strains, cgMLST and cgSNP analysis provided an appropriate and accurate tool for effective traceback analyses for Brucella spp. from Iran. The results of cgSNP and cgMLST helped to understand the geographic distribution and interspecies transmission of Iranian strains and highlight the importance of specific brucellosis control measures in Iran with regard to the One-Health approach.

Public and animal health risks associated with spillover of Brucella melitensis into dairy farms

Brucellosis is a worldwide zoonosis with important public health, animal health and economic implications. Brucella melitensis, commonly associated with small ruminants, is an emerging bovine pathogen in dairy farms. We analysed all B. melitensis outbreaks affecting dairy farms in Israel since 2006, combining traditional and genomic epidemiology to explore the public health implications of this One Health challenge. Whole-genome sequencing was applied to bovine and related human B. melitensis isolates from dairy farm outbreaks. cgMLST-based and SNP-based typing was integrated with epidemiological and investigation data. A secondary analysis combining the bovine-human isolates with endemic human isolates from southern Israel was performed. A total of 92 isolates from dairy cows and related human cases originating from 18 epidemiological clusters were analysed. Most genomic and epi-clusters were congruent, but sequencing showed relatedness between apparently unrelated farm outbreaks. Nine secondary human infections were also genomically confirmed. The bovine-human cohort appeared intermixed with 126 endemic human isolates in southern Israel. We show a persistent and widespread circulation of B. melitensis in dairy farms in Israel with secondary occupational human infection. The genomic epidemiology also uncovered cryptic connections between outbreaks. A regional connection between bovine and endemic human brucellosis cases points to a common reservoir, most probably local small ruminant herds. Control of humans and bovine brucellosis is inseparable. Epidemiological and microbiological surveillance and implementation of control measures across the entire range of farm animals is needed to mitigate this public health challenge.

Seroprevalence of Brucella Infection in Wild Boars (Sus scrofa) of Bavaria, Germany, 2019 to 2021 and Associated Genome Analysis of Five B. suis Biovar 2 Isolates

Brucella species are highly pathogenic zoonotic agents and are found in vertebrates all over the world. To date, Germany is officially declared free from brucellosis and continuous surveillance is currently limited to farm ruminants. However, porcine brucellosis, mostly caused by B. suis biovar 2, is still found in wild boars and hares. In the present study, a three-year monitoring program was conducted focusing on the wild boar population in the state of Bavaria. Serologic screening of 11,956 animals and a direct pathogen detection approach, including a subset of 681 tissue samples, was carried out. The serologic incidence was 17.9%, which is in approximate accordance with previously published results from various European countries. Applying comparative whole genome analysis, five isolated B. suis biovar 2 strains from Bavaria could be assigned to three known European genetic lineages. One isolate was closely related to another strain recovered in Germany in 2006. Concluding, porcine brucellosis is endemic in Bavaria and the wild boar population represents a reservoir for genetically distinct B. suis biovar 2 strains. However, the transmission risk of swine brucellosis to humans and farm animals is still regarded as minor due to low zoonotic potential, awareness, and biosafety measures.

Brucella abortus in Kazakhstan, population structure and comparison with worldwide genetic diversity

Brucella abortus is the main causative agent of brucellosis in cattle, leading to severe economic consequences in agriculture and affecting public health. The zoonotic nature of the infection increases the need to control the spread and dynamics of outbreaks in animals with the incorporation of high resolution genotyping techniques. Based on such methods, B. abortus is currently divided into three clades, A, B, and C. The latter includes subclades C1 and C2. This study presents the results of whole-genome sequencing of 49 B. abortus strains isolated in Kazakhstan between 1947 and 2015 and of 36 B. abortus strains of various geographic origins isolated from 1940 to 2004. In silico Multiple Locus Sequence Typing (MLST) allowed to assign strains from Kazakhstan to subclades C1 and to a much lower extend C2. Whole-genome Single-Nucleotide Polymorphism (wgSNP) analysis of the 46 strains of subclade C1 with strains of worldwide origins showed clustering with strains from neighboring countries, mostly North Caucasia, Western Russia, but also Siberia, China, and Mongolia. One of the three Kazakhstan strains assigned to subclade C2 matched the B. abortus S19 vaccine strain used in cattle, the other two were genetically close to the 104 M vaccine strain. Bayesian phylodynamic analysis dated the introduction of B. abortus subclade C1 into Kazakhstan to the 19th and early 20th centuries. We discuss this observation in view of the history of population migrations from Russia to the Kazakhstan steppes.

Serum PCR Diagnosis of Brucella melitensis Infection in Rev. 1 Vaccinated Sheep

Serological diagnosis provides a robust and effective approach to monitoring and controlling small ruminant brucellosis. Brucella melitensis Rev. 1 is a live vaccine strain used in prophylactic vaccination against small ruminant brucellosis. Because the vaccine strain shares identical serological antigens with the corresponding field strains, differentiating infected from vaccinated animals (DIVA) serological responses hamper surveillance campaigns and interventions that involve vaccination. We have developed a serum PCR-based approach in which we amplify and sequence Brucella omp2a as a DIVA solution and tRNA (uracil-5-)-methyltransferase as a species marker in the serum samples to determine the etiological agent involved in brucellosis field cases. Using this method, we identified the involvement of both the Rev. 1 vaccine strain and a field strain in an outbreak of brucellosis in a flock. This method represents a novel approach in studying the etiology of brucellosis using serum samples as a source of the pathogen’s DNA.

First record of the human infection of Brucella melitensis in Kyrgyzstan: evidence from whole-genome sequencing-based analysis

BACKGROUND

Brucellosis, a zoonosis mainly transmitted by consumption of unpasteurized dairy products as well as direct contact with infected animals, is endemic in Kyrgyzstan. However, Brucella species in humans have not been investigated and the origin of the disease remains poorly known in wide parts of Сentral Asia. Thus, molecular characterization of the circulating strains is a critical first step in understanding Brucella diversity in the country.

METHODS

In this study, isolates were collected from patients with suspected brucellosis from different regions in Kyrgyzstan between 2019 and 2020. The detection and identification of Brucella was carried out by Bruce-ladder PCR. Next generation sequencing was used to sequence the 89 Brucella isolates, which were genotyped by cgSNP and cgMLST to identify epidemiological connection between Brucella isolates as well as placing them in the context of the global Brucella phylogeny.

RESULTS

The Brucella strains isolated from all regions of Kyrgyzstan were identified as B. melitensis. Based on cgSNP analysis, 18 sequence types were differentiated. The highest numbers of different sequence types were found in Batken (n = 8), Osh (n = 8) and Jalal-Abad (n = 6) oblasts. According to cgSNP and cgMLST analyses, different B. melitensis lineages circulate in Kyrgyzstan, all of them belonging to the Eastern Mediterranean group of the global Brucella phylogeny with the highest similarity to strains from Turkmenistan, Iran and Turkey.

CONCLUSION

In the present study, B. melitensis was identified as a causative agent of human brucellosis in Kyrgyzstan and different lineages could be identified. Since this study focused on isolates of human origin, the identity of Brucella species and lineages circulating among animal populations remains elusive. Implementing culture techniques and use of most recent molecular, bioinformatic and epidemiological tools are needed to set up a One Health approach to combat brucellosis in Kyrgyzstan. Further, other Сentral Asian countries need to take part in this effort as brucellosis is a transboundary disease in these regions.

Core Genome Multilocus Sequence Typing Scheme for Improved Characterization and Epidemiological Surveillance of Pathogenic Brucella

Brucellosis poses a significant burden to human and animal health worldwide. Robust and harmonized molecular epidemiological approaches and population studies that include routine disease screening are needed to efficiently track the origin and spread of Brucella strains. Core genome multilocus sequence typing (cgMLST) is a powerful genotyping system commonly used to delineate pathogen transmission routes for disease surveillance and control. Except for Brucella melitensis, cgMLST schemes for Brucella species are currently not established. Here, we describe a novel cgMLST scheme that covers multiple Brucella species. We first determined the phylogenetic breadth of the genus using 612 Brucella genomes. We selected 1,764 genes that were particularly well conserved and typeable in at least 98% of these genomes. We tested the new scheme on 600 genomes and found high agreement with the whole-genome-based single nucleotide polymorphism (SNP) analysis. Next, we applied the scheme to reanalyze the genome of Brucella strains from epidemiologically linked outbreaks. We demonstrated the applicability of the new scheme for high-resolution typing required in outbreak investigations as previously reported with whole-genome SNP methods. We also used the novel scheme to define the global population structure of the genus using 1,322 Brucella genomes. Finally, we demonstrated the possibility of tracing distribution of Brucella strains by performing cluster analysis of cgMLST profiles and found nearly identical cgMLST profiles in different countries. Our results show that sequencing depth of more than 40-fold is optimal for allele calling with this scheme. In summary, this study describes a novel Brucella-wide cgMLST scheme that is applicable in Brucella molecular epidemiology and helps in accurately tracking and thus controlling the sources of infection. The scheme is publicly accessible and should represent a valuable resource for laboratories with limited computational resources and bioinformatics expertise.

Long-Read Sequencing and Hybrid Assembly for Genomic Analysis of Clinical Brucella melitensis Isolates

Brucella melitensis is a key etiological agent of brucellosis and has been increasingly subject to characterization using sequencing methodologies. This study aimed to investigate and compare short-read, long-read, and hybrid assemblies of B. melitensis. Eighteen B. melitensis isolates from Southern Israel were sequenced using Illumina and the Oxford Nanopore (ONP) MinION, and hybrid assemblies were generated with ONP long reads scaffolded on Illumina short reads. Short reads were assembled with INNUca with SPADes, long reads and hybrid with dragonflye. Abricate with the virulence factor database (VFDB) and in silico PCR (for the genes BetB, BPE275, BSPB, manA, mviN, omp19, perA, PrpA, VceC, and ureI) were used for identifying virulence genes, and a total of 61 virulence genes were identified in short-read, long-read, and hybrid assemblies of all 18 isolates. The phylogenetic analysis using long-read assemblies revealed several inconsistencies in cluster assignment as compared to using hybrid and short-read assemblies. Overall, hybrid assembly provided the most comprehensive data, and stand-alone short-read sequencing provided comparable data to stand-alone long-read sequencing regarding virulence genes. For genomic epidemiology studies, stand-alone ONP sequencing may require further refinement in order to be useful in endemic settings.

Genomic Epidemiology of Clinical Brucella melitensis Isolates from Southern Israel

Brucellosis, a zoonosis mainly transmitted by consumption of unpasteurized dairy products, is endemic in Southern Israel, mainly among the Bedouin Arab population. However, the genomic epidemiology of B. melitensis in this region has not yet been elucidated. A cohort of brucellosis cases (n = 118) diagnosed between 2017–2019 was studied using whole-genome sequencing (WGS). Phylogenetic analyses utilized core genome MLST (cgMLST) for all local isolates and core genome SNPs for 347 human-associated B. melitensis genomes, including Israeli and publicly available sequences. Israeli isolates formed two main clusters, presenting a notable diversity, with no clear dominance of a specific strain. On a global scale, the Israeli genomes clustered according to their geographical location, in proximity to genomes originating from the Middle East, and formed the largest cluster in the tree, suggesting relatively high conservation. Our study unveils the genomic epidemiology of B. melitensis in Southern Israel, implicating that rather than a common source, the transmission pattern of brucellosis among Bedouin communities is complex, predominantly local, and household-based. Further, genomic surveillance of B. melitensis is expected to inform future public health and veterinary interventions and clinical care.

Genetic Characterization of Brucella spp.: Whole Genome Sequencing-Based Approach for the Determination of Multiple Locus Variable Number Tandem Repeat Profiles

Brucellosis is an important zoonosis that is emerging in some regions of the world, gaining increased relevance with the inclusion of the causing agent Brucella spp. in the class B bioterrorism group. Until now, multi-locus VNTR Analysis (MLVA) based on 16 loci has been considered as the gold standard for Brucella typing. However, this methodology is laborious, and, with the rampant release of Brucella genomes, the transition from the traditional MLVA to whole genome sequencing (WGS)-based typing is on course. Nevertheless, in order to avoid a disruptive transition with the loss of massive genetic data obtained throughout the last decade and considering that the transition timings will vary considerably among different countries, it is important to determine WGS-based MLVA alleles of the nowadays sequenced genomes. On this regard, we aimed to evaluate the performance of a Python script that had been previously developed for the rapid in silico extraction of the MLVA alleles, by comparing it to the PCR-based MLVA procedure over 83 strains from different Brucella species. The WGS-based MLVA approach detected 95.3% of all possible 1,328 hits (83 strains×16 loci) and showed an agreement rate with the PCR-based MLVA procedure of 96.4% for MLVA-16. According to our dataset, we suggest the use of a minimal depth of coverage of ~50x and a maximum number of ~200 contigs as guiding “boundaries” for the future application of the script. In conclusion, the evaluated script seems to be a very useful and robust tool for the in silico determination of MLVA profiles of Brucella strains, allowing retrospective and prospective molecular epidemiological studies, which are important for maintaining an active epidemiological surveillance of brucellosis.

Comparison of Whole Genome Sequencing and Repetitive Element PCR for Multidrug- Resistant Pseudomonas aeruginosa Strain Typing

Hospital-acquired infections pose significant costly global challenges to patient care. Rapid and sensitive methods to identify potential outbreaks are integral to infection control measures. Whole genome sequencing (WGS)-based bacterial strain typing provides higher discriminatory power over standard nucleotide banding pattern-based methods like repetitive sequence-based PCR (rep-PCR). However, integration of WGS into clinical epidemiology is limited by the lack of consensus in methodology and data analysis/interpretation. In this study, WGS was performed on 22 multidrug-resistant Pseudomonas aeruginosa (MDR-PA) genomic DNA using the Illumina MiSeq platform. Resulting high quality reads were analyzed for phylogenetic relatedness using whole genome multi-locus sequence typing (wgMLST)-based BIOMÉRIEUX EPISEQ CS (EpiSeq; Durham, NC) and Single Nucleotide Variant PHYLogenomics (SNVPhyl). WGS-based results were compared to conventional MLST and archived rep-PCR results. Rep-PCR identified three independent clonal clusters of MDR-PA. Only one clonal cluster identified by rep-PCR, an endemic strain within the pediatric cystic fibrosis population at Texas Children's Hospital, was concordantly identified using wgMLST (EpiSeq) and SNVPhyl. Results were highly consistent between the three sequence-based analyses (conventional MLST, wgMLST, and SNVPhyl), and these results remained consistent with the addition of 74 MDR-PA genomes. These WGS-based methods provided greater resolution for strain discrimination than rep-PCR or standard MLST classification, and the ease of use of EpiSeq renders it clinically viable for analysis, interpretation, and reporting of WGS-based strain typing.

Tracking the Distribution of Brucella abortus in Egypt Based on Core Genome SNP Analysis and In Silico MLVA-16

Brucellosis, caused by the bacteria of the genus Brucella, is one of the most neglected common zoonotic diseases globally with a public health significance and a high economic loss among the livestock industry worldwide. Since little is known about the distribution of B. abortus in Egypt, a total of 46 B. abortus isolates recovered between 2012–2020, plus one animal isolate from 2006, were analyzed by examining the whole core genome single nucleotide polymorphism (cgSNP) in comparison to the in silico multilocus variable number of tandem repeat analysis (MLVA). Both cgSNP analysis and MLVA revealed three clusters and one isolate only was distantly related to the others. One cluster identified a rather widely distributed outbreak strain which is repeatedly occurring for at least 16 years with marginal deviations in cgSNP analysis. The other cluster of isolates represents a rather newly introduced outbreak strain. A separate cluster comprised RB51 vaccine related strains, isolated from aborted material. The comparison with MLVA data sets from public databases reveals one near relative from Argentina to the oldest outbreak strain and a related strain from Spain to a newly introduced outbreak strain in Egypt. The distantly related isolate matches with a strain from Portugal in the MLVA profile. Based on cgSNP analysis the oldest outbreak strain clusters with strains from the UK. Compared to the in silico analysis of MLVA, cgSNP analysis using WGS data provides a much higher resolution of genotypes and, when correlated to the associated epidemiological metadata, cgSNP analysis allows the differentiation of outbreaks by defining different outbreak strains. In this respect, MLVA data are error-prone and can lead to incorrect interpretations of outbreak events.

Whole Genome Sequence Analysis of Brucella melitensis Phylogeny and Virulence Factors

Brucellosis has a wide range of clinical severity in humans that remains poorly understood. Whole genome sequencing (WGS) analysis may be able to detect variation in virulence genes. We used Brucella melitensis sequences in the NCBI Sequence Read Archive (SRA) database to assemble 248 whole genomes, and additionally, assembled 27 B. melitensis genomes from samples of human patients in Southern Israel. We searched the 275 assembled genomes for the 43 B. melitensis virulence genes in the Virulence Factors of Pathogenic Bacteria Database (VFDB) and 10 other published putative virulence genes. We explored pan-genome variation across the genomes and in a pilot analysis, explored single nucleotide polymorphism (SNP) variation among the ten putative virulence genes. More than 99% of the genomes had sequences for all Brucella melitensis virulence genes included in the VFDB. The 10 other virulence genes of interest were present across all the genomes, but three of these genes had SNP variation associated with particular Brucella melitensis genotypes. SNP variation was also seen within the Israeli genomes obtained from a small geographic region. While the Brucella genome is highly conserved, this novel and large whole genome study of Brucella demonstrates the ability of whole genome and pan-genome analysis to screen multiple genomes and identify SNP variation in both known and novel virulence genes that could be associated with differential disease virulence. Further development of whole genome techniques and linkage with clinical metadata on disease outcomes could shed light on whether such variation in the Brucella genome plays a role in pathogenesis.

Campylobacter coli isolated in Brazil typed by core genome Multilocus Sequence Typing shows high genomic diversity in a global context

Campylobacter has been one of the most common causative agent of bacterial food-borne gastroenteritis in humans worldwide. However, in Brazil the campylobacteriosis has been a neglected disease and there is insufficient data to estimate the incidence of this pathogen in the country.

AIMS

The current study aimed to determine the phylogenetic relationships among Campylobacter coli strains isolated in Brazil and to compare them with international Campylobacter isolates available in some public databases.

METHODS AND RESULTS

A total of 63C. coli strains isolated in Brazil were studied. The MLST analysis showed 18 different STs including three STs not yet described in the PubMLST database. The cgMLST allocated the Brazilian strains studied into five main clusters and each cluster comprised groups of strains with nearly identical cgMLST profiles and with significant genetic distance observed among the distinct clusters. The comparison of the Brazilian strains with 3401 isolates from different countries showed a wide distribution of these strains isolated in this country.

CONCLUSIONS

The results showed a high similarity among some strains studied and a wide distribution of the Brazilian strains when compared to isolates from different countries, which is an interesting data set since it showed a high genetic diversity of these strains from Brazil in a global context. This study contributed for a better genomic characterization of C. coli strains isolated in Brazil and provided important information about the diversity of this clinically-relevant pathogen.

Development and Validation of a Burkholderia pseudomallei Core Genome Multilocus Sequence Typing Scheme To Facilitate Molecular Surveillance

Burkholderia pseudomallei causes the severe disease melioidosis. Whole-genome sequencing (WGS)-based typing methods currently offer the highest resolution for molecular investigations of this genetically diverse pathogen. Still, its routine application in diagnostic laboratories is limited by the need for high computing power, bioinformatic skills, and variable bioinformatic approaches, with the latter affecting the results. We therefore aimed to establish and validate a WGS-based core genome multilocus sequence typing (cgMLST) scheme, applicable in routine diagnostic settings. A soft defined core genome was obtained by challenging the B. pseudomallei reference genome K96243 with 469 environmental and clinical genomes, resulting in 4,221 core and 1,351 accessory targets. The scheme was validated with 320 WGS data sets. We compared our novel typing scheme with single nucleotide polymorphism-based approaches investigating closely and distantly related strains. Finally, we applied our scheme for tracking the environmental source of a recent infection. The validation of the scheme detected >95% good cgMLST target genes in 98.4% of the genomes. Comparison with existing typing methods revealed very good concordance. Our scheme proved to be applicable to investigating not only closely related strains but also the global B. pseudomallei population structure. We successfully utilized our scheme to identify a sugarcane field as the presumable source of a recent melioidosis case. In summary, we developed a robust cgMLST scheme that integrates high resolution, maximized standardization, and fast analysis for the nonbioinformatician. Our typing scheme has the potential to serve as a routinely applicable classification system in B. pseudomallei molecular epidemiology.

[Association of sepM gene mutation with mutacin Ⅳ production by Streptococcus mutans]

OBJECTIVE

To investigate the types of sepM gene mutations and their distribution in clinical isolates of Streptococcus mutans (S. mutans) and explore the association of sepM gene mutation with the capacity of mutacin Ⅳ production by S. mutans.

OBJECTIVE

We assessed the capacity of mutacin Ⅳ production in 80 clinical isolates of S. mutans using an inhibition zone assay. The minimum spanning tree and phylogenetic tree of these isolates were constructed using core genome multilocus sequence typing and maximum likelihood method, respectively. GeneMarkS software was used to predict the coding genes of these isolates, and the predicted genes were blasted against the sepM gene sequence of the reference genome UA159 to determine sepM gene mutations and their distribution characteristics in the clinical isolates. The mutation types affecting mutacin Ⅳ production were identified by analyzing the differentially distributed mutations between mutacin Ⅳ-producing isolates and mutacin Ⅳ-free isolates and by comparing the inhibition zones between isolates with sepM gene mutations and those without mutations.

OBJECTIVE

Among the 80 clinical isolates of S. mutans, 25 isolates were capable of mutacin Ⅳ production and 55 did not produce mutacin Ⅳ. The minimum spanning tree showed that the allelic differences were less among the mutacin Ⅳproducing isolates than among the mutacin Ⅳ-free isolates, and the origins of the mutacin Ⅳ-producing isolates were similar. We identified a total of 34 single base mutations in the 80 clinical isolates, and among them, C31T (P=0.001), G533A (P < 0.001), C756T (P=0.025), and C1036T (P=0.003) showed significant differential distributions between the mutacin Ⅳ-producing and mutacin Ⅳ-free isolates. These differentially distributed mutations were positively correlated with the capacity of mutacin Ⅳ production of the bacteria.

OBJECTIVE

sepM gene mutations that affect the capacity of mutacin Ⅳ production are present in the clinical strains of S. mutans.

A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis

Whole-genome sequencing (WGS) has been established for bacterial subtyping and is regularly used to study pathogen transmission, to investigate outbreaks, and to perform routine surveillance. Core-genome multilocus sequence typing (cgMLST) is a bacterial subtyping method that uses WGS data to provide a high-resolution strain characterization. This study aimed at developing a novel cgMLST scheme for Bacillus anthracis, a notorious pathogen that causes anthrax in livestock and humans worldwide. The scheme comprises 3,803 genes that were conserved in 57 B. anthracis genomes spanning the whole phylogeny. The scheme has been evaluated and applied to 584 genomes from 50 countries. On average, 99.5% of the cgMLST targets were detected. The cgMLST results confirmed the classical canonical single-nucleotide-polymorphism (SNP) grouping of B. anthracis into major clades and subclades. Genetic distances calculated based on cgMLST were comparable to distances from whole-genome-based SNP analysis with similar phylogenetic topology and comparable discriminatory power. Additionally, the application of the cgMLST scheme to anthrax outbreaks from Germany and Italy led to a definition of a cutoff threshold of five allele differences to trace epidemiologically linked strains for cluster typing and transmission analysis. Finally, the association of two clusters of B. anthracis with human cases of injectional anthrax in four European countries was confirmed using cgMLST. In summary, this study presents a novel cgMLST scheme that provides high-resolution strain genotyping for B. anthracis. This scheme can be used in parallel with SNP typing methods to facilitate rapid and harmonized interlaboratory comparisons, essential for global surveillance and outbreak analysis. The scheme is publicly available for application by users, including those with little bioinformatics knowledge.

Whole-Genome Sequencing for Tracing the Genetic Diversity of Brucella abortus and Brucella melitensis Isolated from Livestock in Egypt

Brucellosis is a highly contagious zoonosis that occurs worldwide. Whole-genome sequencing (WGS) has become a widely accepted molecular typing method for outbreak tracing and genomic epidemiology of brucellosis. Twenty-nine Brucella spp. (eight B. abortus biovar 1 and 21 B. melitensis biovar 3) were isolated from lymph nodes, milk, and fetal abomasal contents of infected cattle, buffaloes, sheep, and goats originating from nine districts in Egypt. The isolates were identified by microbiological methods and matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS). Differentiation and genotyping were confirmed using multiplex PCR. Illumina MiSeq^® was used to sequence the 29 Brucella isolates. Using MLST typing, ST11 and ST1 were identified among B. melitensis and B. abortus, respectively. Brucella abortus and B. melitensis isolates were divided into two main clusters (clusters 1 and 2) containing two and nine distinct genotypes by core-genome SNP analysis, respectively. The genotypes were irregularly distributed over time and space in the study area. Both Egyptian B. abortus and B. melitensis isolates proved to be genomically unique upon comparison with publicly available sequencing from strains of neighboring Mediterranean, African, and Asian countries. The antimicrobial resistance mechanism caused by mutations in rpoB, gyrA, and gyrB genes associated with rifampicin and ciprofloxacin resistance were identified. To the best of our knowledge, this is the first study investigating the epidemiology of Brucella isolates from livestock belonging to different localities in Egypt based on whole genome analysis.

Direct identification and molecular characterization of zoonotic hazards in raw milk by metagenomics using Brucella as a model pathogen

Metagenomics is a valuable diagnostic tool for enhancing microbial food safety because (i) it enables the untargeted detection of pathogens, (ii) it is fast since primary isolation of micro-organisms is not required, and (iii) it has high discriminatory power allowing for a detailed molecular characterization of pathogens. For shotgun metagenomics, total nucleic acids (NAs) are isolated from complex samples such as foodstuff. Along with microbial NAs, high amounts of matrix NAs are extracted that might outcompete microbial NAs during next-generation sequencing and compromise sensitivity for the detection of low abundance micro-organisms. Sensitive laboratory methods are indispensable for detecting highly pathogenic foodborne bacteria like Brucella spp., because a low infectious dose is sufficient to cause human disease through the consumption of contaminated dairy or meat products. In our study, we applied shotgun metagenomic sequencing for the identification and characterization of Brucella spp. in artificially and naturally contaminated raw milk from various ruminant species. With the depletion of eukaryotic cells prior to DNA extraction, Brucella was detectable at 10 bacterial cells ml⁻¹, while at the same time microbiological culture and isolation of the fastidious bacteria commonly failed. Moreover, we were able to retrieve the genotype of a Brucella isolate from a metagenomic dataset, indicating the potential of metagenomics for outbreak investigations using SNPs and core-genome multilocus sequence typing (cgMLST). To improve diagnostic applications, we developed a new bioinformatics approach for strain prediction based on SNPs to identify the correct species and define a certain strain with only low numbers of genus-specific reads per sample. This pipeline turned out to be more sensitive and specific than Mash Screen. In raw milk samples, we simultaneously detected numerous other zoonotic pathogens, antimicrobial resistance genes and virulence factors. Our study showed that metagenomics is a highly sensitive tool for biological risk assessment of foodstuffs, particularly when pathogen isolation is hazardous or challenging.

Emerging diversity and ongoing expansion of the genus Brucella

Remarkable genetic diversity and breadth of host species has been uncovered in the Brucella genus over the past decade, fundamentally changing our concept of what it means to be a Brucella. From ocean fishes and marine mammals, to pond dwelling amphibians, forest foxes, desert rodents, and cave-dwelling bats, Brucella have revealed a variety of previously unknown niches. Classical microbiological techniques have been able to help us classify many of these new strains but at times have limited our ability to see the true relationships among or within species. The closest relatives of Brucella are soil bacteria and the adaptations of Brucella spp. to live intracellularly suggest that the genus has evolved to live in vertebrate hosts. Several recently discovered species appear to have phenotypes that are intermediate between soil bacteria and core Brucella, suggesting that they may represent ancestral traits that were subsequently lost in the traditional species. Remarkably, the broad relationships among Brucella species using a variety of sequence and fragment-based approaches have been upheld when using comparative genomics with whole genomes. Nonetheless, genomes are required for fine-scale resolution of many of the relationships and for understanding the evolutionary history of the genus. We expect that the coming decades will reveal many more hosts and previously unknown diversity in a wide range of environments.

Molecular characteristics of Brucella melitensis isolates from humans in Qinghai Province, China

BACKGROUND

The prevalence of human brucellosis in Qinghai Province of China has been increasing rapidly, with confirmed cases distributed across 31 counties. However, the epidemiology of brucellosis transmission has not been fully elucidated. To characterize the infecting strains isolated from humans, multiple-locus variable-number tandem repeats analysis (MLVA) and whole-genome single-nucleotide polymorphism (SNP)-based approaches were employed.

METHODS

Strains were isolated from two males blood cultures that were confirmed Brucella melitensis positive following biotyping and MLVA. Genomic DNA was extracted from these two strains, and whole-genome sequencing was performed. Next, SNP-based phylogenetic analysis was performed to compare the two strains to 94 B. melitensis strains (complete genome and draft genome) retrieved from online databases.

RESULTS

The two Brucella isolates were identified as B. melitensis biovar 3 (QH2019001 and QH2019005) following conventional biotyping and were found to have differences in their variable number tandem repeats (VNTRs) using MLVA-16. Phylogenetic examination assigned the 96 strains to five genotype groups, with QH2019001 and QH2019005 assigned to the same group, but different subgroups. Moreover, the QH2019005 strain was assigned to a new subgenotype, IIj, within genotype II. These findings were then combined to determine the geographic origin of the two Brucella strains.

CONCLUSIONS

Utilizing a whole-genome SNP-based approach enabled differences between the two B. melitensis strains to be more clearly resolved, and facilitated the elucidation of their different evolutionary histories. This approach also revealed that QH2019005 is a member of a new subgenotype (IIj) with an ancient origin in the eastern Mediterranean Sea.

Whole Genome Sequence Analysis of Brucella abortus Isolates from Various Regions of South Africa

The availability of whole genome sequences in public databases permits genome-wide comparative studies of various bacterial species. Whole genome sequence-single nucleotide polymorphisms (WGS-SNP) analysis has been used in recent studies and allows the discrimination of various Brucella species and strains. In the present study, 13 Brucella spp. strains from cattle of various locations in provinces of South Africa were typed and discriminated. WGS-SNP analysis indicated a maximum pairwise distance ranging from 4 to 77 single nucleotide polymorphisms (SNPs) between the South African Brucella abortus virulent field strains. Moreover, it was shown that the South African B. abortus strains grouped closely to B. abortus strains from Mozambique and Zimbabwe, as well as other Eurasian countries, such as Portugal and India. WGS-SNP analysis of South African B. abortus strains demonstrated that the same genotype circulated in one farm (Farm 1), whereas another farm (Farm 2) in the same province had two different genotypes. This indicated that brucellosis in South Africa spreads within the herd on some farms, whereas the introduction of infected animals is the mode of transmission on other farms. Three B. abortus vaccine S19 strains isolated from tissue and aborted material were identical, even though they originated from different herds and regions of South Africa. This might be due to the incorrect vaccination of animals older than the recommended age of 4-8 months or might be a problem associated with vaccine production.

Tracking the Origin of Austrian Human Brucellosis Cases Using Whole Genome Sequencing

Brucellosis is a zoonotic disease caused by Brucella spp. and a major concern for livestock. Most human cases are caused by B. melitensis and clinical presentation is usually a mild febrile illness. However, treatment failure is frequent and more severe complications can occur. In Austria, every human brucellosis is investigated to determine whether it was imported from endemic areas or is the sign of an undetected autochthonous transmission. For this study, 21 B. melitensis strains isolated in Austria between 2005 and 2019 were collected, 17 strains from 15 different patients and four strains from cattle. Whole genome sequencing combined with core-genome MLST analysis was used to characterize these strains. A cluster of seven isolates from 2018 (three human and four cattle isolates) was identified, with fewer than two allelic differences. They corresponded to the only Austrian B. melitensis outbreak that happened over the past 15 years. The other 12 Austrian brucellosis cases were single cases, and geographical origins were available for 8/12. Genomic data was used to locate probable geographical origins and compared with the results of the epidemiological investigations. Austrian strains were compared with 67 published B. melitensis sequences available on NCBI. The result of genomic analysis matched for 7/8 cases with documented conclusion of the epidemiological investigation. Genome analysis also pointed to the geographical origin for three of the four cases with missing epidemiological data. Strains from six cases were grouped together (<40 allelic differences) with 4/6 cases imported from the Balkans. Additional B. melitensis isolates from Serbian animals were analyzed and grouped with this branch, suggesting frequent importation from Balkan countries to Austria. Overall, this study highlights the specificities of human brucellosis in Austria. It also underlines the value of whole genome sequencing as a tool to investigate brucellosis cases, allowing to identify and investigate outbreaks but also to support epidemiological investigation of imported cases. However, the reliability of such methods depends on the number of strains for comparison, which can be challenging in low incidence countries. Increasing the availability of published sequences with documented geographical origins would help establishing genomic-based methods for investigating brucellosis cases.

A Vibrio cholerae Core Genome Multilocus Sequence Typing Scheme To Facilitate the Epidemiological Study of Cholera

Core genome multilocus sequence typing (cgMLST) has gained popularity in recent years in epidemiological research and subspecies-level classification. cgMLST retains the intuitive nature of traditional MLST but offers much greater resolution by utilizing significantly larger portions of the genome. Here, we introduce a cgMLST scheme for Vibrio cholerae, a bacterium abundant in marine and freshwater environments and the etiologic agent of cholera. A set of 2,443 core genes ubiquitous in V. cholerae were used to analyze a comprehensive data set of 1,262 clinical and environmental strains collected from 52 countries, including 65 newly sequenced genomes in this study. We established a sublineage threshold based on 133 allelic differences that creates clusters nearly identical to traditional MLST types, providing backwards compatibility to new cgMLST classifications. We also defined an outbreak threshold based on seven allelic differences that is capable of identifying strains from the same outbreak and closely related isolates that could give clues on outbreak origin. Using cgMLST, we confirmed the South Asian origin of modern epidemics and identified clustering affinity among sublineages of environmental isolates from the same geographic origin. Advantages of this method are highlighted by direct comparison with existing classification methods, such as MLST and single-nucleotide polymorphism-based methods. cgMLST outperforms all existing methods in terms of resolution, standardization, and ease of use. We anticipate this scheme will serve as a basis for a universally applicable and standardized classification system for V. cholerae research and epidemiological surveillance in the future. This cgMLST scheme is publicly available on PubMLST (https://pubmlst.org/vcholerae/).IMPORTANCE Toxigenic Vibrio cholerae isolates of the O1 and O139 serogroups are the causative agents of cholera, an acute diarrheal disease that plagued the world for centuries, if not millennia. Here, we introduce a core genome multilocus sequence typing scheme for V. cholerae Using this scheme, we have standardized the definition for subspecies-level classification, facilitating global collaboration in the surveillance of V. cholerae In addition, this typing scheme allows for quick identification of outbreak-related isolates that can guide subsequent analyses, serving as an important first step in epidemiological research. This scheme is also easily scalable to analyze thousands of isolates at various levels of resolution, making it an invaluable tool for large-scale ecological and evolutionary analyses.

Brucella Genomics: Macro and Micro Evolution

Brucella organisms are responsible for one of the most widespread bacterial zoonoses, named brucellosis. The disease affects several species of animals, including humans. One of the most intriguing aspects of the brucellae is that the various species show a ~97% similarity at the genome level. Still, the distinct Brucella species display different host preferences, zoonotic risk, and virulence. After 133 years of research, there are many aspects of the Brucella biology that remain poorly understood, such as host adaptation and virulence mechanisms. A strategy to understand these characteristics focuses on the relationship between the genomic diversity and host preference of the various Brucella species. Pseudogenization, genome reduction, single nucleotide polymorphism variation, number of tandem repeats, and mobile genetic elements are unveiled markers for host adaptation and virulence. Understanding the mechanisms of genome variability in the Brucella genus is relevant to comprehend the emergence of pathogens.

Evolutionary history and current distribution of the West Mediterranean lineage of Brucella melitensis in Italy

Ovine and caprine brucellosis, caused by Brucella melitensis, is one of the world’s most widespread zoonoses and is a major cause of economic losses in domestic ruminant production. In Italy, the disease remains endemic in several southern provinces, despite an ongoing brucellosis eradication programme. In this study, we used whole-genome sequencing to detail the genetic diversity of circulating strains, and to examine the origins of the predominant sub-lineages of B. melitensis in Italy. We reconstructed a global phylogeny of B. melitensis, strengthened by 339 new whole-genome sequences, from Italian isolates collected from 2011 to 2018 as part of a national livestock surveillance programme. All Italian strains belonged to the West Mediterranean lineage, which further divided into two major clades that diverged roughly between the 5th and 7th centuries. We observed that Sicily serves as a brucellosis burden hotspot, giving rise to several distinct sub-lineages. More than 20 putative outbreak clusters of ovine and caprine brucellosis were identified, several of which persisted over the 8 year survey period despite an aggressive brucellosis eradication campaign. While the outbreaks in Central and Northern Italy were generally associated with introductions of single clones of B. melitensis and their subsequent dissemination within neighbouring territories, we observed weak geographical segregation of genotypes in the southern regions. Biovar determination, recommended in routine analysis of all Brucella strains by the World Organisation for Animal Health (OIE), could not discriminate among the four main global clades. This demonstrates a need for updating the guidelines used for monitoring B. melitensis transmission and spread, both at the national and international level, and to include whole-genome-based typing as the principal method for identification and tracing of brucellosis outbreaks.

A Core Genome Multilocus Sequence Typing Scheme for Streptococcus mutans

Streptococcus mutans is one of the primary pathogens responsible for the development of dental caries. Recent whole-genome sequencing (WGS)-based core genome multilocus sequence typing (cgMLST) approaches have been employed in epidemiological studies of specific human pathogens. However, this approach has not been reported in studies of S. mutans Here, we therefore developed a cgMLST scheme for S. mutans We surveyed 199 available S. mutans genomes as a means of identifying cgMLST targets, developing a scheme that incorporated 594 targets from the S. mutans UA159 reference genome. Sixty-eight sequence types (STs) were identified in this cgMLST scheme (cgSTs) in 80 S. mutans isolates from 40 children that were sequenced in this study, compared to 35 STs identified by multilocus sequence typing (MLST). Fifty-six cgSTs (82.35%) were associated with a single isolate based on our cgMLST scheme, which is significantly higher than in the MLST scheme (11.43%). In addition, 58.06% of all MLST profiles with ≥2 isolates were further differentiated by our cgMLST scheme. Topological analyses of the maximum likelihood phylogenetic trees revealed that our cgMLST scheme was more reliable than the MLST scheme. A minimum spanning tree of 145 S. mutans isolates from 10 countries developed based upon the cgMLST scheme highlighted the diverse population structure of S. mutans This cgMLST scheme thus offers a new molecular typing method suitable for evaluating the epidemiological distribution of this pathogen and has the potential to serve as a benchmark for future global studies of the epidemiological nature of dental caries.IMPORTANCE Streptococcus mutans is regarded as a major pathogen responsible for the onset of dental caries. S. mutans can transmit among people, especially within families. In this study, we established a new epidemiological approach to S. mutans classification. This approach can effectively differentiate among closely related isolates and offers superior reliability relative to that of the traditional MLST molecular typing method. As such, it has the potential to better support effective public health strategies centered around this bacterium that are aimed at preventing and treating dental caries.

MLVA-16 Genotyping of Brucella abortus and Brucella melitensis Isolates from Different Animal Species in Egypt: Geographical Relatedness and the Mediterranean Lineage

Brucellosis is a common zoonotic disease in Egypt. However, there are limited data available on the genetic diversity of brucellae circulating in Egypt and other Mediterranean areas. One hundred and nine Brucella (B.) strains were isolated from different animal species in thirteen Egyptian governorates. Multi-locus variable number tandem repeats (VNTRs) analysis (MLVA-16) was employed to determine the geographical relatedness and the genetic diversity of a panel of selected Egyptian strains (n = 69), with strains originating from Italy (n = 49), Portugal (n = 52), Greece (n = 63), and Tunisia (n = 4). Egyptian B. melitensis strains clustered into two main clusters containing 21 genotypes. Egyptian B. abortus strains clustered into three main clusters containing nine genotypes. The genotypes were irregularly distributed over time and space in the study area. Egyptian strains of B. melitensis showed MLVA-16 patterns closer to that of Italian strains. Egyptian B. abortus strains isolated from cattle share the same genotype with strains from Portugal and similar to strains from Italy with low genetic diversity. Strains with similar MLVA patterns isolated from different governorates highlight the movement of the pathogen among governorates. Hence, it may also reflect the long endemicity of brucellosis in Egypt with earlier dispersal of types and great local genetic diversity. Open markets may contribute to cross-species transmission and dissemination of the new types nationwide. The presence of West Mediterranean lineages of B. melitensis and relatedness of B. abortus strains from the studied countries is a result of the socio-historical connections among the Mediterranean countries. Transnational eradication of brucellosis in the Mediterranean basin is highly demanded.

A Core Genome Multilocus Sequence Typing Scheme for Pseudomonas aeruginosa

Pseudomonas aeruginosa is a ubiquitous microorganism and an important opportunistic pathogen responsible for a broad spectrum of infections mainly in immunosuppressed and critically ill patients. Molecular investigations traditionally rely on pulsed field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). In this work we propose a core genome multilocus sequence typing (cgMLST) scheme for P. aeruginosa, a methodology that combines traditional MLST principles with whole genome sequencing data. All publicly available complete P. aeruginosa genomes, representing the diversity of this species, were used to establish a cgMLST scheme targeting 2,653 genes. The scheme was then tested using genomes available at contig, chromosome and scaffold levels. The proposed cgMLST scheme for P. aeruginosa typed over 99% (2,314/2,325) of the genomes available for this study considering at least 95% of the cgMLST target genes present. The absence of a certain number gene targets at the threshold considered for both the creation and validation steps due to low genome sequence quality is possibly the main reason for this result. The cgMLST scheme was compared with previously published whole genome single nucleotide polymorphism analysis for the characterization of the population structure of the epidemic clone ST235 and results were highly similar. In order to evaluate the typing resolution of the proposed scheme, collections of isolates belonging to two important STs associated with cystic fibrosis, ST146 and ST274, were typed using this scheme, and ST235 isolates associated with an outbreak were evaluated. Besides confirming the relatedness of all the isolates, earlier determined by MLST, the higher resolution of cgMLST denotes that it may be suitable for surveillance programs, overcoming possible shortcomings of classical MLST. The proposed scheme is publicly available at: https://github.com/BioinformaticsHIAEMolecularMicrobiology/cgMLST-Pseudomonas-aeruginosa.

Molecular epidemiology of Salmonella Infantis in Europe: insights into the success of the bacterial host and its parasitic pESI-like megaplasmid

Salmonella Infantis is one of the five serovars most frequently causing human salmonellosis in Europe, mainly associated with poultry. A clone harbouring a conjugative plasmid of emerging S. Infantis (pESI)-like megaplasmid, carrying multidrug resistant (MDR) and extended-spectrum beta-lactamases (ESBL) genes, has spread in the Italian broiler chicken industry also causing human illness. This work is aimed at elucidating the molecular epidemiology of S. Infantis and pESI-like in Europe using whole-genome sequencing and bioinformatics analysis, and to investigate the genetic relatedness of S. Infantis clones and pESI-like from animals, meat, feed and humans provided by institutions of nine European countries. Two genotyping approaches were used: chromosome or plasmid SNP-based analysis and the minimum spanning tree (MST) algorithm based on core-genome multilocus sequence typing (cgMLST). The European S. Infantis population appeared heterogeneous, with different genetic clusters defined at core-genome level. However, pESI-like variants present in 64.1 % of the isolates were more genetically homogeneous and capable of infecting different clonal lineages in most of the countries. Two different pESI-like with ESBL genes (n=82) were observed: bla_CTX-M-1-positive in European isolates and bla_CTX-M-65-positive in American isolates (study outgroup). Both variants had toxin-antitoxin systems, resistance genes towards tetracyclines, trimethoprim, sulphonamides and aminoglycosides, heavy metals (merA) and disinfectants (qacEΔ). Worryingly, 66 % of the total isolates studied presented different gyrA chromosomal point mutations associated with (fluoro)quinolone resistance (MIC range 0.125–0.5 mg/L), while 18 % displayed transferable macrolide resistance mediated by mph, mef and erm(B) genes. Proper intervention strategies are needed to prevent further dissemination/transmission of MDR S. Infantis and pESI-like along the food chain in Europe.

Genome-scale approach to study the genetic relatedness among Brucella melitensis strains

Brucellosis is an important zoonotic disease that affects both humans and animals. To date, laboratory surveillance is still essentially based on the traditional MLVA-16 methodology and the associated epidemiological information is frequently scarce. Our goal was to contribute to the improvement of Brucella spp. surveillance through the implementation of a whole genome sequencing (WGS) approach. We created a curated ready-to-use species-specific wgMLST scheme enrolling a panel of 2656 targets (http://doi.org/10.5281/zenodo.3575026) and used this schema to perform a retrospective analysis of the genetic relatedness among B. melitensis strains causing human infection in Portugal (a country where brucellosis is an endemic disease) from 2010 to 2018. The strains showed a phylogenetic clustering within genotype II (25 out of 36) and IV (4 out of 36), and shared clades with strains isolated from countries with which Portugal has intense food trading, tourism and similar eating habits, such as Spain, Italy and Greece. In addition, our results point to the identification of strong associations between B. melitensis strains, likely underlying missed "outbreaks" as 22 out of the 36 strains showed genetic linkage with others. In fact, the applied gene-by-gene approach grouped these strains into six genetic clusters each one containing putative epidemiological links. Nevertheless, more studies will be needed in order to define the appropriate range of cut-offs (probable non-static cut-offs) that best illustrate the association between genetic linkage and epidemiological information and may serve as alerts for the health authorities. The release of this freely available and scalable schema contributes to the required technological transition for laboratorial surveillance of brucellosis and will facilitate the assessment of ongoing and future outbreaks in order to prevent the transmission spread.

Whole-genome sequencing of rough Brucella melitensis in China provides insights into its genetic features

Brucella spp. can cause the zoonosis brucellosis, which affects public health and safety and even economic development. B. melitensis has a smooth phenotype, while 28 B. melitensis isolates had a rough phenotype in 2018. In this study, rough phenotype detection and whole genome sequencing methods were used to analyze the genetic features of rough B. melitensis. A drug susceptibility test was also performed. The results showed that the rough B. melitensis strains originated from strains isolated in China rather than from foreign strains. Furthermore, an MS tree showed that 9 complexes to be epidemic in China. For the rough B. melitensis strains, expression of the metabolic function genes varied in the earlier stages of evolution compared to the cellular process and signalling function genes. Expression of some transcriptional regulatory factors also varied in the later stages of evolution, and compared to MFS transporter genes, ABC transporter genes varied in the earlier stages. Moreover, as there was no significant difference in rifampicin, doxycycline and streptomycin susceptibility between the smooth and rough B. melitensis strains, treatment of brucellosis was not affected by strain type. This study provided important information for understanding the genetics and evolution of rough B. melitensis in China.

Detailed Evaluation of Data Analysis Tools for Subtyping of Bacterial Isolates Based on Whole Genome Sequencing: Neisseria meningitidis as a Proof of Concept

Whole genome sequencing is increasingly recognized as the most informative approach for characterization of bacterial isolates. Success of the routine use of this technology in public health laboratories depends on the availability of well-characterized and verified data analysis methods. However, multiple subtyping workflows are now often being used for a single organism, and differences between them are not always well described. Moreover, methodologies for comparison of subtyping workflows, and assessment of their performance are only beginning to emerge. Current work focuses on the detailed comparison of WGS-based subtyping workflows and evaluation of their suitability for the organism and the research context in question. We evaluated the performance of pipelines used for subtyping of Neisseria meningitidis, including the currently widely applied cgMLST approach and different SNP-based methods. In addition, the impact of the use of different tools for detection and filtering of recombinant regions and of different reference genomes were tested. Our benchmarking analysis included both assessment of technical performance of the pipelines and functional comparison of the generated genetic distance matrices and phylogenetic trees. It was carried out using replicate sequencing datasets of high- and low-coverage, consisting mainly of isolates belonging to the clonal complex 269. We demonstrated that cgMLST and some of the SNP-based subtyping workflows showed very good performance characteristics and highly similar genetic distance matrices and phylogenetic trees with isolates belonging to the same clonal complex. However, only two of the tested workflows demonstrated reproducible results for a group of more closely related isolates. Additionally, results of the SNP-based subtyping workflows were to some level dependent on the reference genome used. Interestingly, the use of recombination-filtering software generally reduced the similarity between the gene-by-gene and SNP-based methodologies for subtyping of N. meningitidis. Our study, where N. meningitidis was taken as an example, clearly highlights the need for more benchmarking comparative studies to eventually contribute to a justified use of a specific WGS data analysis workflow within an international public health laboratory context.

Whole Genome Sequencing for Tracing Geographical Origin of Imported Cases of Human Brucellosis in Sweden

Human infections with Brucella melitensis are occasionally reported in Sweden, despite the fact that the national flocks of sheep and goats are officially free from brucellosis. The aim of our study was to analyze 103 isolates of B. melitensis collected from patients in Sweden between 1994 and 2016 and determine their putative geographic origin using whole genome sequencing (WGS)-based tools. The majority of the strains were assigned to East Mediterranean and African lineages. Both in silico Multiple Loci VNTR (Variable Number of Tandem Repeats) Analysis (MLVA) and core genome Multilocus Sequence Typing (cgMLST) analyses identified countries of the Middle East as the most probable source of origin of the majority of the strains. Isolates collected from patients with travel history to Iraq or Syria were often associated with genotypes from Turkey, as the cgMLST profiles from these countries clustered together. Sixty strains were located within a distance of 20 core genes to related genotypes from the publicly available database, and for eighteen isolates, the closest genotype was different by more than 50 loci. Our study showed that WGS based tools are effective in tracing back the geographic origin of infection of patients with unknown travel status, provided that public sequences from the location of the source are available.

A headache with surprising outcome: first case of brucellosis caused by Brucella suis biovar 1 in Germany

In July 2018, brucellosis was diagnosed in a German patient without a travel history to regions endemic for Brucella. Microbiological analysis, including whole-genome sequencing, revealed Brucella suis biovar 1 as the etiologic agent. Core-genome-based multilocus sequence-typing analysis placed the isolate in close proximity to strains originating from Argentina. Notably, despite a strong IgM response, the patient did not develop Brucella-specific IgG antibodies during infection. Here, we describe the clinical course of infection, the extensive epidemiological investigations, and discuss possible routes of transmission.

Phylogeography and epidemiology of Brucella suis biovar 2 in wildlife and domestic swine

Swine brucellosis due to Brucella suis biovar 2 (bv2) is enzootic in wild boar and hare in continental Europe and may cause major economic losses to the pig industry, mainly in free-ranged pig farms. The high nucleotide identity found among the B. suis biovar 2 isolates has long hindered the full understanding of the epidemiology and the phylogeography of the disease. Here, we used multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) and whole-genome analysis to identify single-nucleotide polymorphisms (SNPs) in order to gain insights from the largest B. suis bv2 dataset analyzed so far composed of domestic pigs and wildlife isolates collected throughout Europe since the 1970s. We found four major clades with a specific phylogeographic pattern. The Iberian clade contains isolates exclusively from the Iberian Peninsula. The Central European clade includes most isolates from France, Northern Italy, Switzerland and an important proportion of those of Northern Spain. The Eastern European clade clustered isolates from Croatia and Hungary mainly but also from areas of France, Germany, Italy and Poland. Finally, a separated Sardinian clade grouped three isolates from this island. At fine scale, MLVA demonstrated an endemic status of the infection in Europe and it allowed tracking a large outbreak formed by different farms from Spain linked to the same infection source. The whole genome SNP analysis showed that the strains form genetically distinct clades, shared between wild boar and pigs, in agreement with the MLVA clades. Interestingly, all hare isolates clustered together within two groups composed exclusively of wildlife isolates. Our results support the hypothesis that maintenance and spread of B. suis bv2 in Europe is a dynamic process linked to the natural expansion of wild boar as the main wild reservoir of the infection, while spread over long distances is found largely dependent on anthropogenic activities.

Epidemiological Typing of Serratia marcescens Isolates by Whole-Genome Multilocus Sequence Typing

Serratia marcescens is an opportunistic bacterial pathogen. It is notorious for its increasing antimicrobial resistance and its potential to cause outbreaks of colonization and infections, predominantly in neonatal intensive care units (NICUs). There, its spread requires rapid infection control response. To understand its spread, detailed molecular typing is key. We present a whole-genome multilocus sequence typing (wgMLST) method for S. marcescens Using a set of 299 publicly available whole-genome sequences (WGS), we developed an initial wgMLST system consisting of 9,377 gene loci. This included 1,455 loci occurring in all reference genomes and 7,922 accessory loci. This closed system was validated using three geographically diverse collections of S. marcescens consisting of 111 clinical isolates implicated in nosocomial dissemination events in three hospitals. The validation procedure showed a full match between epidemiological data and the wgMLST analyses. We set the cutoff value for epidemiological (non)relatedness at 20 different alleles, though for the majority of outbreak-clustered isolates, this difference was limited to 4 alleles. This shows that the wgMLST system for S. marcescens provides prospects for successful future monitoring for the epidemiological containment of this opportunistic pathogen.

Distribution of Brucella field strains isolated from livestock, wildlife populations, and humans in Italy from 2007 to 2015

Brucellosis is a major public health problem still prevalent as a neglected endemic zoonosis requiring proactive attention in many communities worldwide. The present study involved analysis of Brucella field strains submitted for typing to the Italian National Reference Laboratory for Brucellosis from 2007 to 2015. Strains were identified at the species and biovar levels by classic and molecular techniques according to the World Organisation for Animal Health Manual. In total, 5,784 strains were typed: 3,089 Brucella abortus (53.4%), 2,497 B. melitensis (43.2%), 10 B. ovis (0.2%), 181 B. suis (3.1%), and 7 B. ceti (0.1%). The 2,981 strains from cattle were typed as B. abortus biovars 1, 3, and 6 (90.1%) and B. melitensis biovar 3 (9.9%). The 318 strains from water buffalo were typed as B. abortus biovars 1, 3 (95.9%) and B. melitensis biovar 3 (4.1%). The 2,279 strains from sheep and goats were typed as B. abortus biovars 1 and 3 (4.3%); B. melitensis biovars 1, 3, (95.3%); and B. ovis (0.4%). The 173 strains from wild boar were typed as B. suis biovar 2 (98.3%) and B. melitensis biovar 3 (1.7%). The 11 strains from pigs were typed as B. suis biovar 2. The 13 strains from humans were typed as B. melitensis biovar 3. The two strains from horses were typed as B. abortus biovar 1, while the seven strains from dolphins were typed as B. ceti. This additional knowledge on the epidemiology of brucellosis in Italy may be useful to formulate policies and strategies for the control and eradication of the disease in animal populations. The animal species affected, biovars typed, geographical origins, and spatial distributions of isolates are herein analyzed and discussed.