Next Article in Journal
Bibliometric Examination of Global Scientific Research about Carbapenem-Resistant Acinetobacter Baumannii (CRAB)
Next Article in Special Issue
AMPs as Host-Directed Immunomodulatory Agents against Skin Infections Caused by Opportunistic Bacterial Pathogens
Previous Article in Journal
Antimicrobial Activities of Aztreonam-Avibactam and Comparator Agents against Enterobacterales Analyzed by ICU and Non-ICU Wards, Infection Sources, and Geographic Regions: ATLAS Program 2016–2020
Previous Article in Special Issue
Antimicrobial Effects of Tetraspanin CD9 Peptide against Microbiota Causing Armpit Malodour
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 12
Issue 11
10.3390/antibiotics12111592
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome Mining Uncovers NRPS and PKS Clusters in Rothia dentocariosa with Inhibitory Activity against Neisseria Species

by
Elvis Achondou Akomoneh
1,2,
Zina Gestels
1,
Saïd Abdellati
3,
Katleen Vereecken
4,
Koen Bartholomeeusen
4,
Dorien Van den Bossche
3,
Chris Kenyon
1,4,5,* and
Sheeba Santhini Manoharan-Basil
1
1
HIV/STI Unit, Department of Clinical Sciences, Institute of Tropical Medicine Antwerp, 2000 Antwerp, Belgium
2
Department of Microbiology and Parasitology, University of Bamenda, Bambili P.O. Box 39, Cameroon
3
Clinical Reference Laboratory, Department of Clinical Sciences, Institute of Tropical Medicine, 2000 Antwerp, Belgium
4
Virology Unit, Department of Biomedical Sciences, Institute of Tropical Medicine, 2000 Antwerp, Belgium
5
Department of Medicine, University of Cape Town, Cape Town 7700, South Africa
*
Author to whom correspondence should be addressed.
Antibiotics 2023, 12(11), 1592; https://doi.org/10.3390/antibiotics12111592
Submission received: 12 October 2023 / Revised: 31 October 2023 / Accepted: 2 November 2023 / Published: 4 November 2023
(This article belongs to the Special Issue Antimicrobial Peptides: A New Hope in the War against Multi-Drug Resistance)
Browse Figures
Review Reports Versions Notes

Abstract

:
The growing global threat of antimicrobial resistance is reaching a crisis point as common bacterial infections, including those caused by pathogenic Neisseria species, are becoming increasingly untreatable. This is compelling the scientific community to search for new antimicrobial agents, taking advantage of computational mining and using whole genome sequences to discover natural products from the human microbiome with antibiotic effects. In this study, we investigated the crude extract from a Rothia dentocariosa strain with demonstrated antimicrobial activity against pathogenic Neisseria spp. by spot-on-lawn assay. The genomic DNA of the R. dentocariosa strain was sequenced, and bioinformatic evaluation was performed using antiSMASH and PRISM to search for biosynthetic gene clusters (BGCs). The crude extract with potential antimicrobial activity was run on Tricine-SDS-PAGE, and the putative peptides were characterised using liquid chromatography–tandem mass spectrometry (LC-MS). The crude extract inhibited the growth of the pathogenic Neisseria spp. Six BGCs were identified corresponding to non-ribosomal peptide synthases (NRPSs), polyketide synthases (PKSs), and ribosomally synthesised and post-translationally modified peptides. Three peptides were also identified corresponding to Actinorhodin polyketide putative beta-ketoacyl synthase 1. These findings serve as a useful reference to facilitate the research and development of NRPS and PKS as antimicrobial products against multidrug-resistant N. gonorrhoeae.

1. Introduction

The human oral microbiome harbours more than 700 bacteria species within different ecological niches and is considered as one of the most diverse in the human body [1]. As it is consistently in contact with the exterior environment, the oral cavity serves as one of the main portals of entry for pathogenic microorganisms, and therefore, a healthy oral microbiome could be an indispensable portion of the first line of defense against invading pathogens [1]. The ability of the human indigenous microbiota to resist colonisation and invasion by pathogenic bacteria is termed colonisation resistance [2]. Different mechanisms have been postulated to explain colonisation resistance, including the production of antagonistic compounds (e.g., bacteriocins and bacteriolytic enzymes) [3], competition for nutrients and adhesion sites [4], and the stimulation of the host immune response [5]. These activities create a microenvironment that prevents colonisation by and the dissemination of pathogens [6].
Rothia species are among the early colonisers of the oropharyngeal cavity, including the tonsils, tongue, and teeth of healthy humans [7,8] and are important members of the normal human oral microbiome [9,10,11]. Rothia dentocariosa is capable of fermenting a variety of carbohydrates and amino acids and is rarely implicated in human disease [12]. There is evidence that some Rothia species encode novel bioactive compounds in their biosynthetic gene clusters (BGCs) that could be used as a potential source of new antimicrobials [11]. Also, Rothia species have been identified with BGCs for archetype siderophore enterobactin, the strongest iron (III)-binding siderophore [12]. Biosynthetic gene clusters (BGCs) produce three major groups of metabolites: ribosomally synthesised and post-translationally modified peptides (RiPPs), non-ribosomal peptides (NRPs), and polyketides (PKS) [11,13]. RiPPs form a large heterogeneous group of peptides that are usually classified into post-translationally modified peptides and unmodified peptides [14]. NRPS-peptides are relatively small peptides containing 2–30 amino acids with hallmarks of non-proteinogenic amino acids, and their biosynthesis is ribosome independent [15,16]. Polyketides, on the other hand, are synthesised by polyketide synthases through repetitive condensations of acyl CoA esters to ward off competing microbes in response to nutrient limitation [17].
NRPS and PKS are the major enzymes of secondary metabolite synthesis that catalyse the synthesis of oligopeptides and the elongation of polyketides, respectively [18,19]. They can be extensively post-translationally modified, and these modifications lead to products with many features resembling the NRPS peptides [20]. Polyketides are classified into types I, II, and III [21] based on their biochemical mechanisms and enzyme architecture [22]. NRPS and PKS are responsible for the synthesis of a wide array of antimicrobial compounds, siderophores, and toxins, and they are valued clinically as antimicrobial, antifungal, antiparasitic, antitumor, and immunosuppressive agents [23]. Products and corresponding extracts from organisms with NRPS and PKS gene clusters have demonstrated interesting antimicrobial activities in vivo [6,11,13,24].
Previously, we isolated and identified Rothia dentocariosa from oropharyngeal swabs of men who have sex with men (MSM) with inhibitory activities against Neisseria gonorrhoeae and N. meningitidis, the two human pathogenic species of the genus Neisseria [25]. Considering the continuing emergence of antimicrobial resistance reported from both these pathogens, alternative treatments are urgently required [26,27]. In this study, we sought to expand our previous findings by sequencing the whole genome of the R. dentocariosa that we isolated and performed biosynthetic gene cluster (BGC) mining. The putative peptides were characterised using the nanoACQUITY LC system. We identified six BGCs that could possibly explain the inhibitory mechanisms employed by R. dentocariosa against the pathogenic Neisseria spp.

