7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
View/Download PDF

Translate this page into:

Original article
12 2023
:35;
102954
doi:
10.1016/j.jksus.2023.102954

Deciphering the genetics of antagonism and antimicrobial resistance in Bacillus velezensis HU-91 by whole genome analysis

Department of Biosciences, COMSATS University Islamabad, Park Road, Islamabad, Pakistan
Department of Agronomy, Institute of Food and Agriculture Sciences, University of Florida, Gainesville, FL 32611, USA
Department of Clinical Pharmacy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia

⁎Corresponding author. nadeem_hassan@comsats.edu.pk (Muhammad Nadeem Hassan)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Peer review under responsibility of King Saud University

Abstract

Bacillus species are well reported plant growth promoting rhizobacteria which play an important role in improving soil fertility, nutrient recycling, and controlling phytopathogens. Their use as biocontrol agents is an ecofriendly strategy but being a rhizospheric community, they could harbor antimicrobial resistance genes. In current study, rhizobacteria were isolated from the mono cropped wheat. Ten out of 50 isolates showed antifungal activity against Fusarium spp., economically important pathogens of wheat. The strain HU-91 with highest antifungal activity showed the biochemical traits of lipopeptide production. It also resisted multiple antibiotics of different classes (index ≥ 0.2) and was categorized as multidrug resistant bacteria. The strain HU-91 was identified as Bacillus velezensis based on its 16S rRNA homology (99 %). Whole genome analysis of B. velezensis HU-91 revealed 76 contigs with a total length of 3,952,697 bp, GC content 46.3 %, and a total of 3912 genes. Functional classification of genome predicted genes involved in drug resistance, and antagonism. A detailed genome analysis identified CAZymes and 10 biosynthetic gene clusters involved in antagonism. Various antimicrobial resistance genes associated with the resistance to different antimicrobial classes such as β lactams, tetracyclines, aminoglycosides, fluoroquinolones, macrolides, sulfonamides, and lincosamides were also found in the genome of HU-91. Furthermore, HU-91 lacked plasmid but harbored a well-developed CRISPR Cas system. The underlying mechanisms involved in the antimicrobial resistance were depicted as intrinsic based on the genes encoding production of antibiotics degrading enzymes and efflux pumps. These findings depict that HU-91 adapts the intrinsic mechanisms of drug resistivity to tolerate environmental stress. The lack of plasmid and presence of CRISPR Cas system advocate its safety to be used as bioinoculant.

Keywords

Antimicrobial resistance
Antagonism
CRISPR
1

1 Introduction

Phytopathogens, that cause economically important diseases are a widely reported problem in agriculture as they pose deleterious effects on crops (Anand et al., 2023). Fusarium spp. are the most devastating phytopathogen for cereal crops especially wheat grown in arid and semi-arid regions where wheat is mono cropped to meet the needs of growing population (Ullah et al., 2020). Different strategies have been employed to control the deleterious effects of pathogens on crops.

Use of biocontrol agents (Bacillus spp., Pseudomonas spp., etc.) is an interesting recourse to inhibit the growth of phytopathogens. Biocontrol agents adapt different strategies to cope with biotic and abiotic stresses such as production of antibiotics, siderophores, cell wall degrading enzymes, and volatile organic compounds (Zhang et al., 2021). These secondary metabolites not only provide interference competition but usually contribute to the self-resistance against the antimicrobial metabolites especially antibiotics (Xia et al., 2022). This self-protection of antibiotic producers is referred to as innate resistance while the appearance of self-protection in non-antibiotic producers is referred to as acquired resistance. Both types of bacteria, the producers, and non-producers of antibiotics, can achieve innate as well as acquired resistance. Thus, the secondary metabolites producing microbes have been an extensive source of antibiotic resistance genes (Wencewicz, 2019).

The antibiotic resistance genes of bacteria could transfer to other microbes through horizontal gene transfer (HGT) which usually involves mobile genetic elements like transposons, plasmid and integrons. HGT occurs through transformation, transduction, and conjugation (Vinayamohan et al., 2022). Based on resistance pattern, bacteria can be multidrug resistance (MDR), extensive drug resistance (XDR) and pan drug resistance (PDR) (Magiorakos et al., 2012).

The genetic mechanisms adapted by bacteria to resist antimicrobial agents includes inhibition of antibiotics, target modification, enzymatic inactivation of antibiotics, efflux pumps and cell wall permeability. Many of the biocontrol agents such as Pseudomonas spp. and Bacillus spp. are characterized as MDR (Fahsi et al., 2019).

