- 1Department of Biology, College of Sciences, Princess Nourah bint Abdulrahman University, Riyadh 11617, Saudi Arabia
- 2Biological Sciences Department, College of Science and Arts, King Abdulaziz University, Rabigh 21911, Saudi Arabia
- 3Department of Biochemistry, College of Science, University of Jeddah, Jeddah, Saudi Arabia
- 4Department of Chemistry, Al-Leith University College, Umm Al Qura University, Makkah, Saudi Arabia
- 5Department of Chemistry, College of Sciences and Arts, King Abdulaziz University, Rabigh, Saudi Arabia
- 6Department of Biochemistry, Faculty of Science—King Abdulaziz University, Jeddah, Saudi Arabia
- 7Department of Biology, Faculty of Applied Science, Umm Al-Qura University, Makkah, Saudi Arabia
- 8Department of Biology, College of Science, University of Jeddah, Jeddah, Saudi Arabia
Moringa oleifera (or the miracle tree) is a wild plant species widely grown for its seed pods and leaves, and is used in traditional herbal medicine. The metagenomic whole genome shotgun sequencing (mWGS) approach was used to characterize antibiotic resistance genes (ARGs) of the rhizobiomes of this wild plant and surrounding bulk soil microbiomes and to figure out the chance and consequences for highly abundant ARGs, e.g., mtrA, golS, soxR, oleC, novA, kdpE, vanRO, parY, and rbpA, to horizontally transfer to human gut pathogens via mobile genetic elements (MGEs). The results indicated that abundance of these ARGs, except for golS, was higher in rhizosphere of M. oleifera than that in bulk soil microbiome with no signs of emerging new soil ARGs in either soil type. The most highly abundant metabolic processes of the most abundant ARGs were previously detected in members of phyla Actinobacteria, Proteobacteria, Acidobacteria, Chloroflexi, and Firmicutes. These processes refer to three resistance mechanisms namely antibiotic efflux pump, antibiotic target alteration and antibiotic target protection. Antibiotic efflux mechanism included resistance-nodulation-cell division (RND), ATP-binding cassette (ABC), and major facilitator superfamily (MFS) antibiotics pumps as well as the two-component regulatory kdpDE system. Antibiotic target alteration included glycopeptide resistance gene cluster (vanRO), aminocoumarin resistance parY, and aminocoumarin self-resistance parY. While, antibiotic target protection mechanism included RbpA bacterial RNA polymerase (rpoB)-binding protein. The study supports the claim of the possible horizontal transfer of these ARGs to human gut and emergence of new multidrug resistant clinical isolates. Thus, careful agricultural practices are required especially for plants used in circles of human nutrition industry or in traditional medicine.
Introduction
Moringa oleifera (or the miracle tree) is a drought-tolerant wild plant species of the family Moringaceae (Olson and Flora of North America Editorial Committee, 2010). This species is native to several habitats in Southern America, Africa and Asia including Saudi Arabia (Al-Eisawi and Al-Ruzayza, 2015). This plant is commonly used in traditional herbal medicine due to its health benefits as it possesses anti-fungal, anti-viral, anti-depressant as well as anti-inflammatory properties (Olson and Flora of North America Editorial Committee, 2010). Seed oil and leaf powder of M. oleifera contain a variety of useful proteins, vitamins (A, B1, B2, B3, B6, and C) and minerals (e.g., Ca, K, Mg, Fe, P, and Zn) (Gopalakrishnan et al., 2016). The plant is also used in treating edema, asthma and diabetes, in reducing high blood pressure, in treating stomach complaints, in fighting against bacterial diseases and cancer, in making bones healthier, and in healing skin wounds (Gopalakrishnan et al., 2016).
Recent advances in next generation sequencing (NGS) technology and subsequent bioinformatics approaches allow the accurate detection of microbial composition and function in the phyllospheric and rhizospheric regions of the plant (Vorholt, 2012; Bai et al., 2015; Vorholt et al., 2017). The rhizospheric region is regarded as a hotspot of effective microbial activity, recovery of new antibiotics and potential transfer of antibiotic resistance genes (ARGs) (Blau et al., 2019). NGS approaches of amplicon (or 16S rRNA) and metagenomic whole genome shotgun sequencing (mWGS) are used to decipher microbial composition and predicted function. Advantages of the latter approach include the ability to study the genomes of archaea, bacteria, fungi, and viruses and reaching accurate insights into taxonomy, biodiversity, phylogenetic relationships, and functionality of microbiomes (Segata et al., 2011). mWGS also allows studying the influences of environmental conditions and host-microbe interactions in shaping microbiome signatures and allows deciphering new ARGs and antibiotics in the phyllospheric and rhizospheric regions of a given plant especially if this plant is edible by human or livestock (Tringe et al., 2005; Raes et al., 2007).
Identification and development of new antibiotics was for years a top research priority among pharmaceutical companies (Dutescu and Hillier, 2021). However, the treatment of pathogenic bacteria is increasingly hindered by the ability of bacteria to generate resistance against antibiotics due to occurrence of gene mutations, or due to a bacterial taxon being transformed with ARGs of other bacterial taxon; a process called horizontal gene transfer (HGT) (Bäckhed et al., 2015). The latter incidence raises a concern that ARGs integrated in mobile genetic elements (MGEs) of rhizospheric or phyllospheric bacteria can contaminate an edible plant, thus, can eventually be transmitted to human gut microbiome. If a certain microbial pathogen exists in the latter environment, then a risk of acquired resistance against clinical antibiotics is raised. Up to date, no enough attention was paid to the serious incidence of horizontal transfer of ARGs existing in plant microbiomes to human gut or skin microbiome (Chen et al., 2019). In addition, the influence of resistome of native or wild plants was almost neglected (Obermeier et al., 2021). The study of rhizobiome of native plants allows the detection of natural antibiotics along with new versatile ARGs due to the existence of highly diversified microbial community in this dynamic region.
In the present study, we have deciphered abundant ARGs and predicted their metabolic processes that exist in rhizobiome of the native wild plant Moringa oleifera in order to detect the risks on human health due to the possible occurrence of HGT if organs of this plant were used by human either in human nutrition or herbal medicine.
Materials and methods
Sample collection and deoxyribonucleic acid extraction
Microbial samples were collected in three replicates from rhizosphere of M. oleifera plant grown naturally in the North Western region of Mecca district of Saudi Arabia near the red sea coast (21°12′17.8″N 39°31′26.4″E) and from surrounding bulk soil (Al-Eisawi and Al-Ruzayza, 2015). The selected three plants for the experiment are single-grown, have similar size and received no rainfall for > 3 months prior sample collection, while bulk soil samples were concurrently taken ≥ 10 m apart from the plants. For soil rhizospheric samples, plant lateral roots were first cut at ∼10–30 cm depth. Then, soil of a maximum of 1 cm apart from the root, not the soil that physically adheres to the root, was collected and immediately put in liquid nitrogen, transported to the lab, and stored at −20°C until use (Hurt et al., 2001).
Whole genome shotgun sequencing and bioinformatics analysis
Deoxyribonucleic acid (DNAs) of different soil samples were extracted using CTAB/SDS method, and purity and integrity were checked by 1% agarose gel electrophoresis. DNA concentration was adjusted to 10 ng/μl using dsDNA Assay kit (Life Technologies, United States). Then, 1 μg of each DNA sample was shipped to Novogene Co. (Singapore) for whole metagenome shotgun sequencing. Physical fractionation of DNAs followed by data pre-processing were done and low quality bases (Q-value ≤ 38) exceeding 40-bp threshold were trimmed, and reads with N nucleotides over 10-bp threshold were removed. Then, effective clean data were sequenced on Illumina HiSeq 2500 platform after library preparation using an Ultra DNA Library Prep kit for Illumina (NEB, United States). Generated data were assembled to recover scaffolds using MEGAHIT (K-mer = 55) and chimeras were removed as described (Karlsson et al., 2012; Mende et al., 2012; Oh et al., 2014). Unassembled reads of different samples were mixed together to generate NOVO_MIX scaffolds exploring the less-abundant genes of the different samples. Generated scaffolds were cut off at “N” to get N-free fragments (or scaftigs) (Mende et al., 2012; Nielsen et al., 2014). Then, clean data were mapped to assembled scaftigs using Soap 2.21, and effective scaftigs (found in > 2 samples) were further utilized. Gene prediction was carried out by MetaGeneMark (Nielsen et al., 2014) and predicted genes were dereplicated using Cluster Database at High Identity with Tolerance (CD-HIT) (Li and Godzik, 2006; Fu et al., 2012). Then, non-redundant gene catalogs (nrGC) were constructed after removing redundancy using a greedy pairwise comparison (Li et al., 2014) and annotation was done using the binning reference-based classification method MEGAN (Huson et al., 2011, 2016).
All unique open reading frames (ORFs) were subjected to Blastp against CARD (Comprehensive Antibiotic Resistance Database)1 (e-value ≤ 1e–5) (Martínez et al., 2015) in order to identify and characterize ARGs. Then, relative abundance of different ARGs in unit ppm was estimated (Yang et al., 2013; Forsberg et al., 2014). ARGs were categorized to antimicrobial resistance (AMR) families (or metabolic processes) and resistance mechanisms were assigned as described (Liu and Pop, 2009). Based on abundance information of ARGs, heat map of abundance clustering for the top 30 highly abundant ARGs was drawn. Then, two circle charts were drawn to describe the overall proportion and distribution of the most redundant ARGs across the two soil types, on one hand, and the most redundant metabolic processes at the bacterial phylum level in ARO (antibiotic-resistant organisms) datasets, on the other hand. The first circle chart was divided into two parts of which the first refers to soil types (right side) and the second part refers to the most highly abundant ARGs (left side). The second chart was drawn to incorporate the most highly abundant bacterial phyla in terms of ARGs (right side) and the metabolic processes related to these highly abundant ARGs (left side). Stretches of different colors in the right side of inner circle refer to sum of relative abundance of all ARGs per sample in the first circus and sum of relative abundance of all metabolic processes per bacterial phylum in the second circus. While, color stretches in the left side refer to sum of relative abundance of all samples per ARG in the first circus and all bacterial phyla per metabolic process in the second circus. Colors in the right side of outer circle refer to detailed relative abundance of different ARGs per sample in the first circus and detailed relative abundance of different metabolic processes or AMR gene families per bacterial phylum in the second circus. While, colors in the left side refer to detailed relative abundance of different samples per ARG in the first circus and detailed relative abundance different bacterial phyla per metabolic process in the second circus.
Validation of selected highly abundant antibiotic resistance genes via qPCR
Total RNAs were isolated from samples gathered for each soil type using RNA PowerSoil® Total RNA isolation kit (Mo Bio, cat. no. 12866-25) following manufacturer instructions. DNA contamination was removed from RNA samples using RQ1 RNase-free DNase (Promega, Madison, WI, United States) and eventually checked by PCR. Then, primers of four highly abundant ARGs in plant rhizosphere in addition to those of 16S rRNA of Bacillus subtilis, used as a house-keeping gene, were designed using Netprimer software2 following standard criteria (Supplementary Table 1). Expression level of the different ARGs was detected by qPCR using Agilent Mx3000P System (Agilent technology, United States). Maxima™ SYBR Green/ROX qPCR was done as previously described (Bahieldin et al., 2015). Similar amounts of RNAs of the two samples were used and calculations were made to detect the expression level of ARGs relative to that of the housekeeping gene.
Results
Statistics of antibiotic resistance open reading frames
Non-redundant gene catalogs resulting from mWGS were subjected to Blastp with CARD (the Comprehensive Antibiotic Research Database) in order to detect and characterize ORFs referring to the ARGs and metabolic processes existing in the rhizobiome of Moringa oleifera and surrounding bulk soil microbiomes. Unique ORFs described in Supplementary Tables 2, 3 were generated either from assembled reads of a given soil sample or from scaftigs resulted from gathering and re-assembling unassembled reads of all samples (NOVO_MIX) referring to low abundant reads, respectively. The two tables described these ORFs, after being aligned with analog sequences in CARD, in terms of the mismatch and gap sizes (in nt) and nucleotides assigning gene start and end points. A total of 809 and 1,142 ORFs were generated from the two previously mentioned categories of gene assembly (Supplementary Tables 2, 3, respectively). Of which, a number of 1,416 ORFs across the two categories of gene assembly referring to a total of 183 ARGs were detected (Supplementary Table 4). The results in Supplementary Table 4 indicated that a number of 79 of these ARGs was represented by only a single ORF (referring to non-abundant ARGs), while the other ARGs were represented by two or more ORFs (referring to abundant, highly abundant, and the most highly abundant ARGs). Supplementary Figure 1 indicated the 30 most abundant ARGs referring to the top highly abundant ARGs with > 20 ORFs/ARG. Supplementary Table 4 indicates that the most highly abundant ARGs include mtrA (189 ORFs), soxR (104 ORFs), oleC (103 ORFs), golS (88 ORFs), and novA (74 ORFs) genes.
Differential abundance of antibiotic resistance genes
The results in Figure 1 indicated almost no differences in the total number of ARGs in samples of the two soil type (average of ∼125 ARGs/sample), while abundance of these ARGs in terms of ORFs/ARG was extremely higher in rhizosphere of M. oleifera than that in bulk soil microbiome. Supplementary Table 4 indicates abundance of ARGs across types of microbiome surrounding M. oleifera, while Figure 2 refers to the most highly abundant ARGs with > 20 ORFs/ARG. Abundance of these 14 ARGs in samples within each soil type was homogeneous, except for three genes, e.g., myrA, vanSO and dfrA3, that showed non-homogeneous abundance of ORFs among bulk soil microbiome samples (Supplementary Table 5). These three ARGs were not analyzed further.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g001.jpg)
Figure 1. Number (A) and abundance (no. ORFs) (B) of antibiotic resistance genes (ARGs) of rhizobiomes (Rl–R3 in red) and bulk soil (Sl–S3 in blue) microbiomes surrounding Moringa oleifera.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g002.jpg)
Figure 2. Means of the top highly abundant antibiotic resistance genes (>20 ORFs/ARG) across type of microbiome surrounding Moringa oleifera. Patterned columns refer to ARGs whose numbers of ORFs in bulk soil samples are non-homogeneous. Box colors of other ARGs match those of Figure 4.
Estimated abundance (in unit ppm) of ARGs for the different samples of the two types of microbiomes is shown in Supplementary Table 5, while the most highly abundant ARGs (>20 ORFs/ARG) in microbiomes of the two soil types are shown in Figure 3A and described based on CARD information in Supplementary Table 6. The results in the table indicated that abundance of ARGs with > 20 ORFs/ARG was significantly higher in rhizobiome than that of bulk soil microbiome, except for golS gene, that showed no significant difference (P = 0.315) between samples of the two soil types. The results of relative abundance of ARGs supported those of abundance only for mtrA, novA, and parY mutant, as relative abundance of these three ARGs in rhizobiome is higher than that of bulk soil microbiome (Figure 3B). However, the results for golS (Pontel et al., 2007) and rpoB2 (Ishikawa et al., 2006) genes showed higher relative abundance in bulk soil microbiome than that in rhizobiome (Figures 3B, 4). We assume the results of both abundance and relative abundance of ARGs ought to be considered.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g003.jpg)
Figure 3. Abundance in unit ppm (A) and relative abundance (B) of the top highly abundant antibiotic resistance genes (>20 ORFs/ARG) in different samples of rhizobiomes (Rl–R3) and bulk soil (Sl–S3) microbiomes surrounding Moringa oleifera. Unit ppm was calculated by magnifying 106 times of the original abundance data. Box colors of other ARGs match those of Figure 4.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g004.jpg)
Figure 4. Circle chart of the top highly abundant antibiotic resistance genes (>20 ORFs/ARG) in different samples of rhizobiomes (Rl–R3) and bulk soil (S1–S3) microbiomes surrounding Moringa oleifera. The chart is divided into two sides, where the right one refers to samples (S and R) information and the left refers to antibiotic resistance genes (ARGs) information. The wideness of different scales of inner and outer circles refers to relative abundances of ARGs and samples. Colors in inner circle refer to sum of relative abundance of all ARGs per sample (right side) and sum of relative abundance of all samples per ARG (left side). While, colors in outer circle refer to detailed relative abundance of different ARGs per sample (right side) and detailed relative abundance of different samples per ARG (left side).
Differential abundance of bacterial phyla and antibiotic resistance mechanisms
ARGs were further categorized to metabolic processes or AMR families at the bacterial phylum level (Supplementary Table 7). The analysis was focused on 12 bacterial phyla and 55 categories of ARG-related metabolic processes. Based on abundance results, Figure 5 describes the five most common phyla, e.g., Actinobacteria, Proteobacteria, Acidobacteria, Chloroflexi, and Firmicutes, and the 13 most common metabolic processes. Abundance of phyla with > 30 ORFs regardless of metabolic process, abundance of metabolic processes regardless of phylum, as well as abundance of the interaction between phyla and metabolic processes are shown in Supplementary Figures 2, 3, and 5, respectively. Actinobacteria, Proteobacteria, Acidobacteria represent the most abundant bacterial phyla regardless of metabolic processes (Supplementary Figure 2). The eight most abundant metabolic processes shown in Supplementary Table 8 almost refer to all the previously detected highly abundant ARGs in Supplementary Table 6, except for rpoB2 and ileS genes, whose metabolic processes, e.g., rifamycin-resistant beta-subunit of RNA polymerase (rpoB) and Bifidobacterium ileS conferring resistance to mupirocin, respectively, were less abundant (with < 30, Supplementary Figure 3), thus, were not analyzed further. These selected metabolic processes refer to three resistance mechanisms (Supplementary Table 8) namely antibiotic efflux, antibiotic target alteration and antibiotic target protection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g005.jpg)
Figure 5. Circle chart of metabolic processes at the bacterial phylum level of the top highly abundant antibiotic resistance genes (>20 ORFs/ARG) in different samples of rhizosphere and bulk soil microbiome s surrounding Moringa oleifera. The chart is divided into two sides, where the right one refers to bacterial phylum information and the left refers to metabolic processes information. The wideness of different scales in inner and outer circles refers to relative abundances of metabolic processes and bacterial phylum. Colors in inner circle refer to sum of relative abundance of all metabolic processes per bacterial phylum (right side) ai1d sum of relative abundance of all bacterial phyla per metabolic process (left side). While, colors in outer circle refer to detailed relative abundance of different metabolic processes per bacterial phylum (right side) and detailed relative abundance of different bacterial phyla per metabolic process (left side).
AMR families referring to antibiotic efflux in this study include resistance-nodulation-cell division (RND), ATP-binding cassette (ABC), major facilitator superfamily (MFS) antibiotics as well as the two-component regulatory kdpDE (potassium dependent D/E) system. AMR families referring to antibiotic target alteration included glycopeptide resistance gene cluster (vanRO), aminocoumarin resistance parY and aminocoumarin self-resistance parY. While, AMR family referring to antibiotic target protection included RbpA bacterial rpoB-binding protein, only. As indicated in Supplementary Table 6, ARGs of mtrA, soxR, and golS participate in the metabolic process of RND antibiotic efflux, while oleC and novA genes participate in the metabolic process of ABC antibiotic efflux. soxR and kdpE genes participate in the metabolic processes of MFS antibiotic efflux and two-component regulatory kdpDE (potassium dependent D/E) system, respectively. Thus, soxR gene is the only participating in two metabolic processes. vanRO, parY mutant, and rbpA genes participate in the metabolic process of glycopeptide resistance gene cluster, in the metabolic process of aminocoumarin resistance/self-resistance, and in the metabolic process of RbpA bacterial rpoB-binding protein, respectively (Supplementary Table 6). The results of circle chart of Figure 5 referring to the left and right sides of inner circle support those of Supplementary Figures 2, 3 in terms of abundance, respectively. In terms of relative abundance, the results in Figure 5 (outer circle) and Supplementary Figure 4 indicate that Proteobacteria relatively dominates in terms of RND antibiotic efflux pump followed by Actinobacteria. However, the other seven metabolic processes are relatively dominated by Actinobacteria. Proteobacteria also showed high relative abundance of processes MFS antibiotic efflux pump and two-component regulatory kdpDE system. The results involving the most highly abundant ARGs, along with their metabolic processes and resistance mechanisms referring to the five bacterial phyla are summarized in Figure 6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g006.jpg)
Figure 6. Summary of the resistance mechanisms along with the participating AMR gene families and their ARGs of bacterial phyla across soil microbiomes surrounding Moringa oleifera. Different intensities of the red color of the circles refer to the level of differential abundance in bacterial phyla.
Validation of the resistomic data for the two soil types surrounding M. oleifera was done at the metatranscriptomic level via qPCR (Supplementary Figure 5). The experiment involved four highly abundant ARGs namely mtrA (acc. no. ARO:3000816), soxR (acc. no. ARO:3004107), oleC (acc. no. ARO:3003748), and novA (acc. no. ARO:3002522). Metabolic process of the first two ARGs is the RND antibiotic efflux pump, while that of the second two ARGs is the ABC antibiotic efflux pump. The results of qPCR for the metatranscriptomic data in Supplementary Figure 5 align with those of the in silico genomic dataset in Figure 3.
Discussion
Native resistomes of the wild plants can provide unexplored antibiotic resistance mechanisms that add to our understanding of their contribution to the incidence of HGT (Berendonk et al., 2015). Obermeier et al. (2021) claimed that native plants likely harbor substantially diversified secondary metabolic processes that result in the recovery of a more diversified pristine resistome and antibiotics. Therefore, we thought it is important to explore native resistomes of naturally growing plants in the wild. Dissemination of ARGs is not only a concern within clinics, rather, this issue is prevalent in bacterial communities of the open environment (Chen et al., 2019). Initially, soil bacteria pose no direct threat to human health, while horizontal transfer of ARGs across different hosts to human pathogenic bacteria (via MGEs, e.g., plasmids, phages and integrons) can represent a major concern. Therefore, it is a demand to study resistome signatures and related resistance mechanisms in rhizobiomes of edible plants, and the factors that promote dissemination of ARGs.
Antibiotic resistance mechanisms generally include the antibiotic modifications, alteration of antibiotic target (decreasing the affinity of the binding site for the drug), antibiotic efflux, and global cell adaptations (Munita and Arias, 2016). Recent reports assign a new class of resistance mechanism known as antibiotic target protection, where a given resistance protein physically associates with an antibiotic target (or binding site) in the bacteria to hide it from the cognate antibiotic (Wilson et al., 2020). This resistance mechanism does not result in a permanent target modification, or requires a change of the nature of the target protein as the latter can be released again and becomes effective upon disappearance of the concerned antibiotic (Wilson et al., 2020).
The results in the present study indicated that interaction of M. oleifera roots with co-existing soil microbes did not result in the emergence of new soil ARGs (Figure 1), rather it affected differential abundance of existing ARGs as a consequence of differential abundance of microbes harboring these ARGs mostly in their MGEs. In accordance with the results of the present study, rate of resistance of pathogens to existing antibiotics was reported to be higher than that referring to the newly discovered ARGs (Boucher et al., 2009). Except for golS gene, abundance of ARGs in rhizobiome of M. oleifera was significantly higher than that of bulk soil microbiome (Figure 3). The main difference between microbiomes of bulk soil vs. those in the plant rhizosphere is the occurrence of plant exudation that participates in shaping structure and physiology of their respective rhizobiomes (Mendes et al., 2019). Exudation pattern mainly depends on plant genotype and developmental stage (Chaparro et al., 2014; Mönchgesang et al., 2016). Exudates can include minerals (e.g., carbon and nitrogen), simple molecules (e.g., sugars, organic acids, and secondary metabolites), and more complex polymers (e.g., mucilage) (Sasse et al., 2018).
In terms of relative abundance, mtrA, novA, and parY mutant genes showed higher level in rhizobiome ARGs than that in bulk soil microbiome, while those of golS and rpoB2 genes showed opposite results (Figures 3, 4 and Supplementary Table 4). Across the highly abundant metabolic processes (with > 30 ORFs across phyla), Actinobacteria and Proteobacteria were the most abundant phyla followed by Acidobacteria, Chloroflexi, and Firmicutes (Figure 5, Supplementary Figure 2, and Supplementary Table 5). These metabolic processes refer to three resistance mechanisms namely antibiotic efflux, antibiotic target alteration and antibiotic target protection (Figure 6). Processes of antibiotic efflux pump included RND (Eswaran et al., 2004; Murakami et al., 2006; Blair and Piddock, 2009), ABC (Fath and Kolter, 1993), and MFS (Pao et al., 1998; Saier et al., 1999; Li and Nikaido, 2009), as well as the two-component regulatory KdpDE system (Freeman et al., 2013), while processes of antibiotic target alteration included glycopeptide resistance gene cluster for vancomycin resistance (vanRO) (Courvalin, 2006; Gudeta et al., 2014), and aminocoumarin resistance and self-resistant parY (Schmutz et al., 2003, 2004; Li and Heide, 2006), and those of antibiotic target protection included RbpA bacterial rpoB-binding protein, only (Newell et al., 2006; Hu et al., 2012; Figure 6). It was possible in the present study to connect highly abundant ARGs in samples of the two soil types (Figure 4) with their metabolic processes along with the bacteria promoting these processes (Figure 5). These connections are summarized in Figure 6.
Resistance mechanisms in resistome of Moringa oleifera
Antibiotic efflux mechanisms
Efflux pump is an active transport protein with important role in bacterial homeostasis via expulsion of toxic substances or antibiotics (Fernández and Hancock, 2012). The efflux system results in the reduced accumulation of antibiotics and possibly results in multidrug resistance. Induction of efflux pump of bacterial antibiotics is an instant response of bacterial pathogens to potentiate their survival and resist toxic compounds during infectivity period. Efflux of toxins by bacteria ought to be fast in order to promote proper bacterial proliferation way before antibiotic therapy or host defense mechanisms are initiated (Koprivnjak and Peschel, 2011).
There are five bacterial efflux pumps classified into five families or metabolic processes. They are the resistance-nodulation-division (RND) family, the ABC family, the MFS, the small multidrug resistance (SMR) family, and the multidrug and toxic compound extrusion (MATE) family (Piddock, 2006; Poole, 2007; Sun et al., 2014). RND only exists in Gram-negative bacteria, mostly of phylum Proteobacteria, while other families are distributed among both Gram-positive and Gram-negative bacteria (Sun et al., 2014). Export of toxins by the different efflux pumps is energy-dependent, where ABC efflux pump relies on ATP hydrolysis, while the other pumps use proton motive forces for efflux (Chitsaz and Brown, 2017). The three efflux processes of RND, ABC, and MFS followed by the two-component regulatory KdpDE system seem to be dominating in rhizospheric region of M. oleifera due to the high abundance of ARGS of mtrA, soxR, oleC, novA, and kdpE (Figures 5, 6 and Supplementary Figure 3). While, efflux process of RND seems also to be dominating in bulk soil surrounding M. oleifera due to the high abundance of golS.
RND superfamily transporters3 are made of large polypeptide chains located in the cytoplasm. RND efflux is composed of seven families including the hydrophobe/amphiphile efflux (or MtrCDE efflux) and the heavy metal efflux (HME) (Tseng et al., 1999). RND transporters are involved in maintaining cell homeostasis and in exporting a broad spectrum of toxic compounds and virulence determinants (Coyne et al., 2010). In the present study, the most abundant ARG in rhizobiome of M. oleifera, e.g., mtrA (multiple transferable resistance A) gene (Figures 3, 4), is known to encode a transcriptional activator or inducer of the MtrCDE multidrug efflux pump mostly in human pathogenic bacteria, e.g., Neisseria gonorrhoeae (Rouquette et al., 1999; Ohneck et al., 2011; Handing et al., 2018), and Pseudomonas aeruginosa (Du et al., 2018). MtrCDE efflux pump is a member of the highly abundant hydrophobic and amphiphilic efflux RND family. N. gonorrhoeae is the causative pathogen of gonorrhea and can also evade protective human immune responses (Wi et al., 2017; Nudel et al., 2018). MtrA is a gonococcal protein belonging to the AraC/XylS family that participates in the construction of the bacterial multiple transferable resistance (MTR) complex (Nikaido, 1996; Rouquette et al., 1999). The latter promotes an energy-dependent MtrCDE efflux pump comprising three cell envelope proteins namely MtrC, MtrD, and MtrE. These three proteins, in turn, serve in exporting antimicrobial agents outside the bacterial cell. These antimicrobial agents include hydrophobic agents (HAs), such as free fatty acids, bile salts, gonadol steroids, and antibacterial peptides (Handing et al., 2018). Many reports indicated that the MtrCDE efflux system can also act on antibiotics, like penicillin, erythromycin, cephalosporin, and rifampin (Olesky et al., 2006; Handing et al., 2018). On the other hand, HA-sensitive bacteria harbors another gene, e.g., mtrR, that otherwise encodes a transcriptional repressor of MtrCDE multidrug efflux pump system (Lucas et al., 1997). The latter gene does not exist in soil microbiome of M. oleifera, thus, does not seem to be involved in modulating the operon system of MtrCDE efflux pump in the resistome of M. oleifera (Figure 7). Ma et al. (1994) indicated that the loss of repressor activity is correlated with increased bacterial resistance to antimicrobial compounds. Thus, high abundance/relative abundance of mtrA gene (promoting positive regulation) and absence of mtrR gene (promoting negative regulation) in the resistome of M. oleifera can results in higher influence on virulence and the ability of the bacteria to resist antimicrobial agents more efficiently, in accordance with prior reports (Ohneck et al., 2011). This information indicates that horizontal transfer of such an mtrA+/mtrR– resistance determinant to human gut can be dangerous. However, Handing et al. (2018) declared that bacterial acquisition of antibiotic resistance is associated with a fitness cost. Other reports claim that bacteria also need the negative regulation as over-activation of efflux pump can slow down rate of bacterial growth (Nikaido, 1996), thus, modulation of positive/negative regulation of the efflux pump is required for AMR, on one hand, and for optimal growth rate, on the other hand.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g007.jpg)
Figure 7. Demonstration of RND multidrug efflux system contributed by MtrA gonococcal protein (positive regulation) in rhizobiome of M oleifera. This efflux system is effective due to action of the downstream proteins MtrC, MtrD, and MtrE in addition to the absence of the MtrR transcriptional repressor. MtrR contributes to slower growth rate to maintain proper fitness cost.
The golS gene encodes a mercury resistance (MerR)-like sensor, which is highly selective for gold (Au) ions (Brown et al., 2003; Pérez Audero et al., 2010). It controls proteins GolT, which is a P-type ATPase, and GolB, which is a small cytoplasmic metal-binding RND protein located in the inner membrane (Pontel et al., 2007). The latter two proteins participate in RND-dependent CBA efflux system that induces a resistance response against gold salts inside Salmonella cells. Salmonella is a poisonous bacterial disease that affects the intestinal tract and mostly infect human through contaminated water or food. The increase of GolS concentration could also result in direct sequestration of intracellular free Au levels (Brown et al., 2003). However, when high Au levels persists, CBA efflux system will be induced. The CBA efflux system is a tripartite protein complex consisting of subunits C, B, and A that generally direct efflux of metal ions or drugs from either the cytoplasm or the periplasm into the extracellular space. It was speculated that this efflux system likely results in transporting cellular metabolites (either damaged by Au or after forming a complex with this ion), rather than in transporting the toxic metal ion per se (Brown et al., 2003). As the efflux system involving golS gene is Salmonella-specific (Checa et al., 2007), it is expected to see in the present study that this gene is less abundant in the rhizospheric soil of M. oleifera (Figures 3, 4) where this microbe does not naturally interact with plant roots, thus golS gene is unlikely to be horizontally transferred through the plant system.
Recent reports indicated that soxR gene or transcription factor is a global RR that can directly or indirectly modulate bacterial multidrug resistance, growth and fitness (Palma et al., 2005; Li et al., 2017). SoxR promotes several efflux pump genes of the RND family including adeB, adeJ, and adeG (Marchand et al., 2004; Gu and Imlay, 2011). In alignment with the results of the present study, several reports indicated that soxR contributes to multidrug resistance mainly in members of the Gram-negative phylum Proteobacteria (Palma et al., 2005; Sakhtah et al., 2016). As an indirect regulator, SoxR protein is firstly oxidized in response to superoxide-generating agents in order to oxidize another protein namely soxS referring to the SoxR/SoxS paradigm (Pomposiello and Demple, 2001; Li et al., 2017). Although SoxR can be oxidized by extracellular superoxide anion, this anion might not be able to diffuse across the double membrane of the bacteria to oxidize SoxR. Instead, the latter can be rather oxidized by a permeable nitric oxide (NO) radical (Pomposiello and Demple, 2001). SoxS, in turn, promotes overexpression of a multidrug efflux machinery namely AcrAB-TolC in Escherichia coli, a member of the Proteobacterial family Enterobacteriaceae. AcrAB-TolC is a tripartite transporter that effluxes periplasm substrates out of the cell (White et al., 1997; Ruiz and Levy, 2014). This resistance pattern to resistome in the rhizospheric region of M. oleifera as SoxS does not exist. On the other hand, SoxR can directly regulate target genes, e.g., a 6-gene regulon, that promotes efflux pump in members of the Proteobacterial families Pseudomonadaceae (Pseudomonas aeruginosa) and Moraxellaceae (Acinetobacter baumannii) (Bialek-Davenet et al., 2011; Li et al., 2017). In particular, A. baumannii is a human pathogen that causes nosocomial and bloodstream infections (Peleg et al., 2008). SoxR was also reported to facilitate induction of MFS efflux pump genes in a way similar to that of the major facilitator MfsA (Saidijam et al., 2006; Dulyayangkul et al., 2016). The MFS system is composed of membrane transport proteins or facilitators that promote movement of solutes (e.g., drugs, metabolites, oligosaccharides, amino acids, etc.) across membranes as a consequence of chemiosmotic gradients (Marger and Saier, 1993). These facilitators open to either the extracellular or cytoplasmic side and simultaneously seal their opposing face to block bidirectional passage across the membrane, while only allow substrate entrance into the cell (Abramson et al., 2003). As soxR gene that promotes several efflux pump genes of the RND family is highly abundant in rhizobiome of M. oleifera, while soxS gene is absent (Figures 3, 4), then we expect that the soxR+/soxS– resistance determinant can only represent a major risk on human health upon laterally transferred to gut microbiome in case soxR gene directly promoted efflux pump following the resistance pattern in either P. aeruginosa or A. baumannii (Figure 8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g008.jpg)
Figure 8. Demonstration of the RND multidrug efflux system contributed by SoxR response regulator in rhizobiome of M oleifera following two resistance patterns, one in Escherichia coli, and the other in Pseudomonas aeruginosa and Acinetobacter baumannii. The first pattern requires three genes namely adeB, adeJ, and adeG to oxidize SoxR that, in tum, oxidizes SoxS that eventually promotes the efflux system. The second pattern only requires a 6-gene regulon to be activated by SoxR in order to directly promote the efflux system. The first pattern unlikely apply to resistome in the rhizospheric region of M. oleifera as SoxS does not exist.
The ABC transporter superfamily is composed of several membrane and membrane-associated (e.g., AAA ATPases) proteins that transport a variety of substrates (lipids, metabolic products, drugs, etc.) across extra- or intracellular membranes. Some of these proteins are involved in tumor resistance (Dean et al., 2001). During this ABC system, ATPases utilize ATP to provide energy required for the bidirectional translocation (influx/efflux) of substrates across membranes (Sauna et al., 2009). Of which, the OleC transporter of Streptomyces antibioticus is constituted by two proteins, e.g., OleC that recognizes, binds and hydrolyses ATP, and OleC5 that is a hydrophobic membrane protein (Olano et al., 1996). oleC gene, referred to as the orf4 resistance determinant, shows similarity with proteins harboring ATP-binding domains, a feature of the ABC-transporter superfamily. This ARG dominates in the rhizospheric region of M. oleifera (Figure 3). The ABC transporter system allows secretion of oleandomycin out of the cell as an efflux transport mechanism of self-resistance against this antibiotic (Ma Rodriguez et al., 1993). Up to date, genus Streptomyces was reported to produce the majority of the clinical antibiotics (e.g., neomycin, cypemycin, grisemycin, bottromycins, chloramphenicol, aminocoumarin, etc.) mostly by S. antibioticus (de Lima Procópio et al., 2012) and S. lividans (Lee et al., 2007) that still can produce new antibiotics. Interestingly, S. antibioticus harbors another enzyme that catalyzes the last step of oleandomycin biosynthesis to release glucose from the inactivated oleandomycin and eventually generate new active antibiotics (Salas et al., 1994). In addition, novA, a highly abundant ARG in the present study (Figure 3), is an aminocoumarin resistance gene of such a biosynthetic gene cluster in S. lividans (Schmutz et al., 2003). The gene encodes a type III ABC-like transporter that confers a moderate level of resistance as this transporter presumably acts in the sequestration of the antibiotic into the medium (Méndez and Salas, 2001). The gene cluster harboring novA gene also contains two gyrB genes, e.g., gyrBS and gyrBR (Thiara and Cundliffe, 1989). Unlike novA gene, gyrBR presumably provides the principal resistance mechanism in the aminocoumarin antibiotic producer Streptomyces lividans due to the synthesis of coumarin-resistant topoisomerases (Schmutz et al., 2003). As abundance of oleC and novA genes of genus Streptomyces is high in the rhizosphere soil of M. oleifera, then, there is a chance that we can decipher new antibiotics as secondary metabolites generated by members of this genus that can be used further in clinical trials against newly emerging diseases. Fortunately, invasive infections with members of this genus are rarely encountered in clinical practice (Kapadia et al., 2007), thus, if their ARGs, like oleC and novA genes of rhizobiome of M. oleifera, were not horizontally transferred through contaminated food to other clinical species, it is unlikely that these bacteria can pose a risk to human health.
The two-component regulatory system (TCS) is ubiquitous and crucial in bacteria to maintain homeostasis and sense, thus, respond to changes in the surrounding environment (Cardona et al., 2018). TCS also responds directly in bacteria to the presence of antibiotics and acts as a major player in the realm of infectious diseases caused by pathogenic bacteria (Cardona et al., 2018; Tierney and Rather, 2019). Mechanisms of antibiotic resistance via TCS are divided into four categories of which the increased drug efflux dominates. The latter regulatory system comprise a histidine kinase (HK) for sensing environmental signals, and a response regulator (RR) for mediating cellular response and altering expression of target genes (Freeman et al., 2013). Examples include the membrane-bound HK and the cytosolic RR that form KdpD (HK)/KdpE (RR) TCS. This system refers to the regulation of the potassium (K+) dependent- (Kdp−) ATPase pump of operon KdpFABC. The gene encoding KdpE was proven to be highly abundant in the rhizobiome of M. oleifera (Figure 3). It was reported that KdpD/KdpE system can increase the ability of bacteria to cause the disease and/or survive in the host cell (Alegado et al., 2011). In Salmonella typhimurium, KdpD was proven to maintain bacterial growth in macrophage of nematode (Caenorhabditis elegans) cell lines as worms fed on the kdpD mutant strain lived longer than those fed on the wild type bacteria (Alegado et al., 2011). In terms of KdpE, prior results indicated its role as a RR of expression of a range of virulence genes in Staphylococcus aureus and Enterohemorrhagic E. coli (EHEC) (Hughes et al., 2009; Xue et al., 2011). This induced expression coincided with the increased survival rate of the two bacteria in human macrophages. Interestingly, deletion of kdpD/kdpE system in S. aureus resulted in the altered expression levels of > 100 genes regulating synthesis of surface protein and capsular polysaccharide, of alpha- and gamma-hemolysin toxins (offensive capabilities), and of metalloproteinase aureolysin that participates in inhibiting phagocytosis and killing the bacteria by neutrophils (defensive capabilities) (Laarman et al., 2011; Freeman et al., 2013). Moreover, kdpE gene was frequently reported to act as a RR gene in conferring resistance to the aminoglycoside antibiotic kanamycin via efflux pump system in E. coli (Hirakawa et al., 2003; Lv et al., 2021). Interestingly, KdpE is involved in two-compound systems with KdpD, QseC, or Cra, to confer pleotropic effects on the structure, and on the offensive and defensive capabilities of several bacteria including Salmonella typhimurium, Staphylococcus aureus and EHEC. Then, high abundance of kdpE gene in the rhizosphere soil of M. oleifera can pose very high risk to human health via contaminated food even if it is not transferred horizontally to pathogens of human gut.
Antibiotic target alteration and protection
Resistance mechanism of alteration of antibiotic target refers to the ability of the bacteria to decrease the affinity for the drug by modifying its target site. The most common example is the glycopeptides (ex., vancomycin) that kill bacteria by inhibiting cell wall synthesis (Munita and Arias, 2016). Vancomycin resistance is common in Enterococcus faecium, a commensal microbe in the gastrointestinal tract of human that turned to be pathogenic causing serious diseases like neonatal meningitis or endocarditis (Arias and Murray, 2012). Vancomycin resistance occurs due to the action of one or more of the 11 van operons and acquisition of resistance in clinical isolates was previously proven to be due to horizontal transfer via MGEs, e.g., conjugative or non-conjugative plasmids carrying Tn3 transposons (Arias and Murray, 2012). The clinical isolates of E. faecium harbor vanA and vanB operons that are either carried on plasmids or integrated in the bacterial chromosome (Munita and Arias, 2016). Vancomycin resistance is also common in Rhodococcus equi, a soil bacterium that can cause pneumonia in young foals and in immune-deficient adult humans (Franklin, 2016). The latter bacterium harbors vanO operon that includes gene cluster vanSO/vanRO acting as a TCS where VanSO is the HK and VanRO is the RR. vanO operon also includes a vanHOX gene cluster positioned opposite the vanSO/vanRO gene cluster in addition to five ORFs, e.g., ORFs1-3 and ORFsA-B, with no clear mode of action (Gudeta et al., 2014). The VanSO/VanRO system acts in promoting expression of the vanHOX gene cluster where initial sensing of the accumulation of substrates by VanSO inhibits glycosyltransferase activity, which, in turn, results in ATP-dependent phosphorylation of the RR VanRO. The latter, then, induces expression of the vanHOX gene cluster where vanH gene encodes a dehydrogenase enzyme necessary for the synthesis of new peptidoglycan precursors, vanX gene encodes an enzyme that destroys the normal D-Ala-D-Ala-ending precursors, while vanO gene encodes a ligase that synthesizes altered D-Ala-D-Lac substrate of penicillin binding protein (PBP) to reduce binding affinity of vancomycin and confer resistance (Figure 9) (Reynolds, 1989; Gudeta et al., 2014; Munita and Arias, 2016). vanRO gene is highly abundant in rhizobiome of M. oleifera (Figure 3). Being carried on a plasmid, the rhizospheric vanRO gene of M. oleifera along with its operon become a real threat of horizontal transfer to a human pathogen, thus, occurrence of acute diseases.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g009.jpg)
Figure 9. Demonstration of the antibiotic target alteration system existing in rhizobiome of M. oleifera and contributed by vanRO gene encoding the response regulator (RR) and vanSO gene encoding the histidine kinase (HK) of the VanRO/VanSO two-component system. The two genes represent a cluster of the vanO operon of which activated VanRO promotes expression of the other gene cluster of vanO operon namely vanHOX The latter cluster, in turn, acts in reducing binding affinity of vancomycin, thus, confers antibiotic resistance.
Topoisomerase IV (topo IV) is a type II topoisomerase involved in relaxing superhelical DNA prior replication and the decatenation of daughter chromosomes following DNA replication, e.g., two energetically driven processes by ATP hydrolysis (Schmutz et al., 2004). The topo IV enzyme is composed of the subunits ParC that contains the catalytic center for DNA cleavage/rejoining and ParE that contains the catalytic center for ATP hydrolysis. In Actinobacteria, e.g., Streptomyces coelicolor and Bifidobacterium longum, analogs of the parC and parE genes, e.g., parX and parY, exist. The latter type II topoisomerase genes exist in the biosynthetic gene cluster conferring aminocoumarin antibiotic resistance. This antibiotic is a potent inhibitor of the subunit B of the other drug target type II topoisomerase, namely gyrase (Maxwell, 1997). Members of genus Streptomyces harbor a modified de novo synthesized subunit of the drug target gyrase, namely aminocoumarin-resistant gyrase B to protect themselves from the inhibitory effects of aminocoumarin (Schmutz et al., 2004). The gene encoding this subunit, namely gyrBR, exists in the aminocoumarin gene clusters. In addition, there is an additional resistance gene, namely parYR or parY mutant, which is an analog of parY also exists. These modified versions of bacteria harboring GyrBR and ParYR are no longer targets for aminocoumarin antibiotics as they provide these bacteria with self-protection against the toxic effects of their own antibiotic, e.g., aminocoumarin. The mutant version of the parY gene is highly abundant in rhizobiome of M. oleifera (Figures 3, 4). Again, infection with members of the genus Streptomyces is rarely encountered in human clinical practice (Kapadia et al., 2007) and transfer of this ARG to clinical isolates of this genus almost represents no risk to human health, while horizontal transfer to pathogens of other bacterial phyla might pose a great risk.
In terms of antibiotic target protection system, bacteria prevent the antibiotic to reach its binding site (Newell et al., 2006; Wilson et al., 2020). Physical association of a resistance protein with the antibiotic target reduces the ability of the antibiotic to reach receptor of its target molecule. The rbpA (or RNA-polymerase binding protein A) gene is an ARG that is highly abundant in rhizosphere of M. oleifera. The encoded protein in Mycobacterium tuberculosis functions as a transcription activator that binds the core RNA-polymerase (RNAP), e.g., a target of the antibiotic rifampicin, to stabilize the principal promoter-specific housekeeping σA factor, thus, to promote action of RNAP (Hu et al., 2012). Stimulation of σA activity by RbpA is required for the expression of virulence genes and for pathogen proliferation in human macrophages (Hu et al., 2012). M. tuberculosis is the causative pathogen of tuberculosis (TB) in human of which rifampicin is the first-line drug used to cure the disease. It was suggested that RbpA binds cluster I of the β subunit of RNAP to prevent association/induce dissociation of the antibiotic and RNAP (Hu et al., 2012). The rbpA gene was first detected in E. coli (Maeda et al., 2000) and a year after in Streptomyces coelicolor (Paget et al., 2001). This gene was found to be highly upregulated during rifampicin treatment. Several resistance mechanisms were speculated for this gene of which we selected the two most appropriate for its action (Figure 10). They are the antibiotic target protection and antibiotic target modification. However, we assume that the first mechanism is likely to represent the rbpA gene action as the enzyme RNAP is still valid to cause transcription of the virulence genes after disappearance of the antibiotic. Figure 10 indicates that rate of expression by RNAP increases in the presence of RbpA compared with that in the absence of this protein (Hu et al., 2012). We expect that neither E. coli nor Streptomyces coelicolor can pose a serious health problem for human as RbpA is dispensable for survival of S. coelicolor (Newell et al., 2006). Rather, the problem rises if rbpA gene is horizontally transferred to clinical isolates of M. tuberculosis as the gene is, otherwise, essential for its survival (Forti et al., 2011). In addition, the rbpA encoded protein helps M. tuberculosis reaches the maximum expression levels of its virulence genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990169/fmicb-13-990169-HTML-r2/image_m/fmicb-13-990169-g010.jpg)
Figure 10. Demonstration of the two suggested resistance mechanisms existing in rhizobiome of M. oleifera against rifampicin conferred by RbpA to allow RNA polymerase (RNAP) properly conduct transcription of virulence genes. Structure of RNAP is made of the holoenzyme (representing the catalytic core that is made of five subunits of ααββ’ω) and a promoter-specific σ factor (for promoter recognition and conduction of transcription on the DNA antisense strand) located between residues Cys478 and Cys677 of β subunit. (A) No rifampicin, (B) presence of rifampicin and absence of RbpA, (C) resistance mechanism of antibiotic target protection where RbpA occupies rifampicin target site in β subunit of RNAP to make the binding site inaccessible to the antibiotic (D) resistance mechanism of alteration of antibiotic target where RbpA causes allosteric changes in the rifampicin binding site in RNAP.
Plant rhizospheric soil is considered as the reservoir of ARGs due to the highly diversified microorganisms. Antibiotic resistance refers to the increase in the minimum inhibitory concentration of a compound (e.g., antibiotic) for a microbe (pathogenic or non-pathogenic) that was previously known to be sensitive at this antibiotic concentration (World Health Organization [WHO], 2012). Antibiotic resistance was hypothesized to be mainly the results of pre-existing environmental resistome. This hypothesis was proved when antibiotic-sensitive pathogens acquired resistance after introduction of new antibiotics and the gene tetA(P), for tetracyclin resistance, was proven to be a suitable biomarker of soil contamination with ARGs (Blau et al., 2018, 2019). When antibiotic substances exist in the rhizospheric region of a plant root, there might be evolutionary forces that pose a selection pressure on the plant rhizobiome and result in the emergence of microbes harboring ARGs. MGEs, such as viruses (e.g., viromes) or plasmids (e.g., plasmidomes), have uncovered the reason for rapidly disseminating antibiotic resistance within lineage of distant or even closely related microorganisms via HGT (World Health Organization [WHO], 2012). ARGs were questioned whether they are a result of gene shuffling/mutation or a result of resistance acquired through HGT (Nesme and Simonet, 2015). The latter is a genetic transformation process that was proven to promote the emergence of antibiotic resistant pathogenic bacteria (Khan et al., 2018). Recent reports interestingly proved that plant can also acquire exogenous microbial ARGs with the help of MGEs (Chen et al., 2017). Thus, MGEs can promote plant and microbial transformation processes with exogenous soil ARGs.
Plant-contaminating microbes that harbor ARGs and their vectors, e.g., MGEs are unlikely to be removed during food production, storage or processing, thus, can horizontally pass into microbiome of human gut. If the gut microbiome harbors pathogenic bacteria, then, the latter can acquire these ARGs, thus pose a great risk to human health (Chen et al., 2017, 2019). Human skin microbiome can also be exposed to antibiotic resistomes of phyllosphere and receive the same fate as that of gut microbiome (Chen et al., 2019). In addition, dead food-borne bacterial cells might lyse and existing extracellular ARGs consequently be disseminated in the environment (Levy-Booth et al., 2007). Examples of these extracellular DNAs include the self-transmissible plasmids conferring tetracycline resistance that were proven to be captured by Escherichia coli referring to the unintended transfer of this bacteria to human gut microbiome (Blau et al., 2018). Thus, HGT can refer to more than one layer of gene transfer within phyllosphere or rhizosphere microbiome of a plant species, e.g., from microbe in the rhizobiome to other microbe, from microbe to plant, from plant-contaminating microbes with ARGs to human microbiome (e.g., gut or skin), and from non-pathogenic to pathogenic microbes in either human organ. All of these layers of horizontal transfer of ARGs are mediated by the MGEs (Nardelli et al., 2012; Djordjevic et al., 2013).
In addition, concerns of pathogen resistance to antibiotic treatments were proven to be mostly limited to bacterial group namely ESKAPE (Enterococcus, Staphylococcus, Klebsiella, Acinetobacter, Pseudomonas, and Escherichia) (Boucher et al., 2009). Except for Klebsiella, ESKAPE group contains most of the selected highly abundant ARGs in the present study indicating that the transfer of these genes to analogs in the human pathogenic strains is possible. Cheng et al. (2012) indicated that the increase of antibiotic resistance in clinical pathogens is the result of gene mutation or fusion rather than the exchanges between soil and clinical resistomes, where only few resistance genes were shared between human gut and soil bacteria. Forsberg et al. (2014) reached similar conclusions and referred this to the limited rate of HGT between these two niches.
Conclusion
In conclusion, the present study provides information on the structure of rhizospheric resistome of the wild plant M. oleifera and the possible ways and consequences of disseminating the soil ARGs in the open environment in order to develop interventions to prevent transfer of these ARGs to serious human pathogens or clinical isolates.
Data availability statement
The raw data were deposited at the European Nucleotide Archive (ENA) (https://www.ebi.ac.uk/ena/browser/). Condition of deposit indicates that these raw data will be publicly available as soon as the article is published online. This is the standard condition with the deposit of nucleic acid sequences. Accession nos. (ERR10100770-74 and ERR10100781), sample no. (ERS12550318-23), and study no. (ERP139990).
Author contributions
AS, RJ, and AA: conceptualization. HA: methodology. MR: software. NB, HB, and MT: validation. LB: formal analysis. AB and HB: investigation. MT: resources. NB: data curation. AS and MR: writing—original draft preparation. RA: writing—review and editing. AA: visualization. RJ: project administration and supervision. AS: funding acquisition. All authors read and agreed to the published version of the manuscript.
Funding
This work was supported by the Princess Nourah bint Abdulrahman University Researchers Supporting Project No. (PNURSP2022R31), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.990169/full#supplementary-material
Footnotes
- ^ https://card.mcmaster.ca
- ^ http://www.premierbiosoft.com/netprimer/index.html
- ^ https://en.wikipedia.org/wiki/Efflux_(microbiology)
References
Abramson, J., Smirnova, I., Kasho, V., Verner, G., Kaback, H. R., and Iwata, S. (2003). Structure and mechanism of the lactose permease of Escherichia coli. Science 301, 610–615. doi: 10.1126/science.1088196
Alegado, R. A., Chin, C. Y., Monack, D. M., and Tan, M. W. (2011). The two-component sensor kinase KdpD is required for Salmonella typhimurium colonization of Caenorhabditis elegans and survival in macrophages. Cell. Microbiol. 13, 1618–1637. doi: 10.1111/j.1462-5822.2011.01645.x
Al-Eisawi, D. M., and Al-Ruzayza, S. (2015). The flora of holy Mecca district, Saudi Arabia. Int. J. Biodivers. Conserv. 7, 173–189. doi: 10.5897/IJBC2014.0773
Arias, C. A., and Murray, B. E. (2012). The rise of the Enterococcus: Beyond vancomycin resistance. Nat. Rev. Microbiol. 10, 266–278. doi: 10.1038/nrmicro2761
Bäckhed, F., Roswall, J., Peng, Y., Feng, Q., Jia, H., Kovatcheva-Datchary, P., et al. (2015). Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe 17, 690–703. doi: 10.1016/j.chom.2015.04.004
Bahieldin, A., Atef, A., Sabir, J. S., Gadalla, N. O., Edris, S., Alzohairy, A. M., et al. (2015). RNA-Seq analysis of the wild barley (H. spontaneum) leaf transcriptome under salt stress. C. R. Biol. 338, 285–297. doi: 10.1016/j.crvi.2015.03.010
Bai, Y., Muller, D. B., Srinivas, G., Garrido-Oter, R., Potthoff, E., Rott, M., et al. (2015). Functional overlap of the Arabidopsis leaf and root microbiota. Nature 528, 364–369. doi: 10.1038/nature16192
Berendonk, T. U., Manaia, C. M., Merlin, C., Fatta-Kassinos, D., Cytryn, E., Walsh, F., et al. (2015). Tackling antibiotic resistance: The environmental framework. Nat. Rev. Microbiol. 13, 310–317. doi: 10.1038/nrmicro3439
Bialek-Davenet, S., Marcon, E., Leflon-Guibout, V., Lavigne, J.-P., Bert, F., Moreau, R., et al. (2011). In vitro selection of ramR and soxR mutants overexpressing efflux systems by fluoroquinolones as well as cefoxitin in Klebsiella pneumoniae. Antimicrob. Agents Chemother. 55, 2795–2802. doi: 10.1128/AAC.00156-11
Blair, J. M., and Piddock, L. J. (2009). Structure, function and inhibition of RND efflux pumps in Gram-negative bacteria: An update. Curr. Opin. Microbiol. 12, 512–519. doi: 10.1016/j.mib.2009.07.003
Blau, K., Bettermann, A., Jechalke, S., Fornefeld, E., Vanrobaeys, Y., Stalder, T., et al. (2018). The Transferable Resistome of Produce. mBio 9:e01300–e01318. doi: 10.1128/mBio.01300-18
Blau, K., Jacquiod, S., Sorensen, S. J., Su, J. Q., Zhu, Y. G., Smalla, K., et al. (2019). Manure and Doxycycline Affect the Bacterial Community and Its Resistome in Lettuce Rhizosphere and Bulk Soil. Front. Microbiol. 10:725. doi: 10.3389/fmicb.2019.00725
Boucher, H. W., Talbot, G. H., Bradley, J. S., Edwards, J. E., Gilbert, D., Rice, L. B., et al. (2009). Bad bugs, no drugs: No eskape! an update from the infectious diseases society of america. Clin. Infect. Dis. 48, 1–12. doi: 10.1086/595011
Brown, N. L., Stoyanov, J. V., Kidd, S. P., and Hobman, J. L. (2003). The MerR family of transcriptional regulators. FEMS Microbiol. Rev. 27, 145–163. doi: 10.1016/S0168-6445(03)00051-2
Cardona, S. T., Choy, M., and Hogan, A. M. (2018). Essential Two-Component Systems Regulating Cell Envelope Functions: Opportunities for Novel Antibiotic Therapies. J. Membr. Biol. 251, 75–89. doi: 10.1007/s00232-017-9995-5
Chaparro, J. M., Badri, D. V., and Vivanco, J. M. (2014). Rhizosphere microbiome assemblage is affected by plant development. ISME J. 8, 790–803. doi: 10.1038/ismej.2013.196
Checa, S. K., Espariz, M., Audero, M. E., Botta, P. E., Spinelli, S. V., and Soncini, F. C. (2007). Bacterial sensing of and resistance to gold salts. Mol. Microbiol. 63, 1307–1318. doi: 10.1111/j.1365-2958.2007.05590.x
Chen, Q. L., An, X. L., Zhu, Y. G., Su, J. Q., Gillings, M. R., Ye, Z. L., et al. (2017). Application of Struvite Alters the Antibiotic Resistome in Soil, Rhizosphere, and Phyllosphere. Environ. Sci. Technol. 51, 8149–8157. doi: 10.1021/acs.est.7b01420
Chen, Q. L., Cui, H. L., Su, J. Q., Penuelas, J., and Zhu, Y. G. (2019). Antibiotic Resistomes in Plant Microbiomes. Trends Plant Sci. 24, 530–541. doi: 10.1016/j.tplants.2019.02.010
Cheng, G., Hu, Y., Yin, Y., Yang, X., Xiang, C., Wang, B., et al. (2012). Functional screening of antibiotic resistance genes from human gut microbiota reveals a novel gene fusion. FEMS Microbiol. Lett. 336, 11–16. doi: 10.1111/j.1574-6968.2012.02647.x
Chitsaz, M., and Brown, M. H. (2017). The role played by drug efflux pumps in bacterial multidrug resistance. Essays Biochem. 61, 127–139. doi: 10.1042/EBC20160064
Courvalin, P. (2006). Vancomycin resistance in gram-positive cocci. Clin. Infect. Dis. 42:S25–S34. doi: 10.1086/491711
Coyne, S., Rosenfeld, N., Lambert, T., Courvalin, P., and Perichon, B. (2010). Overexpression of resistance-nodulation-cell division pump AdeFGH confers multidrug resistance in Acinetobacter baumannii. Antimicrob. Agents Chemother. 54, 4389–4393. doi: 10.1128/AAC.00155-10
de Lima Procópio, R. E., da Silva, I. R., Martins, M. K., de Azevedo, J. L., and de Araújo, J. M. (2012). Antibiotics produced by Streptomyces. Braz. J. Infect. Dis. 16, 466–471. doi: 10.1016/j.bjid.2012.08.014
Dean, M., Hamon, Y., and Chimini, G. (2001). The human ATP-binding cassette (ABC) transporter superfamily. J. Lipid Res. 42, 1007–1017. doi: 10.1016/S0022-2275(20)31588-1
Djordjevic, S. P., Stokes, H. W., and Chowdhury, P. R. (2013). Mobile elements, zoonotic pathogens and commensal bacteria: Conduits for the delivery of resistance genes into humans, production animals and soil microbiota. Front. Microbiol. 4:86. doi: 10.3389/fmicb.2013.00086
Du, D., Wang-Kan, X., Neuberger, A., Van Veen, H. W., Pos, K. M., Piddock, L. J., et al. (2018). Multidrug efflux pumps: Structure, function and regulation. Nat. Rev. Microbiol. 16, 523–539. doi: 10.1038/s41579-018-0048-6
Dulyayangkul, P., Charoenlap, N., Srijaruskul, K., Mongkolsuk, S., and Vattanaviboon, P. (2016). Major facilitator superfamily MfsA contributes to multidrug resistance in emerging nosocomial pathogen Stenotrophomonas maltophilia. J. Antimicrob. Chemother. 71, 2990–2991. doi: 10.1093/jac/dkw233
Dutescu, I. A., and Hillier, S. A. (2021). Encouraging the Development of New Antibiotics: Are Financial Incentives the Right Way Forward? A Systematic Review and Case Study. Infect. Drug Resist. 14, 415–434. doi: 10.2147/IDR.S287792
Eswaran, J., Koronakis, E., Higgins, M. K., Hughes, C., and Koronakis, V. (2004). Three’s company: Component structures bring a closer view of tripartite drug efflux pumps. Curr. Opin. Struct. Biol. 14, 741–747. doi: 10.1016/j.sbi.2004.10.003
Fath, M. J., and Kolter, R. (1993). ABC transporters: Bacterial exporters. Microbiol. Rev. 57, 995–1017. doi: 10.1128/mr.57.4.995-1017.1993
Fernández, L., and Hancock, R. E. (2012). Adaptive and mutational resistance: Role of porins and efflux pumps in drug resistance. Clin. Microbiol. Rev. 25, 661–681. doi: 10.1128/CMR.00043-12
Forsberg, K. J., Patel, S., Gibson, M. K., Lauber, C. L., Knight, R., Fierer, N., et al. (2014). Bacterial phylogeny structures soil resistomes across habitats. Nature 509, 612–616. doi: 10.1038/nature13377
Forti, F., Mauri, V., Deho, G., and Ghisotti, D. (2011). Isolation of conditional expression mutants in Mycobacterium tuberculosis by transposon mutagenesis. Tuberculosis 91, 569–578. doi: 10.1016/j.tube.2011.07.004
Franklin, R. (2016). “Rhodococcus equi Pneumonia in Foals: An Update on Epidemiology, Diagnosis, Treatment and Prevention,” in Proceedings of the American Association of Equine Practitioners, Lexington: American Association of Equine Practitioners.
Freeman, Z. N., Dorus, S., and Waterfield, N. R. (2013). The KdpD/KdpE two-component system: Integrating K(+) homeostasis and virulence. PLoS Pathog. 9:e1003201. doi: 10.1371/journal.ppat.1003201
Fu, L., Niu, B., Zhu, Z., Wu, S., and Li, W. (2012). CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. doi: 10.1093/bioinformatics/bts565
Gopalakrishnan, L., Doriya, K., and Kumar, D. S. (2016). Moringa oleifera: A review on nutritive importance and its medicinal application. Food Sci. Hum. Wellness 5, 49–56. doi: 10.1016/j.fshw.2016.04.001
Gu, M., and Imlay, J. A. (2011). The SoxRS response of Escherichia coli is directly activated by redox-cycling drugs rather than by superoxide. Mol. Microbiol. 79, 1136–1150. doi: 10.1111/j.1365-2958.2010.07520.x
Gudeta, D. D., Moodley, A., Bortolaia, V., and Guardabassi, L. (2014). van0 a new glycopeptide resistance operon in environmental Rhodococcus equi isolates. Antimicrob. Agents Chemother. 58, 1768–1770. doi: 10.1128/AAC.01880-13
Handing, J. W., Ragland, S. A., Bharathan, U. V., and Criss, A. K. (2018). The MtrCDE efflux pump contributes to survival of Neisseria gonorrhoeae from human neutrophils and their antimicrobial components. Front. Microbiol. 9:2688. doi: 10.3389/fmicb.2018.02688
Hirakawa, H., Nishino, K., Hirata, T., and Yamaguchi, A. (2003). Comprehensive studies of drug resistance mediated by overexpression of response regulators of two-component signal transduction systems in Escherichia coli. J. Bacteriol. 185, 1851–1856. doi: 10.1128/JB.185.6.1851-1856.2003
Hu, Y., Morichaud, Z., Chen, S., Leonetti, J.-P., and Brodolin, K. (2012). Mycobacterium tuberculosis RbpA protein is a new type of transcriptional activator that stabilizes the σ A-containing RNA polymerase holoenzyme. Nucleic Acids Res. 40, 6547–6557. doi: 10.1093/nar/gks346
Hughes, D. T., Clarke, M. B., Yamamoto, K., Rasko, D. A., and Sperandio, V. (2009). The QseC adrenergic signaling cascade in Enterohemorrhagic E. coli (EHEC). PLoS Pathog. 5:e1000553. doi: 10.1371/journal.ppat.1000553
Hurt, R. A., Qiu, X., Wu, L., Roh, Y., Palumbo, A. V., Tiedje, J. M., et al. (2001). Simultaneous recovery of RNA and DNA from soils and sediments. Appl. Environ. Microbiol. 67, 4495–4503. doi: 10.1128/AEM.67.10.4495-4503.2001
Huson, D. H., Beier, S., Flade, I., Górska, A., El-Hadidi, M., Mitra, S., et al. (2016). MEGAN community edition-interactive exploration and analysis of large-scale microbiome sequencing data. PLoS Comput. Biol. 12:e1004957. doi: 10.1371/journal.pcbi.1004957
Huson, D. H., Mitra, S., Ruscheweyh, H.-J., Weber, N., and Schuster, S. C. (2011). Integrative analysis of environmental sequences using MEGAN4. Genome Res. 21, 1552–1560. doi: 10.1101/gr.120618.111
Ishikawa, J., Chiba, K., Kurita, H., and Satoh, H. (2006). Contribution of rpoB2 RNA polymerase beta subunit gene to rifampin resistance in Nocardia species. Antimicrob. Agents Chemother. 50, 1342–1346. doi: 10.1128/AAC.50.4.1342-1346.2006
Kapadia, M., Rolston, K. V., and Han, X. Y. (2007). Invasive Streptomyces infections: Six cases and literature review. Am. J. Clin. Pathol. 127, 619–624. doi: 10.1309/QJEBXP0BCGR54L15
Karlsson, F. H., Fåk, F., Nookaew, I., Tremaroli, V., Fagerberg, B., Petranovic, D., et al. (2012). Symptomatic atherosclerosis is associated with an altered gut metagenome. Nat. Commun. 3:1245. doi: 10.1038/ncomms2266
Khan, T. M., Kok, Y. L., Bukhsh, A., Lee, L. H., Chan, K. G., and Goh, B. H. (2018). Incidence of methicillin resistant Staphylococcus aureus (MRSA) in burn intensive care unit: A systematic review. Germs 8, 113–125. doi: 10.18683/germs.2018.1138
Koprivnjak, T., and Peschel, A. (2011). Bacterial resistance mechanisms against host defense peptides. Cell. Mol. Life Sci. 68, 2243–2254. doi: 10.1007/s00018-011-0716-4
Laarman, A. J., Ruyken, M., Malone, C. L., van Strijp, J. A., Horswill, A. R., and Rooijakkers, S. H. (2011). Staphylococcus aureus metalloprotease aureolysin cleaves complement C3 to mediate immune evasion. J. Immunol. 186, 6445–6453. doi: 10.4049/jimmunol.1002948
Lee, L.-F., Chen, Y.-J., Kirby, R., Chen, C., and Chen, C. W. (2007). A multidrug efflux system is involved in colony growth in Streptomyces lividans. Microbiology 153, 924–934. doi: 10.1099/mic.0.2006/000018-0
Levy-Booth, D. J., Campbell, R. G., Gulden, R. H., Hart, M. M., Powell, J. R., Klironomos, J. N., et al. (2007). Cycling of extracellular DNA in the soil environment. Soil Biol. Biochem. 39, 2977–2991. doi: 10.1016/j.soilbio.2007.06.020
Li, H., Wang, Q., Wang, R., Zhang, Y., Wang, X., and Wang, H. (2017). Global regulator SoxR is a negative regulator of efflux pump gene expression and affects antibiotic resistance and fitness in Acinetobacter baumannii. Medicine 96:e7188. doi: 10.1097/MD.0000000000007188
Li, J., Jia, H., Cai, X., Zhong, H., Feng, Q., Sunagawa, S., et al. (2014). An integrated catalog of reference genes in the human gut microbiome. Nat. Biotechnol. 32, 834–841. doi: 10.1038/nbt.2942
Li, S. M., and Heide, L. (2006). The biosynthetic gene clusters of aminocoumarin antibiotics. Planta Med. 72, 1093–1099. doi: 10.1055/s-2006-946699
Li, W., and Godzik, A. (2006). Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659. doi: 10.1093/bioinformatics/btl158
Li, X. Z., and Nikaido, H. (2009). Efflux-mediated drug resistance in bacteria: An update. Drugs 69, 1555–1623. doi: 10.2165/11317030-000000000-00000
Liu, B., and Pop, M. (2009). ARDB antibiotic resistance genes database. Nucleic Acids Res. 37:D443–D447. doi: 10.1093/nar/gkn656
Lucas, C. E., Balthazar, J. T., Hagman, K. E., and Shafer, W. M. (1997). The MtrR repressor binds the DNA sequence between the mtrR and mtrC genes of Neisseria gonorrhoeae. J. Bacteriol. 179, 4123–4128. doi: 10.1128/jb.179.13.4123-4128.1997
Lv, F., Cai, J., He, Q., Wang, W., Luo, Y., Wang, X., et al. (2021). Overexpression of Efflux Pumps Mediate Pan Resistance of Klebsiella pneumoniae Sequence Type 11. Microb. Drug Resist. 27, 1405–1411. doi: 10.1089/mdr.2020.0395
Ma Rodriguez, A., Olano, C., Vilches, C., Méndez, C., and Salas, J. A. (1993). Streptomyces antibioticus contains at least three oleandomycin-resistance determinants, one of which shows similarity with proteins of the ABC-transporter superfamily. Mol. Microbiol. 8, 571–582. doi: 10.1111/j.1365-2958.1993.tb01601.x
Ma, D., Cook, D. N., Hearst, J. E., and Nikaido, H. (1994). Efflux pumps and drug resistance in gram-negative bacteria. Trends Microbiol. 2, 489–493. doi: 10.1016/0966-842X(94)90654-8
Maeda, H., Fujita, N., and Ishihama, A. (2000). Competition among seven Escherichia coli σ subunits: Relative binding affinities to the core RNA polymerase. Nucleic Acids Res. 28, 3497–3503. doi: 10.1093/nar/28.18.3497
Marchand, I., Damier-Piolle, L., Courvalin, P., and Lambert, T. (2004). Expression of the RND-type efflux pump AdeABC in Acinetobacter baumannii is regulated by the AdeRS two-component system. Antimicrob. Agents Chemother. 48, 3298–3304. doi: 10.1128/AAC.48.9.3298-3304.2004
Marger, M. D., and Saier Jr, M. H. (1993). A major superfamily of transmembrane facilitators that catalyse uniport, symport and antiport. Trends Biochem. Sci. 18, 13–20. doi: 10.1016/0968-0004(93)90081-W
Martínez, J. L., Coque, T. M., and Baquero, F. (2015). What is a resistance gene? Ranking risk in resistomes. Nat. Rev. Microbiol. 13, 116–123. doi: 10.1038/nrmicro3399
Maxwell, A. (1997). DNA gyrase as a drug target. Trends Microbiol. 5, 102–109. doi: 10.1016/S0966-842X(96)10085-8
Mende, D. R., Waller, A. S., Sunagawa, S., Järvelin, A. I., Chan, M. M., Arumugam, M., et al. (2012). Assessment of metagenomic assembly using simulated next generation sequencing data. PLoS One 7:e31386. doi: 10.1371/journal.pone.0031386
Mendes, L. W., de Chaves, M. G., Fonseca, M. C., Mendes, R., Raaijmakers, J. M., and Tsai, S. M. (2019). Resistance Breeding of Common Bean Shapes the Physiology of the Rhizosphere Microbiome. Front. Microbiol. 10:2252. doi: 10.3389/fmicb.2019.02252
Méndez, C., and Salas, J. A. (2001). The role of ABC transporters in antibiotic-producing organisms: Drug secretion and resistance mechanisms. Res. Microbiol. 152, 341–350. doi: 10.1016/S0923-2508(01)01205-0
Mönchgesang, S., Strehmel, N., Schmidt, S., Westphal, L., Taruttis, F., Müller, E., et al. (2016). Natural variation of root exudates in Arabidopsis thaliana-linking metabolomic and genomic data. Sci. Rep. 6:29033. doi: 10.1038/srep29033
Munita, J. M., and Arias, C. A. (2016). Mechanisms of Antibiotic Resistance. Microbiol. Spectr. 4. doi: 10.1128/microbiolspec.VMBF-0016-2015
Murakami, S., Nakashima, R., Yamashita, E., Matsumoto, T., and Yamaguchi, A. (2006). Crystal structures of a multidrug transporter reveal a functionally rotating mechanism. Nature 443, 173–179. doi: 10.1038/nature05076
Nardelli, M., Scalzo, P. M., Ramirez, M. S., Quiroga, M. P., Cassini, M. H., and Centron, D. (2012). Class 1 integrons in environments with different degrees of urbanization. PLoS One 7:e39223. doi: 10.1371/journal.pone.0039223
Nesme, J., and Simonet, P. (2015). The soil resistome: A critical review on antibiotic resistance origins, ecology and dissemination potential in telluric bacteria. Environ. Microbiol. 17, 913–930. doi: 10.1111/1462-2920.12631
Newell, K. V., Thomas, D. P., Brekasis, D., and Paget, M. S. (2006). The RNA polymerase-binding protein RbpA confers basal levels of rifampicin resistance on Streptomyces coelicolor. Mol Microbiol. 60, 687–696. doi: 10.1111/j.1365-2958.2006.05116.x
Nielsen, H. B., Almeida, M., Juncker, A. S., Rasmussen, S., Li, J., Sunagawa, S., et al. (2014). Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat. Biotechnol. 32, 822–828. doi: 10.1038/nbt.2939
Nikaido, H. (1996). Multidrug efflux pumps of gram-negative bacteria. J. Bacteriol. 178, 5853–5859. doi: 10.1128/jb.178.20.5853-5859.1996
Nudel, K., McClure, R., Moreau, M., Briars, E., Abrams, A. J., Tjaden, B., et al. (2018). Transcriptome analysis of Neisseria gonorrhoeae during natural infection reveals differential expression of antibiotic resistance determinants between men and women. Msphere 3:e00312–e00318. doi: 10.1128/mSphereDirect.00312-18
Obermeier, M. M., Wicaksono, W. A., Taffner, J., Bergna, A., Poehlein, A., Cernava, T., et al. (2021). Plant resistome profiling in evolutionary old bog vegetation provides new clues to understand emergence of multi-resistance. ISME J. 15, 921–937. doi: 10.1038/s41396-020-00822-9
Oh, J., Byrd, A. L., Deming, C., Conlan, S., Program, N. C. S., Kong, H. H., et al. (2014). Biogeography and individuality shape function in the human skin metagenome. Nature 514, 59–64. doi: 10.1038/nature13786
Ohneck, E. A., Zalucki, Y. M., Johnson, P. J., Dhulipala, V., Golparian, D., Unemo, M., et al. (2011). A novel mechanism of high-level, broad-spectrum antibiotic resistance caused by a single base pair change in Neisseria gonorrhoeae. mBio 2:e00187–11. doi: 10.1128/mBio.00187-11
Olano, C., Rodriguez, A. M., Mendez, C., and Salas, J. A. (1996). Topological studies of the membrane component of the OleC ABC transporter involved in oleandomycin resistance in Streptomyces antibioticus. FEMS Microbiol. Lett. 143, 133–139. doi: 10.1111/j.1574-6968.1996.tb08472.x
Olesky, M., Zhao, S., Rosenberg, R. L., and Nicholas, R. A. (2006). Porin-mediated antibiotic resistance in Neisseria gonorrhoeae: Ion, solute, and antibiotic permeation through PIB proteins with penB mutations. J. Bacteriol. 188, 2300–2308. doi: 10.1128/JB.188.7.2300-2308.2006
Olson, M., and Flora of North America Editorial Committee (2010). eFlora summary: Moringaceae: Drumstick family. Flora North Am. 7, 167–169.
Paget, M. S., Molle, V., Cohen, G., Aharonowitz, Y., and Buttner, M. J. (2001). Defining the disulphide stress response in Streptomyces coelicolor A3 (2): Identification of the σR regulon. Mol. Microbiol. 42, 1007–1020. doi: 10.1046/j.1365-2958.2001.02675.x
Palma, M., Zurita, J., Ferreras, J. A., Worgall, S., Larone, D. H., Shi, L., et al. (2005). Pseudomonas aeruginosa SoxR does not conform to the archetypal paradigm for SoxR-dependent regulation of the bacterial oxidative stress adaptive response. Infect. Immun. 73, 2958–2966. doi: 10.1128/IAI.73.5.2958-2966.2005
Pao, S. S., Paulsen, I. T., and Saier, M. H. Jr. (1998). Major facilitator superfamily. Microbiol. Mol. Biol. Rev. 62, 1–34. doi: 10.1128/MMBR.62.1.1-34.1998
Peleg, A. Y., Seifert, H., and Paterson, D. L. (2008). Acinetobacter baumannii: Emergence of a successful pathogen. Clin. Microbiol. Rev. 21, 538–582. doi: 10.1128/CMR.00058-07
Pérez Audero, M. E., Podoroska, B. M., Ibáñez, M. M., Cauerhff, A., Checa, S. K., and Soncini, F. C. (2010). Target transcription binding sites differentiate two groups of MerR-monovalent metal ion sensors. Mol. Microbiol. 78, 853–865. doi: 10.1111/j.1365-2958.2010.07370.x
Piddock, L. J. (2006). Multidrug-resistance efflux pumps? not just for resistance. Nat. Rev. Microbiol. 4, 629–636. doi: 10.1038/nrmicro1464
Pomposiello, P. J., and Demple, B. (2001). Redox-operated genetic switches: The SoxR and OxyR transcription factors. Trends Biotechnol. 19, 109–114. doi: 10.1016/S0167-7799(00)01542-0
Pontel, L. B., Audero, M. E. P., Espariz, M., Checa, S. K., and Soncini, F. C. (2007). GolS controls the response to gold by the hierarchical induction of Salmonella-specific genes that include a CBA efflux-coding operon. Mol. Microbiol. 66, 814–825. doi: 10.1111/j.1365-2958.2007.05963.x
Poole, K. (2007). Efflux pumps as antimicrobial resistance mechanisms. Ann. Med. 39, 162–176. doi: 10.1080/07853890701195262
Raes, J., Foerstner, K. U., and Bork, P. (2007). Get the most out of your metagenome: Computational analysis of environmental sequence data. Curr. Opin. Microbiol. 10, 490–498. doi: 10.1016/j.mib.2007.09.001
Reynolds, P. E. (1989). Structure, biochemistry and mechanism of action of glycopeptide antibiotics. Eur. J. Clin. Microbiol. Infect. Dis. 8, 943–950. doi: 10.1007/BF01967563
Rouquette, C., Harmon, J. B., and Shafer, W. M. (1999). Induction of the mtrCDE-encoded efflux pump system of Neisseria gonorrhoeae requires MtrA, an AraC-like protein. Mol. Microbiol. 33, 651–658. doi: 10.1046/j.1365-2958.1999.01517.x
Ruiz, C., and Levy, S. B. (2014). Regulation of acrAB expression by cellular metabolites in Escherichia coli. J. Antimicrob. Chemother. 69, 390–399. doi: 10.1093/jac/dkt352
Saidijam, M., Benedetti, G., Ren, Q., Xu, Z., Hoyle, C. J., Palmer, S. L., et al. (2006). Microbial drug efflux proteins of the major facilitator superfamily. Curr. Drug Targets 7, 793–811. doi: 10.2174/138945006777709575
Saier, M. H. Jr., Beatty, J. T., Goffeau, A., Harley, K. T., Heijne, W., Huang, S.-C., et al. (1999). The major facilitator superfamily. J. Mol. Microbiol. Biotechnol. 1, 257–279.
Sakhtah, H., Koyama, L., Zhang, Y., Morales, D. K., Fields, B. L., Price-Whelan, A., et al. (2016). The Pseudomonas aeruginosa efflux pump MexGHI-OpmD transports a natural phenazine that controls gene expression and biofilm development. Proc. Natl. Acad. Sci.U.S.A. 113:E3538–E3547. doi: 10.1073/pnas.1600424113
Salas, J., Hernandez, C., Mendez, C., Olano, C., Quiros, L., Rodriguez, A., et al. (1994). Intracellular glycosylation and active efflux as mechanisms for resistance to oleandomycin in Streptomyces antibioticus, the producer organism. Microbiologia 10, 37–48.
Sasse, J., Martinoia, E., and Northen, T. (2018). Feed your friends: Do plant exudates shape the root microbiome? Trends Plant Sci. 23, 25–41. doi: 10.1016/j.tplants.2017.09.003
Sauna, Z. E., Ambudkar, S. V., Pradines, V. P., Orlandi-Pradines, E., Sanwal, S. L. P., Prasad, R., et al. (2009). ABC Transporters in Microorganisms. Norfolk: Caister Academic Press.
Schmutz, E., Hennig, S., Li, S. M., and Heide, L. (2004). Identification of a topoisomerase IV in actinobacteria: Purification and characterization of ParYR and GyrBR from the coumermycin A1 producer Streptomyces rishiriensis DSM 40489. Microbiology 150, 641–647. doi: 10.1099/mic.0.26867-0
Schmutz, E., Muhlenweg, A., Li, S. M., and Heide, L. (2003). Resistance genes of aminocoumarin producers: Two type II topoisomerase genes confer resistance against coumermycin A1 and clorobiocin. Antimicrob. Agents Chemother. 47, 869–877. doi: 10.1128/AAC.47.3.869-877.2003
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Sun, J., Deng, Z., and Yan, A. (2014). Bacterial multidrug efflux pumps: Mechanisms, physiology and pharmacological exploitations. Biochem. Biophys. Res. Commun. 453, 254–267. doi: 10.1016/j.bbrc.2014.05.090
Thiara, A. S., and Cundliffe, E. (1989). Interplay of novobiocin-resistant and -sensitive DNA gyrase activities in self-protection of the novobiocin producer, Streptomyces sphaeroides. Gene 81, 65–72. doi: 10.1016/0378-1119(89)90337-5
Tierney, A. R., and Rather, P. N. (2019). Roles of two-component regulatory systems in antibiotic resistance. Future Microbiol. 14, 533–552. doi: 10.2217/fmb-2019-0002
Tringe, S. G., von Mering, C., Kobayashi, A., Salamov, A. A., Chen, K., Chang, H. W., et al. (2005). Comparative metagenomics of microbial communities. Science 308, 554–557. doi: 10.1126/science.1107851
Tseng, T.-T., Gratwick, K. S., Kollman, J., Park, D., Nies, D. H., and Goffeau, A. (1999). The RND permease superfamily: An ancient, ubiquitous and diverse family that includes human disease and development proteins. J. Mol. Microbiol. Biotechnol. 1, 107–125.
Vorholt, J. A. (2012). Microbial life in the phyllosphere. Nat. Rev. Microbiol. 10, 828–840. doi: 10.1038/nrmicro2910
Vorholt, J. A., Vogel, C., Carlstrom, C. I., and Muller, D. B. (2017). Establishing Causality: Opportunities of Synthetic Communities for Plant Microbiome Research. Cell Host Microbe 22, 142–155. doi: 10.1016/j.chom.2017.07.004
White, D. G., Goldman, J. D., Demple, B., and Levy, S. B. (1997). Role of the acrAB locus in organic solvent tolerance mediated by expression of marA, soxS, or robA in Escherichia coli. J. Bacteriol. 179, 6122–6126. doi: 10.1128/jb.179.19.6122-6126.1997
Wi, T., Lahra, M. M., Ndowa, F., Bala, M., Dillon, J. R., Ramon-Pardo, P., et al. (2017). Antimicrobial resistance in Neisseria gonorrhoeae: Global surveillance and a call for international collaborative action. PLoS Med. 14:e1002344. doi: 10.1371/journal.pmed.1002344
Wilson, D. N., Hauryliuk, V., Atkinson, G. C., and O’Neill, A. J. (2020). Target protection as a key antibiotic resistance mechanism. Nat. Rev. Microbiol. 18, 637–648. doi: 10.1038/s41579-020-0386-z
World Health Organization [WHO] (2012). Fact Sheet No. 194: Antimicrobial Resistance. Geneva: World Health Organization.
Xue, T., You, Y., Hong, D., Sun, H., and Sun, B. (2011). The Staphylococcus aureus KdpDE two-component system couples extracellular K+ sensing and Agr signaling to infection programming. Infect. Immun. 79, 2154–2167. doi: 10.1128/IAI.01180-10
Keywords: resistome, horizontal gene transfer, antibiotic target, rhizosphere, microbiome, efflux pump, two-component system
Citation: Shami AY, Abulfaraj AA, Refai MY, Barqawi AA, Binothman N, Tashkandi MA, Baeissa HM, Baz L, Abuauf HW, Ashy RA and Jalal RS (2022) Abundant antibiotic resistance genes in rhizobiome of the human edible Moringa oleifera medicinal plant. Front. Microbiol. 13:990169. doi: 10.3389/fmicb.2022.990169
Received: 12 July 2022; Accepted: 17 August 2022;
Published: 15 September 2022.
Edited by:
Muhammad Waseem, Hainan University, ChinaReviewed by:
Antt Htet Wai, Yangon University of Education (YUOE), MyanmarMehtab Muhammad Aslam, The Chinese University of Hong Kong, Hong Kong SAR, China
Copyright © 2022 Shami, Abulfaraj, Refai, Barqawi, Binothman, Tashkandi, Baeissa, Baz, Abuauf, Ashy and Jalal. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rewaa S. Jalal, Rsjalal@uj.edu.sa