Due to its rapid spread and association with the numerous outbreaks, the global spread of East Asian lineage of Mycobacterium tuberculosis strains presents a global concern. Although there were many attempts to describe its population structure, no consensus has been reached yet. To define unbiased classification that will facilitate future studies of this lineage, we analyzed the performance and congruence of eight different genotyping schemes based on phylogenetic analysis of 1,398 strains from 32 countries using whole-genome sequencing (WGS) data. We confirm that East Asian lineage comprises two major clades, designated proto-Beijing, which harbors unusual 43-signal spoligoprofile, and Beijing, with well-known spoligoprofile (deleted signals from 1 to 34). We show that different genotyping methods give high consistency results in description of ancient Beijing strains while the classification of modern Beijing strains is significantly divergent due to star-shaped phylogeny. Using WGS data we intersect different studies and for the first time provide balanced classification with well-defined major groups and their genetic markers. Our reconstructed phylogenetic tree can also be used for further analysis of epidemiologically important clusters and their ancestors as well as white spots of unclassified strains, which are prospective areas of research.

BackgroundTuberculosis (TB) poses a worldwide threat due to advancing multidrug-resistant strains and deadly co-infections with Human immunodeficiency virus. Today large amounts of Mycobacterium tuberculosis whole genome sequencing data are being assessed broadly and yet there exists no comprehensive online resource that connects M. tuberculosis genome variants with geographic origin, with drug resistance or with clinical outcome.DescriptionHere we describe a broadly inclusive unifying Genome-wide Mycobacterium tuberculosis Variation (GMTV) database, (http://mtb.dobzhanskycenter.org) that catalogues genome variations of M. tuberculosis strains collected across Russia. GMTV contains a broad spectrum of data derived from different sources and related to M. tuberculosis molecular biology, epidemiology, TB clinical outcome, year and place of isolation, drug resistance profiles and displays the variants across the genome using a dedicated genome browser. GMTV database, which includes 1084 genomes and over 69,000 SNP or Indel variants, can be queried about M. tuberculosis genome variation and putative associations with drug resistance, geographical origin, and clinical stages and outcomes.ConclusionsImplementation of GMTV tracks the pattern of changes of M. tuberculosis strains in different geographical areas, facilitates disease gene discoveries associated with drug resistance or different clinical sequelae, and automates comparative genomic analyses among M. tuberculosis strains.

Thirty urogenital Chlamydia trachomatis isolates collected in Moscow in 2005 were typed using newly developed molecular typing approaches: (1) multilocus sequence typing (MLST(7)) based on sequences of seven housekeeping genes (http://pubmlst.org/chlamydiales/), (2) MLST(5) based on the investigation of five target regions of the chlamydial genome and (3) ompA gene sequencing supplemented with three variable number tandem repeat (VNTR) loci of the genome. ompA typing divided all isolates into 11 groups with E serotype dominating, while MLST7, MLST5 and VNTR analysis divided them into eight, 20 and 18 groups, respectively. The discriminatory power of each method calculated using the Hunter-Gaston discriminatory index was found to be 0.83 for the ompA typing scheme, 0.82 for MLST(7) and 0.95 for MLST(5). A novel sequence type combining 13% of all strains was discovered, as well as new alleles of genes. This is the first study characterizing the genetic diversity of the urogenital C. trachomatis population in Central Russia using MLST. We conclude that the MLST(7) scheme is the best possible choice for global epidemiological purposes, whereas MLST(5) is more appropriate for tracing local outbreaks.

Tuberculosis caused by multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mycobacterium tuberculosis (MTB) strains is a growing problem in many countries. The availability of the complete nucleotide sequences of several MTB genomes allows to use the comparative genomics as a tool to study the relationships of strains and differences in their evolutionary history including acquisition of drug-resistance. In our work, we sequenced three genomes of Russian MTB strains of different phenotypes – drug susceptible, MDR and XDR. Of them, MDR and XDR strains were collected in Tomsk (Siberia, Russia) during the local TB outbreak in 1998–1999 and belonged to rare KQ and KY families in accordance with IS6110 typing, which are considered endemic for Russia. Based on phylogenetic analysis, our isolates belonged to different genetic families, Beijing, Ural and LAM, which made the direct comparison of their genomes impossible. For this reason we performed their comparison in the broader context of all M. tuberculosis genomes available in GenBank. The list of unique individual non-synonymous SNPs for each sequenced isolate was formed by comparison with all SNPs detected within the same phylogenetic group. For further functional analysis, all proteins with unique SNPs were ascribed to 20 different functional classes based on Clusters of Orthologous Groups (COG). We have confirmed drug resistant status of our isolates that harbored almost all known drug-resistance associated mutations. Unique SNPs of an XDR isolate CTRI-4XDR, belonging to a Beijing family were compared in more detail with SNPs of additional 14 Russian XDR strains of the same family. Only type specific mutations in genes of repair, replication and recombination system (COG category L) were found common within this group. Probably the other unique SNPs discovered in CTRI-4XDR may have an important role in adaptation of this microorganism to its surrounding and in escape from antituberculosis drugs treatment.

Bacteriophage therapy is considered one of the most promising therapeutic approaches against multi-drug resistant bacterial infections. Infections caused by Staphylococcus aureus are very efficiently controlled with therapeutic bacteriophage cocktails, containing a number of individual phages infecting a majority of known pathogenic S. aureus strains. We assessed the contribution of individual bacteriophages comprising a therapeutic bacteriophage cocktail against S. aureus in order to optimize its composition. Two lytic bacteriophages vB_SauM-515A1 (Myoviridae) and vB_SauP-436A (Podoviridae) were isolated from the commercial therapeutic cocktail produced by Microgen (Russia). Host ranges of the phages were established on the panel of 75 S. aureus strains. Phage vB_SauM-515A1 lysed 85.3% and vB_SauP-436A lysed 68.0% of the strains, however, vB_SauP-436A was active against four strains resistant to vB_SauM-515A1, as well as to the therapeutic cocktail per se. Suboptimal results of the therapeutic cocktail application were due to extremely low vB_SauP-436A1 content in this composition. Optimization of the phage titers led to an increase in overall cocktail efficiency. Thus, one of the effective ways to optimize the phage cocktails design was demonstrated and realized by using bacteriophages of different families and lytic spectra.