Besides the natural producers of antibiotics, human anthropogenic activities such as use of pesticides, and antibiotics could also promote antimicrobial resistance. The presence of these compounds provides a selective pressure on dispersal of antimicrobial resistance genes (ARGs) (Mazhar et al., 2021). The use of bioinoculants in agriculture are usually considered as environment friendly but they could be a potential risk for dissemination of ARGs in the environment. It is therefore necessary to ensure the biosafety, efficacy, and quality of the bio products, especially their antimicrobial resistance profile.

Many approaches like amplification of genes, Sanger sequencing, and metagenomics, have been used to detect the genetics, diversity, and evolution of bacteria. However, the genetics of antimicrobial resistivity and antagonism had not been deciphered yet. The next generation sequencing is a powerful approach for deciphering the underlying genetics of different mechanisms in bacteria through whole genome analysis (WGA). Several studies have mentioned the significance of WGA in prediction of ARGs and antagonism (Wang et al., 2023). WGS has a higher throughput, more flexible and has greater sensitivity over Sanger sequencing. This approach provides new insight into the genetics of unique genes and correlation with their phenotypes (Argimón et al., 2020). The present study aims to isolate bacteria antagonistic to the economically important pathogens (Fusarium spp.), from mono cropping system of wheat, elucidate its antimicrobial resistance phenotype and decipher the underlying genetics of its antagonism and antimicrobial resistance by whole genome analysis.

2

2 Materials and methods

2.1

2.1 Field site and sample collection

A mono cropped field under the cultivation of wheat for 6–7 consecutive years was identified in a rainfed area with subtropic semi-arid conditions of District Chakwal, Pakistan (32.8322° N, 72.6151° E). Plants (n = 5) were sampled at the pre-harvesting stage. Healthy plants showing variable height, color, size, and the number of grains in each panicle were uprooted randomly with bulk rhizospheric soil from different places (corner and center with 2 m distance between plants) of the field. Samples were placed in a sterilized polythene zip bag (60 × 90 cm, hole at one end for aeriation), tagged and stored at 4°C for further processing.

2.2

2.2 Isolation of bacteria

Bacteria were isolated from rhizosphere and endosphere by serial dilution method. The rhizospheric soil tightly adhering to different plant roots was gently removed, whereas for the endospheric sample, roots were pooled to make a composite sample. Root samples were surface sterilized with 0.1 % mercuric chloride (HgCl2) for 1–2 min and washed 3 times with sterilized distilled water. The plant roots were dried with sterile filter paper and crushed in sterile pestle mortar. A serial dilution (10−1 – 10−9) was made from 1 g of rhizospheric soil or crushed plant roots and spread on sterile Luria Bertani agar plates. The plates were incubated until the appearance of the colonies (Ullah et al., 2020).

2.3

2.3 Antagonistic activity

The selected isolates were tested for their antifungal activity against economically important pathogens such as F. oxysporum, F. graminearum, F. solani and F. moniliforme using a dual culture plate method. The fungal disc (5 mm) was placed at the center and the bacterial strain was streaked 4 cm away from the fungal disc. The antagonistic activity was observed after 7 days of inoculation and measured using the following formula. B. velezensis SM39 (Mulk et al., 2022) was used as a positive control. Cell free broth was used as a negative control. I n h i b i t i o n % = C - T / C × 100 ( C = m y c e l i u m g r o w t h under c o n t r o l ; T = m y c e l i u m g r o w t h uunder t r e a t m e n t )

2.4

2.4 Biochemical detection of lipopeptides

Lipopeptides were detected in HU-91 by centrifugation (20 min, 10,000 rpm) of overnight grown culture. The supernatant was used for different biochemical analyses. Bacillus cabrialesii SM 93 (unpublished data) was used as a positive control.

2.4.1

2.4.1 Hydrophobicity assay (BATH)

The harvested cells were washed twice with a 50 mM potassium phosphate buffer and were re-suspended in the same buffer to optical density 600 nm (0.51). The cell suspension was mixed with tested hydrocarbon (p-xylene), and vortexed (3 min). The mixture was allowed to stand for 45 min and the aqueous phase was subjected to spectrophotometer to measure the OD at 610 nm. Cell hydrophobicity was measured by the following formula (Rosenberg et al., 1980). B A T H % = B 0 - B / B 0 × 100

Where B0 and B = Optical density before and after addition of hydrocarbons respectively.

2.4.2

2.4.2 Oil spread assay

For oil spread assay, 10 µL of bacterial suspension, 20 µL crude oil and 25 mL distilled water was mixed in a petri plate. Appearance of a hallow zone indicated the presence of oil displacement (Zargar et al., 2022).

2.4.3

2.4.3 Emulsification assay

For emulsification assay, cell free supernatant was mixed with crude oil. The suspension was vortexed and incubated at room temperature for 24 h. The EA was measured by the following formula (Cooper & Goldenberg, 1987). E A = [ H e i g h t o f t h e e m u l s i f i c a t i o n l a y e r / T o t a l h e i g h t o f l i q u i d ( c m ) ] × 100

