Search

The aim of this work was to study the resistance to macrolides (azithromycin) in the modern Russian population of N. gonorrhoeae with the analysis of genetic resistance determinants. Azithromycin is not used to treat gonococcal infection in Russia. However, among 162 isolates collected in 2020–2021, 22 isolates (13.6%) were phenotypically resistant to azithromycin. Mutations in 23S rRNA genes were found only in two isolates; erm and mefA genes were absent. Azithromycin resistance was shown to be predominantly associated with mutations in the mtrR and mtrD genes of the MtrCDE efflux pump and their mosaic alleles which may have formed due to a horizontal transfer from N. meningitidis. A total of 30 types of mtrR alleles and 10 types of mtrD alleles were identified including mosaic variants. Matching between the mtrR and mtrD alleles was revealed to indicate the cooperative molecular evolution of these genes. A link between the mtrR and mtrD alleles and NG-MAST types was found only for NG-MAST 228 and 807, typical of N. gonorrhoeae in Russia. The high level of resistance to azithromycin in Russia may be related to the spread of multiple transferable resistance to antimicrobials regardless of their use in the treatment of gonococcal infection.

Chromosomal rearrangements in N. gonorrhoeae and N. meningitidis were studied with the determination of mobile elements and their role in rearrangements. The results of whole-genome sequencing and de novo genome assembly for 50 N. gonorrhoeae isolates collected in Russia were compared with 96 genomes of N. gonorrhoeae and 138 genomes of N. meningitidis from the databases. Rearrangement events with the determination of the coordinates of syntenic blocks were analyzed using the SibeliaZ software, the minimum number of events that allow one genome to pass into another was calculated using the DCJ-indel model using the UniMoG program. Population-level analysis revealed a stronger correlation between changes in the gene order and phylogenetic proximity for N. meningitidis in contrast to N. gonorrhoeae. Mobile elements were identified, including Correa elements, Spencer-Smith elements (in N. gonorrhoeae), Neisserial intergenic mosaic elements, IS elements of IS5, IS30, IS110, IS1595 groups, Nf1-Nf3 prophages, NgoФ1-NgoФ9 prophages, Mu-like prophages Pnm1, Pnm2, MuMenB (in N. meningitidis). More than 44% of the observed rearrangements most likely occurred with the participation of mobile elements, including prophages. No differences were found between the Russian and global N. gonorrhoeae population both in terms of rearrangement events and in the number of transposable elements in genomes.

Chromosomal rearrangements in

BackgroundDecreased susceptibility of Neisseria gonorrhoeae to extended-spectrum cephalosporins is a major concern. Elucidation of the phenotypic and genetic characteristics of such isolates is a priority task.MethodsWe developed a method for predicting the N. gonorrhoeae ceftriaxone susceptibility level (MICcro) by identifying genetic determinants of resistance using low-density hydrogel microarrays and a regression equation. A training dataset, containing 5631 isolates from the Pathogenwatch database and 181 isolates obtained in the Russian Federation during 2018–19, was used to build a regression model. The regression equation was tested on 14 WHO reference strains. Ceftriaxone resistance determinants for the 448 evaluated clinical isolates collected in Russia were identified using microarray analysis, and MICcro values were calculated using the regression equation and compared with those measured by the serial dilution method.ResultsThe regression equation for calculating MICcro values included 20 chromosomal resistance determinants. The greatest contributions to the increase in MICcro were shown to be PBP2: Ala-501→Pro, Ala-311→Val, Gly-545→Ser substitutions, Asp(345–346) insertion; and PorB: Gly-120→Arg substitution. The substitutions PBP2: Ala-501→Thr/Val, PorB: Gly-120→Asn/Asp/Lys and PBP1: Leu-421→Pro had weaker effects. For 94.4% of the isolates in the evaluation set, the predicted MICcro was within one doubling dilution of the experimentally determined MICcro. No ceftriaxone-resistant isolates were identified in the analysed samples from Russia, and no interpretative errors were detected in the MICcro calculations.ConclusionsThe developed strategy for predicting ceftriaxone MIC can be used for the continuous surveillance of known and emerging resistant N. gonorrhoeae isolates.

Gonococcal infection caused by the Gram-negative bacteria

