Insight into pathogenomics and phylogeography of hypervirulent and highly-lethal Mycobacterium tuberculosis strain cluster.

BACKGROUND: . The Mycobacterium tuberculosis Beijing genotype is globally spread lineage with important medical properties that however vary among its subtypes. M. tuberculosis Beijing 14717-15-cluster was recently discovered as both multidrug-resistant, hypervirulent, and highly-lethal strain circulating in the Far Eastern region of Russia. Here, we aimed to analyze its pathogenomic features and phylogeographic pattern. RESULTS: . The study collection included M. tuberculosis DNA collected between 1996 and 2020 in different world regions. The bacterial DNA was subjected to genotyping and whole genome sequencing followed by bioinformatics and phylogenetic analysis. The PCR-based assay to detect specific SNPs of the Beijing 14717-15-cluster was developed and used for its screening in the global collections. Phylogenomic and phylogeographic analysis confirmed endemic prevalence of the Beijing 14717-15-cluster in the Asian part of Russia, and distant common ancestor with isolates from Korea (> 115 SNPs). The Beijing 14717-15-cluster isolates had two common resistance mutations RpsL Lys88Arg and KatG Ser315Thr and belonged to spoligotype SIT269. The Russian isolates of this cluster were from the Asian Russia while 4 isolates were from the Netherlands and Spain. The cluster-specific SNPs that significantly affect the protein function were identified in silico in genes within different categories (lipid metabolism, regulatory proteins, intermediary metabolism and respiration, PE/PPE, cell wall and cell processes). CONCLUSIONS: . We developed a simple method based on real-time PCR to detect clinically significant MDR and hypervirulent Beijing 14717-15-cluster. Most of the identified cluster-specific mutations were previously unreported and could potentially be associated with increased pathogenic properties of this hypervirulent M. tuberculosis strain. Further experimental study to assess the pathobiological role of these mutations is warranted.

PCR Assay for Rapid Taxonomic Differentiation of Virulent Staphylococcus aureus and Klebsiella pneumoniae Bacteriophages.

Phage therapy is now seen as a promising way to overcome the current global crisis in the spread of multidrug-resistant bacteria. However, phages are highly strain-specific, and in most cases one will have to isolate a new phage or search for a phage suitable for a therapeutic application in existing libraries. At an early stage of the isolation process, rapid screening techniques are needed to identify and type potential virulent phages. Here, we propose a simple PCR approach to differentiate between two families of virulent Staphylococcus phages (Herelleviridae and Rountreeviridae) and eleven genera of virulent Klebsiella phages (Przondovirus, Taipeivirus, Drulisvirus, Webervirus, Jiaodavirus, Sugarlandvirus, Slopekvirus, Jedunavirus, Marfavirus, Mydovirus and Yonseivirus). This assay includes a thorough search of a dataset comprising S. aureus (n = 269) and K. pneumoniae (n = 480) phage genomes available in the NCBI RefSeq/GenBank database for specific genes that are highly conserved at the taxonomic group level. The selected primers showed high sensitivity and specificity for both isolated DNA and crude phage lysates, which permits circumventing DNA purification protocols. Our approach can be extended and applied to any group of phages, given the large number of available genomes in the databases.

Isolation and Characterization of the First Zobellviridae Family Bacteriophage Infecting Klebsiella pneumoniae.

In order to address the upcoming crisis in the treatment of Klebsiella pneumoniae infections, caused by an increasing proportion of resistant isolates, new approaches to antimicrobial therapy must be developed. One approach would be to use (bacterio)phages and/or phage derivatives for therapy. In this study, we present a description of the first K. pneumoniae phage from the Zobellviridae family. The vB_KpnP_Klyazma podovirus, which forms translucent halos around the plaques, was isolated from river water. The phage genome is composed of 82 open reading frames, which are divided into two clusters located on opposite strands. Phylogenetic analysis revealed that the phage belongs to the Zobellviridae family, although its identity with the closest member of this family was not higher than 5%. The bacteriophage demonstrated lytic activity against all (n = 11) K. pneumoniae strains with the KL20 capsule type, but only the host strain was lysed effectively. The receptor-binding protein of the phage was identified as a polysaccharide depolymerase with a pectate lyase domain. The recombinant depolymerase protein showed concentration-dependent activity against all strains with the KL20 capsule type. The ability of a recombinant depolymerase to cleave bacterial capsular polysaccharides regardless of a phage's ability to successfully infect a particular strain holds promise for the possibility of using depolymerases in antimicrobial therapy, even though they only make bacteria sensitive to environmental factors, rather than killing them directly.