2. Methods

2.1. Samples Used in This Study and the Isolation and Identification of Rothia dentocariosa

Previously, we analysed oropharyngeal samples collected from men who have sex with men available from the Preventing Resistance to Gonorrhoeae (PReGo) study (ClinicalTrials.gov, NCT03881007) [28]. Briefly, bacterial isolates from 118 oropharyngeal swabs were cultured on BDTM Columbia Agar with 5% sheep blood (BA) plates (Becton Dickinson, Heidelberg, Germany). Single colonies were obtained by plating a loopful of culture that was diluted (10−6) on BDTM Columbia nalidixic acid agar (CNA) plates with 5% sheep blood (Becton Dickinson). This was followed by replica plating [29] of the colonies from each plate to gonococcal (GC) medium base (BD Difco™ Dehydrated Culture Media) agar plates (n = 3) that contained a lawn of N. gonorrhoeae (WHO P, reference isolate) or N. meningitidis (M00003, clinical isolate from the ITM collection). The replica plates were incubated overnight at 37 °C with 5% CO2. Inhibitory isolates were sub-cultured from the third set of replica plates, identified using matrix-assisted laser desorption/ionisation-time of flight (MALDI-TOF) Biotyper IVD (Bruker Daltonics, Bremen, Germany) (library updated to v.10.0.0.0_8326-9468 (IVD). Furthermore, growth inhibition of the inhibitory isolates was confirmed by agar overlay assay using the N. gonorrhoeae (WHO P) and N. meningitidis (M00003/1) isolates. This process led to the identification of a strain of Rothia dentocariosa that was inhibitory to the pathogenic Neisseria spp. [25].

2.2. Production of Crude Extracts from the Rothia dentocariosa Isolate

The protocol was adapted from Cheng et al. [30]. R. dentocariosa was grown in GC broth overnight for the production of crude extract. Briefly, 500 mL of GC semi-solid medium (GC broth supplemented with 6 g/L of agar) was prepared and allowed to cool for 24 h. This was stabbed with a sterile glass Pasteur pipet about 100 times, and 4 mL of the overnight culture was poured on the semisolid agar surface. The mixture was incubated at 37 °C with 5% CO2 in a CO2 incubator for 72 h. Following incubation, the semisolid agar was mashed with a spatula, dispensed into 50 mL falcon tubes, and stored overnight at −80 °C. The frozen mixture was then thawed for an hour in a water bath at 65 °C. The freeze/thaw cycle (1X) was used to disrupt the cell membranes of the Rothia bacteria. The thawed mixture was centrifuged at 4000× g for 140 min at 4 °C, and the supernatant was carefully collected and mixed with equal volumes of chloroform, vigorously shaken by hand for 20 min, and stored overnight at 4 °C. The upper layer was then removed, and the white interfacial layer was transferred into a beaker and kept in a fume hood (for approximately 3 days) for any residual chloroform to evaporate. The brownish extract-containing residue was dissolved in 7 mL of 35% acetonitrile-0.1% trifluoroacetic acid and centrifuged at 3200× g for 40 min at 4 °C to remove insoluble materials. The clear supernatant (henceforth referred to as crude extract) was passed through a 0.22 µm syringe filter, concentrated using a SpeedVac vacuum concentrator at room temperature, aliquoted, and stored at −20 °C for further analyses.

2.3. Detection of the Antimicrobial Activity of the Crude Extract

The antimicrobial activity of the crude extract was detected using the spot-on-lawn agar overlay assay protocol [25]. Briefly, 0.5 MacFarland of 24 h culture of target strains (N. gonorrhoeae or N. meningitidis) grown in GC broth was diluted in 4 mL of GC semi-solid medium pre-warmed to 55 °C to obtain a concentration of approximately 106 CFU/mL. This was overlayed onto GC agar plates and allowed to solidify. Twenty microlitre aliquots of the crude extract were spotted on top of the agar overlay and allowed for complete adsorption. The plates were incubated at 37 °C with 5% CO2 in a CO2 incubator for 18 h and observed for the zone of inhibition. Twenty microlitres of 35% acetonitrile-0.1% trifluoroacetic was also spotted on an agar overlay as a negative control.

2.4. Separation of the Crude Extract by Tricine-SDS-PAGE Electrophoresis

The crude extract was run on sodium dodecyl sulphate–polyacrylamide gel electrophoresis (SDS-PAGE) to fractionate putative proteins [31]. Fifty microlitres of crude extract was mixed with 15 µL of sample loading dye and 3 µL of SDS-sample loading buffer (Sigma-Aldrich, Schnelldorf, Germany) and loaded into 4–20% Bis-Tris precast polyacrylamide gel (Bio-Rad, Hercules, CA, USA) alongside a low-range rainbow molecular weight marker (Spectra low molecular range marker, Thermo Scientific, Waltham, MA, USA) of sizes ranging from 3.4 kDa to 250 kDa. The setup was run at 140 V for 50 min in a mini vertical electrophoresis apparatus (Bio-Rad), after which the gel was stained with Coomassie brilliant blue G-250 (Thermo Scientific, Waltham, MA, USA). Following staining, the gel was washed (3×) with distilled water at room temperature to destain and then observed for visible protein bands. These bands were then sequenced for further analysis.

2.5. Whole-Genome Sequencing and Mining for Biosynthetic Gene Clusters (BGCs) in Rothia dentocariosa

The sequencing of Rothia dentocariosa was outsourced to Eurofins Genomics (Konstanz, Germany), where total DNA was isolated. Library preparation was carried out using a Stranded TruSeq DNA library preparation kit. Sequencing was performed on NextSeq6000, v2, 2 × 150 bp (Illumina Inc., San Diego, CA, USA), followed by analysis of the raw reads. For the WGS analysis, initial quality control (QC) of the raw reads was carried out using FASTQC [32]. The raw reads were trimmed using trimmomatic (v0.39) (Phred score ≥ 20 and length of the bases ≥32 bases) [33]. The processed raw reads were de novo assembled using Shovill (v1.0.4) [34], which uses SPAdes (v3.14.0) [35] using the following parameters: trim—depth 150—opts—isolate. The quality of the de novo assembled contigs was evaluated using Quast (v5.0.2) [36] and annotated using Prokka (v1.14.6) [37].
Further computational analysis was carried out as follows: briefly, the contigs were reordered using the Rothia dentocariosa ATCC 17931 genome sequence (GenBank accession CP002280.1), followed by joining the contigs using the ‘join sequences’ option implemented in CLC GenomicsWorkbench V20. The genome was run through the biosynthetic gene cluster (BGC) prediction program antiSMASH v 6.1.1 [38] (antismash --enable-genefunctions --enable-lanthipeptides --enable-lassopeptides --enable-nrps-pks --enable-sactipeptides --enable-t2pks --enable-thiopeptides --enable-tta --enable-html --fullhmmer --genefinding-tool prodigal -c 48 Rothia.gbk) and PRISM (https://prism.adapsyn.com/ (accessed on 4 July 2023)). The simplified molecular input line entry system (SMILES) and structures for the BGC products were predicted using PRISM [39]. The circular genome diagram was generated using CLC GenomicsWorkbench V20.
The raw reads generated were deposited at BioProject ID:PRJNA1013005.

2.6. Characterisation of Putative Peptides by Liquid Chromatography–Tandem Mass Spectrometry (LC-MS)

Bands on the SDS-PAGE gel were excised using a scalpel and collected on low-binding Eppendorf tubes (ThermoFisher scientific). The excised gel plugs were subjected to in-gel digestion following published protocols [40]. Briefly, the bands were washed with 100 µL water and shaken for 5 min. Water was then removed and replaced by 95% ACN to shrink the bands. The procedure was repeated several times to remove most of the staining solution. The bands were then swollen in 6.66 mM dithiothreitol dissolved in 50 mM ammonium bicarbonate (ABC, pH 7.8) and incubated at 56 °C for 45 min. They were then cooled to room temperature and washed with 95% ACN. The plugs were then swollen in 55 mM iodoacetamide, incubated in the dark for 30 min, and then washed with 95% ACN. Approximately 20 µL digestion buffer composed of 1 µg/µL trypsin suspended in 50 mM ABC was added and the tubes were incubated at 37 °C overnight. Digestion was stopped the following day by freezing the tubes. The extracts were then collected before and after the gel bands were shrunk in 95% ACN/H2O. These were dried under vacuum and resuspended in 20 µL 1% ACN/H2O. Approximately 2 µL was used to calculate the total peptide concentration using a NanoDrop® ND-1000 UV-Vis Spectrophotometer (Thermo Scientific). Next, 0.5 µg of each sample (where possible; samples with lower peptide yields were all used) was loaded on a micropillar array (mPAC™) trapping column and injected on a 200 cm C18 mPAC™ column (Pharmafluidics, Zwijnaarde, Belgium) connected to a nanoACQUITY LC system (Waters). Separation was performed in reverse phase using a linear gradient of mobile phase B (0.1% formic acid in 98% acetonitrile) from 1% to 40% in 80 min, followed by a steep increase to 100% mobile phase B in 5 min. After 5 min at 100% mobile phase B, a steep decrease to 1% mobile phase B was achieved in 5 min, and 1% mobile phase B was maintained for 35 min. The flow rate was 750 nL per minute. The LC system was coupled to a Q-Exactive Plus orbitrap mass spectrometer (ThermoFisher Scientific) programmed to acquire in data-dependent mode. The survey scans were acquired in the orbitrap mass analyser operating at 70,000 (FWHM) resolving power at the mass range of 350–1850 m/z, with a target of 3E6 ions and 100 ms injection time. Precursors were selected “on the fly” for high-energy collision-induced dissociation (HCD) fragmentation with an isolation window of 1.6 amu and a normalised collision energy of 28%. A target of 1.7E3 ions and a maximum injection time of 80 ms were used for MS/MS. The method was set to analyse the top 20 most intense ions from the survey scan, and dynamic exclusion was enabled for 20 s.
Tandem mass spectra were processed using PEAKS Online X build 1.7.2022-08-03_160501 (Bioninformatics Solution, Inc., Waterloo, ON, Canada). Database searches were performed with the precursor and fragment tolerance set to 20 ppm and 0.05 Da, respectively, against an in-house database constructed from the protein translations of the whole-genome sequence of R. dentocariosa, as well as antiSMASH [38] identifications derived from the WGS. Trypsin was specified as the cleavage enzyme, allowing for a maximum of 2 miscleavages. Carbamidomethylation (C), and oxidation (M) and deamidation (NQ) were set as fixed and variable modifications, respectively. Peptide-to-spectrum matches (PSMs) were filtered at 1% FDR, with at least one peptide assigned per protein. Peptide de novo search was also performed, and the results were filtered at an average local confidence (ALC) ≥ 50%. PEAKS PTM and SPIDER, integrated with the PEAKS software, were also used to search for potential post-translational modifications and possible mutations using default search parameters.

3. Results

An overview of the study is provided in Figure 1.

3.1. Antibacterial Activity of the Crude Extract and Tricine-SDS-PAGE Analysis of Antibacterial Compounds

The crude extract and the concentrate that was obtained from R. dentocariosa showed inhibition against both the N. gonorrhoeae and N. meningitidis strains used in this study via the spot-on-lawn assay (Figure 2).
A prominent band with an apparent mass of around 80 kDa was identified for the crude extract (Figure 3).

3.2. The Rothia dentocariosa Genome Possesses Multiple Secondary Metabolite Biosynthetic Gene Clusters

The genomic DNA of Rothia dentocariosa was isolated and sequenced. We performed a bioinformatic evaluation of the R. dentocariosa genome using antiSMASH v.6.1.1 and PRISM to search for BGCs such as polyketide synthases (PKSs), non-ribosomal peptide synthases (NRPSs), and ribosomally synthesised and post-translationally modified peptides (RiPPs). Six biosynthetic gene clusters were identified and characterised (Figure 4 and Table 1). The BGCs corresponded to PKSs, NRPSs, and RiPPs and are as follows: one NRPS BGC was predicted to be related to the biosynthesis of enterobactin, EntF (Figure S1a) and two genes belonging to one cluster were identified as type II PKS hybrid BGCs and related to the biosynthesis of Actinorhodin polyketide putative beta-ketoacyl synthase 1 and 2 (Figure S1b). Additionally, two RiPPs, a lantipeptide-class-ii peptide (Figure S1c) and a LAP (Linear azol(in)e-containing peptides) (Figure S1d), along with Butyrolactone (Figure S1e) were identified. Notably, out of the six BGCs, no homologous gene clusters could be identified for one cluster and hence represent an orphan BGCs (Figure 4).

3.3. Identification of Peptides

The crude extract with potential antimicrobial activity was run on Tricine-SDS-PAGE. The fractionated proteins were excised from the gel and the peptides were identified and characterised using the nanoACQUITY LC system for the presence of putative bacteriocins (peptides/proteins). As shown in Supplementary Figure S2, three peptides, Ser-Gly-Phe-Ser-Glu-Glu-Glu-Ile-Ser-Thr-Leu-Asp-Lys, Glu-Ala-Glu-Met-Tyr-Gly-Ser-Ile-Thr-Gly-Tyr-Gly-Ala-Arg, and Ala-Tyr-Asp-Val-Ser-Ser-Leu-Lys (Table 1) were identified. The three peptides corresponded to Actinorhodin polyketide putative beta-ketoacyl synthase 1, identified by antiSMASH 6.0. No peptides were identified from any of the other five BGCs.

4. Discussion

As the search for novel approaches to contain multidrug-resistant bacteria continues, exploration of the mechanisms of colonisation resistance is increasingly gaining attention [41]. Previously, we identified a Rothia dentocariosa strain from oropharyngeal samples from men who have sex with men with inhibitory activity against pathogenic Neisseria spp. [25]. The crude extract derived from R. dentocariosa showed inhibitory activity against two pathogens, Neisseria gonorrhoeae and Neisseria meningitidis, as demonstrated in the spot-on-lawn assay (Figure 2).
In this study, we further characterised R. dentocariosa by sequencing the isolate and conducting in silico analyses of BGCs. More than 5000 BGCs have been identified in oral bacteria including Rothia species [11,42]. We identified six BGCs in this study corresponding to PKSs, NRPSs, and RiPPs whose products could be responsible for the inhibition of the pathogenic Neisseria spp. by R. dentocariosa. An in silico study by Isabela et al. 2022 [11], identified BGCs from 155 Rothia genomes predicted to produce antibiotic non-ribosomal peptides and other secondary metabolites. Similarly, using antiSMASH software, Uranga et al. 2020 [12] predicted that R. dentocariosa draft genomes (n = 12) harbored a butyrolactone BGC, a type 1 polyketide synthase (T1PKS) BGC, a lanthipeptide BGC, and a catechol siderophore-like (cat-sid) BGC. The same study identified an archetype siderophore enterobactin characterised from R. mucilaginosa ATCC 25296 that was shown to be involved in the reduced growth of some strains of S. mutans, oral Streptococci species, and oral Actinomyces timomensis [12]. In the present study, we identified enterobactin, EntF BGC (Figure S1a), type II PKS hybrid BGCs, Actinorhodin polyketide putative beta-ketoacyl synthase 1 and 2 BGCs, two RiPPs, a lantipeptide-class-ii peptide (Figure S1c), and a LAP (Linear azol(in)e-containing peptides) (Figure S1d) BGCs, along with a Butyrolactone BGC (Figure S1e). These peptides could explain some of Rothia’s competitive success in the oral cavity [11,12]. These results are therefore of great significance as PKS and NRPS gene clusters identified elsewhere are responsible for the synthesis of exciting novel antimicrobial compounds with a low risk of antimicrobial resistance, like valinomycin (from R. nasisuis) [11] and lugdinin (from Staphylococcus lugdunensis) [43].
Using the SDS-PAGE gel, a distinct band at approximately 80 kDa was identified that likely represents a bioactive molecule. By combining genome mining and liquid chromatography–tandem mass spectrometry (LC MS) analysis, we identified three peptides corresponding to Actinorhodin polyketide putative beta-ketoacyl synthase 1 and 2 in this band (Figure S2) which catalyse the biosynthesis of natural antimicrobials such as actinorhodin [44]. While these peptides are intriguing, they may represent only a fraction of the bioactive compounds present in the crude extract.
The existence of an orphan BGC within Rothia dentocariosa’s genome presents an opportunity for novel drug discovery. Similar orphan BGCs have been a focus of investigation in other studies [45]. Li et al. (2018) identified orphan BGCs in the Streptomyces atratus ZH16 genome, leading to the production of nocardamine-type (desferrioxamine) compounds via metabolic engineering, highlighting the importance of exploring uncharacterised gene clusters [46].
The caveats of the study include the following: we were unable to synthesise the peptides and prove that they had activity against the pathogenic Neisseria spp. We did however contact three companies/researchers linked to the compound in manuscripts and none could manufacture it for us. Also, we did not screen for antimicrobial activity against other bacteria genera in vitro and the bioactivity of the crude extract in vivo. Likewise, we did not characterise the secreted proteomes of this bacteria.
Nevertheless, this study provides novel evidence of the antibacterial activity of Rothia dentocariosa’s crude extract against Neisseria pathogens, offering a potential source of new antimicrobial agents. The genomic analysis sheds light on the biosynthetic potential of this commensal bacterium, and the discovery of an orphan BGC underscores the importance of further exploration. The identification of specific peptides raises exciting prospects for understanding the mechanisms of antimicrobial action within Rothia dentocariosa and highlights its potential contribution to combating antibiotic-resistant pathogens. Our findings highlight the potential for producing novel antimicrobial compounds from Rothia dentocariosa. Furthermore, we suspect that screening the nasal microbiome for bacteria that produce secondary metabolites will yield more peptides with anti-gonococcal and anti-meningococcal activity that could be developed as therapeutics.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antibiotics12111592/s1. Figure S1: The putative genetic organisation of (a) enterobactin, EntF, (b) actinorhodin polyketide putative beta-ketoacyl synthase 1 and 2, (c) lantipeptide-class-ii peptide, (d) LAP (Linear azol(in)e-containing peptides), and (e) Butyrolactone. Figure S2: Putative genetic organisation of (a) enterobactin, EntF and (b) its corresponding three peptides identified using LC-MS system.

Author Contributions

Conceptualisation, E.A.A., C.K., and S.S.M.-B.; methodology, E.A.A., Z.G., S.A., K.V., K.B., C.K., and S.S.M.-B.; validation, C.K.; formal analysis, E.A.A., Z.G., S.A., K.V., K.B., D.V.d.B., C.K., and S.S.M.-B.; software, E.A.A., C.K., and S.S.M.-B.; investigation, E.A.A., Z.G., S.A., K.V., K.B., D.V.d.B., C.K., and S.S.M.-B.; resources, C.K.; data curation, E.A.A., C.K., and S.S.M.-B.; writing—original draft, E.A.A., C.K., and S.S.M.-B.; writing—review and editing, E.A.A., Z.G., S.A., K.V., K.B., D.V.d.B., C.K., and S.S.M.-B.; visualisation, C.K.; supervision, C.K.; project administration, C.K.; funding acquisition, C.K. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by a SOFI 2021 grant—“PReventing the Emergence of untreatable STIs via radical Prevention” (PRESTIP).

Institutional Review Board Statement

This study being a continuation of the PReGo study was also covered by the ethical approvals from the Institutional Review Board of the Institute of Tropical Medicine, Antwerp (Ref: 1276/18) and the Ethics Committee of the University of Antwerp (Ref: 19/06/058).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the findings of this study are retained at ITM and because of ethical and privacy concerns will not be made openly accessible. ITM adheres to the FAIR data principles (findable, accessible, interoperable, and reusable) and recognises that data should be “as open as possible and as closed as necessary”. Anonymised, individual participant data from the study as well as additional related documents, such as the study protocol, the annotated case report form, the data dictionary, and statistical analysis scripts, can be made available within 12 months of the publication of the study results and without an end date. Data will be retained at the ITM data repository and can be requested via email to ITM’s central point for research data access at [email protected]. A governed data access mechanism applies including: (1) the completion of a data request form, (2) evaluation by a data access committee, (3) the signing of a data sharing agreement, and (4) the secure transfer of data.

Acknowledgments

Center for Proteomics—University of Antwerp for the proteomic analyses.

Conflicts of Interest

All authors declare no competing interests.

References

  1. Human Microbiome Project Consortium. Structure, Function and Diversity of the Healthy Human Microbiome. Nature 2012, 486, 207–214. [Google Scholar] [CrossRef] [PubMed]
  2. Fehervari, Z. Mechanisms of Colonization Resistance. Available online: https://www.nature.com/articles/d42859-019-00018-y#author-0. (accessed on 25 May 2023).
  3. Vonaesch, P.; Anderson, M.; Sansonetti, P.J. Pathogens, Microbiome and the Host: Emergence of the Ecological Koch’s Postulates. FEMS Microbiol. Rev. 2018, 42, 273–292. [Google Scholar] [CrossRef] [PubMed]
  4. Vandenbergh, P.A. Lactic Acid Bacteria, Their Metabolic Products and Interference with Microbial Growth. FEMS Microbiol. Rev. 1993, 12, 221–237. [Google Scholar] [CrossRef]
  5. Corthésy, B.; Gaskins, H.R.; Mercenier, A. Cross-Talk between Probiotic Bacteria and the Host Immune System. J. Nutr. 2007, 137, 781S–790S. [Google Scholar] [CrossRef]
  6. Bernet-Camard, M.F.; Liévin, V.; Brassart, D.; Neeser, J.R.; Servin, A.L.; Hudault, S. The Human Lactobacillus Acidophilus Strain LA1 Secretes a Nonbacteriocin Antibacterial Substance(s) Active in Vitro and in Vivo. Appl. Environ. Microbiol. 1997, 63, 2747–2753. [Google Scholar] [CrossRef] [PubMed]
  7. Sulyanto, R.M.; Thompson, Z.A.; Beall, C.J.; Leys, E.J.; Griffen, A.L. The Predominant Oral Microbiota Is Acquired Early in an Organized Pattern. Sci. Rep. 2019, 9, 10550. [Google Scholar] [CrossRef]
  8. Baker, J.L.; Morton, J.T.; Dinis, M.; Alvarez, R.; Tran, N.C.; Knight, R.; Edlund, A. Deep Metagenomics Examines the Oral Microbiome during Dental Caries, Revealing Novel Taxa and Co-Occurrences with Host Molecules. Genome Res. 2021, 31, 64–74. [Google Scholar] [CrossRef]
  9. Wilbert, S.A.; Mark Welch, J.L.; Borisy, G.G. Spatial Ecology of the Human Tongue Dorsum Microbiome. Cell Rep. 2020, 30, 4003–4015.e3. [Google Scholar] [CrossRef]
  10. Parello, C.S.L. Microbiomics. In Translational Systems Medicine and Oral Disease; Elsevier: Amsterdam, The Netherlands, 2020; pp. 137–162. [Google Scholar]
  11. De Oliveira, I.M.F.; Ng, D.Y.K.; van Baarlen, P.; Stegger, M.; Andersen, P.S.; Wells, J.M. Comparative Genomics of Rothia Species Reveals Diversity in Novel Biosynthetic Gene Clusters and Ecological Adaptation to Different Eukaryotic Hosts and Host Niches. Microb. Genom. 2022, 8, 000854. [Google Scholar] [CrossRef]
  12. Uranga, C.C.; Arroyo, P.; Duggan, B.M.; Gerwick, W.H.; Edlund, A. Commensal Oral Rothia Mucilaginosa Produces Enterobactin, a Metal-Chelating Siderophore. mSystems 2020, 5, e00161-20. [Google Scholar] [CrossRef]
  13. Wang, L.; Ravichandran, V.; Yin, Y.; Yin, J.; Zhang, Y. Natural Products from Mammalian Gut Microbiota. Trends Biotechnol. 2019, 37, 492–504. [Google Scholar] [CrossRef]
  14. Letzel, A.-C.; Pidot, S.J.; Hertweck, C. Genome Mining for Ribosomally Synthesized and Post-Translationally Modified Peptides (RiPPs) in Anaerobic Bacteria. BMC Genom. 2014, 15, 983. [Google Scholar] [CrossRef] [PubMed]
  15. Ahmad, V.; Khan, M.S.; Jamal, Q.M.S.; Alzohairy, M.A.; Al Karaawi, M.A.; Siddiqui, M.U. Antimicrobial Potential of Bacteriocins: In Therapy, Agriculture and Food Preservation. Int. J. Antimicrob. Agents 2017, 49, 1–11. [Google Scholar] [CrossRef]
  16. Miller, B.R.; Gulick, A.M. Structural Biology of Nonribosomal Peptide Synthetases. In Nonribosomal Peptide and Polyketide Biosynthesis; Evans, B., Ed.; Humana Press: New York, NY, USA, 2016; Volume 1401, pp. 3–29. [Google Scholar]
  17. Kealey, J.T.; Liu, L.; Santi, D.V.; Betlach, M.C.; Barr, P.J. Production of a Polyketide Natural Product in Nonpolyketide-Producing Prokaryotic and Eukaryotic Hosts. Proc. Natl. Acad. Sci. USA 1998, 95, 505–509. [Google Scholar] [CrossRef] [PubMed]
  18. Ichikawa, N.; Sasagawa, M.; Yamamoto, M.; Komaki, H.; Yoshida, Y.; Yamazaki, S.; Fujita, N. DoBISCUIT: A Database of Secondary Metabolite Biosynthetic Gene Clusters. Nucleic Acids Res. 2012, 41, D408–D414. [Google Scholar] [CrossRef] [PubMed]
  19. Mach, B.; Reich, E.; Tatum, E.L. Separation of the Biosynthesis of the Antibiotic Polypeptide Tyrocidine from Protein Biosynthesis. Proc. Natl. Acad. Sci. USA 1963, 50, 175–181. [Google Scholar] [CrossRef]
  20. McIntosh, J.A.; Donia, M.S.; Schmidt, E.W. Ribosomal Peptide Natural Products: Bridging the Ribosomal and Nonribosomal Worlds. Nat. Prod. Rep. 2009, 26, 537. [Google Scholar] [CrossRef]
  21. Gokulan, K.; Khare, S.; Cerniglia, C. METABOLIC PATHWAYS|Production of Secondary Metabolites of Bacteria. In Encyclopedia of Food Microbiology; Elsevier: Amsterdam, The Netherlands, 2014; pp. 561–569. [Google Scholar]
  22. Shen, B. Polyketide Biosynthesis beyond the Type I, II and III Polyketide Synthase Paradigms. Curr. Opin. Chem. Biol. 2003, 7, 285–295. [Google Scholar] [CrossRef]
  23. Newman, D.J.; Cragg, G.M. Natural Products as Sources of New Drugs from 1981 to 2014. J. Nat. Prod. 2016, 79, 629–661. [Google Scholar] [CrossRef]
  24. Fritz, S.; Rajaonison, A.; Chabrol, O.; Raoult, D.; Rolain, J.-M.; Merhej, V. Full-Length Title: NRPPUR Database Search and in Vitro Analysis Identify an NRPS-PKS Biosynthetic Gene Cluster with a Potential Antibiotic Effect. BMC Bioinform. 2018, 19, 463. [Google Scholar] [CrossRef]
  25. Akomoneh, E.A.; Laumen, J.G.E.; Abdellati, S.; Van Dijck, C.; Vanbaelen, T.; Britto, X.B.; Manoharan-Basil, S.S.; Kenyon, C. The Discovery of Oropharyngeal Microbiota with Inhibitory Activity against Pathogenic Neisseria Gonorrhoeae and Neisseria Meningitidis: An In Vitro Study of Clinical Isolates. Microorganisms 2022, 10, 2497. [Google Scholar] [CrossRef] [PubMed]
  26. Eyre, D.W.; Sanderson, N.D.; Lord, E.; Regisford-Reimmer, N.; Chau, K.; Barker, L.; Morgan, M.; Newnham, R.; Golparian, D.; Unemo, M.; et al. Gonorrhoea Treatment Failure Caused by a Neisseria gonorrhoeae Strain with Combined Ceftriaxone and High-Level Azithromycin Resistance, England, February 2018. Eurosurveillance 2018, 23, 1800323. [Google Scholar] [CrossRef] [PubMed]
  27. Unemo, M.; Seifert, H.S.; Hook, E.W.; Hawkes, S.; Ndowa, F.; Dillon, J.-A.R. Gonorrhoea. Nat. Rev. Dis. Primers 2019, 5, 79. [Google Scholar] [CrossRef]
  28. van Dijck, C.; Tsoumanis, A.; Rotsaert, A.; Vuylsteke, B.; van den Bossche, D.; Paeleman, E.; de Baetselier, I.; Brosius, I.; Laumen, J.; Buyze, J.; et al. Antibacterial Mouthwash to Prevent Sexually Transmitted Infections in Men Who Have Sex with Men Taking HIV Pre-Exposure Prophylaxis (PReGo): A Randomised, Placebo-Controlled, Crossover Trial. Lancet Infect. Dis. 2021, 21, 657–667. [Google Scholar] [CrossRef] [PubMed]
  29. Lederberg, J.; Lederberg, E.M. Replica Plating and Indirect Selection of Bacterial Mutants. J. Bacteriol. 1952, 63, 399–406. [Google Scholar] [CrossRef]
  30. Cheng, M.; Gong, S.-G.; Lévesque, C. Rapid Isolation and Purification of Secreted Bacteriocins from Streptococcus Mutans and Other Lactic Acid Bacteria. Bio Protoc. 2020, 10, e3824. [Google Scholar] [CrossRef]
  31. Schägger, H.; von Jagow, G. Tricine-Sodium Dodecyl Sulfate-Polyacrylamide Gel Electrophoresis for the Separation of Proteins in the Range from 1 to 100 KDa. Anal. Biochem. 1987, 166, 368–379. [Google Scholar] [CrossRef]
  32. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 11 July 2023).
  33. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  34. Seemann, T. Shovill. Available online: https://github.com/tseemann/shovill (accessed on 15 July 2023).
  35. Prjibelski, A.; Antipov, D.; Meleshko, D.; Lapidus, A.; Korobeynikov, A. Using SPAdes De Novo Assembler. Curr. Protoc. Bioinform. 2020, 70, e102. [Google Scholar] [CrossRef]
  36. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality Assessment Tool for Genome Assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  37. Seemann, T. Prokka: Rapid Prokaryotic Genome Annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef] [PubMed]
  38. Blin, K.; Shaw, S.; Kloosterman, A.M.; Charlop-Powers, Z.; van Wezel, G.P.; Medema, M.H.; Weber, T. AntiSMASH 6.0: Improving Cluster Detection and Comparison Capabilities. Nucleic Acids Res. 2021, 49, W29–W35. [Google Scholar] [CrossRef] [PubMed]
  39. Skinnider, M.A.; Merwin, N.J.; Johnston, C.W.; Magarvey, N.A. PRISM 3: Expanded Prediction of Natural Product Chemical Structures from Microbial Genomes. Nucleic Acids Res. 2017, 45, W49–W54. [Google Scholar] [CrossRef]
  40. Jafari, M.; Primo, V.; Smejkal, G.B.; Moskovets, E.V.; Kuo, W.P.; Ivanov, A.R. Comparison of In-Gel Protein Separation Techniques Commonly Used for Fractionation in Mass Spectrometry-Based Proteomic Profiling. Electrophoresis 2012, 33, 2516–2526. [Google Scholar] [CrossRef] [PubMed]
  41. Lewis, K. Recover the Lost Art of Drug Discovery. Nature 2012, 485, 439–440. [Google Scholar] [CrossRef] [PubMed]
  42. Aleti, G.; Baker, J.L.; Tang, X.; Alvarez, R.; Dinis, M.; Tran, N.C.; Melnik, A.V.; Zhong, C.; Ernst, M.; Dorrestein, P.C.; et al. Identification of the Bacterial Biosynthetic Gene Clusters of the Oral Microbiome Illuminates the Unexplored Social Language of Bacteria during Health and Disease. mBio 2019, 10, e00321-19. [Google Scholar] [CrossRef]
  43. Zipperer, A.; Konnerth, M.C.; Laux, C.; Berscheid, A.; Janek, D.; Weidenmaier, C.; Burian, M.; Schilling, N.A.; Slavetinsky, C.; Marschal, M.; et al. Human Commensals Producing a Novel Antibiotic Impair Pathogen Colonization. Nature 2016, 535, 511–516. [Google Scholar] [CrossRef]
  44. Malpartida, F.; Hopwood, D.A. Molecular Cloning of the Whole Biosynthetic Pathway of a Streptomyces Antibiotic and Its Expression in a Heterologous Host. Nature 1984, 309, 462–464. [Google Scholar] [CrossRef]
  45. Machado, H.; Tuttle, R.N.; Jensen, P.R. Omics-Based Natural Product Discovery and the Lexicon of Genome Mining. Curr. Opin. Microbiol. 2017, 39, 136–142. [Google Scholar] [CrossRef]
  46. Li, Y.; Zhang, C.; Liu, C.; Ju, J.; Ma, J. Genome Sequencing of Streptomyces Atratus SCSIOZH16 and Activation Production of Nocardamine via Metabolic Engineering. Front. Microbiol. 2018, 9, 1269. [Google Scholar] [CrossRef]
