{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,1]],"date-time":"2026-03-01T06:12:16Z","timestamp":1772345536019,"version":"3.50.1"},"reference-count":52,"publisher":"MDPI AG","issue":"5","license":[{"start":{"date-parts":[[2026,2,27]],"date-time":"2026-02-27T00:00:00Z","timestamp":1772150400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["IJMS"],"abstract":"<jats:p>Azithromycin remains an important agent in gonorrhea treatment, yet resistance is a growing global threat. To comprehensively define its genetic basis, we performed a large-scale genome-wide association study of 14,727 Neisseria gonorrhoeae genomes with linked azithromycin MICs from 66 countries. We identified 113 genetic variants significantly associated with elevated MICs. Beyond well-known mutations in 23S rRNA (A2059G, C2611T) and mtrCDE operon, we uncovered a broad repertoire of potential resistance determinants, including multiple amino acid substitutions in 16 ribosomal proteins (e.g., L2, L4, L13, L23) forming the nascent peptide exit tunnel (NPET), and porin PorB alterations (G120K, A121D\/N). Systematic pairwise analysis revealed extensive synergistic interactions, particularly between variants affecting drug influx\/efflux (PorB, MtrCDE) and ribosomal target affinity. Phylogenetic analysis identified successful, globally circulating lineages employing distinct resistance strategies: NPET-dominated, 23S rRNA-associated, and porin\/efflux-mediated. Our findings demonstrate that azithromycin resistance is a polygenic trait shaped by functional complementarity and epistasis between target modification, membrane permeability, and efflux. This integrated model is essential for accurate resistance prediction from genomic data and highlights key lineages for focused surveillance.<\/jats:p>","DOI":"10.3390\/ijms27052258","type":"journal-article","created":{"date-parts":[[2026,2,27]],"date-time":"2026-02-27T15:52:38Z","timestamp":1772207558000},"page":"2258","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["The Synergistic Armory: A Global Genome-Wide Association Study Reveals the Integrated Mechanisms of Azithromycin Resistance in Neisseria gonorrhoeae"],"prefix":"10.3390","volume":"27","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-0316-2262","authenticated-orcid":false,"given":"Boris","family":"Shaskolskiy","sequence":"first","affiliation":[{"name":"Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0009-0003-1647-1018","authenticated-orcid":false,"given":"Konstantin","family":"Tutaev","sequence":"additional","affiliation":[{"name":"Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0009-0009-8225-6870","authenticated-orcid":false,"given":"Dmitry","family":"Kravtsov","sequence":"additional","affiliation":[{"name":"Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-9416-875X","authenticated-orcid":false,"given":"Ilya","family":"Kandinov","sequence":"additional","affiliation":[{"name":"Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-3183-318X","authenticated-orcid":false,"given":"Dmitry","family":"Gryadunov","sequence":"additional","affiliation":[{"name":"Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2026,2,27]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"1003","DOI":"10.1134\/S0026893324700535","article-title":"Unveiling Neisseria gonorrhoeae Survival: Genetic Variability, Pathogenesis, and Antimicrobial Drug Resistance","volume":"58","author":"Shaskolskiy","year":"2024","journal-title":"Mol. Biol."},{"key":"ref_2","unstructured":"WHO (2024). Implementing the Global Health Sector Strategies on HIV, Viral Hepatitis and Sexually Transmitted Infections, 2022\u20132030: Report on Progress and Gaps 2024."},{"key":"ref_3","unstructured":"WHO (2024). Updated Recommendations for the Treatment of Neisseria gonorrhoeae, Chlamydia trachomatis, and Treponema pallidum (Syphilis) and New Recommendations on Syphilis Testing and Partner Services."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"73","DOI":"10.1093\/jac\/25.suppl_A.73","article-title":"The pharmacokinetics of azithromycin in human serum and tissues","volume":"25","author":"Foulds","year":"1990","journal-title":"J. Antimicrob. Chemother."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"159","DOI":"10.2165\/00044011-199815020-00009","article-title":"Pharmacokinetics of Oral Azithromycin in Serum, Urine, Polymorphonuclear Leucocytes and Inflammatory vs Non-Inflammatory Skin Blisters in Healthy Volunteers","volume":"15","author":"Ballow","year":"1998","journal-title":"Clin. Drug Investig."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"17158","DOI":"10.1073\/pnas.1008685107","article-title":"Revisiting the structures of several antibiotics bound to the bacterial ribosome","volume":"107","author":"Bulkley","year":"2010","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"814","DOI":"10.1038\/35101544","article-title":"Structural basis for the interaction of antibiotics with the peptidyl transferase centre in eubacteria","volume":"413","author":"Zarivach","year":"2001","journal-title":"Nature"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"508","DOI":"10.1016\/j.cell.2012.09.018","article-title":"Selective Protein Synthesis by Ribosomes with a Drug-Obstructed Exit Tunnel","volume":"151","author":"Kannan","year":"2012","journal-title":"Cell"},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"668","DOI":"10.1016\/j.tibs.2018.06.011","article-title":"How Macrolide Antibiotics Work","volume":"43","author":"Mankin","year":"2018","journal-title":"Trends Biochem. Sci."},{"key":"ref_10","doi-asserted-by":"crossref","unstructured":"Sergiev, P.V., Kabilov, M.R., Komarova, E.S., Serebryakova, M.V., Lazarev, V.N., Manakova, T.E., Khabibullina, N.F., Golovin, A.V., and Osterman, I.A. (2024). Context specificity of translation inhibitors revealed by Toe-seq. Res. Sq.","DOI":"10.21203\/rs.3.rs-3832009\/v1"},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"12026","DOI":"10.1038\/ncomms12026","article-title":"A combined cryo-EM and molecular dynamics approach reveals the mechanism of ErmBL-mediated translation arrest","volume":"7","author":"Arenz","year":"2016","journal-title":"Nat. Commun."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"787","DOI":"10.1038\/nsmb.2871","article-title":"A proton wire to couple aminoacyl-tRNA accommodation and peptide-bond formation on the ribosome","volume":"21","author":"Polikanov","year":"2014","journal-title":"Nat. Struct. Mol. Biol."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"3020","DOI":"10.1128\/AAC.46.9.3020-3025.2002","article-title":"Mutation in 23S rRNA associated with macrolide resistance in Neisseria gonorrhoeae","volume":"46","author":"Ng","year":"2002","journal-title":"Antimicrob. Agents Chemother."},{"key":"ref_14","unstructured":"CLSI (2025). Performance Standards for Antimicrobial Susceptibility Testing, Clinical and Laboratory Standards Institute. [35th ed.]. CLSI Supplement M100."},{"key":"ref_15","unstructured":"(2025, May 01). The European Committee on Antimicrobial Susceptibility Testing (EUCAST). Clinical Breakpoint Tables. Available online: https:\/\/www.eucast.org\/clinical_breakpoints\/."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"624","DOI":"10.1128\/AAC.02093-13","article-title":"First three Neisseria gonorrhoeae isolates with high-level resistance to azithromycin in Sweden: A threat to currently available dual-antimicrobial regimens for treatment of gonorrhea?","volume":"58","author":"Unemo","year":"2014","journal-title":"Antimicrob. Agents Chemother."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"3812","DOI":"10.1128\/AAC.00309-10","article-title":"High-level azithromycin resistance occurs in Neisseria gonorrhoeae as a result of a single point mutation in the 23S rRNA genes","volume":"54","author":"Chisholm","year":"2010","journal-title":"Antimicrob. Agents Chemother."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"001231","DOI":"10.1099\/mic.0.001231","article-title":"Transcriptional regulation of the mtrCDE efflux pump operon: Importance for Neisseria gonorrhoeae antimicrobial resistance","volume":"168","author":"Ayala","year":"2022","journal-title":"Microbiology"},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"e02281-18","DOI":"10.1128\/mBio.02281-18","article-title":"Mechanistic Basis for Decreased Antimicrobial Susceptibility in a Clinical Isolate of Neisseria gonorrhoeae Possessing a Mosaic-Like mtr Efflux Pump Locus","volume":"9","author":"Reimche","year":"2018","journal-title":"mBio"},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"621","DOI":"10.1099\/00221287-144-3-621","article-title":"Loss-of-function mutations in the mtr efflux system of Neisseria gonorrhoeae","volume":"144","author":"Veal","year":"1998","journal-title":"Microbiology"},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"5374","DOI":"10.1038\/s41467-020-19250-6","article-title":"Increased power from conditional bacterial genome-wide association identifies macrolide resistance mutations in Neisseria gonorrhoeae","volume":"11","author":"Ma","year":"2020","journal-title":"Nat. Commun."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"1804","DOI":"10.1086\/522964","article-title":"Regulation of the MtrC-MtrD-MtrE efflux-pump system modulates the in vivo fitness of Neisseria gonorrhoeae","volume":"196","author":"Warner","year":"2007","journal-title":"J. Infect. Dis."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"462","DOI":"10.1111\/j.1365-2958.2008.06424.x","article-title":"Clinically relevant mutations that cause derepression of the Neisseria gonorrhoeae MtrC-MtrD-MtrE Efflux pump system confer different levels of antimicrobial resistance and in vivo fitness","volume":"70","author":"Warner","year":"2008","journal-title":"Mol. Microbiol."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"e01419-18","DOI":"10.1128\/mBio.01419-18","article-title":"Azithromycin Resistance through Interspecific Acquisition of an Epistasis-Dependent Efflux Pump Component and Transcriptional Regulator in Neisseria gonorrhoeae","volume":"9","author":"Wadsworth","year":"2018","journal-title":"mBio"},{"key":"ref_25","doi-asserted-by":"crossref","unstructured":"Kandinov, I., Shaskolskiy, B., Kravtsov, D., Vinokurova, A., Gorshkova, S., Kubanov, A., Solomka, V., Shagabieva, J., Deryabin, D., and Dementieva, E. (2023). Azithromycin Susceptibility Testing and Molecular Investigation of Neisseria gonorrhoeae Isolates Collected in Russia, 2020\u20132021. Antibiotics, 12.","DOI":"10.3390\/antibiotics12010170"},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"e452","DOI":"10.1016\/S2666-5247(22)00044-1","article-title":"Europe-wide expansion and eradication of multidrug-resistant Neisseria gonorrhoeae lineages: A genomic surveillance study","volume":"3","author":"Cole","year":"2022","journal-title":"Lancet Microbe"},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"e01509-20","DOI":"10.1128\/mBio.01509-20","article-title":"Efflux Pump Antibiotic Binding Site Mutations Are Associated with Azithromycin Nonsusceptibility in Clinical Neisseria gonorrhoeae Isolates","volume":"11","author":"Ma","year":"2020","journal-title":"mBio"},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"1752","DOI":"10.1093\/jac\/dkab084","article-title":"Molecular pathways to high-level azithromycin resistance in Neisseria gonorrhoeae","volume":"76","author":"Laumen","year":"2021","journal-title":"J. Antimicrob. Chemother."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"1579","DOI":"10.1093\/infdis\/jiw420","article-title":"Genomic Epidemiology of Gonococcal Resistance to Extended-Spectrum Cephalosporins, Macrolides, and Fluoroquinolones in the United States, 2000\u20132013","volume":"214","author":"Grad","year":"2016","journal-title":"J. Infect. Dis."},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"3877","DOI":"10.1128\/AAC.47.12.3877-3880.2003","article-title":"Acquired macrolide resistance genes in pathogenic Neisseria spp. isolated between 1940 and 1987","volume":"47","author":"Cousin","year":"2003","journal-title":"Antimicrob. Agents Chemother."},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"412","DOI":"10.1038\/s41589-020-00715-0","article-title":"Structure of Erm-modified 70S ribosome reveals the mechanism of macrolide resistance","volume":"17","author":"Svetlov","year":"2021","journal-title":"Nat. Chem. Biol."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"2471","DOI":"10.1093\/jac\/dkw182","article-title":"Molecular epidemiology and mechanisms of resistance of azithromycin-resistant Neisseria gonorrhoeae isolated in France during 2013\u201314","volume":"71","author":"Belkacem","year":"2016","journal-title":"J. Antimicrob. Chemother."},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"272","DOI":"10.2174\/1389450115666141120110724","article-title":"History and epidemiology of antibiotic susceptibilities of Neisseria gonorrhoeae","volume":"16","author":"Shigemura","year":"2015","journal-title":"Curr. Drug Targets"},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"e478","DOI":"10.1016\/S2666-5247(23)00370-1","article-title":"Antimicrobial-resistant Neisseria gonorrhoeae in Europe in 2020 compared with in 2013 and 2018: A retrospective genomic surveillance study","volume":"5","author":"Golparian","year":"2024","journal-title":"Lancet Microbe"},{"key":"ref_35","doi-asserted-by":"crossref","unstructured":"Mauffrey, F., Poncet, F., Jacot, D., Greub, G., Nordmann, P., and Blanc, D.S. (2024). Impact of mutations in the mtrR, rpdlVD and rrl genes on azithromycin resistance in Neisseria gonorrhoeae. PLoS ONE, 19.","DOI":"10.1371\/journal.pone.0306695"},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"354","DOI":"10.1128\/AAC.37.2.354","article-title":"Outer membrane permeability barrier to azithromycin, clarithromycin, and roxithromycin in gram-negative enteric bacteria","volume":"37","author":"Vaara","year":"1993","journal-title":"Antimicrob. Agents Chemother."},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"18923","DOI":"10.3390\/molecules201018923","article-title":"Membrane Interactions of Phytochemicals as Their Molecular Mechanism Applicable to the Discovery of Drug Leads from Plants","volume":"20","author":"Tsuchiya","year":"2015","journal-title":"Molecules"},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"eadc9566","DOI":"10.1126\/sciadv.adc9566","article-title":"Lipids mediate supramolecular outer membrane protein assembly in bacteria","volume":"8","author":"Webby","year":"2022","journal-title":"Sci. Adv."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"236","DOI":"10.1038\/s41579-021-00638-0","article-title":"Physical properties of the bacterial outer membrane","volume":"20","author":"Sun","year":"2022","journal-title":"Nat. Rev. Microbiol."},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"569","DOI":"10.1021\/acsmedchemlett.7b00499","article-title":"Azithromycin-Induced Changes to Bacterial Membrane Properties Monitored in Vitro by Second-Harmonic Light Scattering","volume":"9","author":"Wilhelm","year":"2018","journal-title":"ACS Med. Chem. Lett."},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"203","DOI":"10.1007\/s00232-002-1076-7","article-title":"The macrolide antibiotic azithromycin interacts with lipids and affects membrane organization and fluidity: Studies on Langmuir-Blodgett monolayers, liposomes and J774 macrophages","volume":"192","author":"Tyteca","year":"2003","journal-title":"J. Membr. Biol."},{"key":"ref_42","doi-asserted-by":"crossref","first-page":"61","DOI":"10.1186\/s13073-021-00858-2","article-title":"A community-driven resource for genomic epidemiology and antimicrobial resistance prediction of Neisseria gonorrhoeae at Pathogenwatch","volume":"13","author":"Yeats","year":"2021","journal-title":"Genome Med."},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"124","DOI":"10.12688\/wellcomeopenres.14826.1","article-title":"Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications","volume":"3","author":"Jolley","year":"2018","journal-title":"Wellcome Open Res."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"12728","DOI":"10.1038\/s41598-021-91456-0","article-title":"AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence","volume":"11","author":"Feldgarden","year":"2021","journal-title":"Sci. Rep."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"772","DOI":"10.1093\/molbev\/mst010","article-title":"MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability","volume":"30","author":"Katoh","year":"2013","journal-title":"Mol. Biol. Evol."},{"key":"ref_46","doi-asserted-by":"crossref","unstructured":"Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., and Madden, T.L. (2009). BLAST+: Architecture and applications. BMC Bioinform., 10.","DOI":"10.1186\/1471-2105-10-421"},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"e15","DOI":"10.1093\/nar\/gku1196","article-title":"Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins","volume":"43","author":"Croucher","year":"2015","journal-title":"Nucleic Acids Res."},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"1530","DOI":"10.1093\/molbev\/msaa015","article-title":"IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era","volume":"37","author":"Minh","year":"2020","journal-title":"Mol. Biol. Evol."},{"key":"ref_49","doi-asserted-by":"crossref","first-page":"307","DOI":"10.1093\/sysbio\/syq010","article-title":"New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0","volume":"59","author":"Guindon","year":"2010","journal-title":"Syst. Biol."},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"685","DOI":"10.1093\/sysbio\/syr041","article-title":"Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes","volume":"60","author":"Anisimova","year":"2011","journal-title":"Syst. Biol."},{"key":"ref_51","unstructured":"R Core Team (2025). R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing. Available online: https:\/\/www.R-project.org\/."},{"key":"ref_52","doi-asserted-by":"crossref","first-page":"28","DOI":"10.1111\/2041-210X.12628","article-title":"ggtree: An r package for visualization and annotation of phylogenetic trees with their covariates and other associated data","volume":"8","author":"Yu","year":"2017","journal-title":"Methods Ecol. Evol."}],"container-title":["International Journal of Molecular Sciences"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/1422-0067\/27\/5\/2258\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,3,1]],"date-time":"2026-03-01T05:18:17Z","timestamp":1772342297000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/1422-0067\/27\/5\/2258"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,2,27]]},"references-count":52,"journal-issue":{"issue":"5","published-online":{"date-parts":[[2026,3]]}},"alternative-id":["ijms27052258"],"URL":"https:\/\/doi.org\/10.3390\/ijms27052258","relation":{},"ISSN":["1422-0067"],"issn-type":[{"value":"1422-0067","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,2,27]]}}}