Phenoxazine-based scaffold for designing G4-interacting agents.

G-quadruplexes (G4) represent one class of non-canonical secondary nucleic acid structures that are currently regarded as promising and attractive targets for anti-cancer, anti-viral and antibacterial therapy. Herein, we probe a new i-clamp-inspired phenoxazine scaffold for designing G4-stabilizing ligands. The length of the protonated aminoalkyl tethers ('arms') of the phenoxazine-based ligand was optimized in silico. Two double-armed ligands differing in the relative orientation of their arms and one single-armed ligand were synthesized. The two-armed ligands significantly enhanced the thermal stability of the G-quadruplex structures (increasing the melting temperature by up to 20 °C) and displayed G4 selectivity over duplex DNA. The ligands look promising for biological studies and the phenoxazine scaffold could be a starting point for designing new G4-interacting compounds.

Simple Assay for Detection of the Central Asia Outbreak Clade of the Mycobacterium tuberculosis Beijing Genotype.

The Central Asia outbreak (CAO) clade is a branch of the Mycobacterium tuberculosis Beijing genotype that is associated with multidrug resistance, increased transmissibility, and epidemic spread in parts of the former Soviet Union. Furthermore, migration flows bring these strains far beyond their areas of origin. We aimed to find a specific molecular marker of the Beijing CAO clade and develop a simple and affordable method for its detection. Based on the bioinformatics analysis of the large M. tuberculosis whole-genome sequencing (WGS) data set (n = 1,398), we identified an IS6110 insertion in the Rv1359-Rv1360 intergenic region as a specific molecular marker of the CAO clade. We further designed and optimized a multiplex PCR method to detect this insertion. The method was validated in silico with the recently published WGS data set from Central Asia (n = 277) and experimentally with M. tuberculosis isolates from European and Asian parts of Russia, the former Soviet Union, and East Asia (n = 319). The developed molecular assay may be recommended for rapid screening of retrospective collections and for prospective surveillance when comprehensive but expensive WGS is not available or practical. The assay may be especially useful in high multidrug-resistant tuberculosis (MDR-TB) burden countries of the former Soviet Union and in countries with respective immigrant communities.

The role of IS6110 in micro- and macroevolution of Mycobacterium tuberculosis lineage 2.

The insertion sequence 6110 (IS6110) is the most studied transposable element in the Mycobacterium tuberculosis complex species. The element plays a significant role in genome plasticity of this important human pathogen, but still many causes and consequences of its transposition have not been fully studied. Here, we analyzed insertion sites for 902 Mycobacterium tuberculosis lineage 2 strains using whole-genome sequencing data. In total, 17,972 insertions were found, corresponding to 827 independent positions in the genome of the reference strain H37Rv. To trace the history of IS6110 expansion since proto-Beijing strains up to modern sublineages, we looked at the distribution of IS6110 across the genome-wide SNP-based phylogenetic tree. This analysis demonstrated a stepwise transposition of IS6110 that occurs by «copy-and-paste¼ mechanism. Additionally, we detected evolutionary-scale and sublineage-specific integration sites, which can be used for typing and for understanding the reasons for the success of the lineage. A significant part of such insertions affected the genes that are essential for the pathogen. Finally, we identified and confirmed deletions that occurred between differently oriented elements, which is uncommon for this family of insertion elements and appears to be another mechanism of genome variability.

Emerging peak on the phylogeographic landscape of Mycobacterium tuberculosis in West Asia: Definitely smoke, likely fire.

To date, a major attention was justly given to the global lineages of Mycobacterium tuberculosis. Here, we demonstrated an importance of the minor ones, on an example of intriguing and underestimated NEW-1 family that belongs to Euro-American lineage (lineage 4). Analysis of the global WGS/NGS datasets (5715 strains) identified 2235 strains of Lineage 4 and 66 strains of sublineage L4.5. This latter is marked with RD122 genomic deletion and includes NEW-1 family. Phylogenomic analysis confirmed a separate position of the NEW-1 family that we tentatively designate L4.5.1/Iran. We propose an evolution/migration scenario starting with origin of L4.5 1000-1300 ya in China, subsequent origin of the pre-NEW-1 intermediate genotype in Tibet, further migration to Xinjiang/Uyghur, and finally to Iran since 800 ya (origin of NEW-1), possibly, via expansion of the Mongol Yuan empire. Analysis of longitudinal phylogeographic datasets revealed a sharp increase in prevalence of NEW-1 strains in Iran and its eastwards neighbors in the last 20years; most alarmingly, it is accompanied with significant association with multidrug resistance (MDR). Ongoing migration, especially, Afghan refugees flows to developed countries emphasize a risk of the wider spread of the epidemic MDR subtype within NEW-1 family that we coin as emerging resistant clone of M. tuberculosis in West Asia.