N. gonorrhoeae strains producing plasmid-mediated broad- and extended-spectrum β-lactamases (ESBLs) have been obtained in vitro and studied for their viability and susceptibility to β-lactam antibiotics. The artificial pblaTEM-1 and pblaTEM-20 plasmids were constructed by site-directed mutagenesis from a pblaTEM-135 plasmid of the Toronto/Rio type, which was extracted from the clinical isolate. MIC values were determined for a series of β-lactam antibiotics, including benzylpenicillin, ampicillin, cefuroxime, ceftriaxone, cefixime, cefotaxime, cefepime, meropenem, imipenem, and doripenem. The N. gonorrhoeae strain carrying the pblaTEM-20 plasmid exhibited a high level of resistance to penicillins and II-IV generation cephalosporins (MIC ≥ 2 mg/L) but not to carbapenems (MIC ≤ 0.008 mg/L). However, this strain stopped growing after 6 hours of cultivation. The reduced viability was not associated with the plasmid loss but can be explained by the following factors: the presence of the plasmid itself, which requires additional costs for its reproduction, and the expression of ESBL, which can affect the structure of the peptidoglycan layer in the cell membrane. The cell growth was mathematically modeled using the generalized Verhulst equation. The parameter, characterizing the maximum possible number of cells grown under given conditions, decreased within the wild type (without plasmids) - pblaTEM-135 - pblaTEM-1 - pblaTEM-20 series, i.e., the plasmid-bearing strains had reduced viability compared to that of the wild-type strain. The kinetics for the cell death of N. gonorrhoeae strains without the pblaTEM-20 plasmid in the presence of ceftriaxone can be described by a modified Chick-Watson law. For the N. gonorrhoeae strain that harbors the pblaTEM-20 plasmid, the kinetics for cell number changes with time reflected several processes: the hydrolysis of ceftriaxone by TEM-20 β-lactamase, growth and gradual death of cells. The demonstrated reduced viability of N. gonorrhoeae strains with the pblaTEM-20 plasmid probably explains the absence of clinical isolates of N. gonorrhoeae that produce ESBL.

The goal of this work was to determine the factors affecting the emergence of azithromycin-resistant Neisseria gonorrhoeae isolates in Russia, where azithromycin was never recommended for the treatment of gonococcal infections. Clinical N. gonorrhoeae isolates collected in 2018–2021 (428 isolates) were analyzed. No azithromycin-resistant isolates were found in 2018–2019, but in 2020–2021, a significant increase in the ratio of azithromycin-resistant isolates was observed: 16.8% and 9.3%, respectively. A hydrogel DNA microarray was developed for the analysis of resistance determinants: mutations in the genes encoding the mtrCDE efflux system and in all four copies of the 23S rRNA gene (position 2611). A majority of the azithromycin-resistant Russian isolates belonged to the NG-MAST G12302 genogroup, and the resistance was associated with the presence of a mosaic structure of the mtrR gene promoter region with the −35 delA deletion, an Ala86Thr mutation in the mtrR gene, and a mosaic structure of the mtrD gene. A comparative phylogenetic study of modern Russian and European N. gonorrhoeae populations allowed us to conclude that the emergence of azithromycin resistance in Russia in 2020 was the result of the appearance and spread of European N. gonorrhoeae strains belonging to the G12302 genogroup due to possible cross-border transfer.

Comparative whole-genome analysis was performed for Neisseria gonorrhoeae isolates belonging to the Neisseria gonorrhoeae multiantigen sequence typing (NG-MAST) types predominant worldwide — 225, 1407, 2400, 2992, and 4186 — and to genogroup 807, the most common genogroup in the Russian Federation. Here, for the first time, the complete genomes of 25 N. gonorrhoeae isolates from genogroup 807 were obtained. For NG-MAST types 225, 1407, 2400, 2992, and 4186, genomes from the Pathogenwatch database were used. The phylogenetic network constructed for 150 genomes showed that the clustering according to NG-MAST type corresponded to the clustering according to genome. Comparisons of genomes of the six sequence types revealed 8-20 genes specific to each sequence type, including the loci for phase variations and genetic components of the gonococcal genetic island (GGI). NG-MAST type 2992 and 4186 isolates either lacked the GGI or carried critical mutations in genes essential for DNA secretion. In all analyzed genogroup 807 isolates, substitution of the essential atlA gene with the eppA gene was found, accompanied by a change in the traG allele, replacement of the ych gene with ych1, and the absence of the exp1 gene, which is likely to result in loss of GGI functionality. For the NG-MAST type 225, 1407 and 2400 isolates, no premature stop codons or reading frameshifts were found in the genes essential for GGI function. A relationship between isolate susceptibility to ciprofloxacin, penicillin, tetracycline and the presence of lesions in GGI genes necessary for DNA secretion was established. The N. gonorrhoeae evolutionary pathways, which allow a particular sequence type to maintain long-term predominance in the population, may include changes in genes responsible for adhesion and virulence, changes in the GGI structure, preservation of genes carrying drug resistance determinants, and changes in genes associated with host adaptation or encoding enzymes of biochemical pathways.