2.5

2.5 Antimicrobial susceptibility testing

The overnight-grown culture of HU-91 was subjected to antimicrobial susceptibility testing by using a disc diffusion assay following Clinical Standards Laboratory Institute (CLSI) and European Committee on Antimicrobial Susceptibility Testing (EUCAST) guidelines (Wash et al., 2022). Briefly, the antimicrobial discs were placed on the lawn of strain on agar plates and incubated overnight at 37°C. The inhibition zone was measured, and the strain was categorized as MDR, XDR or PDR. Multiple antimicrobial resistance index (MAR) was calculated by using the following formula (Dela Peña et al., 2022). M A R i n d e x = a / b ( a = n u m b e r o f a n t i m i c r o b i a l a g e n t s a g a i n s t w h i c h H U - 91 s h o w s r e s i s t a n c e ; b = t o t a l n u m b e r o f a n t i m i c r o b i a l a g e n t s t e s t e d ) .

2.6

2.6 Whole genome sequencing and analysis

Genomic DNA was isolated from the fresh bacterial culture by DNA extraction kit (Thermo Fisher Scientific, USA). Quantification was done by Nanodrop spectrophotometer. The genome of bacteria was sequenced commercially by Azenta life sciences USA C/O Alpha Genomics (Pvt) Ltd, Islamabad, Pakistan. Briefly, sequencing was performed using Illumina HiSeq 4000 platform. The quality of raw reads were filtered by FastQC version 0.11.9 (Chandra et al., 2021). De novo assembly of the FastQC filtered genome was done by Velevet assembler version 1.2.10 using Bacillus velezensis FZB42 as a reference strain.

2.6.1

2.6.1 Gene prediction and annotation

The genome was annotated using PROKA which is used for the rapid annotation of genome using Circular Genome Viewer (CG viewer). The comparison of the sequences of genome with neighbor strain was also uploaded by this tool. The protein coding genes were predicted by Rapid Annotation Subsystem Technology (RAST) which annotates the genome built upon SEED. This system identified protein-encoding, rRNA and tRNA genes, assigned specific functions to each gene and categorized them into subsystems. The genome was constructed in the form of a circular map by using CG viewer tool using default parameters (Narayanasamy et al., 2023).

2.6.2

2.6.2 Taxonomic classification of HU-91

The average nucleotide identity (ANI) was used for the phylogenetic analysis and similarity of nucleotides of HU-91 was calculated using OAT database (orthoANI) (Han et al., 2016).

2.6.3

2.6.3 Genome mining to elucidate the genetics of antagonism

The non-ribosomal peptides and antibiotics were predicted by antiSMASH 3.0 and compared with that of 15 closest strains belonging to same taxonomic group. The carbohydrate related enzymes involved in antifungal activities were detected by using dbCAN3 sever with a combination of DIAMOND and HMMER databases (Xu et al., 2022). The percent identity obtained from biosynthetic gene cluster was plotted as heatmap by using Originpro (version 2021).

2.6.4

2.6.4 Genome mining to elucidate the genetics of antimicrobial resistance

Antimicrobial resistance genes in the assembled genome were detected through comprehensive antibiotic resistance gene database (CARD) and homolog detection was done through loose, strict, and perfect algorithm.

2.6.5

2.6.5 Detection of CRISPR and plasmid

The genome of HU-91 was analyzed for the detection of CRISPR Cas system and plasmid by using CRISPR Cas Finder and Plasmid Finder on default settings (Kamruzzaman & Iredell, 2020). All the tools used for genome analysis were used in default settings.

3

3 Results

3.1

3.1 Antagonistic bacteria associated with wheat

Ten out of 50 bacterial isolates showed antifungal activity against one or more phytopathogens of wheat such as F. graminearum, F. moniliforme, F. oxysporum and F. solani. The antagonistic bacteria were found to be rhizospheric (6), root endophytic (2) and both rhizospheric as well as root endophytic (2) (Table 1). They had a variable colony morphology such as round, oval with irregular and smooth margins. Most of the isolates (n = 7) were Gram-positive while a few (n = 3) were Gram-negative (Table 1). The antagonistic bacteria showed antifungal activity against more than one phytopathogen (Table 1). The maximum number of strains antagonized F. moniliforme, followed by F. oxysporum, F. graminearum and F. solani. HU-91 showing maximum antagonistic activity was selected for further analysis.