Antibiotics 12 01592 g001
Figure 1. Overview of the study. The figure was generated using BioRender.
Figure 1. Overview of the study. The figure was generated using BioRender.
Antibiotics 12 01592 g001
Antibiotics 12 01592 g002
Figure 2. Spot-on-lawn assay showing the inhibitory effect (marked with *) of crude extract on N. gonorrhoeae growing on GC semi-solid medium.
Figure 2. Spot-on-lawn assay showing the inhibitory effect (marked with *) of crude extract on N. gonorrhoeae growing on GC semi-solid medium.
Antibiotics 12 01592 g002
Antibiotics 12 01592 g003
Figure 3. SDS–PAGE analysis of secretory products from Streptococcus parasanguinis and Rothia dentocariosa. Marker (lane 1), S. parasanguinis (undiluted extract) (lanes 2–4), R. dentocariosa (undiluted extract) (lanes 5–6), S. parasanguinis (diluted 1:1) (lanes 7–9), and R. dentocariosa (diluted 1:1) (lanes 10–11).
Figure 3. SDS–PAGE analysis of secretory products from Streptococcus parasanguinis and Rothia dentocariosa. Marker (lane 1), S. parasanguinis (undiluted extract) (lanes 2–4), R. dentocariosa (undiluted extract) (lanes 5–6), S. parasanguinis (diluted 1:1) (lanes 7–9), and R. dentocariosa (diluted 1:1) (lanes 10–11).
Antibiotics 12 01592 g003
Antibiotics 12 01592 g004
Figure 4. The Rothia dentocariosa (57_M6) genome possesses multiple biosynthetic gene clusters (BGCs). The BGCs were identified using antiSMASH v6.0 and the structures were predicted using PRISM. * No predicted structure available. The grey colored BGC depicts orphan the BGC.
Figure 4. The Rothia dentocariosa (57_M6) genome possesses multiple biosynthetic gene clusters (BGCs). The BGCs were identified using antiSMASH v6.0 and the structures were predicted using PRISM. * No predicted structure available. The grey colored BGC depicts orphan the BGC.
Antibiotics 12 01592 g004
Table 1. Putative secondary metabolites producing biosynthetic clusters of Rothia dentocariosa with homology to the most similar known biosynthetic clusters as predicted by antiSMASH 6.1.1.
Table 1. Putative secondary metabolites producing biosynthetic clusters of Rothia dentocariosa with homology to the most similar known biosynthetic clusters as predicted by antiSMASH 6.1.1.
CategoryTypeContigRegionStart (bp)End (bp)Size (bp)Core Biosynthetic Genes/ProductLocus TagMass
Spectrometry
Identification
(Separated Protein on Tricine SDS-PAGE Gel)
NRPSNRPS112.1381,301425,53344,232Enterobactin non-ribosomal peptide synthetase EntFINOPIIDP_01825Not detected
RiPPlantipeptide-class-ii24.1130,371178,14547,774Type 2 lanthipeptide synthetase LanMINOPIIDP_00698Not detected
YcaO-like family proteinINOPIIDP_00709Not detected
OtherButyrolactone32.1259,764275,62615,862AfsA-relatedINOPIIDP_00236Not detected
PKSHybrid (type II PKS)53.1128,334200,79572,461Actinorhodin polyketide putative beta-ketoacyl synthase 2INOPIIDP_00466Not detected
Actinorhodin polyketide putative beta-ketoacyl synthase 1INOPIIDP_00467Detected
RiPPLAP65.162,13385,66123,528YcaO-like family protein INOPIIDP_00975Not detected
SagB family peptide dehydrogenaseINOPIIDP_00976Not detected
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Akomoneh, E.A.; Gestels, Z.; Abdellati, S.; Vereecken, K.; Bartholomeeusen, K.; Van den Bossche, D.; Kenyon, C.; Manoharan-Basil, S.S. Genome Mining Uncovers NRPS and PKS Clusters in Rothia dentocariosa with Inhibitory Activity against Neisseria Species. Antibiotics 2023, 12, 1592. https://doi.org/10.3390/antibiotics12111592