Comparative whole-genome analysis was performed for

Chromosomal rearrangements in N. gonorrhoeae and N. meningitidis were studied with the determination of mobile elements and their role in rearrangements. The results of whole-genome sequencing and de novo genome assembly for 50 N. gonorrhoeae isolates collected in Russia were compared with 96 genomes of N. gonorrhoeae and 138 genomes of N. meningitidis from the databases. Rearrangement events with the determination of the coordinates of syntenic blocks were analyzed using the SibeliaZ software v.1.2.5, the minimum number of events that allow one genome to pass into another was calculated using the DCJ–indel model using the UniMoG program v.1.0. Population-level analysis revealed a stronger correlation between changes in the gene order and phylogenetic proximity for N. meningitidis in contrast to N. gonorrhoeae. Mobile elements were identified, including Correa elements; Spencer-Smith elements (in N. gonorrhoeae); Neisserial intergenic mosaic elements; IS elements of IS5, IS30, IS110, IS1595 groups; Nf1–Nf3 prophages; NgoФ1–NgoФ9 prophages; and Mu-like prophages Pnm1, Pnm2, MuMenB (in N. meningitidis). More than 44% of the observed rearrangements most likely occurred with the participation of mobile elements, including prophages. No differences were found between the Russian and global N. gonorrhoeae population both in terms of rearrangement events and in the number of transposable elements in genomes.

A multiplex assay based on a low-density hydrogel microarray was developed to identify genomic substitutions in N. gonorrhoeae that determine resistance to the currently recommended treatment agents ceftriaxone and azithromycin and the previously used drugs penicillin, tetracycline, and ciprofloxacin. The microarray identifies 74 drug resistance determinants in the N. gonorrhoeae&nbsp, penA, ponA, porB, gyrA, parC, rpsJ, mtrR, blaTEM,&nbsp, tetM, and 23S rRNA genes. The hydrogel elements were formed by automated dispensing of nanoliter-volume droplets followed by UV-induced copolymerization of NH2-containing oligonucleotides with gel-forming monomers. Polybutylene terephthalate plates without special modifications were used as microarray substrates. Sequences and concentrations of immobilized oligonucleotides, gel composition, and hybridization conditions were carefully selected, and the median discrimination ratio ranged from 2.8 to 29.4, allowing unambiguous identification of single-nucleotide substitutions. The mutation identification results in a control sample of 180 N. gonorrhoeae isolates were completely consistent with the Sanger sequencing results. In total, 648 clinical N. gonorrhoeae isolates obtained in Russia during the last 5 years were analyzed and genotyped using these microarrays. The results allowed us to draw conclusions about the present situation with antimicrobial susceptibility of N. gonorrhoeae in Russia and demonstrated the possibility of using hydrogel microarrays to control the spread of antibiotic resistance.

Surveillance of multidrug-resistant infections is a priority task for contemporary epidemiology. The aim of this study was to genotype modern clinical isolates of N. gonorrhoeae using the NG-MAST technique (Neisseria gonorrhoeae multi-antigen sequence typing) and to compare the phylogeny of the gonococcal pathogens coming from Russia, European Union and Japan. We studied a total of 822 isolates collected in Russia from 2013 through 2018. We also used NG-MAST data from the following databases: PathogenWatch (European Union, 1,071 isolates) and PubMLST (Japan, 206 isolates). Russian isolates represented 301 different NG-MAST types. The most common were types 807, 228, 1993, 5714, and 9476 (8.3%, 3.3%, 3.2%, 3.2%, and 2.7%, respectively). There were only 3 isolates (0.4%) from Russia that represented the epidemiologically significant sequence type 1407 prevailing in many countries and characterized by multiple determinants of antimicrobial resistance. A phylogenetic tree for the NG-MAST types found in Russia and European countries was constructed. The cluster analysis of the proportion of isolates belonging to unique sequence types and the country population size allowed us to identify 2 clusters (significance level — 0.01): the first cluster included Russia and Japan, the second, European countries. A distribution pattern was identified for unique sequence types: the greater is the population size, the higher is their proportion. The phylogenetic analysis demonstrated a genetic distance between the most common Russian, European and Japanese sequence types, suggesting that the Russian population of N. gonorrhoeae has been evolving relatively locally.

