BackgroundDiphtheria remains a major public health concern with multiple recent outbreaks around the world. Moreover, invasive non-toxigenic strains have emerged globally causing severe infections. A diphtheria epidemic in the former Soviet Union in the 1990s resulted in ~5000 deaths. In this study, we analysed the genome sequences of a collection of 93 C. diphtheriae strains collected during and after this outbreak (1996 – 2014) in a former Soviet State, Belarus to understand the evolutionary dynamics and virulence capacities of these strains.Results C. diphtheriae strains from Belarus belong to ten sequence types (STs). Two major clones, non-toxigenic ST5 and toxigenic ST8, encompassed 76% of the isolates that are associated with sore throat and diphtheria in patients, respectively. Core genomic diversity is limited within outbreak-associated ST8 with relatively higher mutation rates (8.9 × 10−7 substitutions per strain per year) than ST5 (5.6 × 10−7 substitutions per strain per year) where most of the diversity was introduced by recombination. A variation in the virulence gene repertoire including the presence of tox gene is likely responsible for pathogenic differences between different strains. However, strains with similar virulence potential can cause disease in some individuals and remain asymptomatic in others. Eight synonymous single nucleotide polymorphisms were observed between the tox genes of the vaccine strain PW8 and other toxigenic strains of ST8, ST25, ST28, ST41 and non-toxigenic tox gene-bearing (NTTB) ST40 strains. A single nucleotide deletion at position 52 in the tox gene resulted in the frameshift in ST40 isolates, converting them into NTTB strains.ConclusionsNon-toxigenic C. diphtheriae ST5 and toxigenic ST8 strains have been endemic in Belarus both during and after the epidemic in 1990s. A high vaccine coverage has effectively controlled diphtheria in Belarus; however, non-toxigenic strains continue to circulate in the population. Recombination is an important evolutionary force in shaping the genomic diversity in C. diphtheriae. However, the relative role of recombination and mutations in diversification varies between different clones.Electronic supplementary materialThe online version of this article (10.1186/s12864-017-4276-3) contains supplementary material, which is available to authorized users.

Corynebacterium ulcerans is an important zoonotic pathogen which is causing diphtheria-like disease in humans globally. In this study, the genomes of three recently isolated C. ulcerans strains, 4940, 2590 and BR-AD 2649, respectively from an asymptomatic carrier, a patient with pharyngitis and a canine host, were sequenced to investigate their virulence potential. A comparative analysis was performed including the published genome sequences of 16 other C. ulcerans isolates. C. ulcerans strains belong to two lineages; 13 strains are grouped together in lineage 1, and six strains comprise lineage 2. Consistent with the zoonotic nature of C. ulcerans infections, isolates from both the human and canine hosts clustered in both the lineages. Most of the strains possessed spaDEF and spaBC gene clusters along with the virulence genes cpp, pld, cwlH, nanH, rpfI, tspA and vsp1. The gene encoding Shiga-like toxin was only present in one strain, and 11 strains carried the tox gene encoding the diphtheria-like toxin. However, none of strains 4940, 2590 and BR-AD 2649 carried any toxin genes. These strains varied in the number of prophages in their genomes, which suggests that they play an important role in introducing diversity in C. ulcerans. The pan-genomic analyses revealed a variation in the number of membrane-associated and secreted proteins that may contribute to the variation in pathogenicity among different strains.

In this study, we evaluated a novel macroarray-based spoligotyping method for Corynebacterium diphtheriae strain typing. A total of 20 C. diphtheriae biotype gravis toxigenic isolates collected in Belarus from suspected foci of diphtheria infection (diphtheria cases, carriers, or contacts) were subjected to DNA fingerprinting. All strains had an identical ribotyping profile that was identified as ribotype 'Rossija' by comparison with the international ribotype database at the Institut Pasteur of Paris. A spoligotyping method based on simultaneous reverse-hybridization analysis of two CRISPR (clustered, regularly interspaced short palindromic repeats) loci differentiated these strains into three spoligotypes. Comparison of the spoligotyping results with the epidemiological linkage network helped us to resolve suspected links in the chains of transmission. To conclude, the C. diphtheriae spoligotyping method demonstrated its utility in the field study, in particular, underlining the importance of the use of both CRISPR loci. The generated discrete data can be presented in digital binary format and be easily exchanged between laboratories and stored in local and global databases.

There is a need for the method which helps to choose the less mutable immunogenic determinant for the design of recombinant or synthetic vaccines and ELISA test-systems. Our method based on the directional mutational pressure theory includes two steps: estimation of symmetric and asymmetric mutational pressure directions in a gene coding for a protein of interest; and selection of regions coding for its immunogenic determinants which are less prone to missense mutation occurrence and so, to immune escaping. Three original computer algorithms (“VVK Sliding Window”, “VVK VarInvar” and “VVK Protective Buffer” available via www.barkovsky.hotmail.ru) have been created to perform all the necessary calculations and tests. “VVK Sliding Window” calculates nucleotide usage in fourfold and twofold degenerated sites, as well as usage of missense, nonsense and synonymous sites for each kind of nucleotide mutation along the length of a coding region, while “VVK Protective Buffer” calculates those indexes in a set of sequences. “VVK VarInvar” calculates percentage of variable sites in a set of aligned sequences, as well as nucleotide usage in invariable sites. Our method has been tested on HIV1 gp120 protein and on diphtheria toxin. The less mutable epitopes have been found for both proteins. Finally, it has been shown that antibodies recognizing the less mutable epitope of gp120 can be found in 80.22% of HIV1-infected persons

Background: The reemergence of epidemic diphtheria in Belarus in 1990s has provided us with important information on the biology of the disease and the diversity of the causative agent Corynebacterium diphtheriae. Molecular investigations were conducted with the aim to analyze the genetic variability of C diphtheriae during the post-epidemic period.