Table 1 Morphological characteristics and antagonistic activity of rhizobacteria.
Isolate code Endosphere/Rhizosphere Colony morphology 1Gram reaction 2Antagonistic activity (%)
Shape Margin F. oxysporum F. moniliforme F. graminearum F. solani
HU-91 Endosphere/Rhizosphere Round Irregular + 51.66a ± 0.65 52.85a ± 0.71 48.51a ± 1.18 55.93a ± 0.91
HU-92 Endosphere Oval Smooth + 33.80d ± 0.72 48.59b ± 0.83
HU-93 Rhizosphere Round Irregular + 29.93e ± 0.50 25.88e ± 0.83 25.21e ± 0.93
HU-94 Rhizosphere Round Irregular + 25.59f ± 0.81 36.99d ± 0.40 52.89b ± 0.72
HU-95 Endosphere/Rhizosphere Round Irregular + 21.74 g ± 0.84 30.22e ± 0.60 21.85f ± 0.81 20.40f ± 0.74
HU-96 Rhizosphere Oval Irregular 45.63b ± 0.72 32.77d ± 0.73
HU-97 Endosphere Round Irregular 37.25c ± 0.82 41.11c ± 0.56 28.55de ± 0.56 33.22c ± 0.77
HU-98 Rhizosphere Oval Smooth 26.83f ± 0.60 36.83c ± 0.60
HU-99 Rhizosphere Round Irregular + 16.92 g ± 0.83 31.15d ± 0.55 28.07d ± 0.59
HU-100 Rhizosphere Round Irregular + 42.70b ± 1.08
+ (Gram positive) - (Gram negative).
- (no antagonism) F (Fusarium).

3.2

3.2 Screening of HU-91 for lipopeptides production

HU-91 showed maximum adherence to hydrocarbons (80 %) which is close to the positive control Bacillus cabrialesii SM 93 (87 %) while the oil spreading, and emulsification activity of HU-91 was 2.1 cm and 73 % respectively (Table 2).

Table 2 Biochemical and genetic profile of lipopeptides in B. velezensis HU-91.
Biochemical attributes of lipopeptides
Biochemical trait Activity
BATH activity 80 % ±3
Oil spread activity 2.12 cm ± 0.56
Emulsification activity 73 % ± 2
Genetics of lipopeptides
Lipopeptide name Genes
Surfactin srfAB, srfAA, srfAC, srfAD
Fengycin fenABCDE
Bacillolysin npr
Plipastatin ppsB-ppsE

3.3

3.3 Antimicrobial resistivity of b. Velezensis HU-91

HU-91 exhibited a varied pattern of resistance against different antimicrobial classes. HU-91 exhibited highest resistance to β lactams and least resistance was shown against lincosamides, and macrolides (Fig. 1). HU-91 showed resistance to rifampin (RIF, 5 µg), clindamycin (CLI, 2 µg), aztronam (AZM, 30 µg), cefepime (FEP, 30 µg), ampicillin (AMP, 10 µg), amoxicillin (AML, 30 µg), oxacillin (OXA, 5 µg) and penicillin (PEN, 10 µg) (Table 2). The tested strain showed intermediate resistance against piperacillin (PIP, 100 µg). HU-91 exhibited susceptibility to all the antibiotics belonging to class nitrofurantoins, phenicols, glycopeptides, fluoroquinolones, aminoglycosides, oxazolidinones, tetracyclines and sulfonamides. HU-91 was categorized as multi drug resistance. The cumulative MAR index of HU-91 was higher than 0.2 while the MAR index based on different classes was variable i.e. 0.25 for macrolides, 1 for lincosamide, 0.7 for β lactams and 0 for nitrofurantoins, phenicol, glycopeptides, fluoroquinolones, aminoglycosides, oxazolidinones, tetracyclines and sulfonamides.

Heat map indicating the antimicrobial resistance of B. velezensis HU-91. Values in the box against each category represent the number of antimicrobial agents.
Fig. 1
Heat map indicating the antimicrobial resistance of B. velezensis HU-91. Values in the box against each category represent the number of antimicrobial agents.

3.4

3.4 Whole genome analysis of b. Velezensis HU-91

The genome of HU-91 was assembled in 76 contigs (3,952,697 bp) with a mean GC content of 46.3 %. The total number of genes predicted was 3912, out of which, 3729 were protein-coding genes (118 total RNAs, 84 tRNA, 29 rRNA, and 5 nc RNA) (Fig. 2). BLAST similarity search indicates sequence similarity with Bacillus velezensis. The genome sequence of HU-91 is available at the Gene Bank under accession number JAAHBR000000000.1. The genome was classified into 27 subsystems with genes encoding different functions (Fig. 4).

Circular genome map of B. velezensis HU91. From outer to inner circle [tRNA, rRNA and open reading frame (green, black, blue)). Middle circle (red) (Blast sequence similarity with reference genome (B. velezensis FZB42)]. Inner circle (GC skew and GC content).
Fig. 2
Circular genome map of B. velezensis HU91. From outer to inner circle [tRNA, rRNA and open reading frame (green, black, blue)). Middle circle (red) (Blast sequence similarity with reference genome (B. velezensis FZB42)]. Inner circle (GC skew and GC content).

3.5

3.5 Phylogenomic analysis

Pairwise genome comparison of HU-91 showed sequence similarity with Bacillus velezensis FZB42 based on the phylogenetic analysis of whole genome with average nucleotide identity of 98.49 (Fig. 3).

Phylogenetic tree showing similarity of the B. velezensis HU-91 with other Bacillus spp. The tree was constructed based on average nucleotide identity (ANI) which was calculated by using Orthologous Average Nucleotide Identity Tool (OAT). Numeric values on the branch represent branch length and values in colored boxes represent percentage sequence similarity.
Fig. 3
Phylogenetic tree showing similarity of the B. velezensis HU-91 with other Bacillus spp. The tree was constructed based on average nucleotide identity (ANI) which was calculated by using Orthologous Average Nucleotide Identity Tool (OAT). Numeric values on the branch represent branch length and values in colored boxes represent percentage sequence similarity.
Subsystem category distribution of B. velezensis HU-91 genome. Numeric values in brackets represent number of genes.
Fig. 4
Subsystem category distribution of B. velezensis HU-91 genome. Numeric values in brackets represent number of genes.

3.6

3.6 Antifungal metabolites detected in genome of b. Velezensis HU-91

A total of ten secondary metabolites genes cluster comprising six non ribosomal peptide synthesis (NRPS), one for each metabolite (bacillibactin, locillomycin, fengycin, bacillaene) and two NRPS for surfactin biosynthesis was identified (Fig. 5). The genome of HU-91 showed similarity in their biosynthetic gene cluster with its closely related strains such as B. velezensis, B. amyloliquefaciens and B. subtilis (Fig. 6a). HU-91 showed variation in its biosynthetic gene clusters with that of distinctly related genera (Fig. 6b). The genome of HU-91 carried a unique metabolite locillomycin, comprising 21 % identity with B. subtilis and butirosin, comprising 7 % identity with B. circulans.

Secondary metabolite gene clusters encoding antibiotics. a) Bacilysin b) Macrolactin, c) Bacillibactin d) Surfactin, and e) Amyloliqucifidin f) Locillomycin g) Difficidine h) Bacillaene. Numeric values show the position of genes.
Fig. 5
Secondary metabolite gene clusters encoding antibiotics. a) Bacilysin b) Macrolactin, c) Bacillibactin d) Surfactin, and e) Amyloliqucifidin f) Locillomycin g) Difficidine h) Bacillaene. Numeric values show the position of genes.
Comparison of antibiotics biosynthetic gene clusters of B. velezensis HU-91. a) Comparison among different taxonomic groups (Bacillus spp.) b) Comparison among same taxonomic group (B. amyloliquefaciens).
Fig. 6
Comparison of antibiotics biosynthetic gene clusters of B. velezensis HU-91. a) Comparison among different taxonomic groups (Bacillus spp.) b) Comparison among same taxonomic group (B. amyloliquefaciens).

3.7

3.7 Analysis of CAZymes in genome of HU-91

A total of 205 CAZymes were detected in the whole genome of HU-91. The CAZymes (205) were classified into 5 families such as glucoside hydrolases (68), glycosyltransferases (85), carbohydrate esterases (19), carbohydrate binding modules (24), auxiliary activities (6) and polysaccharide lyases (3). The genes belonging to glucoside hydrolase family were involved in production of enzymes such as chitinases, pectinases, and beta glucanases respectively (Fig. 7).

Distribution of CAZymes in the genome of B. velezensis HU-91. a) Number of CAZymes families in HU-91b) Number of genes present in subfamilies of glycoside hydrolase involved in biosynthesis of major hydrolytic enzymes.
Fig. 7
Distribution of CAZymes in the genome of B. velezensis HU-91. a) Number of CAZymes families in HU-91b) Number of genes present in subfamilies of glycoside hydrolase involved in biosynthesis of major hydrolytic enzymes.

3.8

3.8 Detection of antimicrobial resistance genes in b. Velezensis HU-91

The B. velezensis HU-91 showed genes responsible for the antimicrobial resistance (Table 3). The antimicrobial resistance genes identified belonged to different antimicrobial classes. Furthermore, a CRISPR Cas system (136 bp) was identified while no plasmid was detected in this strain (Fig. 4).

Table 3 Biochemical and genetic profile of antimicrobial resistance in B. velezensis HU-91.
Antibiotic class Resistance phenotype MAR index Resistance genotype
β lactams ATM, FEP, AMP, AMX, OXA PEN 0.75 Bcl, penP, bla2, blm, ccrA, blaB
Fosfomycins FosB
Glycopeptides vanT
Lincosamides Streptogramins Oxazolidinones Phenicols
Pleuromutilins
CLI 1 ClbA, ImrB
Macrolides RIF 0.25 msr, vmIR
Tetracyclines tetA
Fluoroquinolones gyrA-B

ATM (Aztreonam), FEP (Cefepime), AMP (Ampicillin), AMX (Amoxicillin), OXA (Oxacillin), PEN (Penicillin), CLI (Clindamycin), RIF (Rifampin), MAR (multiple antimicrobial resistance), “-” (no resistance phenotype).

4

4 Discussion

Rhizobacteria are an important part of soil ecosystem playing role in soil fertility, nutrient recycling, and control phytopathogens. The use of rhizobacteria as biocontrol agent is an ecofriendly strategy. However, many rhizobacteria species isolated from soil have antimicrobial resistance genes (Zhang et al., 2021). In the current study, rhizobacteria were isolated from rhizosphere and endosphere of wheat which was grown in monocropping system under semi-arid region. In mono cropped wheat, the rhizospheric bacteria were higher in number than that of the endophytic bacteria. This difference may be due to the highly nutritious environment of rhizosphere which makes it a hotspot for the microbial activity, habitat specificity, and passive bacterial transport through roots (Vandana et al., 2021).

These rhizobacteria showed variable antagonistic activity against more than one Fusarium pathogens. The antagonistic property of a bacteria depends upon the type of pathogen, their genetic makeup, production of secondary metabolites, and hydrolytic enzymes. The antagonistic activity of rhizobacteria could also be due to the repeated mono culturing of wheat which governs the rhizospheric microbial community and determines the nature of soil like suppressive or non-suppressive (Ullah et al., 2020).

HU-91 showed maximum antagonistic activity against Fusarium spp. A tremendous amount of literature has been documented on the antagonistic behaviors of B. velezensis against Fusarium (Diabankana et al., 2022; Moreno-Velandia et al., 2021). The strain also produces lipopeptides which are broad spectrum antifungal metabolites.

Biocontrol agents could be a carrier of antimicrobial resistance genes. Therefore, HU-91 was screened for drug resistance traits. In the current study, HU-91 was categorized as MDR because it resisted three antimicrobial classes. The highest resistance was shown against β lactam while lowest was shown against macrolides, and lincosamides. The MAR index of HU-91 was higher than 0.2 which exhibited a potential risk associated with this beneficial bacterium. The resistance/susceptibility of bacteria to these antibiotics have been reported previously by Fraccalvieri et al., (2022). The resistance/susceptibility of bacteria to antimicrobial agents depends upon several factors such as continuous exposure to antimicrobial contaminants, production of antibiotic degrading enzymes by bacteria and presence of efflux pumps (Bianco et al., 2021). The bio control agents could use the same mechanisms of resistance for the antibiotics production and resistance (multi drug efflux pumps) (Zhang et al., 2021).

These traits of HU-91 were further validated by the whole genome analysis. The genomic sequencing of HU-91 revealed 67 contigs and it was identified as a member of Bacillus velezensis. The size of the HU-91 genome was somewhat comparable with its closest strain Bacillus velezensis FZB42 (3918596 bp) as reported earlier (Zaid et al., 2022). The inclusion of some genes involved in other unique functions may account for a slight variation in genome size. The phylogenetic analysis based on the whole genome revealed the association of HU-91 with B. velezensis FZB42. ANI a new metric for pairwise genome comparison also confirmed the 98 % similarity of HU-91 with Bacillus velezensis (Han et al., 2016).

WGS predicted ten secondary metabolites genes clusters in HU-91 genome. The presence of antimicrobial functions such as surfactin, fengycin, bacilleane, difficidin, macrolactin and bacillibactin (catechol iron siderophore) endows the potential to HU-91 as effective biocontrol agent (Nifakos et al., 2021; Zhang et al., 2022). Biosynthetic gene clusters for butirosin, locillomycin and amyloliquecidin showed 7 %, 21 % and 93 % homology which showed a sequence similarity with B. amyloliquefecians genome (Samaras et al., 2021). Genome mining also revealed the presence of carbohydrate active enzymes which cleave the polysaccharides in fungal cell wall into simple monomers (Tsalgatidou et al., 2022). CAZymes analysis of HU-91 showed genes involved in the degradation of fungal cell wall such as chitinases and β glucanase. The presence of these enzymes serves as the inhibitory action of HU-91 which indicates that this strain could play an important role in suppressing the phytopathogens.

The annotated genome carried antimicrobial resistance genes by following the default settings in CARD database. HU-91 genome carried a CRISPR Cas system. CRISPR Cas system plays an important role in bacterial immunity and protects the bacterial genome from invading plasmids (Tao et al., 2022). The presence of CRISPR Cas system could be a possible reason for the absence of plasmid in HU-91 genome. These findings are in line with previous studies (Wheatley & MacLean, 2021). However, the presence of antimicrobial resistant genes in chromosomal DNA could lead to the accumulation of AMRs in the soil upon repeated inoculation. Hence the presence of PGP traits, antibiotic resistance genes, and antagonism provide a new insight into the adaptation of this strain as a potential new candidate for plant growth promotion, antagonism, and drug resistance.

5

5 Conclusion

The present study depicts that B. velezensis HU-91 harbors antifungal and antimicrobial resistance traits but the presence of CRISPR Cas system and absence of plasmid make it potentially safe. However, it is necessary to screen these biocontrol agents before applying them on a commercial scale.

Consent to participate

All authors consent to participate in the manuscript publication.

Consent for publication

All authors approved the manuscript to be published.

Availability of data and material

The data supporting the conclusions of this article are included within the article.

Funding

These experiments were funded by Higher Education Commission (HEC) Pakistan, grant number (NRPU-8062).

Author contributions

PW and MNH have contributed in conceptualization, acquisition of data, and drafting. HU and HY have helped in sampling, research coordination, editing and analysis of data. WH, NK, MNH and AA have contributed in the analysis and interpretation of data, revising and editing the manuscript, and acquisition of funding.

Acknowledgements

The authors would like to thank Higher Education Commission (HEC) Pakistan for providing funds for conducting these experiments under the research grant (NRPU-8062). The authors would like to extend their sincere appreciation to the Researchers Supporting Project number (RSP2023R350), King Saud University, Riyadh, Saudi Arabia.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. , , , , , , , , , , , , , . Current scenario and future prospects of endophytic microbes: promising candidates for abiotic and biotic stress management for agricultural and environmental sustainability. Microb. Ecol. 2023:1-32.
    [CrossRef] [Google Scholar]
  2. , , , , , , , , , , , , , , , , , , , , . Integrating whole-genome sequencing within the National Antimicrobial Resistance Surveillance Program in the Philippines. Nat. Commun.. 2020;11(1):2719.
    [CrossRef] [Google Scholar]
  3. , , , , , , , , . Characterization of Bacillus cereus group isolates from human bacteremia by whole-genome sequencing. Front. Microbiol.. 2021;11:599524
    [CrossRef] [Google Scholar]
  4. , , , . Whole genome sequence insight of two plant growth-promoting bacteria (B. subtilis BS87 and B. megaterium BM89) isolated and characterized from sugarcane rhizosphere depicting better crop yield potentiality. Microbiol. Res.. 2021;247:126733
    [CrossRef] [Google Scholar]
  5. , , . Surface-active agents from two Bacillus species. Appl. Environ. Microbiol.. 1987;53(2):224-229.
    [CrossRef] [Google Scholar]
  6. , , , , . Multiple antibiotic resistance in Escherichia coli isolates from fecal and water sources in Laguna lake, Philippines. Water. 2022;14(9):1517.
    [CrossRef] [Google Scholar]
  7. , , , , . Genetic characteristics and enzymatic activities of Bacillus velezensis KS04AU as a stable biocontrol agent against phytopathogens. Int. J. Plant Biol.. 2022;13(3):201-222.
    [CrossRef] [Google Scholar]
  8. , , , , , , . Preliminary characterization of the first phosphate solubilizing bacteria isolated from the rhizospheric soil of the Jujube Plant Ziziphus. SSRN Electron. J. 2019
    [CrossRef] [Google Scholar]
  9. , , , , , , , , , , , , , . Toxigenic genes, pathogenic potential and antimicrobial resistance of Bacillus cereus group isolated from ice cream and characterized by whole genome sequencing. Foods. 2022;11(16):2480.
    [CrossRef] [Google Scholar]
  10. , , , . ANItools web: a web tool for fast genome comparison within multiple bacterial strains. Database. 2016;2016:1-5.
    [CrossRef] [Google Scholar]
  11. , , . CRISPR-Cas system in antibiotic resistance plasmids in Klebsiella pneumoniae. Front. Microbiol.. 2020;10:2934.
    [CrossRef] [Google Scholar]
  12. , , , , , , , , , , , , , , , , , . Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect.. 2012;18(3):268-281.
    [CrossRef] [Google Scholar]
  13. , , , , , , , , , , , , . Co-selection of antibiotic resistance genes, and mobile genetic elements in the presence of heavy metals in poultry farm environments. Sci. Total Environ.. 2021;755:142702
    [CrossRef] [Google Scholar]
  14. , , , , . Biosynthesis of cyclic lipopeptides by Bacillus velezensis Bs006 and its antagonistic activity are modulated by the temperature and culture media conditions. Curr. Microbiol.. 2021;78(9):3505-3515.
    [CrossRef] [Google Scholar]
  15. , , , , , , , . Prevalence of wheat associated Bacillus spp. and their bio-control efficacy against Fusarium root rot. Front. Microbiol.. 2022;12:798619
    [CrossRef] [Google Scholar]
  16. , , , , , . Complete genome sequence analysis of a plant growth-promoting phylloplane Bacillus altitudinis FD48 offers mechanistic insights into priming drought stress tolerance in rice. Genomics. 2023;115(1):110550
    [CrossRef] [Google Scholar]
  17. , , , , , , , , , , , . Genomic analysis and secondary metabolites production of the endophytic Bacillus velezensis bvel1: a biocontrol agent against Botrytis cinerea causing bunch rot in post-harvest table grapes. Plants. 2021;10(8):1716.
    [CrossRef] [Google Scholar]
  18. , , , . Adherence of bacteria to hydrocarbons: a simple method for measuring cell-surface hydrophobicity. FEMS Microbiol. Lett.. 1980;9(1):29-33.
    [CrossRef] [Google Scholar]
  19. , , , , , , , , . Whole genome sequencing and root colonization studies reveal novel insights in the biocontrol potential and growth promotion by Bacillus subtilis MBI 600 on cucumber. Front. Microbiol.. 2021;11:600393
    [CrossRef] [Google Scholar]
  20. , , , , . The application of the CRISPR-Cas system in antibiotic resistance. Infection Drug Resistance. 2022;15:4155-4168.
    [CrossRef] [Google Scholar]
  21. , , , , , , , . Integrated genomic and metabolomic analysis illuminates key secreted metabolites produced by the novel endophyte Bacillus halotolerans Cal.l.30 involved in diverse biological control activities. Microorganisms. 2022;10(2):399.
    [CrossRef] [Google Scholar]
  22. , , , , , , , . Multitrait Pseudomonas spp. isolated from monocropped wheat (Triticum aestivum) suppress Fusarium root and crown rot. Phytopathology. 2020;110(3):582-592.
    [CrossRef] [Google Scholar]
  23. , , , , , , , , , , , . The endophytic microbiome as a hotspot of synergistic interactions, with prospects of plant growth promotion. Biology. 2021;10(2):101.
    [CrossRef] [Google Scholar]
  24. , , , . Role of horizontal gene transfer in the dissemination of antimicrobial resistance in food animal production. Curr. Opin. Food Sci.. 2022;47:100882
    [CrossRef] [Google Scholar]
  25. , , , , , . Genome sequence and comparative analysis of fungal antagonistic strain Bacillus velezensis LJBV19. Folia Microbiol.. 2023;68(1):73-86.
    [CrossRef] [Google Scholar]
  26. , , , , , , , , , . Prevalence of antimicrobial resistance and respective genes among Bacillus spp., a versatile bio-fungicide. Int. J. Environ. Res. Public Health. 2022;19(22):14997.
    [CrossRef] [Google Scholar]
  27. , . Crossroads of antibiotic resistance and biosynthesis. J. Mol. Biol.. 2019;431(18):3370-3399.
    [CrossRef] [Google Scholar]
  28. , , . CRISPR-Cas systems restrict horizontal gene transfer in Pseudomonas aeruginosa. ISME J.. 2021;15(5):1420-1433.
    [CrossRef] [Google Scholar]
  29. , , , , , , , , , , , , . Biosynthetic gene cluster profiling predicts the positive association between antagonism and phylogeny in Bacillus. Nat. Commun.. 2022;13(1):1023.
    [CrossRef] [Google Scholar]
  30. , , , , , , , , , , , . Genomic and phenotypic insights into the potential of Bacillus subtilis YB-15 isolated from rhizosphere to biocontrol against crown rot and promote growth of wheat. Biology. 2022;11(5):778.
    [CrossRef] [Google Scholar]
  31. , , , , , . Comparative genome analysis reveals phylogenetic identity of Bacillus velezensis HNA3 and genomic insights into its plant growth promotion and biocontrol effects. Microbiol. Spectr.. 2022;10(1) e02169-21
    [CrossRef] [Google Scholar]
  32. , , , , . Isolation and chemical characterization of the biosurfactant produced by Gordonia sp. IITR100. Plos One. 2022;17(4):e0264202.
    [CrossRef] [Google Scholar]
  33. , , , , , , , . Coexistence between antibiotic resistance genes and metal resistance genes in manure-fertilized soils. Geoderma. 2021;382(114760):114760
    [CrossRef] [Google Scholar]
  34. , , , , , , , , , . Mechanism of a volatile organic compound (6-methyl-2-heptanone) emitted from Bacillus subtilis ZD01 against Alternaria solani in potato. Front. Microbiol.. 2022;12:808337
    [CrossRef] [Google Scholar]
Show Sections