Gonococcal infection caused by the Gram-negative bacteria Neisseria gonorrhoeae is one of the most common sexually transmitted infections (STIs) worldwide [...]

This work aimed to study penA gene polymorphisms in clinical isolates of Neisseria gonorrhoeae collected in Russia in 2018–2019 and the contribution of the penA allele type to susceptibility to β-lactam antibiotics. A total of 182 isolates were analyzed. penA allele types were determined by sequencing, and the minimum inhibitory concentrations (MICs) of benzylpenicillin and ceftriaxone were measured. The influence of genetic factors on MICs was evaluated by regression analysis. All isolates were susceptible to ceftriaxone, and 40.1% of isolates were susceptible to penicillin. Eleven penA allele types were identified. The mosaic type XXXIV penA allele and the Gly120Lys substitution in PorB made the greatest contributions to increasing the ceftriaxone MIC; the presence of the blaTEM plasmid, Gly120Asp, Ala121Gly/Asn substitutions in PorB, and the adenine deletion in the promoter region of the mtrR gene caused an increase in the penicillin MIC. Among 61 NG-MAST types identified, the most frequent were types 228, 807, 9486, 1993, and 6226. A link between penA alleles and Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) types was established. Resistance to two groups of β-lactam antibiotics was associated with non-identical changes in penA alleles. To prevent the emergence of ceftriaxone resistance in Russia, NG-MAST genotyping must be supplemented with penA allele analysis.

Objectives: The goal of this work was to assess the genetic diversity of Neisseria gonorrhoeae isolates in Russia and Europe and to compare the distribution of the N. gonorrhoeae multi-antigen sequencing types (NG-MAST) of Russian isolates with that of isolates from European countries. Methods: NG-MAST typing was performed for 804 N. gonorrhoeae isolates collected in Russia in 2013–2018. For isolates from European countries, data from the https://pathogen.watch/collection/eurogasp2013 database were used. Results: Among the isolates from Russia, 296 NG-MAST types were found. A maximum likelihood phylogenetic tree was constructed. Phylogenetic analysis revealed seven major genogroups uniting the most frequent Russian sequence types: G807, G1993, G9476, G14942, G1152, G9486, and G12531. Conclusions: The NG-MAST type distribution in Russia differed from that in European countries. Most of the Russian isolates had sequence types that were not found in Europe. Only 33% of the Russian isolates belonged to genogroups established for European countries, and the widespread European genogroup G1407 was represented by only nine isolates. Analysis of the Russian isolates belonging to phylogenetically close European genogroups indicated similarities in drug resistance, although no epidemically dangerous drug-resistant clones were found among the Russian isolates. Keywords: Neisseria gonorrhoeae, N. gonorrhoeae multi-antigen sequence typing (NG-MAST), Phylogenetic tree, Phylogenetic analysis

This work aimed to study penA gene polymorphisms in clinical isolates of Neisseria gonorrhoeae collected in Russia in 2018–2019 and the contribution of the penA allele type to susceptibility to β-lactam antibiotics. A total of 182 isolates were analyzed. penA allele types were determined by sequencing, and the minimum inhibitory concentrations (MICs) of benzylpenicillin and ceftriaxone were measured. The influence of genetic factors on MICs was evaluated by regression analysis. All isolates were susceptible to ceftriaxone, and 40.1% of isolates were susceptible to penicillin. Eleven penA allele types were identified. The mosaic type XXXIV penA allele and the Gly120Lys substitution in PorB made the greatest contributions to increasing the ceftriaxone MIC; the presence of the blaTEM plasmid, Gly120Asp, Ala121Gly/Asn substitutions in PorB, and the adenine deletion in the promoter region of the mtrR gene caused an increase in the penicillin MIC. Among 61 NG-MAST types identified, the most frequent were types 228, 807, 9486, 1993, and 6226. A link between penA alleles and Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) types was established. Resistance to two groups of β-lactam antibiotics was associated with non-identical changes in penA alleles. To prevent the emergence of ceftriaxone resistance in Russia, NG-MAST genotyping must be supplemented with penA allele analysis.