Large scale analysis of amino acid substitutions in bacterial proteomics.

BACKGROUND: Proteomics of bacterial pathogens is a developing field exploring microbial physiology, gene expression and the complex interactions between bacteria and their hosts. One of the complications in proteomic approach is micro- and macro-heterogeneity of bacterial species, which makes it impossible to build a comprehensive database of bacterial genomes for identification, while most of the existing algorithms rely largely on genomic data. RESULTS: Here we present a large scale study of identification of single amino acid polymorphisms between bacterial strains. An ad hoc method was developed based on MS/MS spectra comparison without the support of a genomic database. Whole-genome sequencing was used to validate the accuracy of polymorphism detection. Several approaches presented earlier to the proteomics community as useful for polymorphism detection were tested on isolates of Helicobacter pylori, Neisseria gonorrhoeae and Escherichia coli. CONCLUSION: The developed method represents a perspective approach in the field of bacterial proteomics allowing to identify hundreds of peptides with novel SAPs from a single proteome.

Genome-wide Mycobacterium tuberculosis variation (GMTV) database: a new tool for integrating sequence variations and epidemiology.

BACKGROUND: Tuberculosis (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. DESCRIPTION: Here 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. CONCLUSIONS: Implementation 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.

A revised SNP-based barcoding scheme for typing Mycobacterium tuberculosis complex isolates.

The development of whole-genome sequencing technologies is gradually leading to a more detailed description of the population structure of the Mycobacterium tuberculosis complex (MTBC). In this study, we correlated previously published classifications on a collection of more than 10,000 genomes and proposed a new, comprehensive nomenclature that unifies the existing ones. In total, we identified 169 lineages and sublineages of M. tuberculosis/M. africanum and 9 animal-adapted species. For the purpose of organizing these genotypes in a more streamlined manner, we stratified them into five hierarchical levels. To represent the classification and compare it with the reference, we compiled a confirmatory data set of 670 high-quality isolates, which includes all genotypes and species of MTBC, and this confirmatory data set can serve as a basis for further studies. We proposed a set of 213 robust barcoding single-nucleotide polymorphisms and a suitable workflow for reliable differentiation of genotypes and species within the complex. This work integrates the results of all the major systematized studies to date to provide an understanding of the global diversity of the MTBC population structure. The results of this work may ultimately help to reliably determine the pathogen genotype and associate it with traits that reflect its prevalence, virulence, vaccination, and treatment efficiency, as well as to reliably find natural features revealed during its spread. IMPORTANCE Through years of research into the Mycobacterium tuberculosis complex (MTBC), a number of ambiguous phylogenetic classifications have emerged, which often overlap with one another. In the present study, we have combined all major studies on MTBC classification and inferred a unified, most complete to date classification and accompanying SNP barcodes.

Transcriptional Landscapes of Herelleviridae Bacteriophages and Staphylococcus aureus during Phage Infection: An Overview.

The issue of antibiotic resistance in healthcare worldwide has led to a pressing need to explore and develop alternative approaches to combat infectious diseases. Among these methods, phage therapy has emerged as a potential solution to tackle this growing challenge. Virulent phages of the Herelleviridae family, known for their ability to cause lysis of Staphylococcus aureus, a clinically significant pathogen frequently associated with multidrug resistance, have proven to be one of the most effective viruses utilized in phage therapy. In order to utilize phages for therapeutic purposes effectively, a thorough investigation into their physiology and mechanisms of action on infected cells is essential. The use of omics technologies, particularly total RNA sequencing, is a promising approach for analyzing the interaction between phages and their hosts, allowing for the assessment of both the behavior of the phage during infection and the cell's response. This review aims to provide a comprehensive overview of the physiology of the Herelleviridae family, utilizing existing analyses of their total phage transcriptomes. Additionally, it sheds light on the changes that occur in the metabolism of S. aureus when infected with virulent bacteriophages, contributing to a deeper understanding of the phage-host interaction.

Global Transcriptomic Response of Staphylococcus aureus to Virulent Bacteriophage Infection.

In light of the ever-increasing number of multidrug-resistant bacteria worldwide, bacteriophages are becoming a valid alternative to antibiotics; therefore, their interactions with host bacteria must be thoroughly investigated. Here, we report genome-wide transcriptional changes in a clinical Staphylococcus aureus SA515 strain for three time points after infection with the vB_SauM-515A1 kayvirus. Using an RNA sequencing approach, we identify 263 genes that were differentially expressed (DEGs) between phage-infected and uninfected host samples. Most of the DEGs were identified at an early stage of phage infection and were mainly involved in nucleotide and amino acid metabolism, as well as in cell death prevention. At the subsequent infection stages, the vast majority of DEGs were upregulated. Interestingly, 39 upregulated DEGs were common between the 15th and 30th minutes post-infection, and a substantial number of them belonged to the prophages. Furthermore, some virulence factors were overexpressed at the late infection stage, which necessitates more stringent host strain selection requirements for further use of bacteriophages for therapeutic purposes. Thus, this work allows us to better understand the influence of kayviruses on the metabolic systems of S. aureus and contributes to a better comprehension of phage therapy.

Genome-Wide Transcriptional Response of Mycobacterium smegmatis MC2155 to G-Quadruplex Ligands BRACO-19 and TMPyP4.

G-quadruplexes (G4s) are non-canonical DNA structures that could be considered as potential therapeutic targets for antimicrobial compounds, also known as G4-stabilizing ligands. While some of these ligands are shown in vitro to have a stabilizing effect, the precise mechanism of antibacterial action has not been fully investigated. Here, we employed genome-wide RNA-sequencing to analyze the response of Mycobacterium smegmatis to inhibitory concentrations of BRACO-19 and TMPyP4 G4 ligands. The expression profile changed (FDR < 0.05, log2FC > |1|) for 822 (515↑; 307↓) genes in M. smegmatis in response to BRACO-19 and for 680 (339↑; 341↓) genes in response to TMPyP4. However, the analysis revealed no significant ligand-induced changes in the expression levels of G4-harboring genes, genes under G4-harboring promoters, or intergenic regions located on mRNA-like or template strands. Meanwhile, for the BRACO-19 ligand, we found significant changes in the replication and repair system genes, as well as in iron metabolism genes which is, undoubtedly, evidence of the induced stress. For the TMPyP4 compound, substantial changes were found in transcription factors and the arginine biosynthesis system, which may indicate multiple biological targets for this compound.

The Gene Expression Profile Differs in Growth Phases of the Bifidobacterium Longum Culture.

To date, transcriptomics have been widely and successfully employed to study gene expression in different cell growth phases of bacteria. Since bifidobacteria represent a major component of the gut microbiota of a healthy human that is associated with numerous health benefits for the host, it is important to study them using transcriptomics. In this study, we applied the RNA-Seq technique to study global gene expression of B. longum at different growth phases in order to better understand the response of bifidobacterial cells to the specific conditions of the human gut. We have shown that in the lag phase, ABC transporters, whose function may be linked to active substrate utilization, are increasingly expressed due to preparation for cell division. In the exponential phase, the functions of activated genes include synthesis of amino acids (alanine and arginine), energy metabolism (glycolysis/gluconeogenesis and nitrogen metabolism), and translation, all of which promote active cell division, leading to exponential growth of the culture. In the stationary phase, we observed a decrease in the expression of genes involved in the control of the rate of cell division and an increase in the expression of genes involved in defense-related metabolic pathways. We surmise that the latter ensures cell survival in the nutrient-deprived conditions of the stationary growth phase.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis.

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a "perfect storm" that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

Proteogenomic Approach for Mycobacterium tuberculosis Investigation.

Recent advances in MS/MS technology have made it possible to use proteomic data to predict protein-coding sequences. This approach is called proteogenomics, and it allows to correctly translate start and stop sites and to reveal new open reading frames. Here, we focus on using proteogenomics to improve the annotation of Mycobacteriumtuberculosis strains. We also describe detail procedures of the extraction of proteins and their further preparation for LC-MS/MS analysis and outline the main steps of data analysis.

Substitutions in SurA and BamA Lead to Reduced Susceptibility to Broad Range Antibiotics in Gonococci.

There is growing concern about the emergence and spread of multidrug-resistant Neisseria gonorrhoeae. To effectively control antibiotic-resistant bacterial pathogens, it is necessary to develop new antimicrobials and to understand the resistance mechanisms to existing antibiotics. In this study, we discovered the unexpected onset of drug resistance in N. gonorrhoeae caused by amino acid substitutions in the periplasmic chaperone SurA and the ß-barrel assembly machinery component BamA. Here, we investigated the i19.05 clinical isolate with mutations in corresponding genes along with reduced susceptibility to penicillin, tetracycline, and azithromycin. The mutant strain NG05 (surAmut bamAmut, and penAmut) was obtained using the pan-susceptible n01.08 clinical isolate as a recipient in the transformation procedure. Comparative proteomic analysis of NG05 and n01.08 strains revealed significantly increased levels of other chaperones, Skp and FkpA, and some transport proteins. Efflux pump inhibition experiments demonstrated that the reduction in sensitivity was achieved due to the activity of efflux pumps. We hypothesize that the described mutations in the surA and bamA genes cause the qualitative and quantitative changes of periplasmic chaperones, which in turn alters the function of synthesized cell envelope proteins.

Transcriptomic Profile of Mycobacterium smegmatis in Response to an Imidazo[1,2-b][1,2,4,5]tetrazine Reveals Its Possible Impact on Iron Metabolism.

Tuberculosis (TB), caused by the Mycobacterium tuberculosis complex bacteria, is one of the most pressing health problems. The development of new drugs and new therapeutic regimens effective against the pathogen is one of the greatest challenges in the way of tuberculosis control. Imidazo[1,2-b][1,2,4,5]tetrazines have shown promising activity against M. tuberculosis and M. smegmatis strains. Mutations in MSMEG_1380 lead to mmpS5-mmpL5 operon overexpression, which provides M. smegmatis with efflux-mediated resistance to imidazo[1,2-b][1,2,4,5]tetrazines, but the exact mechanism of action of these compounds remains unknown. To assess the mode of action of imidazo[1,2-b][1,2,4,5]tetrazines, we analyzed the transcriptomic response of M. smegmatis to three different concentrations of 3a compound: 1/8×, 1/4×, and 1/2× MIC. Six groups of genes responsible for siderophore synthesis and transport were upregulated in a dose-dependent manner, while virtual docking revealed proteins involved in siderophore synthesis as possible targets for 3a. Thus, we suggest that imidazo[1,2-b][1,2,4,5]tetrazines may affect mycobacterial iron metabolism.

Probable long-term prevalence for a predominant Mycobacterium tuberculosis clone of a Beijing genotype in Colon, Panama.

Beijing genotype Mycobacterium tuberculosis strains associate with increased virulence, resistance and/or higher transmission rates. This study describes a specific Beijing strain predominantly identified in the Panamanian province of Colon with one of the highest incidences of tuberculosis in the country. Retrospective mycobacterial interspersed repetitive unit/variable number of tandem repeats analysis of 42 isolates collected between January and August 2018 allowed to identify a cluster (Beijing A) with 17 (40.5%) Beijing isolates. Subsequent prospective strain-specific PCR-based surveillance from September 2019 to March 2020 confirmed the predominance of the Beijing A strain (44.1%) in this province. Whole-genome sequencing revealed higher-than-expected diversity within the cluster, suggesting long-term prevalence of this strain and low number of cases caused by recent transmission. The Beijing A strain belongs to the Asian African 3 (Bmyc13, L2.2.5) branch of the modern Beijing sublineage, with their closest isolates corresponding to cases from Vietnam, probably introduced in Panama between 2000 and 2012.

Novel Klebsiella pneumoniae K23-Specific Bacteriophages From Different Families: Similarity of Depolymerases and Their Therapeutic Potential.

Antibiotic resistance is a major public health concern in many countries worldwide. The rapid spread of multidrug-resistant (MDR) bacteria is the main driving force for the development of novel non-antibiotic antimicrobials as a therapeutic alternative. Here, we isolated and characterized three virulent bacteriophages that specifically infect and lyse MDR Klebsiella pneumoniae with K23 capsule type. The phages belonged to the Autographiviridae (vB_KpnP_Dlv622) and Myoviridae (vB_KpnM_Seu621, KpS8) families and contained highly similar receptor-binding proteins (RBPs) with polysaccharide depolymerase enzymatic activity. Based on phylogenetic analysis, a similar pattern was also noted for five other groups of depolymerases, specific against capsule types K1, K30/K69, K57, K63, and KN2. The resulting recombinant depolymerases Dep622 (phage vB_KpnP_Dlv622) and DepS8 (phage KpS8) demonstrated narrow specificity against K. pneumoniae with capsule type K23 and were able to protect Galleria mellonella larvae in a model infection with a K. pneumoniae multidrug-resistant strain. These findings expand our knowledge of the diversity of phage depolymerases and provide further evidence that bacteriophages and phage polysaccharide depolymerases represent a promising tool for antimicrobial therapy.