{"id":134835,"date":"2025-09-13T15:32:08","date_gmt":"2025-09-13T15:32:08","guid":{"rendered":"https:\/\/www.newsbeep.com\/uk\/134835\/"},"modified":"2025-09-13T15:32:08","modified_gmt":"2025-09-13T15:32:08","slug":"bioactive-molecules-unearthed-by-terabase-scale-long-read-sequencing-of-a-soil-metagenome","status":"publish","type":"post","link":"https:\/\/www.newsbeep.com\/uk\/134835\/","title":{"rendered":"Bioactive molecules unearthed by terabase-scale long-read sequencing of a soil metagenome"},"content":{"rendered":"<p>Bacterial strains and growth conditions<\/p>\n<p>All strains, unless otherwise specified, were grown in Luria\u2013Bertani (LB) at 30\u2009\u00b0C with 200-rpm shaking. Antibiotic activity test strains and ESKAPE pathogens were: B.\u2009subtilis 168, S.\u2009aureus SH1000, E.\u2009coli ATCC 25922, E.\u2009faecium COM15, K.\u2009pneumoniae ATCC 10031, A.\u2009baumannii ATCC 17978, P.\u2009aeruginosa PAO1 and E.\u2009cloacae ATCC 13047.<\/p>\n<p>Sample collection<\/p>\n<p>In late May 2023, we collected samples from Rockefeller University\u2019s 1,200-acre Center for Field Research in Ethology located in Dutchess County, New York. In order of collection, these comprised soil from a healthy forest (forest 1), sediment from a muddy bog, soil from the edge of a lake (forest 2), soil from a grove of dead trees (forest 3), decaying material from a dead tree and soil from a dry forest area (forest 4). Soil from a maximum depth of 10\u2009cm was collected from an area of 50\u2009cm2 into a gallon-sized Ziploc bag. A single collected bag of forest soil 1 was used for testing and scaling. After collection, bags were stored at 4\u2009\u00b0C.<\/p>\n<p>Nycodenz isolation of soil microorganisms<\/p>\n<p>Gradient centrifugation for separation of cells from the soil matrix and testing washes for removal of copurified contaminants were motivated by previous approaches in the field<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Lindahl, V. &amp; Bakken, L. R. Evaluation of methods for extraction of bacteria from soil. FEMS Microbiol. Ecol. 16, 135&#x2013;142 (1995).\" href=\"#ref-CR42\" id=\"ref-link-section-d217346668e1340\">42<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Cheng, F. et al. Soil pretreatment and fast cell lysis for direct polymerase chain reaction from forest soils for terminal restriction fragment length polymorphism analysis of fungal communities. Bras. J. Microbiol. 47, 817&#x2013;827 (2016).\" href=\"#ref-CR43\" id=\"ref-link-section-d217346668e1340_1\">43<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 44\" title=\"Amalfitano, S. &amp; Fazi, S. Recovery and quantification of bacterial cells associated with streambed sediments. J. Microbiol. Methods 75, 237&#x2013;243 (2008).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR44\" id=\"ref-link-section-d217346668e1343\" rel=\"nofollow noopener\" target=\"_blank\">44<\/a>. Soil (62.5\u2009g) was added to 200\u2009ml of prechilled 10\u2009mM sodium pyrophosphate (Sigma) and blended (Hamilton Beach 10 speed) for 1\u2009min using the LO\/batter setting. The container was then placed on ice for 5\u2009min. The 1-min blending and 5-min ice break were repeated two more times. Blended soil was then filtered through a cheesecloth (grade 90, bleached) to remove large insoluble particles and the filtrate was gently layered onto 8\u2009ml of 36% w\/v nycodenz (Accurate Chemical and Scientific Corporation) to the top of 50-ml conicals (Sarstedt). The tubes were centrifuged at 4\u2009\u00b0C, 5,000g for 1.5\u2009h in a swing bucket rotor yielding a fuzzy beige cell layer at the interface of the top (aqueous) and bottom (nycodenz) layer. The cell layer was gently aspirated and combined in a fresh 50-ml conical. The cell suspension was mixed 1:1 with sterile water and the mixture centrifuged at 4\u2009\u00b0C, 5,000g for 1\u2009h. The supernatant was discarded and the cell pellet was washed twice with 40\u2009ml of JSM buffer (0.9% w\/v NaCl and 1% w\/v instant nonfat dry milk (Nestle)), followed by two washes with 40\u2009ml of TE lysis buffer (25\u2009mM Tris, 25\u2009mM EDTA and 50\u2009mM glucose). Lastly, the cell pellet was resuspended in a small volume (1\u20135\u2009ml) of lysis buffer, the optical density at 600\u2009nm (OD600) was checked and the bacterial suspension diluted as necessary for the desired DNA extraction protocol.<\/p>\n<p>DNA extractionsDirect soil extraction<\/p>\n<p>DNA extraction directly from soil was performed as previously described<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 18\" title=\"Brady, S. F. Construction of soil environmental DNA cosmid libraries and screening for clones that produce biologically active small molecules. Nat. Protoc. 2, 1297&#x2013;1305 (2007).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR18\" id=\"ref-link-section-d217346668e1367\" rel=\"nofollow noopener\" target=\"_blank\">18<\/a>. Briefly, 2.4\u2009ml of CTAB lysis buffer (100\u2009mM Tris-HCl, 100\u2009mM Na EDTA, 1.5\u2009M NaCl, 1% w\/v CTAB and 2% w\/v SDS, pH 8.0) was added to 2\u2009g of soil in a 5-ml tube and the slurry was incubated at 70\u2009\u00b0C for 2\u2009h. For the nycodenz-isolated bacteria, 3\u2009ml of a 30-OD600 bacterial suspension was pelleted, then resuspended in 4\u2009ml of CTAB lysis buffer and incubated at 70\u2009\u00b0C for 2\u2009h. After incubation, the suspensions were centrifuged at 20,000g for 10\u2009min and the supernatant was moved to new tubes. Supernatants were combined with 0.7 volumes of isopropanol and incubated at room temperature for 10\u2009min. Samples were centrifuged at 20,000g for 30\u2009min and the resulting DNA pellets was washed twice with 1\u2009ml of 70% ethanol. The pellets were then air-dried for 1\u20132\u2009h at room temperature and resuspended in 100\u2013200\u2009\u00b5l of 5\u2009mM Tris pH 8 overnight at 4\u2009\u00b0C.<\/p>\n<p>Gel plug<\/p>\n<p>This protocol is an adaptation of gel plug lysis described by Walter et al.<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 45\" title=\"Walter, J., Mangold, M. &amp; Tannock, G. W. Construction, analysis, and &#x3B2;-glucanase screening of a bacterial artificial chromosome library from the large-bowel microbiota of mice. Appl. Environ. Microbiol. 71, 2347&#x2013;2354 (2005).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR45\" id=\"ref-link-section-d217346668e1387\" rel=\"nofollow noopener\" target=\"_blank\">45<\/a> Nycodenz-isolated bacteria were resuspended in lysis buffer to 10 OD600. The bacterial suspension was mixed 1:1 with 1.5% UltraPure low-melting-point agarose (Thermo Fisher Scientific) prepared in lysis buffer and cooled to 50\u2009\u00b0C. The mixture was used to cast gel plugs in CHEF Mapper XA system 50-well plug molds (BioRad, 1703713). Gel plugs were then suspended in a solution of 5\u2009mg\u2009ml\u22121 lysozyme and 1\u2009mg\u2009ml\u22121 achromopeptidase prepared in lysis buffer and gently rotated at 37\u2009\u00b0C overnight using a Roto-therm plus (Benchmark Scientific). The buffer was decanted and the plugs were topped up with TESP buffer (50\u2009mM Tris, 100\u2009mM EDTA, 1% w\/v SDS and 1\u2009mg\u2009ml\u22121 proteinase K) and incubated overnight with gentle rotation at 50\u2009\u00b0C. The gel plug suspension was decanted and the plugs were topped up with fresh TESP and once again incubated overnight with gentle rotation at 50\u2009\u00b0C. TESP was then decanted and the gel plugs were rinsed several times with milli-Q water before being suspended in T10E50 (10\u2009mM Tris and 50\u2009mM EDTA, pH 8) with 1\u2009mM PMSF and incubated at room temperature for 2\u2009h. The liquid was decanted and the plugs washed four times with T10E50 with each wash lasting 1\u2009h on ice. Gel plugs were then stored suspended in T10E50 in the fridge.<\/p>\n<p>Direct lysis<\/p>\n<p>Cell lysis was performed similarly to previously established approaches<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 46\" title=\"Morita, H. et al. An improved DNA isolation method for metagenomic analysis of the microbial flora of the human intestine. Microbes Environ. 22, 214&#x2013;222 (2007).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR46\" id=\"ref-link-section-d217346668e1407\" rel=\"nofollow noopener\" target=\"_blank\">46<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 47\" title=\"Barsotti, O. et al. Rapid isolation of DNA from Actinomyces. Ann. Inst. Pasteur Microbiol. 138, 529&#x2013;536 (1987).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR47\" id=\"ref-link-section-d217346668e1410\" rel=\"nofollow noopener\" target=\"_blank\">47<\/a>. Nycodenz-isolated bacteria were resuspended in lysis buffer to 10 OD600 and 500\u2009\u00b5l of the suspension was pelleted by centrifugation for 1\u2009min at 20,000g. The cell pellet was resuspended in 800\u2009\u00b5l of lysis buffer containing 5\u2009mg\u2009ml\u22121 lysozyme and 1\u2009mg\u2009ml\u22121 achromopeptidase and the suspension was incubated at 37\u2009\u00b0C for 2\u2009h. After incubation, 1\u2009mg of proteinase K powder was added, followed by 200\u2009\u00b5l of 5% w\/v SDS. The sample was very gently inverted to mix and incubated for 3\u2009h at 55\u2009\u00b0C. The viscous solution was then layered onto a sucrose gradient for size selection.<\/p>\n<p>Monarch HMW DNA extraction kit for tissue<\/p>\n<p>About 1\u2009ml of 10-OD600 nycodenz-isolated bacteria was centrifuged at 20,000g for 1\u2009min and the cell pellet was resuspended in 300\u2009\u00b5l of TE lysis buffer containing 5\u2009mg\u2009ml\u22121 lysozyme and 1\u2009mg\u2009ml\u22121 achromopeptidase. The sample was incubated for 2\u2009h at 37\u2009\u00b0C and then carried forward using the New England Biolabs protocol for HMW DNA extraction from bacteria using the Monarch HMW DNA extraction kit for tissue (New England Biolabs, T3060). Briefly, 20\u2009\u00b5l of proteinase K was added and the sample was incubated at 56\u2009\u00b0C with varied shaking speeds (1,900, 1,700, 1,400 or 650\u2009rpm) using a Fisherbrand Heat\/Cool Thermal Mixer II (Fisher Scientific). After 30\u2009min, 10\u2009\u00b5l of RNase A was added and the sample was incubated a further 10\u2009min with shaking at the same speed. Then, 300\u2009\u00b5l of protein separation solution was added and the sample was mixed by inversion for 1\u2009min. The sample was centrifuged for 10\u2009min at 16,000g and the upper phase was then transferred to a Monarch 2-ml tube using a wide-bore pipette tip. Two DNA capture beads were added, followed by 550\u2009\u00b5l of isopropanol. The sample was mixed by gentle rotation using the Roto-therm Plus on the lowest speed setting for 5\u2009min. The liquid was then discarded and the beads were washed twice with 500\u2009\u00b5l of genomic DNA wash buffer. Finally, DNA was eluted by addition of 100\u2009\u00b5l of elution buffer II and incubation at 56\u2009\u00b0C for 10\u2009min with 300-rpm agitation. The DNA solution was left to resuspend at 4\u2009\u00b0C overnight.<\/p>\n<p>DNA size-selection methodsElectroelution<\/p>\n<p>The electroelution protocol was performed as previously described<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 18\" title=\"Brady, S. F. Construction of soil environmental DNA cosmid libraries and screening for clones that produce biologically active small molecules. Nat. Protoc. 2, 1297&#x2013;1305 (2007).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR18\" id=\"ref-link-section-d217346668e1457\" rel=\"nofollow noopener\" target=\"_blank\">18<\/a> with minor modifications. Briefly, a 0.9% agarose gel was run in 0.5\u00d7 TBE (VWR) at 6\u2009V\u2009cm\u22121 for 2\u2009h to size-select crude DNA extracts. The compression band containing HMW DNA was then excised, the gel slice was placed in dialysis tubing (Spectra\/Por; 12,000\u201314,000-kDa molecular weight cutoff (MWCO)) and topped up with 0.5\u00d7 TBE. DNA was then electroeluted in 0.5\u00d7 TBE at 6\u2009V\u2009cm\u22121 for 2\u2009h. The gel slice was removed from the dialysis bag and the remaining liquid was concentrated to dryness against 25% PEG 8000 with 5\u2009mM Tris pH 8. DNA was then rehydrated in the dialysis bag against 5\u2009mM Tris pH 8 for isolation.<\/p>\n<p>Sucrose gradient<\/p>\n<p>Linear 10\u201340% sucrose gradients, with 10\u2009mM Tris pH 8, were prepared using the freeze\u2013thaw method<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Luthe, D. S. A simple technique for the preparation and storage of sucrose gradients. Anal. Biochem. 135, 230&#x2013;232 (1983).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR48\" id=\"ref-link-section-d217346668e1473\" rel=\"nofollow noopener\" target=\"_blank\">48<\/a>. Gradients were prepared by liquid nitrogen flash-freezing at the end of the day and immediately placed to thaw at 4\u2009\u00b0C until use 22\u201324\u2009h later. Two sizes of polypropylene tubes were used, 25\u2009\u00d7\u200989\u2009mm and 14\u2009\u00d7\u200995\u2009mm (Beckman), with 8-ml and 2.9-ml sucrose layers, respectively. Larger tubes were centrifuged using an SW28 rotor, while the smaller tubes were centrifuged using an SW40 Ti rotor. In both cases, centrifugation was at 4\u2009\u00b0C, 40,000g for 16\u2009h. After centrifugation, fractions were gently pipetted from the top of the tubes using wide-bore tips. Pulse-field gel electrophoresis was used to analyze the fractions. For shorter-fragment purification (that is, DNA from hot CTAB and Monarch 1,900\u2009rpm), analysis was conducted using 0.9% agarose with 1\u00d7 TAE at 6\u2009V\u2009cm\u22121 for 4\u2009h with a 5-s switch time. For ultralarge DNA, analysis was conducted using 0.9% agarose with 1\u00d7 TAE at 6\u2009V\u2009cm\u22121 for 16\u2009h with a linear 3\u20137-s switch time. Ultralarge DNA from direct lysis was isolated using large tubes from fractions at 14\u2009ml to 32\u2009ml. Hot CTAB DNA was purified using small tubes from fractions at 2.5\u2009ml to 4.5\u2009ml. DNA from Monarch 1,900\u2009rpm was purified using large tubes from fractions at 6\u2009ml to 8.5\u2009ml, aiming to isolate DNA ~40\u2009kb in size.<\/p>\n<p>Short fragment eliminator (SFE) kit<\/p>\n<p>DNA size selection with the SFE kit (Oxford Nanopore EXP-SFE001) was performed according to the manufacturer\u2019s instructions.<\/p>\n<p>A-domain survey<\/p>\n<p>DNA from gel-plug-extracted metagenomic DNA was isolated using the NucleoSpin gel and PCR cleanup kit (Macherey-Nagel) according to the manufacturer\u2019s instructions. A-domain amplicon sequencing was performed as previously described<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 49\" title=\"Libis, V. et al. Uncovering the biosynthetic potential of rare metagenomic DNA using co-occurrence network analysis of targeted sequences. Nat. Commun. 10, 3848 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR49\" id=\"ref-link-section-d217346668e1501\" rel=\"nofollow noopener\" target=\"_blank\">49<\/a>. Read counts were 1,424,301 for forest soil 1, 923,413 for forest soil 2, 551,873 for forest soil 3, 1,965,137 for forest soil 4, 717,416 for the muddy bog and 467,043 for the decaying tree. A-domain amplicons were trimmed to remove primer sequences using vsearch (version 2.28.1)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 50\" title=\"Rognes, T., Flouri, T., Nichols, B., Quince, C. &amp; Mahe, F. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4, e2584 (2016).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR50\" id=\"ref-link-section-d217346668e1505\" rel=\"nofollow noopener\" target=\"_blank\">50<\/a> fastx_filter and then assembled into operational taxonomic units (OTUs) clustered at 95% nucleotide similarity using vsearch cluster_size. A-domain OTUs for each soil were then compared across all samples at 95% similarity. Parallel processing was achieved using GNU parallel (version 20230922)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 51\" title=\"Tange, O. GNU Parallel 2018. Zenodo &#010;                https:\/\/doi.org\/10.5281\/zenodo.1146014&#010;                &#010;               (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR51\" id=\"ref-link-section-d217346668e1509\" rel=\"nofollow noopener\" target=\"_blank\">51<\/a>.<\/p>\n<p>Nanopore ultralong DNA sequencing<\/p>\n<p>Ultralong DNA sequencing kit V14 (SQK-ULK114) preparation was performed according to the manufacturer\u2019s instructions. Three DNA inputs were attempted, one following Oxford Nanopore directions on the Monarch HMW DNA extraction kit for tissue extraction, the manufacturer-recommended protocol for Monarch HMW DNA extraction kit for tissue and sucrose gradient purification. An additional attempt using the Monarch extracted DNA input was tried with a slightly altered protocol; instead of heat-killing the fragmentation mix, proteinase K was added to 1\u2009mg\u2009ml\u22121 for 1\u2009h at room temperature, followed by the addition of 3\u2009mM PMSF for 1\u2009h at room temperature to inhibit the proteinase K before continuing the protocol.<\/p>\n<p>Nanopore ligation kit sequencing<\/p>\n<p>Oxford Nanopore\u2019s ligation sequencing kit V14 (SQK-LSK114) was used to prepare sequencing samples, with DNA extraction and size selection as indicated, according to the manufacturer\u2019s instructions and minor modifications. Input DNA was increased to 3\u2009\u00b5g, blunting and end-repair mix reaction incubation times were extended to 15\u2009min at 20\u2009\u00b0C and 10\u2009min at 65\u2009\u00b0C and adaptor ligation incubation time was increased to 30\u2009min. Library solution was used for all samples and the DNA concentrations loaded are indicated in Supplementary Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>. On the basis of our optimization experiments, the 1,700-rpm Monarch isolation and SFE size-selection protocol applied to 62.5\u2009g of the forest soil would be sufficient to generate &gt;2\u2009Tbp of sequence.<\/p>\n<p>Bioinformatic analysisSequencing<\/p>\n<p>PromethION R10.4.1 flow cells were used to generate all of the sequencing data. Nanopore sequencing was conducted using a P2 solo with MinKnow 23.04.6 (ultra kits) or 23.11.4 (Ligation kits) and a P24 with MinKnow 23.07.12 (scaled ligation). All raw data were duplex-basecalled with standalone Dorado 0.4.1\u2009+\u20096c4c636 for subsequent analysis.<\/p>\n<p>Global assembly<\/p>\n<p>The longest, highest-quality reads (&gt;20\u2009kbp, Q20+) were extracted from the 2.5-Tbp sequence dataset using Chopper (version 0.6.0)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 52\" title=\"De Coster, W. &amp; Rademakers, R.NanoPack2: population-scale evaluation of long-read sequencing data. Bioinformatics 39, btad311 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR52\" id=\"ref-link-section-d217346668e1554\" rel=\"nofollow noopener\" target=\"_blank\">52<\/a>. The resulting 528.5-Gbp dataset was assembled using Flye (version 2.9.3-b1797) in metagenomic mode (metaFlye) with parameters \u2018&#8211;nano-hq &#8211;i 0 &#8211;threads 64 &#8211;meta\u2019 (ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 19\" title=\"Kolmogorov, M. et al. metaFlye: scalable long-read metagenome assembly using repeat graphs. Nat. Methods 17, 1103&#x2013;1110 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR19\" id=\"ref-link-section-d217346668e1558\" rel=\"nofollow noopener\" target=\"_blank\">19<\/a>). The assembly was run on Rockefeller University\u2019s large-memory high-performance computing node (64 cores, 3\u2009TB of RAM) and took just under 8\u2009days to complete with RAM usage peaking at 845\u2009GB. The assembly was polished once (the emerging best practice<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 53\" title=\"Luan, T. et al. Benchmarking short and long read polishing tools for nanopore assemblies: achieving near-perfect genomes for outbreak isolates. BMC Genomics 25, 679 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR53\" id=\"ref-link-section-d217346668e1562\" rel=\"nofollow noopener\" target=\"_blank\">53<\/a>) with Medaka (version 1.11.3; <a href=\"https:\/\/github.com\/nanoporetech\/medaka\" rel=\"nofollow noopener\" target=\"_blank\">https:\/\/github.com\/nanoporetech\/medaka<\/a>) and annotated by Bakta (version 1.10.4)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 54\" title=\"Schwengers, O. et al. Bakta: rapid and standardized annotation of bacterial genomes via alignment-free sequence identification. Microb. Genom. 7, 000685 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR54\" id=\"ref-link-section-d217346668e1573\" rel=\"nofollow noopener\" target=\"_blank\">54<\/a>. Assembled contigs\u2009&gt;\u20091\u2009Mbp were extracted and their taxonomy assigned using the GTDB toolkit (GTDB-tk version 2.4.0)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 55\" title=\"Chaumeil, P. A., Mussig, A. J., Hugenholtz, P. &amp; Parks, D. H. GTDB-Tk v2: memory friendly classification with the genome taxonomy database. Bioinformatics 38, 5315&#x2013;5316 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR55\" id=\"ref-link-section-d217346668e1578\" rel=\"nofollow noopener\" target=\"_blank\">55<\/a> classify workflow with GTDB release 09-RS220. The presence of rRNA was detected within the Bakta annotation using a GNU AWK script (version 5.3.0), while tRNAscan-SE (version 2.0.9) was used to quantify total unique tRNAs. A set of 563 contigs\u2009&gt;\u20091\u2009Mbp were identified containing all rRNAs and at least 18 unique tRNAs to establish the complete or near-complete metagenomic genome dataset. The dataset was then frame-shift-corrected using proovframe<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 56\" title=\"Hackl, T. et al. proovframe: frameshift-correction for long-read (meta)genomics. Preprint at bioRxiv &#010;                https:\/\/doi.org\/10.1101\/2021.08.23.457338&#010;                &#010;               (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR56\" id=\"ref-link-section-d217346668e1582\" rel=\"nofollow noopener\" target=\"_blank\">56<\/a> and completeness was estimated by CheckM (version 1.2.3)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 57\" title=\"Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. &amp; Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043&#x2013;1055 (2015).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR57\" id=\"ref-link-section-d217346668e1586\" rel=\"nofollow noopener\" target=\"_blank\">57<\/a>. The taxonomic classification of these contigs was plotted as a Sankey diagram using the networkD3 R package 0.4 (Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2a<\/a>). AntiSMASH (version 6.1.1)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 58\" title=\"Blin, K. et al. antiSMASH 6.0: improving cluster detection and comparison capabilities. Nucleic Acids Res. 49, W29&#x2013;W35 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR58\" id=\"ref-link-section-d217346668e1593\" rel=\"nofollow noopener\" target=\"_blank\">58<\/a> results from the output JSON files were used to generate the appended BGC summaries. Subtrees expanding previously underexplored clades were generated by extracting the alignment information from the GTDB-tk classify workflow for desired sequences and plotting with FastTree (version 2.1.11)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 59\" title=\"Price, M. N., Dehal, P. S. &amp; Arkin, A. P. FastTree 2&#x2014;approximately maximum-likelihood trees for large alignments. PLoS ONE 5, e9490 (2010).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR59\" id=\"ref-link-section-d217346668e1597\" rel=\"nofollow noopener\" target=\"_blank\">59<\/a> using default settings. Additional information was then manually added.<\/p>\n<p>Targeted NRPS assembly<\/p>\n<p>HMMER (version 3.1b2)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 60\" title=\"Eddy, S. R. Accelerated profile HMM searches. PLoS Comput. Biol. 7, e1002195 (2011).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR60\" id=\"ref-link-section-d217346668e1609\" rel=\"nofollow noopener\" target=\"_blank\">60<\/a> was used to scan the total sequencing dataset for AMP-binding domains using an available HMM (AMP-binding, Pfam <a href=\"https:\/\/www.ebi.ac.uk\/interpro\/entry\/pfam\/PF00501\/\" rel=\"nofollow noopener\" target=\"_blank\">PF00501.23<\/a>). Reads containing AMP-binding domains were assembled (metaFlye), polished (Medaka) and annotated (Bakta), as described above. The AMP-binding domain search using HMMER was repeated on the assemblies to extract the subset containing at least four AMP-binding domains. Similar to the complete and near-complete genome assembly dataset, frameshifts were minimized using proovframe and then the subset was analyzed by antiSMASH (version 6.1.1). For analysis of NRPS content, the antiSMASH results were analyzed by BiG-SCAPE (version 1.1.5) to generate Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3a<\/a>. The phylogeny of NRPS containing BGCs was determined using the mmseqs2 (ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 61\" title=\"Steinegger, M. &amp; Soding, J. MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets. Nat. Biotechnol. 35, 1026&#x2013;1028 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR61\" id=\"ref-link-section-d217346668e1623\" rel=\"nofollow noopener\" target=\"_blank\">61<\/a>) toolkit\u2019s easy-taxonomy workflow and the National Center for Biotechnology Information nonredundant protein database.<\/p>\n<p>16S rRNA analysis<\/p>\n<p>Reads containing 16S rRNA sequences were extracted from the 2.5\u2009Tbp of sequencing data using barrnap 0.9 (<a href=\"https:\/\/github.com\/tseemann\/barrnap\" rel=\"nofollow noopener\" target=\"_blank\">https:\/\/github.com\/tseemann\/barrnap<\/a>). The 16S rRNA read subset was then assembled (metaFlye), polished (Medaka) and annotated (Bakta) as described above. Annotated 16S rRNA gene sequences were then extracted and, using vsearch, dereplicated and clustered at 99%, the recommended threshold for clustering full-length 16S sequences<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 62\" title=\"Edgar, R. C. Updating the 97% identity threshold for 16S ribosomal RNA OTUs. Bioinformatics 34, 2371&#x2013;2375 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR62\" id=\"ref-link-section-d217346668e1642\" rel=\"nofollow noopener\" target=\"_blank\">62<\/a>. Lastly, to ensure full-length sequences, only OTUs between 1,400 and 1,700\u2009bp were retained. The same extraction, dereplication and clustering process was used to identify 16S rRNA genes within the global assembly dataset. Phylogeny was assigned using mmseqs2 and the SILVA database (release 138.2)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 63\" title=\"Yilmaz, P. et al. The SILVA and &#x2018;All-Species Living Tree Project (LTP)&#x2019; taxonomic frameworks. Nucleic Acids Res. 42, D643&#x2013;D648 (2014).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR63\" id=\"ref-link-section-d217346668e1646\" rel=\"nofollow noopener\" target=\"_blank\">63<\/a>.<\/p>\n<p>Forest soil 1 cultured isolatesCulturing<\/p>\n<p>Forest soil 1 microorganisms were isolated using the nycodenz method with the last two washes replaced with Dulbecco\u2019s PBS (DPBS; Gibco). The purified suspension was diluted to 0.05 OD600 in DPBS and 200\u2009\u00b5l was spread on individual agar plates (150\u2009mm\u2009\u00d7\u200915\u2009mm). Ten plates each were used with three different media: R2A agar (Fisher Scientific), which was previously used to isolate the Aquella species<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 27\" title=\"Chan, K. G. et al. Aquella oligotrophica gen. nov. sp. nov.: a new member of the family Neisseriaceae isolated from laboratory tap water. Microbiologyopen 8, e00793 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR27\" id=\"ref-link-section-d217346668e1668\" rel=\"nofollow noopener\" target=\"_blank\">27<\/a>, and two versions of dilute nutrient broth (Criterion). In one case, nutrient broth was diluted 200-fold compared to the manufacturer\u2019s recommendation; in the second case, dilute nutrient broth was with supplemented 1:1 (v\/v) with soil extract<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 64\" title=\"George, I. F., Hartmann, M., Liles, M. R. &amp; Agathos, S. N. Recovery of as-yet-uncultured soil acidobacteria on dilute solid media. Appl. Environ. Microbiol. 77, 8184&#x2013;8188 (2011).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR64\" id=\"ref-link-section-d217346668e1672\" rel=\"nofollow noopener\" target=\"_blank\">64<\/a>. In all cases, nystatin (50\u2009\u00b5g\u2009ml\u22121) was added to limit fungal growth. Plates were incubated at 30\u2009\u00b0C with R2A colonies picked after 3\u2009days and nutrient broth colonies picked after 10\u2009days. Attempting to maximize isolate diversity, colonies were picked to maximize different morphologies. Liquid R2A or unmodified nutrient broth were used to culture isolates for sequencing (3\u2009ml, 30\u2009\u00b0C, 200-rpm shaking).<\/p>\n<p>Genomic DNA extraction and sequencing<\/p>\n<p>Bacterial cells were pelleted from turbid cultures (1\u20131.5\u2009ml) by centrifugation for 1\u2009min at 20,000g. The supernatant was removed and the cell pellets were caried forward for genomic DNA isolation using the PureLink microbiome DNA purification kit (Thermo Fisher Scientific) as per manufacturer\u2019s instructions for soil samples with elution volume lowered to 50\u2009\u00b5l. Oxford Nanopore\u2019s native barcoding kit 24 V14 (SQK-NBD114.24) was used as per the manufacturer\u2019s instructions to prepare sequencing samples. From R2A medium culturing, 22 genomes were individually barcoded and sequenced. From nutrient medium culturing, 20 genomes were individually barcoded and 12 were barcoded as pools of three. Each medium set was sequenced using a PromethION flow cell.<\/p>\n<p>Bioinformatics<\/p>\n<p>Isolate genomes were assembled using Flye with Medaka polishing and Bakta annotation. Genome assemblies were dereplicated by extracting 16S rRNA genes, clustering sequences at 99% identity and then selecting a single assembly per cluster. This yielded 30 genomes (Supplementary Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">7<\/a>). BGCs were identified in each dereplicated cultured genome using antiSMASH. Large NRPS BGCs (that is, \u22654 A domains) were compared to those from the targeted metagenomic assembly using mmseqs2. Matches are shown in Supplementary Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>.<\/p>\n<p>SynBNP predictions<\/p>\n<p>SynBNP prediction was carried out according to our standard workflow<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 65\" title=\"Wang, Z., Koirala, B., Hernandez, Y., Zimmerman, M. &amp; Brady, S. F. Bioinformatic prospecting and synthesis of a bifunctional lipopeptide antibiotic that evades resistance. Science 376, 991&#x2013;996 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR65\" id=\"ref-link-section-d217346668e1713\" rel=\"nofollow noopener\" target=\"_blank\">65<\/a>. Briefly, \u2018bioinformatically tractable clusters\u2019 were identified and their products were predicted as follows. BGCs were scanned for NRPS domains and modules. The order in which these domains are likely to function was predicted on the basis of canonical NRPS module organization combined with the principle of collinearity. The following BGCs were then removed: (1) BGCs lacking predicted initiation or termination domains; (2) BGCs with modules that do not follow canonical NRPS domain organization (that is, condensation domain, A-domain and thiolation domain); (3) BGCs rich in tailoring genes; and (4) BGCs having fewer than four A domains, as these make poor synBNP targets. Structure predictions were then generated on the basis of A-domain specificity and derived structural modifications predicted on the basis of an analysis of additional domains present in the BGC (for example, epimerization and N-methylation). These data were combined on the basis of the initial domain order analysis to generate the linear structure encoded by a BGC. Specificity codes were compared to both the antiSMASH output and an in-house specificity code table collated from previously characterized natural products. The highest A-domain matches were solely considered for predictions. In rare cases where multiple amino acids share the same code, one was selected from the consensus of predictions generated by NRPSPredictor2 within antiSMASH. If multiple predictions were tied after the consensus stage, one was selected at random. The set of bioinformatically tractable clusters considered and their corresponding predictions are given in Supplementary Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>. This represented roughly 19% of detected NRPS BGCs. On average, A domains matched at 90.1%\u2009\u00b1\u200913.9%. The average A-domain match per peptide was 89.9%\u2009\u00b1\u20098.0%. From the subset of predictions that contained at most one &lt;70% A-domain specificity code (~83% of the predictions passed this cutoff), we randomly selected 20for synthesis. These covered a range of phylogenetic classifications (Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3b<\/a> and Supplementary Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">9<\/a>). Contig_44957 had two mismatched codes but was retained as the codes were closely related and predicted to specify the same amino acid. While we used a cutoff of NRPs containing at most a single A-domain that matched &lt;70%, this can be loosened or more restrictive if desired. A breakdown of predictions that passed different cutoff levels is provided in Supplementary Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>.<\/p>\n<p>SynBNP synthesis<\/p>\n<p>Reagents and solvents used for synthesis were obtained from commercial sources and used without further purification. Chromatography solvents were of high-performance liquid chromatography (HPLC) grade or higher. Peptides were purified using a CombiFlash EZ Prep purification system with ultraviolet detection and equipped with a RediSep gold HP C18 column (30\u2009g), RediSep gold HP C8 column or Phenomenex Luna 5\u2009\u00b5m C18 prepHPLC column using a 5\u201395% acetonitrile\u2013water gradient supplemented with 0.1% formic acid. High-resolution mass spectrometry (HRMS) data were acquired on a SCIEX ExionLC ultra-HPLC (UPLC) instrument coupled to an X500R quadrupole time-of-flight mass spectrometer with a Phenomenex Kinetex PS C18 100-\u00c5 column (2.1\u2009mm\u2009\u00d7\u200950\u2009mm, 2.6\u2009\u00b5m) and analyzed using SCIEXOS software. 1H and 13C nuclear magnetic resonance (NMR) spectra were obtained on a Bruker Avance DMX 600-MHz spectrometer equipped with cryogenic probe (Rockefeller University) and all spectra were analyzed using MestRe Nova software.<\/p>\n<p>Linear peptides were synthesized using standard Fmoc-based solid-phase peptide synthesis (SPPS) methods. 2-Chlorotritylchloride resin (0.1\u2009mmol) was swollen in dichloromethane (DCM) at room temperature for 30\u2009min. After draining, the resin was washed with DCM (three washes, 5\u2009ml each). The C-terminal Fmoc-protected amino acid (four equivalents, 0.4\u2009mmol) was dissolved in 5\u2009ml of DCM and treated with N,N-diisopropylethylamine (DIPEA; five equivalents, 0.5\u2009mmol). This mixture was added to the swollen resin and agitated for 1\u2009h. The solution was drained and the resin was washed with DCM (three washes, 5\u2009ml each). The resin was then capped by treatment with a 17:2:1 mixture of DCM, methanol and DIPEA. The resin was agitated for 30\u2009min, drained and washed with DCM (three washes, 5\u2009ml each) and dimethylformamide (DMF; three washes, 5\u2009ml each). For peptide couplings, the Fmoc group was removed by treating the resin with 20% piperidine in 5\u2009ml of DMF for 7\u2009min. The solution was drained, the deprotection reaction was repeated and the resin was washed with DMF (five washes, 5\u2009ml each). The next Fmoc amino acid to be coupled (three equivalents, 0.3\u2009mmol) and hexafluorophosphate azabenzotriazole tetramethyl uronium (HATU; three equivalents, 0.3\u2009mmol) was dissolved in 4\u2009ml of DMF and treated with DIPEA (three equivalents, 0.3\u2009mmol). This solution was added to the resin and agitated for 45\u2009min. The solution was drained and the resin was washed with DMF (three washes, 5\u2009ml each). This deprotection and coupling process was repeated for each subsequent amino acid. Finally, the N-terminal Fmoc amino acid was deprotected with 20% piperidine in DMF (two cycles of 5\u2009ml for 7\u2009min) and reacted with myristic acid (three equivalents, 0.3\u2009mmol), HATU (three equivalents, 0.3\u2009mmol) and DIPEA (three equivalents, 0.3\u2009mmol) in DMF (4\u2009ml). The reaction was agitated for 45\u2009min, at which point the resin was washed with DMF (three washes, 5\u2009ml each) and DCM (three washes, 5\u2009ml each). The peptide was cleaved from resin and the side-chain-protecting groups were removed by treating the resin with a trifluoroacetic acid (TFA) cleavage cocktail containing 95% TFA, 2.5% water and 2.5% triisopropylsilane for 2\u2009h. The solution was evaporated under air flow and the resulting residue was purified on a CombiFlash EZ Prep system (30-g HP C18 column, water\u2013acetonitrile with 0.1% v\/v formic acid 5%\u201395%).<\/p>\n<p>Peptides cyclized through the side chain of an amino acid were synthesized starting from the penultimate amino acid following the same Fmoc SPPS loading, Fmoc deprotection and coupling procedures described above. Side-chain serine and threonine residues at which cyclization was predicted to occur were added to the peptide with the side chain unprotected. After the N-terminal amino acid was deprotected, myristic acid was coupled to the terminal amine as described above (for peptides not predicted to have a lipid tail, the Boc-protected N-terminal amino acid was coupled). Next, the C-terminal amino acid was esterified onto the free serine or threonine side chain by adding a solution of Fmoc amino acid (15 equivalents, 1.5\u2009mmol), N,N\u2032-diisopropylcarbodiimide (DIC; 15 equivalents, 1.5\u2009mmol) and 4-dimethylaminopyridine (DMAP; 0.5 equivalents, 0.05\u2009mmol) in 5\u2009ml of DMF to the resin and agitating for 16\u2009h. The solution was drained and the resin was washed with DMF (three washes, 5\u2009ml each) and DCM (three washes, 5\u2009ml each), followed by deprotection of the Fmoc group using 20% piperidine in DMF (two cycles of 5\u2009ml for 7\u2009min). The linear precursor peptide was gently cleaved from resin by treating with a solution of 20% hexafluoroisopropanol in 6\u2009ml of DCM for 45\u2009min. The solution was collected and the process repeated once more. After removal of solvent under air flow, the resulting peptide was dissolved in DMF (50\u2009ml, 0.002\u2009M) and treated with (7-azabenzotriazol-1-yloxy)tripyrrolidinophosphonium hexafluorophosphate (seven equivalents, 0.7\u2009mmol) and DIPEA (30 equivalents, 3\u2009mmol). After stirring for 2\u2009h, the solution was diluted with ethyl acetate and washed with brine (three washes, 100\u2009ml each). The organic layer was dried (Na2SO4), filtered and concentrated in vacuo. The resulting cyclic peptide was globally deprotected by treatment with TFA cleavage cocktail for 2\u2009h. After evaporation under air flow, the crude cyclic peptide was purified on a CombiFlash EZ Prep system as described above.<\/p>\n<p>Peptides cyclized through the fatty acid were synthesized starting from the penultimate amino acid residue following the Fmoc SPPS method described above. After synthesis of the linear precursor was completed, the N-terminal Fmoc group was removed using 20% piperidine in DMF (two cycles of 5\u2009ml for 7\u2009min) and reacted with (\u00b1)-3-hydroxymyristic acid (three equivalents, 0.3\u2009mmol), HATU (three equivalents, 0.3\u2009mmol) and DIPEA (three equivalents, 0.3\u2009mmol) in DMF (4\u2009ml). The reaction was agitated for 45\u2009min, at which point the resin was washed with DMF (three washes, 5\u2009ml each) and DCM (three washes, 5\u2009ml each). Next, the C-terminal amino acid was esterified onto the hydroxyl group of the fatty acid by adding a solution of Fmoc amino acid (15 equivalents, 1.5\u2009mmol), DIC (15 equivalents, 1.5\u2009mmol) and DMAP (0.5 equivalents, 0.05\u2009mmol) in 5\u2009ml of DMF to the resin and agitating for 16\u2009h. The solution was drained and the resin was washed with DMF (three washes, 5\u2009ml each) and DCM (three washes, 5\u2009ml each), followed by deprotection of the Fmoc group using 20% piperidine in DMF (two cycles of 5\u2009ml for 7\u2009min). The linear precursor was cleaved from resin, cyclized, deprotected and purified as described above.<\/p>\n<p>Successful peptide synthesis was confirmed by UPLC, HRMS and\/or NMR (Supplementary Table <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">15<\/a> and Supplementary Figs. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a>\u2013<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">9<\/a>). All compounds that displayed antibacterial activity were confirmed to have &gt;95% purity.<\/p>\n<p>MIC assay<\/p>\n<p>To avoid solubility issues, each synBNP concentration was prepared individually in DMSO before addition to LB. A series of eight concentrations were tested, with the last one being a DMSO only control. First, the synBNP was diluted to 3.2\u2009mg\u2009ml\u22121 in DMSO and then serially diluted twofold six times in DMSO. Then, 4\u2009\u00b5l of each concentration was added to 100\u2009\u00b5l of LB dispensed across a 96-well plate (Costar, 3370). Test bacteria were grown overnight shaking at 200\u2009rpm at 30\u2009\u00b0C and the culture was diluted to an OD600 of 0.005. Next, 100\u2009\u00b5l of the dilution was added to the antibiotic dilution series, yielding a test range of 64 to 1\u2009\u00b5g\u2009ml\u22121 along with a control (0\u2009\u00b5g\u2009ml\u22121). The plate was then incubated statically at 30\u2009\u00b0C for 18\u2009h, followed by visual inspection. The MIC was determined as the lowest concentration within the dilution series that did not exhibit turbidity. MICs were determined by two biologically independent replicates for the initial synBNP screen and three independent replicates for all other determinations.<\/p>\n<p>Human cell toxicity assay<\/p>\n<p>The inhibitory activity of synBNPs was assessed using a 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 66\" title=\"Heo, D. S. et al. Evaluation of tetrazolium-based semiautomatic colorimetric assay for measurement of human antitumor cytotoxicity. Cancer Res. 50, 3681&#x2013;3690 (1990).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR66\" id=\"ref-link-section-d217346668e1807\" rel=\"nofollow noopener\" target=\"_blank\">66<\/a>. HEK239 and HeLa cells were cultured in DMEM supplemented with 10% FBS, 10\u2009\u00b5g\u2009ml\u22121 penicillin\u2013streptomycin, 2\u2009\u00b5M glutamine and nonessential amino acids. Cells were passaged and seeded in 96-well plates at the exponential growth phase, inoculated at ~1\u2009\u00d7\u2009105 cells per ml and maintaining a cell density below 1\u2009\u00d7\u2009106 cells per ml. The plates were incubated at 37\u2009\u00b0C with 5% CO2. All assays were performed using 5,000 cells per well, with 50\u2009\u00b5l of a homogeneous cell suspension at 1\u2009\u00d7\u2009105 cells per ml seeded in each well. For screening, synBNPs were diluted to 6.4\u2009\u00b5g\u2009ml\u22121 in DMSO to prepare a 200\u00d7 stock solution. For the high-concentration assays (32\u2009\u00b5g\u2009ml\u22121), 1.5\u2009\u00b5l of synBNP stocks were used to prepare a master solution by dilution in 150\u2009\u00b5l of DMEM without phenol red. Two 50-\u00b5l aliquots of the master solution were added to cell-seeded plates, obtaining duplicate 100-\u00b5l wells with synBNPs at 32\u2009\u00b5g\u2009ml\u22121 (0.5% DMSO). For the low concentrations (3.2\u2009\u00b5g\u2009ml\u22121), 15\u2009\u00b5l of the master solution was diluted to 150\u2009\u00b5l with DMEM and 50\u2009\u00b5l was added to cell-seeded plates to obtain duplicate wells of 100\u2009\u00b5l at 3.2\u2009\u00b5g\u2009ml\u22121 (0.05% DMSO). DMEM solutions containing 0.5% and 0.05% DMSO were used as vehicle controls to define 100% cell viability. Camptothecin at 40\u2009\u00b5g\u2009ml\u22121 was used as a positive control to define 0% cell viability. After treatment, plates were incubated for 48\u2009h followed by addition of 10\u2009\u00b5l of 5\u2009mg\u2009ml\u22121 MTT solution in DPBS and the plate was incubated for an additional 4\u2009h (SK-N-SH was incubated for 5\u20136\u2009h). After incubation, 90\u2009\u00b5l of solubilization solution (40% DMF, 16% SDS and 2% acetic acid) was added to each well and formazan crystals were allowed to dissolve overnight. The absorbance of each well was measured at 570\u2009nm using a TECAN Infinite M NANO+. Absorbance values were normalized to the positive and negative controls and percentage survival values were calculated.<\/p>\n<p>Kill curve assay<\/p>\n<p>S.\u2009aureus SH1000 was streaked onto a LB agar plate and incubated overnight at 30\u2009\u00b0C. Three independent colonies were then inoculated into 3\u2009ml of LB and grown overnight at 30\u2009\u00b0C with shaking at 200\u2009rpm. The cultures were diluted 1:100 and grown at 30\u2009\u00b0C and 200\u2009rpm until they reached an OD600 between 0.4 and 0.6. The growing cultures were then diluted into a preprepared series of 3\u2009ml of LB containing the indicated antibiotic concentrations to ~6\u2009\u00d7\u2009106 colony-forming units (cfu) per ml using the expectation that 1 OD600 of S. aureus contained 5\u2009\u00d7\u2009108 cfu per ml. When a series was inoculated, a 100-\u00b5l sample was taken, representing 1\u2009min of treatment. The cultures were then placed shaking at 200\u2009rpm and 30\u2009\u00b0C and samples were taken at 10\u2009min, 30\u2009min, 1\u2009h and then every subsequent hour for a total of 6\u2009h. For cfu quantification, 100-\u00b5l samples were placed immediately into the top well of a 96-well plate and serially diluted tenfold (20\u2009\u00b5l into 180\u2009\u00b5l of LB) down to 10\u22127. Then, 6\u2009\u00b5l of each dilution from the series was spotted onto LB agar, the spots were air-dried and the plate was incubated overnight at 30\u2009\u00b0C. Resulting colonies were enumerated to calculate the cfu per ml.<\/p>\n<p>Mutant isolationDirect plating<\/p>\n<p>S.\u2009aureus SH1000 was streaked onto an LB agar plate and incubated overnight at 30\u2009\u00b0C. A single colony was inoculated into 3\u2009ml of LB and grown to confluence overnight at 30\u2009\u00b0C. The culture was then diluted to roughly 1\u2009\u00d7\u2009107 cells per ml in LB containing antibiotic at 8\u00d7 MIC. Five 96-well plates were then seeded with 100\u2009\u00b5l of the dilution and grown stationary overnight at 30\u2009\u00b0C. Mutants appeared as small colonies at the bottom of the well and were picked into fresh 8\u00d7 MIC medium to ensure resistance. The confirmation culture was then streaked onto LB agar and incubated overnight at 30\u2009\u00b0C; two colonies per mutant were inoculated for an overnight culture in LB at 30\u2009\u00b0C. The MIC of each mutant was checked to confirm resistance. For erutamicin, no resistant colonies appeared and the direct plating was repeated at 4\u00d7 and 2\u00d7 MIC concentrations, which also did not produce resistant mutants.<\/p>\n<p>Passaging<\/p>\n<p>S.\u2009aureus SH1000 and A.\u2009baumannii ATCC 17978 were streaked onto LB agar and incubated overnight at 30\u2009\u00b0C. Three independent colonies for each of the bacteria were then inoculated into 3\u2009ml of LB and grown overnight at 30\u2009\u00b0C with shaking at 200\u2009rpm. The cultures were used to prepare a MIC assay as described above with the indicated antibiotics. The following day, 1\u2009\u00b5l was taken from the well containing the highest concentration of antibiotic exhibiting turbid growth and mixed into 1\u2009ml of fresh LB. This dilution was used to seed a new MIC assay. This was repeated for each subsequent day and the change in MIC was plotted.<\/p>\n<p>Mutant sequencing and identification<\/p>\n<p>Genomic DNA from resistant mutants was extracted for sequencing using the PureLink microbiome DNA purification kit (Thermo Fisher Scientific) per the manufacturer\u2019s instructions. Isolated DNA was then Illumina-sequenced using standard Nextera preparation and the MiSeq reagent kit v3 with 2\u00d7 300-bp reads. Reads for both parent and mutants were mapped to a reference genome using minimap2 (version 2.28-r1209)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 67\" title=\"Li, H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094&#x2013;3100 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#ref-CR67\" id=\"ref-link-section-d217346668e1912\" rel=\"nofollow noopener\" target=\"_blank\">67<\/a>. Changes detected in the parent strain were eliminated from consideration and mutations were uncovered using snippy (<a href=\"https:\/\/github.com\/tseemann\/snippy\" rel=\"nofollow noopener\" target=\"_blank\">https:\/\/github.com\/tseemann\/snippy<\/a>).<\/p>\n<p>ClpX expression vector<\/p>\n<p>The clpX gene was amplified from S.\u2009aureus SH1000 genomic DNA by PCR using the primers ClpXF (TATTACTCGAGATGTTTAAATTCAATGAAGATGAAG) and ClpXR (TAATAGGATCCACATCAATGATTAAGCTGATG). The pET19b (Sigma) vector backbone was amplified by 19BF (TAATAGGATCCGAAAGGAAGCTGAGTTGG) and 19BR (TATTACTCGAGCATATGCTTGTCGTCGTCGTC). The two PCR products were then digested with XhoI and BamHI and ligated together; the ligation reaction was transformed into E.\u2009coli EPI300 for colony screening on LB agar plates containing 100\u2009\u00b5g\u2009ml\u22121 carbapenem. The successful construction of the ClpX expression vector pHisClpX (which added a 10\u00d7 His tag to the N terminus of ClpX) was confirmed by Sanger sequencing performed by Genewiz using the universal T7 primer.<\/p>\n<p>ClpX expression and purification<\/p>\n<p>The pHisClpX vector was transformed into E.\u2009coli Rosetta 2(DE3) (Sigma) and a single colony was inoculated into 5\u2009ml of LB containing carbapenem to grow overnight at 37\u2009\u00b0C and 200-rpm shaking. The overnight culture was then added to 1\u2009L of LB\u2009+\u2009carbapenem prewarmed to 37\u2009\u00b0C and grown at 37\u2009\u00b0C and 200-rpm shaking until an OD600 of 0.2\u20130.4. Once the desired OD range was achieved, the culture was transferred to 16\u2009\u00b0C and 200-rpm shaking for 30\u2009min, after which IPTG was added to 1\u2009mM to induce protein expression, which continued overnight at 16\u2009\u00b0C and 200\u2009rpm. The culture was pelleted by centrifugation at 4,000g for 45\u2009min and the cell pellet was resuspended in 10\u2009ml of protein lysis buffer (50\u2009mM Tris pH 7.5, 100\u2009mM NaCl, 5\u2009mM imidazole, 0.1\u2009mM EDTA and 1\u2009mM \u03b2-mercaptoethanol). The suspension was sonicated using a Fisherbrand sonicator for a total of 7.5\u2009min of sonication time, with 45\u2009s on and 45\u2009s off at 45% amplitude. The sonicated suspension was then centrifuged at 4,000g for 30\u2009min. The supernatant was then centrifuged at 23,000g for an additional 30\u2009min. The His-tagged ClpX was isolated from the cleared supernatant using Ni-NTA agarose (Quiagen) according to the manufacturer\u2019s instructions. Briefly, a column containing 2\u2009ml of Ni-NTA was washed with 20\u2009ml of protein lysis buffer, followed by passage of the cleared lysate. Bound ClpX was then washed with 20\u2009ml of increasing concentrations of imidazole in protein lysis buffer (10\u2009mM, 20\u2009mM and 30\u2009mM), followed by elution with protein lysis buffer containing 250\u2009mM imidazole, while collecting 0.5-ml fractions. The protein content in the eluted fractions was checked by absorbance at 280\u2009nm, with any fractions at an OD\u2009&gt;\u20090.6 pooled together. About 1\u2009ml of the pooled fractions were then dialyzed (MWCO: 12,000\u201314,000\u2009kDa; Spectra\/Por) against 1\u2009L of protein storage buffer (50\u2009mM Tris pH 7.5, 200\u2009mM KCl, 25\u2009mM MgCl2, 1\u2009mM \u03b2-mercaptoethanol and 10% glycerol) for 3\u2009h at 4\u2009\u00b0C, followed by an additional dialysis against 3\u2009L of protein storage buffer overnight. The protein concentration was then determined to be 5.83\u2009mg\u2009ml\u22121 (117\u2009\u00b5M) by spectroscopy at 280\u2009nm with the calculation of 1 OD corresponding to 3.69\u2009mg\u2009ml\u22121 of His-ClpX. A sample (50\u2009ml) of the 3-L storage buffer after dialysis was saved for resuspending trigintamicin and its isomer for ITC.<\/p>\n<p>ITC<\/p>\n<p>ITC measurements were carried out at 25\u2009\u00b0C using a Malvern ITC PEAQ system. The ligands (trigintamicin and its biologically inactive isoform) were dissolved in storage buffer taken post dialysis of ClpX. The concentration of ClpX was adjusted to ~30\u2009\u03bcM and the ligand concentrations were adjusted to ~300\u2009\u03bcM using the postdialysis storage buffer batch. Then, 300\u2009\u03bcl of ClpX solution was placed in the sample cell and titrated by ligand solution through 19 successive injections of 2\u2009\u03bcl with a 2.5-min interval between injections and a reference power of 5\u2009\u03bccal\u2009s\u22121. Titration data were analyzed using the PEAQ-ITC analysis software (version 1.41). The isotherm curves were best fitted to a single-site binding model.<\/p>\n<p>Depolarization assay<\/p>\n<p>An overnight culture of S.\u2009aureus SH1000 was pelleted, washed twice in equal volume of DPBS (Thermo Fisher Scientific) and diluted to 0.35 OD600 in DPBS. Then, 100\u2009\u00b5l of the dilution was added to 300\u2009\u00b5l of DPBS, followed by the addition of 50\u2009\u00b5l of 20\u2009\u00b5M DISC3(5) dye in DPBS. The mixture was incubated at room temperature in the dark for 15\u2009min, followed by the addition of 50\u2009\u00b5l of 2\u2009M KCl and another 15-min incubation. The samples were then transferred into a 384-well flat, clear-bottom black microtiter plate with 30\u2009\u00b5l per well. The initial fluorescence intensity of each well was recorded using a TECAN Infinite M NANO+ (excitation, 620\u2009nm; emission, 675\u2009nm) at 10-s intervals until the baseline stabilized. After signal stabilization, 30\u2009\u00b5l of test compounds diluted in DPBS were added to the indicated concentrations and readings resumed. Gramicidin at 16\u2009\u00b5g\u2009ml\u22121 (2\u00d7 MIC) was used as a positive control for depolarization.<\/p>\n<p>Lipid antagonism assay<\/p>\n<p>Test lipids were dissolved in the manufacturer\u2019s recommended solvent to 0.5\u20131\u2009mg\u2009ml\u22121. An aliquot containing the appropriate amount of lipid for the desired molar ratio to a target antibiotic was transferred to a 1.5-ml Eppendorf and dried under vacuum for 1\u2009h. The dried pellet was then resuspended by pipetting with 6\u2009\u00b5l of the test antibiotic dissolved in ethanol at 1.28\u2009\u00b5g\u2009ml\u22121. The suspension was then topped up with 114\u2009\u00b5l of LB and the sample was mixed by pipetting. Then, 100\u2009\u00b5l was transferred to a 96-well assay plate and 50\u2009\u00b5l was transferred down a column to generate a 1:2 dilution series. An overnight culture of S.\u2009aureus SH1000 was diluted 1:1,000 into LB and 50\u2009\u00b5l was added to the dilution series to generate a final antibiotic concentration range of 0.25 to 32\u2009\u00b5g\u2009ml\u22121. The plate was then incubated overnight at 30\u2009\u00b0C for 18\u2009h, after which growth was visually inspected. For experiments with telomycin and daptomycin, all LB was supplemented with Ca2+ to 50\u2009mg\u2009L\u22121 (183\u2009mg\u2009L\u22121 CaCl2).<\/p>\n<p>Reporting summary<\/p>\n<p>Further information on research design is available in the <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41587-025-02810-w#MOESM2\" rel=\"nofollow noopener\" target=\"_blank\">Nature Portfolio Reporting Summary<\/a> linked to this article.<\/p>\n","protected":false},"excerpt":{"rendered":"Bacterial strains and growth conditions All strains, unless otherwise specified, were grown in Luria\u2013Bertani (LB) at 30\u2009\u00b0C with&hellip;\n","protected":false},"author":2,"featured_media":134836,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[7],"tags":[5083,4414,10401,18924,3251,3252,59,3250,6194,90,63080,56,54,55],"class_list":{"0":"post-134835","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","7":"category-science","8":"tag-agriculture","9":"tag-antibiotics","10":"tag-bioinformatics","11":"tag-biomedical-engineering-biotechnology","12":"tag-biomedicine","13":"tag-biotechnology","14":"tag-gb","15":"tag-general","16":"tag-life-sciences","17":"tag-science","18":"tag-soil-microbiology","19":"tag-uk","20":"tag-united-kingdom","21":"tag-unitedkingdom"},"_links":{"self":[{"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/posts\/134835","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/comments?post=134835"}],"version-history":[{"count":0,"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/posts\/134835\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/media\/134836"}],"wp:attachment":[{"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/media?parent=134835"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/categories?post=134835"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.newsbeep.com\/uk\/wp-json\/wp\/v2\/tags?post=134835"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}