AMA Style

Akomoneh EA, Gestels Z, Abdellati S, Vereecken K, Bartholomeeusen K, Van den Bossche D, Kenyon C, Manoharan-Basil SS. Genome Mining Uncovers NRPS and PKS Clusters in Rothia dentocariosa with Inhibitory Activity against Neisseria Species. Antibiotics. 2023; 12(11):1592. https://doi.org/10.3390/antibiotics12111592

Chicago/Turabian Style

Akomoneh, Elvis Achondou, Zina Gestels, Saïd Abdellati, Katleen Vereecken, Koen Bartholomeeusen, Dorien Van den Bossche, Chris Kenyon, and Sheeba Santhini Manoharan-Basil. 2023. "Genome Mining Uncovers NRPS and PKS Clusters in Rothia dentocariosa with Inhibitory Activity against Neisseria Species" Antibiotics 12, no. 11: 1592. https://doi.org/10.3390/antibiotics12111592

APA Style

Akomoneh, E. A., Gestels, Z., Abdellati, S., Vereecken, K., Bartholomeeusen, K., Van den Bossche, D., Kenyon, C., & Manoharan-Basil, S. S. (2023). Genome Mining Uncovers NRPS and PKS Clusters in Rothia dentocariosa with Inhibitory Activity against Neisseria Species. Antibiotics, 12(11), 1592. https://doi.org/10.3390/antibiotics12111592

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop