Analysis of Biosynthetic Gene Clusters, Secretory, and Antimicrobial Peptides Reveals Environmental Suitability of Exiguobacterium profundum PHM11

Halotolerant bacteria produce a wide range of bioactive compounds with important applications in agriculture for abiotic stress amelioration and plant growth promotion. In the present study, 17 biosynthetic gene clusters (BGCs) were identified in Exiguobacterium profundum PHM11 belonging to saccharides, desmotamide, pseudaminic acid, dipeptide aldehydes, and terpene biosynthetic pathways representing approximately one-sixth of genomes. The terpene biosynthetic pathway was conserved in Exiguobacterium spp. while the E. profundum PHM11 genome confirms the presence of the 1-deoxy-d-xylulose 5-phosphate (DXP) pathway for the isopentenyl diphosphate (IPP) synthesis. Further, 2,877 signal peptides (SPs) were identified using the PrediSi server, out of which 592 proteins were prophesied for the secretion having a transmembrane helix (TMH). In addition, antimicrobial peptides (AMPs) were also identified using BAGEL4. The transcriptome analysis of PHM11 under salt stress reveals the differential expression of putative secretion and transporter genes having SPs and TMH. Priming of the rice, wheat and maize seeds with PHM11 under salt stress led to improvement in the root length, root diameters, surface area, number of links and forks, and shoot length. The study shows that the presence of BGCs, SPs, and secretion proteins constituting TMH and AMPs provides superior competitiveness in the environment and make E. profundum PHM11 a suitable candidate for plant growth promotion under salt stress.


INTRODUCTION
The studies of extremophiles like Exiguobacterium are of high importance to address the potential of agricultural growth and production under abiotic stresses using these microorganisms for plant growth promotion (PGP) and stress amelioration (Hayat et al., 2012;Srivastava et al., 2019a,b). The genus Exiguobacterium possesses a variety of traits to combat environmental challenges and have been isolated from diverse environments, including permafrost, hot springs, oceans, the rhizosphere of plants, and food-processing plants as well as from alimentary canals (Vishnivetskaya et al., 2009). The presence of antimicrobial peptides (AMPs), secondary metabolite production, non-ribosomal peptide synthesis (NRPS), and polyketide (PK) synthesis provides extra-environmental traits to make it a suitable candidate for application in a variety of pharmaceutical, agricultural, and industrial applications (Kulshreshtha et al., 2013;Singh et al., 2013;Dastager et al., 2015;Meng et al., 2020). The comparative genome studies and advancements in many technologies have been exploited to identify NRPS, PK, antibiotic resistance, and salt tolerance in addition to PGP activities (Hugenholtz and Tyson, 2008;Menzel et al., 2016;Patel et al., 2018;Bharati et al., 2020).
In the last few decades, many secondary metabolite pathways have been identified in the microbes, which contribute toward its potential and environmental suitability (Pichersky et al., 2006;Dickschat, 2016). With the development of bioinformatics and biotechnologies, many secondary metabolite biosynthetic gene clusters (BGCs) have been characterized in bacteria. These clusters commonly contain a gene encoding one of several key signature enzymes (Yamada et al., 2015). The finding of BGCs in bacteria shows that the secondary metabolites are also widely distributed in microbes (Pichersky et al., 2006;Dickschat, 2016). These secondary metabolites are not necessary for the growth or development of the host but sometimes confer pathogenicity, virulence, and adaptation to the environment and play an important role in defense mechanisms (Osbourn, 2010;Wang and Lu, 2017).
The genome of Exiguobacterium sp. 8-11-1 has 134 genes for metabolism and inorganic ion transport, including five NhaC-type Na + /H + antiporters and one MnhB-type Na + /H + antiporter which allow bacteria to survive and grow at alkaline pH (Jiang et al., 2013). Most of the transporters have been identified with the signal peptides (SPs) with a secretory nature. These SPs have a role in the targeting of the protein, but the mechanism is still unknown (Galán et al., 2014). These secretion proteins possess transmembrane helices (TMHs) which have a role in the transportation of the substrate (Krampen et al., 2018). Further studies have shown the presence of PGP and biocontrol genes for auxin biosynthesis and siderophore and iron acquisition in different Exiguobacterium spp. (Tang et al., 2013;White et al., 2013). E. chiriqhucha N139 isolated from cold lakes is a recently reported species, and its genome was sequenced due to its stress defense ability, like growth in high concentrations of different metals especially arsenic, high ultraviolet B radiation, and salinity (Gutiérrez-Preciado et al., 2017). The genome of N139 contains arsenite efflux pump, ArsB, and ATPase. An ATPase that provides energy to ArsB for extrusion of arsenite and antimonite, ArsA, cotranscribed with ArsD, an arsenic chaperone for the ArsAB pump (Páez-Espino et al., 2009).
In previous studies, we have described the halotolerance potential of E. profundum PHM11 and its pan genome comparison with other genomes of Exiguobacterium. PHM11 exhibits a gene for shape determination (MreB), mannitol/Lproline pathways, carotenoid biosynthesis, and shikimate pathways (Patel et al., 2018). These studies have also described the impact of salinity on the expression of these genes. Further, we reported the presence of several pathways and its comparison with the already sequenced genome of Exiguobacterium . In the present study, different BGCs, secretory peptides (SPs), TMHs, and AMPs were mined from the genome of PHM11 with the aim to decipher its environmental suitability for PGP applications under salt stress. Transcriptome analysis and validation of different PGP traits in a pot trial were also conducted in rice, maize, and wheat.

Phylogenetic Tree Analysis
The 16S rRNA gene sequences of Exiguobacterium sp. were obtained from the NCBI, and alignment was implemented by ClustalW. The aligned file was subsequently used in neighbor joining tree construction by MEGA-X (Tamura et al., 2011). For the pairwise alignment of E. profundum, the sequence was aligned using the ClustalW server 1 .

Comparative Genome Analysis Using the Genome-to-Genome Calculator
Genomes of all type strains of Exiguobacterium spp. were retrieved from NCBI. Only 12 complete or partial whole-genome sequences of the Exiguobacterium were found in the database. The genome of E. profundum PHM 11 was compared with eleven other type strains of Exiguobacterium using the genometo-genome distance (GGD) distance, dDDH prediction, and model confidence interval between pairs of entirely or partially sequenced genomes. The values were calculated in silico by using the genome-to-genome distance calculator (GGDC) tool 2 . GGDC was performed for type strains of Exiguobacterium strains considering E. profundum PHM11 as a standard reference. GGDC functions on the basis of the distances from the set of HSPs or MUMs obtained by comparing each pair of genomes, transformed to values analogous to DDH.

Datasets and Secondary Metabolite Biosynthetic Gene Clusters Finding
We obtained all finished Exiguobacterium genome sequences from EBI at http://www.ebi.ac.uk/genomes/archaea.html. The number and types of secondary metabolite BGCs in these sequences were identified by antiSMASH version 5.1.2, and the integrated ClusterFinder algorithm, a hidden Markov modelbased probabilistic algorithm, was used to detect BGC-like regions in genomes (Cimermancic et al., 2014;Weber et al., 2015;Blin et al., 2017). The characterized and unknown secondary metabolite biosynthesis gene clusters in the E. profundum PHM11 were identified. Additionally, both the Known Cluster Blast and Cluster Blast modules were selected to identify similar clusters in sequenced genomes by genome comparisons. Further, domain functions and genetic similarities with known BGCs in studied gene clusters were predicted and annotated using antiSMASH 5.1.2.

Signal Peptides and Transmembrane Prediction
The SP and its prediction for secretion were calculated using the PrediSi server 3 . All the protein sequences present in PHM11 were given as input file in FASTA format. The outputs represent a single numeric score, predicted cleavage site, and Boolean flag denoting CTT (co-translationally translocated) SP presence or absence. PrediSi is designed for extremely fast analysis and is well suited for high-throughput processing (Hiller et al., 2004). The transmembrane domain prediction was done using the TMHMM server 4 . The FASTA format of all the proteins present in PHM11 was uploaded on the server which uses cyclic hidden Markov models to predict membrane helices from the proteins.
Mining of Antimicrobial Peptides From the Genome of Exiguobacterium profundum PHM 11 BAGEL 4 was used for the analysis of AMPs in genomic data. BAGEL is a web server that identifies putative bacteriocin ORFs in a DNA sequence using novel, knowledge-based bacteriocin databases and motif databases. Gene clusters of interest are discovered using the core-peptide database and/or through HMM motifs that are present in associated context genes.

Validation for Antimicrobial Activity and Environmental Fitness of Exiguobacterium profundum
Validation for antimicrobial activity and environmental fitness of E. profundum PHM11 was done by an antagonism experiment against 63 bacterial and 31 fungal strains isolated from different environmental samples (Supplementary Table S1). A loop full of bacterial culture was inoculated in nutrient broth (NB) and incubated at 28 • C for 72 h. Cell-free supernatant from NB was collected by centrifugation at 10,000 g for 10 min. The supernatant was sterilized by filtration through a 0.45-µm syringe filter. The extract obtained was screened for antibacterial activity by the well diffusion method. The test bacterial cultures were spread on a set nutrient agar plate, and immediately four wells of 5 mm were punched with a sterile borer, and the wells were filled with 20 µl of cell-free extract or NB (2 wells each). The plates were then incubated at 28-30 • C for 18 h in the upright position. Inhibition was detected by a zone of clearing around the extract well (Baskaran et al., 2012). Effects of volatile organic compounds (VOCs) of E. profundum PHM11 on growth inhibition of thirty-one selected fungi were recorded. Seventytwo hour-old NA plates of E. profundum PHM11 and 48-h-old inoculated plates of test fungal isolates were used for inhibition assay. The sandwich plates were made under laminar airflow, and edges were sealed with parafilm to avoid any loss of VOCs. The

Transcriptome Analysis of Exiguobacterium profundum PHM11
To investigate the transcriptional response of E. profundum PHM11 to osmotic stress, two RNA-seq libraries were generated in two different conditions, viz., 0 mM NaCl and 100 mM NaCl. Cell samples were collected during the exponential phase when cultures reached enough biomass to assure the isolation of 100-500 ng of mRNA for RNA-seq library construction. Total RNA was isolated using the manufacturer's protocol of the GeneJET RNA Purification Kit with slight modification, as described earlier (Bharati et al., 2016). Two biological replicates were used for transcriptome analysis. The de novo transcriptome sequencing was performed using the Illumina HiSeq platform in the paired-end module (Liu et al., 2014).
The raw fastq files were processed before performing assembly. Prior to the assembly, base trimming, removal of adapter sequences, and filtering out of reads with an average quality score less than 30 were performed in every paired end read. Further, the rRNAs were removed from the generated sequences based on the reference ribosomal RNA sequences available in the SILVA database (Quast et al., 2013). The cleaned reads were aligned (length ≥ 200 bp) using the Bowtie2 program (Langmead and Salzberg, 2012), normalized, and finally subjected to de novo transcriptome assembly using Trinity (Grabherr et al., 2011;Huang et al., 2016). Normalization was performed using the variance analysis package EdgeR program (Robinson et al., 2010). The assembled transcripts were annotated using BLASTX against the non-redundant nucleotide databases. The gene ontology (GO) annotation for Molecular Function (MF), Cellular Component (CC), and Biological Process (BP) for upregulated transcripts was performed using the gene ontology tool (Ashburner et al., 2000). DESeq was used for the differential gene expression analysis (Anders and Huber, 2010). The gene considered differentially expressed if the p-value for differential expression was <0.05 and the absolute log (base 2)-fold change was ≥1. The final figure was represented in the form of heatmap.

Evaluation of Plant Growth-Promoting Attributes of Exiguobacterium profundum PHM 11
In order to evaluate PGP attributes of PHM11, a pot experiment was designed using rice (cv. Rajendra sweta), maize (cv. Bio9637), and wheat (cv. HD2967). The seeds were surface sterilized with 0.5% (v/v) NaOCl for 10 min and were subsequently washed with sterile deionized water. Seeds were air dried and coated with active culture (10 8 cfu ml −1 ) of E. profundum PHM11 suspended in 1% CMC. Approximately 15 seeds were sown per pot, and 90% germination was observed. Plant growth characteristics like total root length, root diameter, projected area, surface area, number of forks, and number of links were determined using the root scanner (Epson Expression 12000XL scanner with WinRHIZO Pro software) for the treated (+) and non-treated (-) in both conditions of salt absence (-) and presence (+) (100 mM NaCl) after 28 days of sowing. The stem length and total chlorophyll content (using a chlorophyll meter, Apogee Instruments, Logan, UT, United States) were also measured. All the laboratory and greenhouse experiments were executed in three replicates.

RESULTS AND DISCUSSION
Phylogenetic Affiliation of Exiguobacterium profundum PHM11 The identity of PHM11 isolated from the soil collected from Mau, Uttar Pradesh, was established as Exiguobacterium profundum on the basis of 16S rRNA sequencing followed by wholegenome sequencing (NCBI accession no. MRSV00000000) (Patel et al., 2018). E. profundum PHM11 was characterized as the bacilli forming a light-orange colony with optimal growth at 37 • C ( Figure 1A). PHM11 was tested positive for catalase and oxidase enzymatic activities. The PGP activities like phosphate solubilization and zinc solubilization as well as siderophore, ammonia, and HCN production were checked and found positive ( Table 1). Similar PGP activities have been already reported in many Exiguobacterium spp. (Tang et al., 2013;White et al., 2013). The comparison of PHM11 with other eighteen type species of Exiguobacterium on the basis of the 16S rRNA sequence put it with the E. profundum ( Figure 1B). Phylogenetic tree analysis showed the two major branching ( Figure 1B) and the FIGURE 1 | Genome comparison of type strains of Exiguobacterium spp. (A) Colony shape and pigmentation of the strain. (B) Phylogenetic tree analysis of type strains of Exiguobacterium spp. using MEGA-X using the 16S rRNA sequence. The sequences were aligned by ClustalW, and the neighbor-joining method was used for the phylogenetic tree preparation. The stars represent branching of the strains in two major groups. (C) The pairwise sequence alignment of the 16S rRNA sequence of E. profundum type strain (10C) and the E. profundum PHM11. The alignment was done using ClustalW, then the alignment file was submitted to ESPript, and a final figure was generated. (D) Represents the graph between GGD and DDH of twelve type strains of Exiguobacterium species with reference to the E. profundum PHM11, and (E) represents AAI vs. ANI graph of 12 type strains of Exiguobacterium species with reference to the E. profundum PHM11.
Frontiers in Microbiology | www.frontiersin.org  Figure 1C) revealed that 58 nucleotides were missing (gap) in the type species.
Furthermore, out of 18 type species of Exiguobacterium, the whole genomes of only 12 species were available in the NCBI database. Comparison of these twelve genomes with E. profundum PHM11 on the basis of the GGD and DDH showed the closeness of E. marinum with PHM11 having the 80.6% DDH value, 0.1296 GGD, and model CI value [76.6 -84%] ( Figure 1D and Table 2). E. chiriqhucha, E. aurantiacum, E. mexicanum, and E. alkaliphilum were found the second closest to E. profundum PHM11 which showed ∼40% DDH value, ∼0.4 GGD, and model CI value [10.7 -16.8%]. The rest of the strains E. sibricum, E. undae, E. antarcticum, E. indicum, E. oxidotolerans, E. acetylicum, and E. enclense showed less than 14% DDH, more than 0.9 GGD, and a model CI value in the range of 10-17%. The result was further confirmed with the comparison of the ANI value and AAI value with reference to PHM11 ( Figure 1E). The graph between ANI and AAI showed a similar pattern to that of the DDH and GGD graph. The graph indicated the genomic differences between the species as well. Because of the genomic differences, they possess a variety of habitat, PGP traits, and environmental suitability (Pichersky et al., 2006;Vishnivetskaya et al., 2009;Dickschat, 2016).

Metabolite Analysis of PHM11 Identified Putative Biosynthetic Gene Clusters
Our earlier reports identified E. profundum PHM11 a rich source of the carotenoids, and the production was enhanced under salt stress (Patel et al., 2018). This led us to mine additional secondary metabolite pathways. In order to mine the secondary metabolite pathways, the E. profundum PHM11 genome was analyzed through antiSMASH version 5.1.2 and identified seventeen regions in the genome having putative BGCs with 514 genes (Figure 2A and Supplementary Table S2). Most of the BGCs (15) corresponded to saccharide type; one BGC corresponded to each fatty acid and terpene biosynthesis. Regions 1.3, 1.5, 1.6, 1.7, 1.9, and 1.14 belonged to most similar known clusters of desotamide, capsular polysaccharides, O-antigen, emulsan, pseudaminic acid, and dipeptide aldehydes, respectively. All these BGCs comprised approximately 16.6% of the genome. Region 1.6 is the largest one with ∼63 kb having more than 157 genes which belonged to the o-antigen and polysaccharide production, while region 1.1 was the smallest BGC with ∼19 kb nucleotides only having 19 genes in the cluster. Regions 1.3, 1.9, and 1.14 represented the presence of the clusters of gene for the biosynthesis of desotamides, dipeptide aldehydes, PK, and NRPS. Desotamides are a natural cyclohexapeptide antibiotic which has a role in the inhibition of certain microbes (Chen et al., 2018). Dipeptide aldehydes and the non-ribosomal peptides are a diverse family of natural products that fall into the class of secondary metabolites with diverse properties as toxins, siderophores, pigments, antibiotics, immune-suppressants, and anticancer agents (Bushley and Turgeon, 2010;Wang et al., 2014;Martínez-Núñez and López, 2016). Similarly, PKs are polycyclic aromatic natural products having anticancerous, antibacterial, antifungal, antiviral, antiparasitic, and other medicinally significant properties synthesized by PK synthases (PKSs) in soil-borne bacteria (Das and Khosla, 2009).
In PHM11, we observed a number of traits for PK as well as non-ribosomal protein synthesis. Other closely related Firmicutes are also known to abundantly produce various NRPS and PKs. Genome mining studies have indicated that 31% of the Firmicutes harbor NRPS and PK secondary metabolite gene clusters (Wang et al., 2014). Various members of the genus Bacillus like B. subtilis, B. thuringiensis, B. cereus, B. velezensis, and B. licheniformis are well known for their ability to produce an array of non-ribosomal and PK antimicrobial compounds. It has been reported that almost 10% of the genome of B. velezensis FZB42 is for the synthesis of different antimicrobial compounds (Borriss et al., 2019). Non-ribosomal peptides (NRPs) are synthesized through enzyme-mediated condensation of amino acid residues where > 300 different precursor molecules help in the assemblage of NRPs (Saxena et al., 2020). A number of researchers indicated that cyclic lipopeptides are one of the most abundant NRPs in Bacillus (Mora et al., 2015;Solanki et al., 2015;Saxena et al., 2017Saxena et al., , 2020. The genomes of B. subtilis, B. thuringiensis, B. velezensis, and B. cereus present high abundance of NRP gene clusters (Zhao and Kuipers, 2016). Among the lipopeptides, iturin and surfactin are the most abundant in Bacillus (Mora et al., 2011;Kushwaha et al., 2020;Saxena et al., 2020). Bacillaene, Difficidin, and Macrolactin are some of the important PK antimicrobial compounds produced by various species of Bacillus (Chen et al., 2006;Wang et al., 2015). B. amyloliquefaciens and B. subtilis were identified as prolific producers of PKs. However, other species like B. atrophaeus and B. mojavensis are also known to produce PKs (Wang et al., 2015). The genome (3.92 Mb) of B. amyloliquefaciens FZB42 dedicated ∼200 kb toward production of PKs. The genomic and transcriptomic data of E. profundum PHM11 presented here provide insight for the need of quantitative shotgun proteomics to further validate and improve the protein identification and integration to their functional context. It has been suggested that changes in gene expression do not necessarily correlate with the adaptive manifestations, while proteins are rather more directly related (Evans, 2015). Further, proteomic analysis are capable of estimating absolute protein amounts and help in the analysis of low-abundance and modified proteins which can be of high biological relevance (Semanjski and Macek, 2016).

Exiguobacterium profundum PHM11's Biosynthetic Gene Clusters Contain Genes Related to Environmental Fitness
The third cluster (Region 1.3) contained 42 genes in ∼38-kb size (Figure 2A). This cluster contained genes for the desotamides and NRPs, which includes diguanylate cyclase, sugar ABC transporter ATP-binding protein, tryptophan 2,3-dioxygenase, and kynureninase (Supplementary Table S2). The NRPs generally contain amino acids like ornithine or amino acids, and their structures are macrocyclic, or branched macrocyclic, dimers or trimers of identical structural elements (Singh et al., 2012). The presence of arginine-ornithine antiporter participates in the direct supply of the ornithine, while the presence of arginine deiminase and ornithine carbamoyltransferase genes supports the biosynthesis (Bushley and Turgeon, 2010;Wang et al., 2014). Another enzyme UDP-N-acetylmuramoyl-Lalanyl-D-glutamate-2,6-diaminopimelate ligase (MurE) was also present which catalyzed a reaction in the cytoplasmic step of peptidoglycan biosynthesis by adding the third amino acid residue to the peptide stem (McGroty et al., 2013). The 9 th cluster (Region 1.9) contained 59 genes in 49 kb of stretch along with the pseudaminic acid synthase (Figure 2A). Pseudaminic acid synthase catalyzes the condensation of phosphoenolpyruvate (PEP) with the hexose, 2,4-diacetamido-2,4,6-trideoxy-L-altrose, to form pseudaminic acid and phosphate (Chou et al., 2005). The isomers of the pseudaminic acid are part of the lipopolysaccharide, and their production is important for the invasiveness of pathogenic organisms (Knirel et al., 1997). Besides these proteins, this cluster contained genes for antibiotic biosynthesis monooxygenase (ABM) and  Table S2). The ABM has the fidelity control during aromatic PK biosynthesis (Qin et al., 2019). The 14 th cluster (Region 1.14) contained 30 genes spanning ∼29 kb (Figure 2A). These genes included asparagine-tRNA ligase, pyridoxal phosphate-dependent aminotransferase, aspartate 1-decarboxylase, pantoate-beta-alanine ligase, CCA tRNA nucleotidyltransferase, N-acetyl-alpha-D-glucosaminyl Lmalate synthase (BshA), bacillithiol biosynthesis deacetylase (BshB1), tetratricopeptide repeat protein, and 3-dehydroquinate synthase which participated in the PK and NRP synthesis (Supplementary Table S2). The cluster also possessed N-acetyl-alpha-D-glucosaminyl L-malate synthase (BshA) and bacillithiol biosynthesis deacetylase (BshB1) which had key roles in the biosynthesis of bacillithiol (BSH) and lowmolecular-weight (LMW) thiols (Gaballa et al., 2010). LMW thiols play critical roles in cell physiology during oxidative stress. In most cells, the cytosol is reducing and protein thiols are largely in the reduced state (Dalle-Donne et al., 2009). During the oxidative stress, the disulfide bond formation between cysteine thiols provides stability and determines the structure of extra-cytoplasmic proteins (Meyer et al., 2009). Thus, the BSH provides protection to the cell against oxidative stress. Furthermore, this cluster possessed 3-dehydroquinate synthase which catalyzed the second step in the shikimate pathway, which was essential for the production of aromatic amino acids in bacteria, plants, and fungi, but not mammals (Liu et al., 2008). This makes it an ideal target for new antimicrobial agents, antiparasitic agents, and herbicides. Previously, our group reported that the genes related to the shikimate pathway get upregulated upon the salinity stress (Patel et al., 2018). In order to further investigate the AMPs, BAGEL 4 was used for the analysis of AMPs in genomic data. A total 181 bacteriocins were identified in the E. profundum PHM11 genome out of which eleven putative bacteriocins and nine cerecidins were identified ( Table 2). We also observed the presence of cerecidin which was a lantibiotic, ribosomally synthesized and posttranslationally modified AMPs ( Table 2). Validation for antimicrobial activity and environmental fitness of E. profundum PHM11 was done by antagonism experiment against a wide range of fungal and bacterial strains isolated from environmental samples. Out of 94 test organisms (63 bacteria and 31 fungi), E. profundum PHM11 was found highly competitive against 43 bacteria and 26 fungi (Supplementary Table S1 and Supplementary Figure S1).

Terpene Biosynthesis Cluster Is
Conserved in Exiguobacterium spp.
The 17 th BGC represented the genes belonging to the saccharide as well as terpene biosynthesis. This region was covered between 2,521,526 and 2,557,456 bp having 31 genes (Supplementary Table S2). The terpene cluster was present between 2,536,720 and 2,556,936 bp having 17 genes ( Figure 2B and Table 3). These 17 genes could be categorized in the core biosynthetic gene, additional biosynthetic genes, transport-related genes, regulatory genes, and other genes. Squalene/phytoene synthase (crtB) is the core biosynthetic gene which is responsible for the converting phytoene into lycopene, an important stage in carotene biosynthesis (Nakashima et al., 1995;Pan et al., 2015). Some additional biosynthetic genes, namely, glycoside hydrolase family 13 protein (crtGH), enoyl-CoA hydratase/isomerase family protein, glycotransferase (crtGT), phytoene desaturase (crtIa, crtIb, and crtIc, three copies at different locations), and class I SAM-dependent methyltransferase (crtMT), were also present. Two transportrelated genes, namely, MATE family efflux transporter (crtT) and ATP-binding cassette domain-containing proteins, were present (crtATP). Two regulatory proteins, GAF domain-containing protein and some other genes, cys-t-RNA (Pro) deacylase, CDPdiacylglycerol-serine O-phosphatidyl transferase, carotenoid biosynthesis protein, and 1-acyl-sn-glycerol-3-phosphate acyltransferase gene, were present in this cluster ( Figure 2B). This cluster had 100% similarity with that of the other Exiguobacterium spp. JLM2, HVEsp1, AT1b, ZOR0005, and NG55; E. marinum DSM 16307; and E. chiriqhucha RW-2 ( Figure 2C). The results indicated that all the genes present in this cluster are ∼100% similar in the Exiguobacterium species which has been characterized for the isolates closely related with E. profundum PHM11 in our previous studies .

Exiguobacterium profundum PHM11 Possesses the 1-Deoxy-D-xylulose 5-phosphate Pathway for IIP and Carotene Biosynthesis
With the rapid development of biocatalysis, β-carotene production through biosynthetic methods has become an active field, and several studies regarding genetic modification to enhance microbial production of β-carotene or carotenoids have been reported (Yang and Guo, 2014). E. profundum PHM11 showed an increased total carotenoid concentration in 100-mM salinity-affected cells (Patel et al., 2018). Cluster 17 shows high similarity to terpene BGCs. Carotenoids are terpenoids found in all photosynthetic organisms and also in some non-phototrophic    Table 4). The first step of the synthesis of the IPP units is the condensation and decarboxylation reaction between the D-glyceraldehyde phosphate and the pyruvate carried out by DXPS (1-deoxy-D-xylulose-5-phosphate synthase) and the production of DXP. The formation of DXP has been observed in many bacteria, and due to the formation of DXP, the pathway is also known as DXP pathway (Rohmer, 1999). These are responsible for the synthesis of the carotenoids as well.
In the second step, DXP is converted to MEP by the enzyme DXR (1-deoxy-D-xylulose 5-phosphate reductoisomerase). The third step represents the formation of 4-diphosphocytidyl-2-C-methyl-D-erythritol from MEP, which is catalyzed by ispD (2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase). Similarly, the 4 th , 5 th , and 6 th steps are catalyzed by ispE, ispF, and ispH, respectively, and after the 6 th step HMBPP was produced which is the substrate for the formation of IPP and DAMPP by the enzyme ispH (4-hydroxy-3-methylbut-2-enyl diphosphate-reductase). Cluster 17 represented the presence of the squalene/phytoene synthase (crtB) and phytoene desaturase (crtI) which is responsible for the conversion of the phytoene to lycopene and then lycopene into carotene, respectively ( Figure 2B). We were unable to trace the enzyme geranylgeranyl diphosphate synthase (crtE) in genomes as well as in cluster analysis which is responsible for the formation of GGDP from the condensation IIP and DMAPP (Nakashima et al., 1995;Pan et al., 2015). Furthermore, we observed the presence of 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase (crtG) enzyme which could catalyze the efficient conversion from β-carotene to 2,2-dihydroxy β-carotene ( Table 4).

Identification of Signal Peptides and Transmembrane Helix Supports the Targeted Substrate Transport
PrediSi predicted 2,877 proteins with SPs having a cutoff score of 0.5. Out of these, 592 were predicted for the secretion and among them, 160 were hypothetical proteins (Supplementary Table S3).
The major secretory proteins with the highest scores were rodshaped-determining protein (MreD), Na/H + antiporter, alphaamylase, nitrous oxidase accessory protein, iron-siderophore ABC transporter substrate-binding protein, and transporter SBD (substrate binding domain-containing) protein with a score more than 0.9 ( Table 5). Along with these, the multidrug efflux SMR transporter, peptidoglycan glycosyltransferase, type I pullulanase, zinc metalloprotease HtpX, zinc ABC transporter solute-binding protein, flagellar type III secretion system pore protein (FliP), and type IV secretory system conjugative DNA transfer family protein were predicted for SPs with secretion having TMH. The result indicated that most of the transporters related to the PGP activity and heavy metal resistance were predicted to have TMH for secretion (Diepold et al., 2015;Krampen et al., 2018). The transmembrane domain of some of the top score-predicted secretion protein was shown (Figure 3). MreD had six TMHs (Figure 3A), the potassium/proton antiporter had twelve TMHs (Figure 3B), ABC transporter permease subunit had seven TMHs (Figure 3C), multidrug efflux SMR transporter had four TMHs (Figure 3D), isoprenylcysteine carboxyl methyltransferase family protein had four TMHs ( Figure 3E) and Na+/H+ antiporter NhaC had ten TMHs ( Figure 3F). Apart from PGP traits, some proteins with the secondary metabolite pathway were also identified which were predicted to have SPs. These include isoprenylcysteine carboxyl methyltransferase family protein and phytoene desaturase ( Table 5). The former is involved in the methyl-transferase activity of isoprenylcysteine (Clarke et al., 1988;Winter-Vann et al., 2005), and the latter one is involved in the terpene biosynthesis pathways by converting phytoene to carotene by desaturase activity (Cunningham and Gantt, 1998;Grünewald et al., 2000). Beside the secondary metabolite proteins, the pullulanase type I enzyme with predicted secretion was also identified ( Table 5). Pullulanase type I (PulA) is a debranching enzyme that specifically cleaves α-1,6-glycosidic linkages in pullulan. Pullulan has not only diverse applications in the food industry but also immune-stimulatory effects on B and T cells and is found to enhance the production of various anti-inflammatory cytokines in human (Coltelli et al., 2020).

Differential Expression of Gene Analysis Validates the Presence of Secretary and Some Unique Hypothetical Proteins
In order to further look into the transcriptomic changes upon 100-mM salinity stress in PHM11, whole-cell transcriptomics was performed. A total of 2,320 transcripts were found, out of which 21 were downregulated ( Figure 4A) and 33 were upregulated (Figure 4B and Supplementary Table S5). Most of the downregulated genes belonged to hypothetical genes whose functions could not be assigned. Some of the proteins identified belonged to extracellular solute-binding protein family 5, flagelline, and apolipoprotein A-I-2-like proteins. The upregulated gene belonged to amino acid carrier protein, cytochrome c oxidase, drug resistance transporter, Na + /H + antiporter (NhaC), multicopper oxidase, acriflavin resistance protein, ferritin, bifunctional purine biosynthesis protein PurH, transport system permease protein, and penicillinbinding protein transpeptidase. Two hypothetical genes were also identified. The drug resistance transporter belonged to the region 9 of the biosynthetic cluster genes belonging to    The left side of the figure represents the root of rice plant grown till 28 days in normal soil (no salt) grown from the untreated seed while the right side represents the root of wheat plant grown till 28 days in normal soil (no salt) grown from the treated seed (PMH11). (B) Represents the appearance of the whole plant shoot grown till 28 days in normal soil (no salt) grown from non-treated (left side) and treated (right side) seed. Similarly, (C) the left side of the figure represents the root of rice plant grown till 28 days in soil supplemented with 100 mM NaCl grown from the untreated seed while the right side represents the root of rice plant grown till 28 days in soil supplemented with 100 mM NaCl grown from the treated seed. (D) Represents the appearance of the whole plant shoot grown till 28 days in soil supplemented with 100 mM NaCl grown from non-treated (left side) and treated (right side) seed.
pseudaminic acid biosynthesis pathways (Knirel et al., 1997;Chou et al., 2005). The upregulated genes were further classified through GO analysis for their functions (Supplementary Table S4). The amino acid carrier protein was predicted with SP with a 0.72 score (WP_074037007.1, Supplementary Table S3) which had four TMHs spanning to the plasma membrane (Supplementary Figure S2). The cytochrome c oxidase was related to the region 14 of the BGCs related to dipeptide aldehyde biosynthesis. These genes were related to the stress response (Dalle-Donne et al., 2009;Meyer et al., 2009;Gaballa et al., 2010). The Na + /H + antiporter (NhaC) was also predicted with the SP having 10 TMH spanning regions ( Figure 3F). The TMHs of the antiporter have a role in substrate transport (Diepold et al., 2015;Krampen et al., 2018).

Exiguobacterium profundum PHM11 Inoculation Promotes the Plant Growth
Exiguobacterium profundum was found positive for phosphate solubilization and zinc solubilization as well as siderophore, ammonia, and HCN production (Páez-Espino et al., 2009;Tang et al., 2013;White et al., 2013;Gutiérrez-Preciado et al., 2017). In addition to these PGP traits, we also investigated the presence of PK, dipeptide aldehydes, and NRP providing the environment suitability to the organism (Das and Khosla, 2009;Bushley and Turgeon, 2010;Wang et al., 2014;Martínez-Núñez and López, 2016;Chen et al., 2018). Because of the presence of all these traits, PHM11 was thought to promote plant growth. Its effect on the growth of paddy (Oryza sativa var. Rajendra sweta) was studied in pot trials under the influence of 100 mM salt (NaCl) ( Table 6). The growth parameters were recorded after 28 days of sowing (Table 6 and Figure 5). The results showed that the roots of the treated plants had multiple branching and were healthier as compared to the non-treated plants ( Figure 5A). The plant height was also improved in treated seeds ( Figure 5B). Similar observations were found when the plants were grown in the soil having 100 mM of salt (NaCl). The plants without treatment had the compromised appearance of root (Figure 5C), plant height, and appearance ( Figure 5D) as compared to that treated with PHM11. Furthermore, in absence of the salt when the seeds were pretreated with PHM11, the root length, root diameter, projected area, surface area, number of forks, number of links, and total chlorophyll content increased significantly ( Table 6). The presence of salt leads to the compromised growth in the absence of PHM11 treatment, which was represented by the quantitative reduction in all measured parameters. Furthermore, the pretreatment of seeds with PHM11 was grown in 100 mM NaCl condition, PHM11 significantly enhanced the root length (Figure 6A), root diameter ( Figure 6B), projected area (Figure 6C), surface area (Figure 6C), number of forks (Figure 6D), number of links (Figure 6D), shoot length ( Figure 6E) and chlorophyll content ( Figure 6F). Similar experiments were performed with wheat ( Table 7) and maize (Table 8) as well. In case of wheat, we observed similar results. The growth of root and shoot was enhanced upon PHM11 treatment whether in the absence (Table 7 and Supplementary Figure S3A) or in the presence (Supplementary Figures S3A,B) of salt. Statistical analyses showed similar results for the root length, root diameter, projected area, surface area, number of forks, shoot length, and total chlorophyll content (Supplementary Figure S4). Similar results were obtained in the case of maize as well (Table 8 and Supplementary Figures S5, S6).

CONCLUSION
The study exhibits the presence of many BGCs in the genome of E. profundum PHM11 which contributes toward environmental suitability through different metabolic pathways. PHM11 possesses biosynthetic pathways for several secondary metabolites like NRP, PK, dipeptide aldehydes, AMP, and carotenoid. Some putative bacteriocins with antibacterial activities were also identified. Besides, these varieties of proteins with the SPs as well as TMH which have a role in the transport and signaling were identified in the genome. The analysis of differential expression of genes (DEGs) showed the upregulation of some transporters like Na + /H + antiporter (NhaC). These transporters and antiporters help them to sustain the high salt conditions. These transporters possess the TMH which form the channel through the membrane to transport the solutes across the membrane. In addition to the environmental suitability, the strain has PGP traits that remain functional under salt stress and make the strain a potential agent for the formulation for crops under salt-affected soil.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm. nih.gov/, MRSV01000001.

AUTHOR CONTRIBUTIONS
AlSr participated in conceptualization, data curation, analysis, investigation, methodology, data validation, and writing of the original draft. RS and AB have done the formal analysis, validation of the data, and writing of the draft. AlSi performed the antagonistic activity of the strain. AS isolated the strain and performed genomic DNA isolation and basic characterization. SD performed the growth promotion experiments. PT and AnSr performed the preliminary analysis of results. PK, HC, and AnSa gave critical inputs. All authors contributed to the article and approved the submitted version.