Skip to main content

Genomic and physiological traits provide insights into ecological niche adaptations of mangrove endophytic Streptomyces parvulus VCCM 22513

Abstract

Purpose

Endophytic Streptomyces parvulus VCCM 22513 isolated from Bruguiera gymnorrhiza in Quang Ninh mangrove forest, northern Vietnam showed abiotic stress tolerance consisting of antioxidant, salt-tolerant, and aromatic-compound degrading activities. The goal of this study was to shed light on genomic bases rendering mangrove endophytic S. parvulus more resilient to environmental stressors.

Methods

Phenotypic analysis including antioxidant activities, hydrogen peroxide and sodium chloride resistance, and aromatic compound utilization were evaluated. The genome of strain VCCM 22513 was sequenced using Illumina Miseq sequencing platform and assembled using SPAdes.

Results

Out of 15 endophytic actinomycetes associated with B. gymnorrhiza in Quang Ninh mangrove, northern Vietnam, VCCM 22513 extract showed remarkable antioxidant activities through (1,1-diphenyl-2-picrylhydrazyl) DPPH and superoxide radical scavenging assays of 72.1 ± 0.04% and 38.3 ± 0.16% at 1.6 mg/ml, respectively. The genome consists of a 7,688,855 bp linear chromosome, 6782 protein-coding sequences, and 68 tRNAs. Genomic analysis identified strain VCCM 22513 as Streptomyces parvulus and confirmed a highly conserved core genome and stability of S. parvulus under natural selection. Genome mining revealed the presence of genetic determinants involved in mycothiol and ergothioneine biosynthesis (26 genes), oxidative stress resistance (43 genes), osmoadaptation (87 genes), heat and cold stress (34 genes), aromatic compound degradation (55 genes). Further genome-wide comparison between S. parvulus VCCM 22513 and 11 Streptomyces genomes showed that VCCM 22513 possesses significantly higher copies of genes involved in mycothiol and ergothioneine biosynthesis. In support of this finding, the strain exhibited much resistance to 0.6–1.0 M H2O2 and 6% (w/v) NaCl as compared to Streptomyces cavourensis YBQ59 isolated from Cinnamomum cassia Prels. In addition, the complete pathways for degradation of aromatic compounds including protocatechuate, gentisate, 4-hydroxyphenylpyruvate, cinnamate, 3-phenylpropionate, and styrene were only identified in the genome of VCCM 22513.

Conclusions

The present study revealed for the first time adaptive responses of mangrove endophytic S. parvulus VCCM 22513 to survive in hostile environment. The information shown here provided better understanding of underlying mechanisms related to adaptation and partially plant-microbe interaction of Streptomyces associated with mangrove plants.

Introduction

Existing along the intertidal zone between terrestrial and the sea, mangroves represent distinct ecological niches, mainly found in Southeast Asia. The natural environment for mangrove plants is considered a complex set of abiotic and biotic stress factors such as salinity, oxygen deficiency, flooding, UV-B, and water logging, which trigger the excessive production of reactive oxygen species (ROS) leading to irreversible cellular damages (Ali et al. 2017; Asaeda and Barnuevo 2019). ROS are constantly generated as toxic by-products of aerobic metabolism including superoxide anion, hydroxyl radical, and hydrogen peroxide (Tan et al. 2017; Kemung et al. 2020). In this context, mangrove plants such as Bruguiera gymnorrhiza have developed mutualistic relationship with endophytic bacteria that provide various enzymatic and non-enzymatic antioxidants to withstand unfavorable environmental conditions (Afzal et al. 2019). Previous studies proved that the utilization of endophytic bacteria can favor host plants by contributing to salt tolerance and ROS scavenging (Alkio et al. 2005; Ali et al. 2017). Endophytic bacteria isolated from mangrove plant Avicennia marina were shown to neutralize oxidative burst under salinity stress by modulating antioxidants including superoxide dismutase, catalase, polyphenol oxidases, and peroxidase (Ali et al. 2017). Therefore, it is valuable to explore strategies for environmental niche adaptations of endophytic bacteria present in mangrove plants.

Actinobacteria stand out as the most ecological and biotechnological important prokaryotes. They are able to synthesize various bioactive secondary metabolites such as antiviral, antibacterial, antioxidant, and anticancer compounds, accounting for two-third of all known antibiotics (Afzal et al. 2019; Quach et al. 2021). Of these actinobacteria, members of the genus Streptomyces produced 80% of the total actinobacterial metabolites, leading to an increasing interest in exploitation of Streptomyces, especially those from terrestrial sources in the last decades (Maiti et al. 2020). However, the ability to find novel metabolites with biological activities recently reaches a stagnant point, thus attention has been shifted towards Streptomyces from unexplored area such as the mangrove environment. Upon adaptations to diverse environmental stressors, changes in the metabolic pathways of mangrove Streptomyces species might lead to the capability to synthesize diverse and novel secondary metabolites. A notable illustration of this potential in mangrove endophytic Streptomyces is the discovery of a novel metabolite, xiamycin A, exhibiting selective anti-HIV activity (Ding et al. 2010). Given that oxidative stress is a condition characterized by elevated amounts of free radicals contributing to cancer initiation and progression (Reuter et al. 2010), natural antioxidants produced by Streptomyces effectively scavenge excessive free radicals and protect cells from severe damages during oxidative burst (Law et al. 2017; Kemung et al. 2020). Most Streptomyces spp. capable of producing antioxidant have mainly been isolated from mangrove soils and sediments such as Streptomyces mangrovisoli, Streptomyces antioxidans, Streptomyces colonosanans, and Streptomyces parvulus (Law et al. 2017; Hu et al. 2018). This readily prompted the study of endophytic Streptomyces derived from mangrove plants.

To date, there are growing interests in exploring bacterial genomes to decipher the genetic bases for environmental niche adaptations and secondary metabolite biosynthesis potential. The rapid development of whole genome sequencing technology has improved data quality and analysis of the bacterial genomes, leading to an increasing number of bacterial genomes available publicly (Tian et al. 2016; Sun et al. 2018). Genome-wide comparison revealed that despite being closely-related to their terrestrial counterparts including Streptomyces albus J1074 and Streptomyces flavovirens ATCC 33331, marine sponge-associated Streptomyces sp. SM17 and SM18 were shown to produce different secondary metabolites and comprise 29 unique environmental niche adaptation genes (Almeida et al. 2019). For adaptation to an endophytic lifestyle, Kitasatospora sp. SUK42 has a reduced genome size due to the burden of encoding unnecessary metabolic pathways (Zin et al. 2021). Besides, a previous study proved that presence of genetic determinants involved in horizontal gene transfer, mobile genetic elements, and aromatic compound degradation was recognised as an adaptation mechanism of Azoarcus sp. CIB to its endophytic lifestyles (Martín-Moldes et al. 2015).

Quang Ninh mangrove forests locating in northern Vietnam is a rich mangrove flora consisting of 5 plant species that have been unexplored for its potent endophytes (Hanh et al. 2018; Dang et al. 2020). An earlier study using genome mining demonstrated S. parvulus 03 genome derived from mangrove plant Kandelia candel was able to secrete melanin and desferrioxamine B (Hu et al. 2018), however environmental niche adaptation of this species remained unknown. This study aims to explore the phylogenomic positions of S. parvulus VCCM 22513 and shed light on the genetic determinants rendering it more resilient to model environmental stressors (salt, H2O2, and aromatic compounds). This is the first report deciphering underlying mechanisms related to multiple stress resistance and aromatic catabolism of mangrove endophytic S. parvulus VCCM 22513 at the genomic and phenotypic level.

Materials and methods

Isolation and preservation of endophytic actinomycetes

Actinomycetes were isolated from healthy mangrove plant Bruguiera gymnorrhiza collected at Quang Ninh mangrove forest, Northern Vietnam (21.0064° N, 107.2925° E). After sampling, the roots, stems, and leaves were placed in sterile plastic bags, transported to the laboratory within 48 h. The samples were washed with running tap water and distilled water. To eliminate unwanted microorganisms, surface sterilization procedure was carried out as described previously (Vu et al. 2020). In brief, the leaves, stems and roots were cut into small segments (1–2 cm) followed by surface sterilization via successive immersion into 15% (v/v) NaClO for 5 min, 2.0% (w/v) Na2S2O3 for 2 min, washed three times with sterile water. After that, the resulting samples were treated with 70% ethanol for 7 min, washed three times with sterile water, then dried under laminar flow conditions. After surface sterilization, the samples were frozen at – 80 °C for 2 weeks, and then spread onto 8 selective isolation media including humic acid-vitamin B, raffinose-histidine, tap water-yeast extract, International Streptomyces Project 5 (ISP5), trehalose-proline, sodium succinate-asparagine, starch, and citric acid agars supplemented with 50 mg/ml nystatin, 25 mg/ml K2Cr2O7, and 25 mg/ml nalidixic acid as described previously (Musa et al. 2020; Vu et al. 2020). Once observed, actinomycete colonies were streaked out several times on International Streptomyces Project (ISP) 2 medium to get pure isolates (Shirling and Gottlieb 1966). The pure cultures were stored in 20% glycerol at – 80 °C.

The 16S rRNA gene of each isolate was amplified using the universal primers 27F and 1492R as described previously (Vu et al. 2020). The resulting PCR products were sent to First BASE Laboratories Sdn. Bhd. (Malaysia) for Sanger sequencing. The obtained sequences were trimmed and compared with representative 16S rRNA gene sequences of related type strains retrieved from the GenBank database (NCBI).

Screening of antioxidant activities

All isolates were inoculated on YIM38 medium (glucose 4 g; yeast extract powder 4 g; malt extract powder 5 g; thiamine-HCl, riboflavin, niacin, pyridoxine-HCl, inositol, calcium pantothenate, and p-aminoben-zoic acid, each 0.5 mg; biotin 0.25 mg; pH 7.0 per liter) for 10 days at 30 °C with shaking at 200rpm (Nguyen et al. 2019). The biomass was removed by centrifugation while the supernatant was subjected to freeze-drying process. Extraction of freeze-dried samples was carried out using absolute methanol that underwent subsequent evaporation at 40 °C (Tan et al. 2017). The resulting extract was dissolved in dimethyl sulfoxide (DMSO) for antioxidant assays.

In the superoxide radical scavenging assay, 4-ml mixtures containing 20 μM phenazine methosulfate, 160 μM NADH, 55 μM nitro blue tetrazolium (NBT), the extract (0.2–1.6 mg/ml), and 0.2 M sodium phosphate buffer (pH 7.3) were incubated at room temperature for 10 min. Superoxide radical scavenging activity was calculated based on the absorbance measured at 560 nm as reported previously (Rajoka et al. 2019). For the 1,1-diphenyl-2-picrylhydrazyl (DPPH) assay, about 0.2 ml of extract solution (0.2; 0.4; 0.8; 1.6 mg/ml) in methanol was reacted with 0.2 ml of 0.1 mM DPPH solution followed by the addition of 2.0 ml of deionized water. The reaction was incubated in the dark for 30 min at room temperature. DPPH radical scavenging activity of VCCM 22513 extract was determined using absorbance measured at 517 nm as previously described (Liu et al. 2010; Vu et al. 2021). In the reducing power assay, the mixtures of 1.5 ml of the crude extract (0.2, 0.4, 0.8, 1.6 mg/ml), 1.5 ml of 0.2 M sodium phosphate buffer (pH 7.3), and 1.5 ml of 1% (w/v) K3Fe(CN)6 were incubated at 50 °C for 25 min. The reactions were terminated by adding 1.5 ml of 12% (w/v) trichloroacetic acid, followed by centrifugation. The supernatants were reacted with 0.5 ml of 0.2% (w/v) FeCl3 and absorbance was measured at 700 nm (Rajoka et al. 2019). Ascorbic acid was used as the positive control for all assays.

Morphological and physiological characteristics of VCCM 22513

The bioactive isolate VCCM 22513 was cultivated on ISP2-7 agar media at 30 °C for 10 days to evaluate morphological characteristics such as aerial and substrate mycelium color as well as pigment production as described previously (Quach et al. 2021). Microscopic cell morphology was observed using a conventional microscope (Carl Zeiss) with 100× oil immersion objective. Physiological properties of VCCM 22513 were assessed under different conditions of pH (2.0-9.0) and temperature (16-42°C). The utilization of different carbon sources (1%, w/v) was examined on ISP9 medium. In addition, the ability to produce exoenzymes such as amylase, cellulase, chitinase, protease, and xylanase was determined following previous protocol (Kemung et al. 2020).

Growth assessment of VCCM 22513 in response to H2O2, NaCl, and aromatic compounds

In the H2O2 tolerance test, the strain VCCM 22513 was grown until exponential phase and then spread on ISP-2 agar plates. Sterile paper disks soaked in 10 μl of hydrogen peroxide (H2O2; 0.6, 0.8, 1.0, and 2.0 M) were deposited on the prepared plates (Undabarrena et al. 2017). Growth inhibition diameters were determined after 4 days of incubation at 30 °C. For the NaCl tolerance test, the strain VCCM 22513 was plated on ISP-2 agar plates with 0, 2, 4, and 6% (w/v) NaCl (Zin et al. 2021). Bacterial growth was monitored after 16 days of incubation at 30 °C. For aromatic compound degradation experiments, the strain VCCM 22513 was inoculated in CM basal medium (l-asparagine 0.5 g; K2HPO4 0.5 g; MgSO4·7H2O 0.2 g; FeSO4·7H2O 0.01 g; casein hydrolysate 0.05 g; pH 7.0 per liter) supplemented with gentisate (60 μM), cinnamate (80 μM), styrene (80 μM), or protocatechuate (60 μM) (Khalil et al. 2019). The CM basal medium without aromatic compounds was used as a negative control. The growth of strain VCCM 22513 was evaluated after 8 days at 30 °C based on the comparison to the control. All experiments were performed in triplicate.

Genome sequencing, assembly, and annotation of the S. parvulus VCCM 22513

Genomic DNA of strain VCCM 22513 was extracted using G-spin™ Total DNA Extraction Mini Kit (Intron Bio, Korea) according to the manufacturer’s protocol. Whole genome of VCCM 22513 was sequenced with Illumina Miseq sequencing platform (Illumina, California, USA). Quality control and read trimming were performed by FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc) and Trimmomatic 3.0, respectively (Bolger et al. 2014). De novo assembly with default parameters and k-mer = 21, 33, 55, 77 was performed using SPAdes 3.15 (Bankevich et al. 2012). Benchmarking Universal Single-Copy Orthologous (BUSCO) 3 (https://gitlab.com/ezlab/busco) was used to evaluate the completeness of the assembled genome. Genome features were analyzed using 2 pipelines including Prokka v1.12 and Prokaryotic Genomes Annotation Pipeline (PGAP) at NCBI (Seemann 2014, Tatusova et al. 2016). Functional annotation of the predicted coding genes was performed based on EggNOG database (Huerta-Cepas et al. 2019). Default parameters were applied for all software tools.

Genome-wide comparison

The assembled sequence of VCCM 22513 was uploaded to Type Strain Genome Server TYGS server (https://tygs.dsmz.de) for whole genome-based taxonomic analysis. Taxogenomic analyses with closely related species were implemented by calculating the overall genome-related indexes including digital DNA–DNA hybridization (dDDH), G + C difference, and distance formula d5. The 16S rRNA gene sequence was retrieved from VCCM 22513 genome and aligned with related type Streptomyces strains accessible on the GenBank databases using CLUSTAL-X software. Sequence similarities were calculated based on the EzBiocloud database (Kim et al. 2012). The phylogenetic tree was built by the neighbor-joining method using Molecular Evolutionary Genetics Analysis (MEGA) software version 7 with Kimura-2-parameter distances (Kumar et al. 2016). Streptomyces inhibens NEAU-D10T was used as an outgroup branch.

The genome of VCCM 22513 was compared to 11 Streptomyces strains available on GenBank: S. parvulus 2297 (CP015866); S. parvulus JCM 4068T (BMRX00000000); S. parvulus LP03 (JAIWPL000000000); Streptomyces cavourensis YBQ59 (QLNH00000000); Streptomyces venezuelae ATCC 10712T (CP029197); Streptomyces albus J1074 (CP004370); Streptomyces flavovirens ATCC 33331 (CP002475), Streptomyces olivaceus B-3009T (JADOEJ000000000), Streptomyces gilvigriseus MUSC26T (MLCF00000000), Streptomyces mangrovisoli 149T (LAVA00000000), Streptomyces monashensis 1JT (MLYO00000000), Streptomyces endocoffecius CA3R110T (JAERRG000000000), Streptomyces kebangsaanensis SUK12T (JUJA00000000). Genes organized into sets of logically related functional roles were clarified based on the SEED-viewer analysis of genome sequences by RASTtk (Brettin et al. 2015). The predicted protein sequences were used to predict the orthologous gene clusters among S. parvulus strains using OrthoVenn2 web server (Xu et al. 2019).

Identification of genes involved in adaptive mechanisms and accession number

Secondary metabolite biosynthetic gene clusters (BGCs) of VCCM 22513 were predicted by antiSMASH 5.1.2 with default parameters (Blin et al. 2021). Potential adaptation mechanisms of strain VCCM 22513 were determined via an extensive literature search for mechanisms involved in redox balance, oxidative stress response, osmoadaptation, heat and cold response, and aromatic compounds degradation that have been demonstrated in microorganisms, especially Streptomyces. Genes encoding for proteins of interest were predicted using the KEGG, RAST, Mapper tools and annotated using the BLAST tool of NCBI, GenBank database, and Conserved Domain Database (CDD). The questionable open reading frames were searched through BLASTp and TBLASTN tools (e-value cut-off = 1e−5) against a database of enzymes, transcriptional regulators reported in UniProtKB-Swissprot (Camacho et al. 2009; Brettin et al. 2015). An e-value cutoff < 10−10, identity > 25% and coverage > 50% were used to filter the outcome.

The 16S rRNA gene and draft genome of S. parvulus VCCM 22513 were deposited in GenBank under accession number ON026158 and JAJVLA000000000, respectively.

Results

Isolation and phenotypic characterization of VCCM 22513

Of the 15 endophytic actinomycetes associated with mangrove plant B. gymnorrhiza in Northern Vietnam, isolate MX9 showed the strongest antioxidant activities (Table S1). Indeed, MX9 extract displayed significant antioxidant activity from 10.9 ± 0.12% to 72.1 ± 0.04% of DPPH radical reduction at doses ranging from 0.2 to 1.6 mg/ml, which could be due to its hydrogen donating ability (Fig. 1A). The antioxidant activity of the extract against superoxide radical was 38.3 ± 0.16% at 1.6 mg/ml, while the reducing power was observed at a low level (Fig. 1B). The strain grew well on all ISP agar and produced whitish-grey aerial mycelium and pale-yellow substrate mycelium after 6 days at 30 °C (Table S2). Physiological analysis revealed that the growth was observed at 16–42 °C (optimum 30 °C) and at pH 2.0–9.0 (optimum pH 6.0). The cells were able to produce broad range of extracellular enzymes such as amylase, cellulase, xylanase, chitinase, protease, and grow on media supplemented with gentisate, cinnamate, styrene, and protocatechuate (Table S2). The strain MX9 was then renamed as Streptomyces sp. VCCM 22513 for further investigation based on the deposition at VAST-Culture Collection of Microorganisms (VCCM, www.vccm.vast.vn).

Fig. 1
figure 1

Antioxidant activities including DPPH and superoxide radical scavenging activities (A) and reducing power (B) determined in the VCCM 22513 extract

Genome sequence and general features of the S. parvulus VCCM 22513

Using Illumina platform, a total of 5,283,212 paired-end reads were obtained and then assembled in 62 contigs with N50 of 246,939 bp. The total size of assembled genome was 7,688,855 bp with 72.1% G + C encoding 6,782 protein-coding sequences (CDSs) and 68 tRNAs (Table S3). Streptomyces chromosomes vary from 6.28 Mb for Streptomyces cattleya NRRL 8057T to 11.94 Mb for Streptomyces bingchenggensis BCW-1. Among them, the genome size of Streptomyces sp. VCCM 22513 felt in the middle size range of Streptomyces and was quite comparable to three selected S. parvulus genomes (Fig. 2A).

Fig. 2
figure 2

Identification and genome-wide comparison between S. parvulus VCCM 22513 and its relatives within the genus Streptomyces. A Comparison of genome size and number of CDS of 18 Streptomyces genomes with the focus on S. parvulus. B Whole genome-based phylogenetic classification indicating Streptomyces sp. VCCM 22513 as a member of S. parvulus. C SEED Sub-systems predicted in S. parvulus VCCM 22513 and 3 other S. parvulus genomes using RASTtk. D Venn diagram showing the number of shared and unique proteins between S. parvulus VCCM 22513 and other S. parvulus strains

To identify VCCM 22513 at species level, the genome-to-genome distance calculator web server was used to calculate the dDDH between VCCM 22513 genome and Streptomyces genomes. Streptomyces sp. VCCM 22513 was found to display the highest similarity to S. parvulus JCM 4068T represented by 94.7% identity dDDH and 0.02% G + C difference values (Fig. 2B). Furthermore, strain VCCM 22513 exhibited 16S rRNA gene sequence similarity values of 100% and 99.9% with S. parvulus NBRC 13193T and S. parvulus JCM 4068T, respectively. Phylogenetic analysis showed that the strain VCCM 22513 only formed a well-delineated subclade with S. parvulus species (Fig. S1). Taken together, the strain was identified as S. parvulus VCCM 22513.

With the aim to study the functional categories of 6,782 coding sequences of S. parvulus VCCM 22513, this strain and 3 other S. parvulus including 2297, JCM 4068T from soils and LP03 from lichen were subjected to SEED using RASTtk where genes were organized into sets of logically related functional roles (Overbeek et al. 2005). The three most abundant categories in S. parvulus VCCM 22513 genome were amino acids and derivatives (396 genes), carbohydrates (323 genes), and protein metabolism (225 genes). Of note, VCCM 22513 encoded 97.7% of the average number of subsystems per S. parvulus genome (Fig. 2C). Subsystems involved in Nitrogen Metabolism (101.3%) and Metabolism of Aromatic Compounds (105.6%) were over-represented with respect to the average. The distribution of subsystems among S. parvulus species is highly similar implying that isolation source does not lead to the loss of the genetic materials. In support of this result, all four species shared a large core-genome of 6162 COGs and only 2 unique COGs annotated as putative proteins were found in the VCCM 22513 genome using OrthoVenn (Fig. 2D).

Maintenance of redox balance

Low-molecular weight (LMW) thiols are required for maintaining the reduced cytoplasm of bacteria, which further protect bacteria from the hostile oxidative environments (Tung et al. 2018). In the VCCM 22513 genome, 18 genes involved in the biosynthesis of low-molecular weight thiols such as mycothiol were identified (Table S4). In the first step, the glycosyltransferase MshA (orf_549, orf_4674, orf_5830, orf_6700) conjugates myo-inositol and N-acetyl glucosamine to N-acetyl glucosamine myo-inositol, followed by deacetylation by metal-dependent deacetylase MshB (orf_1123, orf_6108) to yield glucosamine inositol (Fig. 3A). In the third step, Cys ligase MshC (orf_1437, orf_1645) conjugates Cys to produce Cys-GlcN-Ins. The final step involved acetylation of the Cys amino group by 10 copies of mycothiol acetyltransferase MshD to generate mycothiol. Similar to mycothiol, the complete biosynthetic pathway of ergothioneine required the presence of 5 main genes including egtA, egtB, egtC, egtD, egtE which were found in the VCCM 22513 genome (Fig. 3B). The four genes egtA, egtB, egtC, and egtD are clustered in an operon, and egtE is separately encoded in a distant locus, which are different to Mycobacterium smegmatis MC2 155 and S. coelicolor A3(2) (Seebeck 2010; Nakajima et al. 2015). Surprisingly, S. parvulus VCCM 22513 also contained two copies of gshA (orf_1595, orf_4097) encoding glutamate-cysteine ligases that might promote ergothioneine biosynthesis through the production of precursor γ-glutamylcysteine.

Fig. 3
figure 3

The biosynthesis of mycothiol and ergothioneine in S. parvulus VCCM 22513. The proposed biosynthetic pathways of mycothiol (A) and ergothioneine (B) in S. parvulus VCCM 22513. Abbreviations: GlcNAc, N-acetyl glucosamine; Ins, myo-inositol; His, histidine; Cys, cysteine. C Comparative analysis of genetic determinants dedicating to mycothiol and ergothioneine biosynthesis between mangrove endophytic S. parvulus VCCM 22513 and 11 Streptomyces

Comparative genomic analysis of VCCM 22513 with 11 mangrove-, plant- and soil-derived Streptomyces spp. indicated that genetic determinants related to the biosynthesis of low-molecular weight thiols are present in all genomes but with a different number of copies (Fig. 3C). Despite sharing a highly conserved core genome, mshD genes in VCCM 22513 (10 copies) was higher than S. parvulus JCM 4068T (9 copies) isolated from soil. Gene egtC was duplicated in the VCCM 22513 genome, which was only similar to S. olivaceus B-3009T, S. venezuelae ATCC 10712T, and S. parvulus JCM 4068T. Of note, S. parvulus VCCM 22513 and S. olivaceus B-3009T recovered from mangrove plant Kandelia candel (L.) Druce showed a similar profile of genes involved in mycothiol and ergothioneine biosynthesis, which were particularly high in gene copies as compared to 10 mangrove and terrestrial genomes selected. Interestingly, mangrove soil-derived S. gilvigriseus MUSC26T lacked the redox protective ergothioneine (Fig. 3C).

Primary responses of S. parvulus VCCM 22513 to oxidative stress

Given that the presence of oxygen, cycles of inundation, salinity fluctuations, and aromatic compounds can indirectly cause oxidative stress such as ROS and reactive electrophilic species (RES) to endophytes (Ali et al. 2017; Asaeda and Barnuevo 2019), mechanisms to prevent cellular damages might be important components of mangrove endophytic S. parvulus VCCM 22513. Investigation of genes involved in primary oxidative stress response resulted in the prediction of 43 corresponding genes (Fig. 4A, Table S5). Under oxidative stress, mycothiol is oxidized to mycothiol disulfide, which is recycled by the mycothiol disulfide reductase Mtr (orf_5308). Parallelly, mycothiol forms mixed disulfides with redox proteins to protect them from irreversible over-oxidation, which subsequently is catalyzed by mycoredoxin Mrx1 (orf_1066) to release reduced proteins. In addition, thioredoxin system found in the VCCM 22513 genome was encoded by 10 putative thioredoxin and thioredoxin reductase genes which play an important role in reducing disulfide bonds of cysteine-containing proteins. Furthermore, genes categorized as superoxide dismutase (2 genes), alkyl hydroperoxide reductase (7 genes), peroxidase and peroxiredoxin (6 genes), catalase (4 genes), and ROS resistance protein (5 genes) were predicted, in which several are duplicated and multiplicated (Fig. 4A). In line with these results, global oxidative stress responses of strain VCCM 22513 might be regulated by 6 universal transcription factors such as redox-sensitive transcriptional regulator SoxR, regulatory protein SoxS, hydrogen peroxide-sensing transcriptional regulator OxyR, peroxide-responsive repressor PerR, organic hydroperoxide resistance transcriptional regulator OhrR, and ferric uptake regulator FurA (Table S5). Further comparative genomics showed that the number of genes encoding for enzymatic scavengers including ROS resistance protein, redox transcriptional regulator, catalase, peroxidase, and peroxiredoxin found in S. parvulus VCCM 22513 were similar to those in the genomes of terrestrial and mangrove Streptomyces strains (Fig. 4A). The differences were that S. parvulus VCCM 22513 possessed an extra copy of alkyl hydroperoxide reductase ahp and bromoperoxidase-catalase bca in comparison to other mangrove genomes including S. olivaceus B-3009T, S. gilvigriseus MUSC26T, S. mangrovisoli 149T, and S. monashensis 1JT.

Fig. 4
figure 4

Functional analysis of oxidative and osmotic stress tolerance of S. parvulus VCCM 22513 at genomic and phenotypic levels. Comparative analysis of genes involved in oxidative stress resistance (A) and osmotic stress (C) present in S. parvulus VCCM 22513 and 11 Streptomyces genomes. Functional response of S. parvulus VCCM 22513 and S. cavourensis YBQ59 to increasing concentrations of hydrogen peroxide (H2O2) (B) and 0–6% (w/v) NaCl (D)

Besides enzymatic scavengers, S. parvulus VCCM 22513 also utilizes secondary metabolites responsible for neutralizing oxidative damages. Using AntiSMASH, in silico prediction of secondary metabolite biosynthetic gene clusters identified a total of 25 clusters with similarity to known clusters of ≥ 16%. Of those, 4 biosynthetic gene clusters including desferrioxamin, melanin, coelichelin, and isorenieratene were reported to be antioxidant and iron chelating agents (Table S6) (Hu et al. 2018; Chen et al. 2019; Williams et al. 2019; Wibowo et al. 2022). This was in agreement with the DPPH, reducing power, superoxide radical scavenging activities observed, indicating the ability of strain VCCM 22513 to produce 4 antioxidant compounds.

To correlate with the genetic content of redox balance and oxidative stress properties, the response of strain VCCM 22513 to H2O2 concentrations was tested and compared to S. cavourensis YBQ59 from Cinnamomum cassia Prels (Tung et al. 2018). Strain VCCM 22513 showed smaller growth inhibition zones in response to 0.6, 0.8, and 1.0 M H2O2 as compared to S. cavourensis YBQ59 (Fig. 4B), indicating great capability of mangrove endophytic VCCM 22513 in neutralizing H2O2 toxicity. At 2.0 M H2O2, no significant difference was observed (Fig. S2).

Distribution of genes involved in osmoadaptation

To protect themselves from variations in external osmotic pressure, actinomycetes have to accumulate and release solutes by regulating water fluxes and maintaining cellular homeostasis (Yaakop et al. 2016; Undabarrena et al. 2017). In the present study, 87 genes involved in the defensive responses to osmotic stress were identified in the VCCM 22513 genome (Fig. 4C). Given that proton pump is crucial to produce ATP via proton-motive force and to maintain a proton gradient under high-salt stress (Sun et al. 2018), NADH quinone oxidoreductase genes such as nuo and ndh encoding proton pumps were found in 2 operon-like structures including nuo (14 genes) and ndhC-nuoN (10 genes), making them the most abundant groups identified in the genome (Table S7). In contrast, the partial nuo operon (8 genes) was present in the genome of marine actinomycetes such as Streptomyces sp. SM17 and Streptomyces sp. SM18 (Almeida et al. 2019).

Additionally, S. parvulus VCCM 22513 might also be able to synthesize ectoine and its derivative, 5-hydroxyectoine, as compatible solutes to maintain the osmotic equilibrium in response to salt stress. In line with having the ectABCD operon, transporters for ectoine encoded by ehuABCD were revealed, supporting the ability to tolerate saline stress of the mangrove endophytic strain VCCM 22513 (Table S7). Conversely, in a low-salt stress scenario, biosynthesis and uptake of compatible solutes including proline (15 genes), glutamate (17 genes), glutamine (14 genes), and glycine betaine/choline (8 genes) are likely employed as another essential defense mechanism of VCCM 22513. In halophilic bacteria, proline is thought to be the abundant compatible solute upon increasing osmotic stress. Meanwhile, the proline/betaine transporter ProP is not specific to its substrates leading to its versatility to transport various solutes such as proline, glycine betaine, proline betaine, ectoine, and other compounds (Saum and Müller 2007; Burg and Ferraris 2008).

Comparative genomics analysis revealed that 5 copies of glutamine transport ATP-binding glnQ related to glutamine synthesis were determined in 3 soil-derived genomes, while only 3–4 copies were present in the mangrove genomes (Fig. 4C). For glutamate synthesis, the genetic determinants found in the VCCM 22513 genome were quite comparable with mangrove endophytic S. olivaceus B-3009T, except for gdhI encoding dehydrogenase reductase. It is worthy to note that this gene was present in the terrestrial genomes, but not in the mangrove-derived genomes. In addition, further analysis revealed 5 copies of a potential proline/betaine transporter proP in the VCCM 22513 genome, unlike terrestrial genomes such as S. albus J1074, S. flavovirens ATCC 33331, and S. venezuelae ATCC 10712T in which proP exists in only three copies (Fig. 4C). In support of these genetic findings and to explore the potential to resist osmotic stress, the strain VCCM 22513 and S. cavourensis YBQ59 were plated on ISP-2 agar supplemented with 0–6% NaCl (w/v). The presence of 2% NaCl did not affect the growth and sporulation of strains VCCM 22513 and YBQ59. However, exposure to 6% NaCl resulted in a significant growth delay of YBQ59, but not VCCM 22513 after 16 days (Fig. 4D). These results indicated that mangrove endophytic VCCM 22513 was notably more tolerant to NaCl than the terrestrial endophytic YBQ59.

Mechanisms of dealing with heat and cold stress

Despite being protected by the host B. gymnorrhiza, mangrove endophytic S. parvulus VCCM 22513 might equip with the genetic determinants to adapt to sudden temperature shift. Genomic analysis confirmed the presence of 34 predicted genes may play a role in such tolerance (Table S8). For heat shock response, the dnaK operon containing dnaK, grpE, dnaJ, and hspR was identified, which is conserved across Actinobacteria such as S. coelicolor, Corynebacterium glutamicum, and Mycobacterium bovis (Segal and Ron 1996; Ventura et al. 2005). Heat shock repressor HrcA may regulate the expression of dnaJ, clpB, groSL, grpE, htpG and itself upon exposure to changes in temperatures. The predicted proteases include 2 Lon (cytoplasmic), 3 HtpX (membrane), and 1 DegP (periplasmic) which play a vital part in eliminating protein aggregates or damaged proteins as a consequence of heat shock. Regarding cold-induced stress, 3 major cold-shock proteins Csp, Cps, and ScoF are present in the VCCM 22513 genome (Table S8). In addition, the cells can also be protected by transcription termination proteins NusAB and DEAD-box ATP-dependent RNA helicase CshA that are expressed during temperature downshift. These results are in agreement with the wide range of temperatures in Northern Vietnam which varies from 5 to 38 °C.

Prediction of metabolic pathways for aromatic compound utilization

Being able to resist aromatic compound stress is important in helping bacteria survive within host plants (Martín-Moldes et al. 2015; de Melo et al. 2018). Phenotypic analysis displayed that S. parvulus VCCM 22513 was able to grow in the presence of gentisate, cinnamate, styrene, and protocatechuate (Fig. S3). In agreement with these results, genome mining revealed the gene cluster encoding protocatechuate meta-cleavage pathway distributed along the chromosome of VCCM 22513, including β-ketoadipate succinyl-CoA transferase subunit α and β (pcaI, pcaJ), protocatechuate 3,4-dioxygenase subunit α and β (pcaH, pcaG), 3-carboxy-cis,cis-muconate cycloisomerase (pcaB), muconolactone D-isomerase (catC) and 3-oxoadipate enol-lactonase (catD), and 2 regulatory proteins (nodD, hosA) (Fig. 5A). It should be noted that the gene catD encoding 3-oxoadipate enol-lactonase that converts β-ketoadipate-enol-lactone to β-ketoadipate was found in 4 copies along the genome of VCCM 22513 and terrestrial endophytic S. endocoffecius CA3R110T and S. kebangsaanensis SUK12T. Surprisingly, all reference genomes contained incomplete protocatechuate degradation pathway due to the lack of some important genes such as catC, pcaHG, pacB, and pcaR (Fig. 5B).

Fig. 5
figure 5

Complete pathways of aromatic compound degradation identified in S. parvulus VCCM 22513. A Six peripheral catabolic pathways for the aromatic compound degradation and their corresponding metabolites predicted in the S. parvulus VCCM 22513 genome. Ring-cleavage product structures of central aromatic catabolic pathways are shown inside the circle, which are fueled into tricarboxylic acid (TCA) cycle. B Distribution of genes involved aromatic compound degradation in S. parvulus VCCM 22513 and 11 Streptomyces strains

In contrast to protocatechuate, genes involved in gentisate catabolism are not clustered in an operon. The identified gentisate pathway consists of 4 main genes including transcriptional regulator (kdgR), gentisate 1,2-dioxygenase (sdgD), fumarylacetoacetate hydrolase (nagK), and maleylpyruvate isomerase (nagL) (Table S9). The SdgD is a key component in oxygenating the benzen ring of gentisate to form 3-maleylpyruvate, which is further metabolized by 6 NagL and 4 NagK enzymes into TCA cycle intermediates such as fumarate and pyruvate. Additionally, the homogentisic pathway, a central pathway for catabolism of phenylalanine and tyrosine, was identified with 4-hydroxyphenylpyruvate as the precursor of homogentisic acid which is converted into the final product maleacetoacetate by 4-hydroxyphenylpyruvate dioxygenases (hpd) and homogentisate 1,2-dioxygenase (hmgA).

Furthermore, genome annotation showed 15 genes related to peripheral pathways which diverged into the styrene upper and lower catabolic pathways (Fig. 5A). Moreover, mhp and hca genes encoding the putative peripheral pathways for the degradation of 3-phenylpropionate as well as cinnamate were found distributed along the VCCM 22513 genome. The styrene upper catabolic pathway consisted of styA and styD genes, leading to formation of phenylacetate which is then converted into metabolites needed for the TCA cycle by proteins encoded by paaABCDEFGHIJKX genes present in the styrene lower catabolic pathway (Table S9). Out of 12 compared genomes, the complete styrene degradation pathway was only found in the genome of strains VCCM 22513, B-3009T, 149T, and JCM 4068T (Fig. 5B). Different to the protocatechuate and styrene pathways, catabolic pathways for gentisate, 4-hydroxyphenylpyruvate, 3-phenylpropionate, and cinnamate were reported to be complete in most compared genomes.

Discussion

The mangroves are a unique ecosystem and characterized by the presence of a few plant species such as B. gymnorrhiza that is adapted to higher-than-normal levels of salinity, fluctuations in tidal gradients, temperature and other oxidative stresses. To survive in harsh environments, mangrove endophytic Streptomyces might be well-equipped with a number of genetic determinants encoding enzymatic, non-enzymatic, and secondary metabolite agents (Tan et al. 2017; Jiang et al. 2018). In this study, a robust adaptive response to environmental stressors of S. parvulus VCCM 22513 associated with B. gymnorrhiza was deciphered at the phenotypic and genetic levels to provide better understanding of adaptation of actinobacteria to mangrove endophytic lifestyle.

In natural habitats, actinobacteria have evolved to adapt to the environmental niches and compete with multispecies communities for resources and space, leading to metabolic and physiological differences (Almeida et al. 2019; Zerouki et al. 2021; Zin et al. 2021). For example, the growth of marine sponge-derived S. albus SM17 was improved in high salinity medium as opposed to soil-derived S. albus J1074, mainly due to a pool of 29 potential environmental niche adaptation genes (Almeida et al. 2019). However, the genome size of mangrove endophytic S. parvulus VCCM 22513 was quite similar to its terrestrial counterparts including JCM 4068T, 2297, and LP03 (~ 7.69 Mb compared to average of ~ 7.52 Mb). Moreover, the subsystems of VCCM 22513 were nearly identical with S. parvulus reference sequences in the functional categories. In addition, the mangrove endophytic VCCM 22513 shared with other strains a large core-genome that harbored all protein families encoding functions related to the basic biology and phenotypes of S. parvulus species. This was in contrast to actinobacterial strain Kitasatospora sp. SUK42 from Antidesma neurocarpum Miq that adapted to endophytic life-style through genomic reduction (Zin et al. 2021). Hence, it seems that environmental niche adaptation of S. parvulus may not be under selective pressures. As this is the first mangrove endophytic S. parvulus sequenced so far, more genome sequences of S. parvulus from different environmental niches are further required to confirm whether S. parvulus has open genome as a result of natural selection.

Genome mining demonstrated that S. parvulus VCCM 22513 effectively utilized both mycothiol and ergothioneine to keep the cytoplasm in a highly reducing state. Under oxidative stress conditions, LMW thiols form mixed disulfides with redox regulation and thiol-protection of proteins termed as S-thiolation, which plays an important role in protecting active site cysteine residues from irreversible oxidation (Tung et al. 2019). In parallel, oxidative stress is neutralized by the oxidation of LMW thiols that is later converted to their reduced states by specific enzymes (Tung et al. 2018). In soil-derived strains S. coelicolor A3(2) and Nocardia asteroides, one copy of mshA, mshB, mshC, and mshD resulted in the production of 3–6 μM mycothiol (Park et al. 2006; Newton et al. 2008). Interestingly, the marine actinomycete Streptomycetaceae CNQ530 was shown to produce approximately 50 μM mycothiol but its genome has not been sequenced yet (Newton et al. 2008). Using comparative genomics, mangrove endophytic VCCM 22513 and S. olivaceus B-3009T were found to have a higher number of msh and egt genes than those of mangrove and terrestrial genomes. This inferred that maintenance of a highly reducing cytoplasm could be an important protective strategy of actinobacterial strains against the hostile environment of mangrove plants. Of note, some plants might also uptake extracellular ergothioneine as an additional antioxidant through symbiotic relationship with beneficial microbes (Samples and Balunas 2020). In support of the genomic findings, strain VCCM 22513 was notably more resistant to H2O2-induced oxidative stress as compared to endophytic S. cavourensis YBQ59. This led to a possibility that a highly reducing cytoplasm may be amongst the main factors that makes S. parvulus VCCM 22513 more resilient to oxidative stress triggerred by B. gymnorrhiza. To our knowledge, this is the first genome mining report of S. parvulus describing biosynthetic pathways of mycothiol and ergothioneine.

Genes encoding ROS–scavenging antioxidant enzymes also contributed to oxidative stress resistance of strain VCCM 22513. Hydrogen peroxide-induced toxicity could be neutralized by VCCM 22513 due to the presence of catalases, thioredoxin system, superoxide dismutases, peroxidases, peroxiredoxins, and alkyl hydroperoxide reductases. All obtained genes are conserved across Streptomyces genera, indicating that enzymatic antioxidants are the driving force allowing Streptomyces to inhabit diverse environmental niches such as soils, lichens, and plants. The differences between VCCM 22513 and all compared genomes were the presence of additional copies of ahp and bca that catalyze the reduction of hydrogen peroxide and organic hydroperoxides to water (Jiang et al. 2019). Resistance to ROS in strain VCCM 22513 may likely correlate with single nucleotide polymorphism and other type of mutations that can improve positively the functioning of ROS-scavenging antioxidant enzymes, which is an interesting subject for future studies.

BGCs of desferrioxamin, coelichelin, melanin, and isorenieratene were present in the VCCM 22513 genome, further supporting the antioxidant properties of this strain. These compounds were previously found in the genome of Streptomyces sp. Babs14 and Osf17 isolated from the sand dunes of Sahara as an antioxidant strategy (Zerouki et al. 2021). Of note, desferrioxamin secreted by S. parvulus 03 from a mangrove plant is a specific iron-complexing agent completely preventing ROS production (Hu et al. 2018) which could contribute to the reducing power ability found in the VCCM 22513 extract. Also, a recent study proved that endophytic bacterial derived from mangrove plant Avicennia marina could protect tomato plant against the salinity-stress-induced oxidative stress via the production of peroxidases, superoxide dismutases, and catalases (Ali et al. 2017). However, the expression of BGCs under abiotic stress remains unknown. Future studies will be carried out to investigate the underlying mechanisms involved in non-enzymatic and enzymatic antioxidants of S. parvulus VCCM 22513 that help B. gymnorrhiza overcome oxidative stress under in vivo conditions.

An interesting outcome of S. parvulus VCCM 22513 genome analysis was the evidence of various genes participating in the catabolism of aromatic compounds. Since mangrove B. gymnorrhiza is a rich source of aromatic compounds due to its own production and pollutants uptake from the surrounding environments (Han et al. 2005; Naidoo and Naidoo 2016), S. parvulus VCCM 22513 must develop prorective mechanisms to resist and utilize these toxic substrates. As expected, strain VCCM 22513 grew well on the mineral medium supplemented with gentisate, cinnamate, styrene, and protocatechuate. It was in agreement with previous studies demonstrating that Streptomyces spp. were able to utilize benzaldehyde, catechol, phenylacetic acid, protocatechuic acid, and even polycyclic aromatic hydrocarbons as growth substrates (Khalil et al. 2019). Using genome-wide comparison, the complete pathways for degradation of aromatic compounds including protocatechuate, gentisate, 4-hydroxyphenylpyruvate, cinnamate, 3-phenylpropionate, and styrene were only identified in S. parvulus genomes including VCCM 22513 and JCM 4068T. In mangrove-derived S. mangrovisoli 149T and S. olivaceus sp. B-3009T, the protocatechuate degradation pathway was predicted to be incomplete due to the absence of catC. Given that ROS could be generated by aromatic compounds (Alkio et al. 2005; Kim et al. 2007), the stress proteins of Streptomyces scabies 87.22 were highly upregulated in response to an aromatic compound such as ferulic acid (Khalil et al. 2019). In addition, the absence of mycothiol resulted in an inability to degrade gentisate and 3-hydroxybenzoate (Liu et al. 2013), indicating the importance of mycothiol as a cofactor of mycothiol-dependent maleylpyruvate isomerase NagL. Thus, the presence of genes encoding catabolic pathways along with antioxidants might render S. parvulus VCCM 22513 more adaptive to mangrove endophytic lifestyle.

Conclusion

Overall, the present study provided novel insights into adaptive mechanisms of mangrove endophytic VCCM 22513 to survive under oxidative hostile environment encountered while residing inside the host plant B. gymnorrhiza. Of the 15 endophytic actinomycetes, the isolate VCCM 22513 exhibited potent antioxidant activities against DPPH, superoxide radicals, and reducing power. Moreover, the strain was strongly resistant to H2O2, NaCl and able to grow in the presence of aromatic compounds. Based on the high-quality genomic information and comparative genome analysis, the endophytic strain VCCM 22513 was identified as a member of S. parvulus species. Genome-wide comparison revealed a highly conserved core genome of S. parvulus that might not be undergoing streamlining selection to adapt to environmental niches. Notably, S. parvulus VCCM 22513 harbored a wide range of genes encoding enzymatic, non-enzymatic, and secondary metabolites to overcome multiple stressors, including oxidative stress, high salinity, dramatic temperature changes, and toxic aromatic compounds. Thus, this work paves the way for a better understanding of underlying genetic mechanisms of S. parvulus adaptation to mangrove endophytic lifestyles, while enabling the discovery of potential antioxidant agents that can possibly be engineered for their future production.

Availability of data and materials

Not applicable.

Abbreviations

DPPH:

1,1-diphenyl-2-picrylhydrazyl

ROS:

Reactive oxygen species

RES:

Reactive electrophilic species

CDSs:

Protein coding sequences

LMW:

Low-molecular weight

CDD:

Conserved Domain Database

dDDH:

Digital DNA–DNA hybridization

H2O2 :

Hydrogen peroxide

DMSO:

Dimethyl sulfoxide

NBT:

Nitro blue tetrazolium

References

  • Afzal I, Shinwari ZK, Sikandar S, Shahzad S (2019) Plant beneficial endophytic bacteria: mechanisms, diversity, host range and genetic determinants. Microbiol Res 221:36–49. https://doi.org/10.1016/j.micres.2019.02.001

    Article  CAS  PubMed  Google Scholar 

  • Ali A, Shahzad R, Khan AL, Halo BA, Al-Yahyai R, Al-Harrasi A, Al-Rawahi A, Lee I-J (2017) Endophytic bacterial diversity of Avicennia marina helps to confer resistance against salinity stress in Solanum lycopersicum. J Plant Interact 12(1):312–322. https://doi.org/10.1080/17429145.2017.1362051

  • Alkio M, Tabuchi TM, Wang X, Colón-Carmona A (2005) Stress responses to polycyclic aromatic hydrocarbons in Arabidopsis include growth inhibition and hypersensitive response-like symptoms. J Exp Bot 56(421):2983–2994. https://doi.org/10.1093/jxb/eri295

    Article  CAS  PubMed  Google Scholar 

  • Almeida EL, Carrillo Rincón AF, Jackson SA, Dobson ADW (2019) Comparative genomics of marine sponge-derived Streptomyces spp. isolates SM17 and SM18 with their closest terrestrial relatives provides novel insights into environmental niche adaptations and secondary metabolite biosynthesis potential. Front Microbiol 10:1713. https://doi.org/10.3389/fmicb.2019.01713

    Article  PubMed  PubMed Central  Google Scholar 

  • Asaeda T, Barnuevo A (2019) Oxidative stress as an indicator of niche-width preference of mangrove Rhizophora stylosa. For Ecol Manag 432:73–82. https://doi.org/10.1016/j.foreco.2018.09.015

  • Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA (2012) SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19(5):455–477. https://doi.org/10.1089/cmb.2012.0021

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Blin K, Shaw S, Kloosterman AM, Charlop-Powers Z, van Wezel GP, Medema Marnix H, Weber T (2021) antiSMASH 6.0: improving cluster detection and comparison capabilities. Nucleic Acids Res 49(W1):W29–W35. https://doi.org/10.1093/nar/gkab335

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30(15):2114–2120. https://doi.org/10.1093/bioinformatics/btu170

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Brettin T, Davis JJ, Disz T, Edwards RA, Gerdes S, Olsen GJ, Olson R, Overbeek R, Parrello B, Pusch GD, Shukla M, Thomason JA 3rd, Stevens R, Vonstein V, Wattam AR, Xia F (2015) RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci Rep 5:8365. https://doi.org/10.1038/srep08365

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Burg MB, Ferraris JD (2008) Intracellular organic osmolytes: function and regulation. J Biol Chem 283(12):7309–7313. https://doi.org/10.1074/jbc.R700042200

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL (2009) BLAST+: architecture and applications. BMC Bioinformatics 10:421. https://doi.org/10.1186/1471-2105-10-421

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Chen Y, Guo M, Yang J, Chen J, Xie B, Sun Z (2019) Potential TSPO ligand and photooxidation quencher isorenieratene from Arctic ocean Rhodococcus sp. B7740. Mar Drugs 17(6):316

    Article  CAS  Google Scholar 

  • Dang KB, Nguyen MH, Nguyen DA, Phan TTH, Giang TL, Pham HH, Nguyen TN, Tran TTV, Bui DT (2020) Coastal wetland classification with deep u-net convolutional networks and sentinel-2 imagery: a case study at the Tien Yen Estuary of Vietnam. Remote Sens 12(19):3270

    Article  Google Scholar 

  • de Melo RR, Tomazetto G, Persinoti GF, Sato HH, Ruller R, Squina FM (2018) Unraveling the cellulolytic and hemicellulolytic potential of two novel Streptomyces strains. Ann Microbiol 68(10):677–688. https://doi.org/10.1007/s13213-018-1374-7

    Article  CAS  Google Scholar 

  • Ding L, Münch J, Goerls H, Maier A, Fiebig HH, Lin WH, Hertweck C (2010) Xiamycin, a pentacyclic indolosesquiterpene with selective anti-HIV activity from a bacterial mangrove endophyte. Bioorg Med Chem Lett 20(22):6685–6687. https://doi.org/10.1016/j.bmcl.2010.09.010

    Article  CAS  PubMed  Google Scholar 

  • Han L, Huang X, Sattler I, Moellmann U, Fu H, Lin W, Grabley S (2005) New aromatic compounds from the marine mangrove Bruguiera gymnorrhiza. Planta Med 71(2):160–164. https://doi.org/10.1055/s-2005-837784

  • Hanh NH, Tuan MS, Tinh PH, Hanh NTH (2018) Characteristics of gap regeneration in mangrove forests of Dong Rui, Tien Yen, Quang Ninh. Acad J Biol 40(2):129–137. https://doi.org/10.15625/0866-7160/v40n2.13243

    Article  Google Scholar 

  • Hu D, Chen Y, Sun C, Jin T, Fan G, Liao Q, Mok KM, Lee M-YS (2018) Genome guided investigation of antibiotics producing actinomycetales strain isolated from a Macau mangrove ecosystem. Sci Rep 8(1):14271. https://doi.org/10.1038/s41598-018-32076-z

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H, Mende DR, Letunic I, Rattei T, Jensen LJ, von Mering C, Bork P (2019) eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res 47(D1):D309–d314. https://doi.org/10.1093/nar/gky1085

    Article  CAS  PubMed  Google Scholar 

  • Jiang G, Yang J, Li X, Cao Y, Liu X, Ling J, Wang H, Zhong Z, Zhu J (2019) Alkyl hydroperoxide reductase is important for oxidative stress resistance and symbiosis in Azorhizobium caulinodans. FEMS Microbiol Lett 366(3). https://doi.org/10.1093/femsle/fnz014

  • Jiang Z-K, Tuo L, Huang D-L, Osterman IA, Tyurin AP, Liu S-W, Lukyanov DA, Sergiev PV, Dontsova OA, Korshun VA, Li F-N, Sun C-H (2018) Diversity, novelty, and antimicrobial activity of endophytic actinobacteria from mangrove plants in Beilun Estuary national nature reserve of Guangxi, China. Front Microbiol 9:868–868. https://doi.org/10.3389/fmicb.2018.00868

    Article  PubMed  PubMed Central  Google Scholar 

  • Kemung HM, Tan LT-H, Chan K-G, Ser H-L, Law JW-F, Lee L-H, Goh B-H (2020) Antioxidant activities of Streptomyces sp. strain MUSC 14 from mangrove forest soil in Malaysia. Biomed Res Int 2020:6402607. https://doi.org/10.1155/2020/6402607

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Khalil M, Lerat S, Beaudoin N, Beaulieu C (2019) The plant pathogenic bacterium Streptomyces scabies degrades the aromatic components of potato periderm via the β-Ketoadipate pathway. Front Microbiol 10:2795. https://doi.org/10.3389/fmicb.2019.02795

    Article  PubMed  PubMed Central  Google Scholar 

  • Kim OS, Cho YJ, Lee K, Yoon SH, Kim M, Na H, Park SC, Jeon YS, Lee JH, Yi H, Won S, Chun J (2012) Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species. Int J Syst Evol Microbiol 62(Pt 3):716–721. https://doi.org/10.1099/ijs.0.038075-0

    Article  CAS  PubMed  Google Scholar 

  • Kim YS, Min J, Hong HN, Park JH, Park KS, Gu MB (2007) Gene expression analysis and classification of mode of toxicity of polycyclic aromatic hydrocarbons (PAHs) in Escherichia coli. Chemosphere 66(7):1243–1248. https://doi.org/10.1016/j.chemosphere.2006.07.040

  • Kumar S, Stecher G, Tamura K (2016) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol 33(7):1870–1874. https://doi.org/10.1093/molbev/msw054

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Law JW-F, Ser H-L, Duangjai A, Saokaew S, Bukhari SI, Khan TM, Ab Mutalib N-S, Chan K-G, Goh B-H, Lee L-H (2017) Streptomyces colonosanans sp. nov., a novel actinobacterium isolated from Malaysia mangrove soil exhibiting antioxidative activity and cytotoxic potential against human colon cancer cell lines. Front Microbiol 8:877. https://doi.org/10.3389/fmicb.2017.00877

  • Liu J, Luo J, Ye H, Sun Y, Lu Z, Zeng X (2010) In vitro and in vivo antioxidant activity of exopolysaccharides from endophytic bacterium Paenibacillus polymyxa EJS-3. Carbohydr Polym 82(4):1278–1283. https://doi.org/10.1016/j.carbpol.2010.07.008

  • Liu YB, Long MX, Yin YJ, Si MR, Zhang L, Lu ZQ, Wang Y, Shen XH (2013) Physiological roles of mycothiol in detoxification and tolerance to multiple poisonous chemicals in Corynebacterium glutamicum. Arch Microbiol 195(6):419–429. https://doi.org/10.1007/s00203-013-0889-3

  • Maiti PK, Das S, Sahoo P, Mandal S (2020) Streptomyces sp. SM01 isolated from Indian soil produces a novel antibiotic picolinamycin effective against multi drug resistant bacterial strains. Sci Rep 10(1):10092. https://doi.org/10.1038/s41598-020-66984-w

  • Martín-Moldes Z, Zamarro MT, Del Cerro C, Valencia A, Gómez MJ, Arcas A, Udaondo Z, García JL, Nogales J, Carmona M, Díaz E (2015) Whole-genome analysis of Azoarcus sp. strain CIB provides genetic insights to its different lifestyles and predicts novel metabolic features. Syst Appl Microbiol 38(7):462–471. https://doi.org/10.1016/j.syapm.2015.07.002

  • Musa Z, Ma J, Egamberdieva D, Abdelshafy Mohamad OA, Abaydulla G, Liu Y, Li W-J, Li L (2020) Diversity and antimicrobial potential of cultivable endophytic actinobacteria associated with the medicinal plant Thymus roseus. Front Microbiol 11:191. https://doi.org/10.3389/fmicb.2020.00191

  • Naidoo G, Naidoo K (2016) Uptake of polycyclic aromatic hydrocarbons and their cellular effects in the mangrove Bruguiera gymnorrhiza. Mar Pollut Bull 113(1):193–199. https://doi.org/10.1016/j.marpolbul.2016.09.012

  • Nakajima S, Satoh Y, Yanashima K, Matsui T, Dairi T (2015) Ergothioneine protects Streptomyces coelicolor A3(2) from oxidative stresses. J Biosci Bioeng 120(3):294–298. https://doi.org/10.1016/j.jbiosc.2015.01.013

    Article  CAS  PubMed  Google Scholar 

  • Newton GL, Jensen PR, Macmillan JB, Fenical W, Fahey RC (2008) An N-acyl homolog of mycothiol is produced in marine actinomycetes. Arch Microbiol 190(5):547–557. https://doi.org/10.1007/s00203-008-0405-3

  • Nguyen QH, Nguyen HV, Vu TH-N, Chu-Ky S, Vu TT, Hoang H, Quach NT, Bui TL, Chu HH, Khieu TN, Sarter S, Li W-J, Phi Q-T (2019) Characterization of endophytic Streptomyces griseorubens MPT42 and assessment of antimicrobial synergistic interactions of its extract and essential oil from host plant Litsea cubeba. Antibiotics 8(4):197

  • Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang H-Y, Cohoon M, de Crécy-Lagard V, Diaz N, Disz T, Edwards R, Fonstein M, Frank ED, Gerdes S, Glass EM, Goesmann A, Hanson A, Iwata-Reuyl D, Jensen R, Jamshidi N, Krause L, Kubal M, Larsen N, Linke B, McHardy AC, Meyer F, Neuweger H, Olsen G, Olson R, Osterman A, Portnoy V, Pusch GD, Rodionov DA, Rückert C, Steiner J, Stevens R, Thiele I, Vassieva O, Ye Y, Zagnitko O, Vonstein V (2005) The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res 33(17):5691–5702. https://doi.org/10.1093/nar/gki866

  • Park JH, Cha CJ, Roe JH (2006) Identification of genes for mycothiol biosynthesis in Streptomyces coelicolor A3(2). J Microbiol 44(1):121–125

    CAS  PubMed  Google Scholar 

  • Quach NT, Nguyen QH, Vu THN, Le TTH, Ta TTT, Nguyen TD, Van Doan T, Van Nguyen T, Dang TT, Nguyen XC, Chu HH, Phi QT (2021) Plant-derived bioactive compounds produced by Streptomyces variabilis LCP18 associated with Litsea cubeba (Lour.) Pers as potential target to combat human pathogenic bacteria and human cancer cell lines. Braz J Microbiol 52(3):1215–1224. https://doi.org/10.1007/s42770-021-00510-6

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Rajoka MSR, Mehwish HM, Hayat HF, Hussain N, Sarwar S, Aslam H, Nadeem A, Shi J (2019) Characterization, the antioxidant and antimicrobial activity of exopolysaccharide isolated from poultry origin Lactobacilli. Probiotics Antimicrob Proteins 11(4):1132–1142. https://doi.org/10.1007/s12602-018-9494-8

    Article  CAS  PubMed  Google Scholar 

  • Reuter S, Gupta SC, Chaturvedi MM, Aggarwal BB (2010) Oxidative stress, inflammation, and cancer: how are they linked? Free Radic Biol Med 49(11):1603–1616. https://doi.org/10.1016/j.freeradbiomed.2010.09.006

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Samples RM, Balunas MJ (2020) Bridging the gap: plant-endophyte interactions as a roadmap to understanding small-molecule communication in marine microbiomes. Chembiochem 21(19):2708–2721. https://doi.org/10.1002/cbic.202000064

    Article  CAS  PubMed  Google Scholar 

  • Saum SH, Müller V (2007) Salinity-dependent switching of osmolyte strategies in a moderately halophilic bacterium: glutamate induces proline biosynthesis in Halobacillus halophilus. J Bacteriol 189(19):6968–6975. https://doi.org/10.1128/jb.00775-07

  • Seebeck FP (2010) In vitro reconstitution of mycobacterial ergothioneine biosynthesis. J Am Chem Soc 132(19):6632–6633. https://doi.org/10.1021/ja101721e

  • Seemann T (2014) Prokka: rapid prokaryotic genome annotation. Bioinformatics 30(14):2068–2069. https://doi.org/10.1093/bioinformatics/btu153

  • Segal G, Ron EZ (1996) Regulation and organization of the groE and dnaK operons in Eubacteria. FEMS Microbiol Lett 138(1):1–10. https://doi.org/10.1111/j.1574-6968.1996.tb08126.x

  • Shirling EB, Gottlieb D (1966) Methods for characterization of Streptomyces species. Int J Syst Evol Microbiol 16(3):313–340. https://doi.org/10.1099/00207713-16-3-313

  • Sun W, Liu C, Zhang F, Zhao M, Li Z (2018) Comparative genomics provides insights into the marine adaptation in sponge-derived Kocuriaflava S43. Front Microbiol 9:1257. https://doi.org/10.3389/fmicb.2018.01257

  • Tan LT-H, Chan K-G, Khan TM, Bukhari SI, Saokaew S, Duangjai A, Pusparajah P, Lee L-H, Goh B-H (2017) Streptomyces sp. MUM212 as a source of antioxidants with radical scavenging and metal chelating properties. Front Pharmacol 8:276. https://doi.org/10.3389/fphar.2017.00276

  • Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, Ostell J (2016) NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res 44(14):6614–6624. https://doi.org/10.1093/nar/gkw569

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Tian X, Zhang Z, Yang T, Chen M, Li J, Chen F, Yang J, Li W, Zhang B, Zhang Z, Wu J, Zhang C, Long L, Xiao J (2016) Comparative genomics analysis of Streptomyces species reveals their adaptation to the marine environment and their diversity at the genomic level. Front Microbiol 7:998. https://doi.org/10.3389/fmicb.2016.00998

  • Tung QN, Linzner N, Loi VV, Antelmann H (2018) Application of genetically encoded redox biosensors to measure dynamic changes in the glutathione, bacillithiol and mycothiol redox potentials in pathogenic bacteria. Free Radic Biol Med 128:84–96. https://doi.org/10.1016/j.freeradbiomed.2018.02.018

    Article  CAS  PubMed  Google Scholar 

  • Tung QN, Loi VV, Busche T, Nerlich A, Mieth M, Milse J, Kalinowski J, Hocke AC, Antelmann H (2019) Stable integration of the Mrx1-roGFP2 biosensor to monitor dynamic changes of the mycothiol redox potential in Corynebacterium glutamicum. Redox Biol 20:514–525. https://doi.org/10.1016/j.redox.2018.11.012

  • Undabarrena A, Ugalde JA, Seeger M, Cámara B (2017) Genomic data mining of the marine actinobacteria Streptomyces sp. H-KF8 unveils insights into multi-stress related genes and metabolic pathways involved in antimicrobial synthesis. PeerJ 5:e2912. https://doi.org/10.7717/peerj.2912

  • Ventura M, Zink R, Fitzgerald GF, van Sinderen D (2005) Gene structure and transcriptional organization of the dnaK operon of Bifidobacterium breve UCC 2003 and application of the operon in bifidobacterial tracing. Appl Environ Microbiol 71(1):487–500. https://doi.org/10.1128/AEM.71.1.487-500.2005

  • Vu THN, Nguyen QH, Dinh TML, Quach NT, Khieu TN, Hoang H, Chu-Ky S, Vu TT, Chu HH, Lee J, Kang H, Li WJ, Phi QT (2020) Endophytic actinomycetes associated with Cinnamomum cassia Presl in Hoa Binh province, Vietnam: distribution, antimicrobial activity and, genetic features. J Gen Appl Microbiol 66(1):24–31. https://doi.org/10.2323/jgam.2019.04.004

    Article  CAS  PubMed  Google Scholar 

  • Vu THN, Quach NT, Nguyen NA, Nguyen HT, Ngo CC, Nguyen TD, Ho P-H, Hoang H, Chu HH, Phi Q-T (2021) Genome mining associated with analysis of structure, antioxidant activity reveals the potential production of levan-rich exopolysaccharides by food-derived Bacillus velezensis VTX20. Appl Sci 11(15):7055

    Article  CAS  Google Scholar 

  • Wibowo JT, Kellermann MY, Petersen L-E, Alfiansah YR, Lattyak C, Schupp PJ (2022) Characterization of an insoluble and soluble form of melanin produced by Streptomyces cavourensis SV 21, a sea cucumber associated bacterium. Mar Drugs 20(1):54

    Article  CAS  Google Scholar 

  • Williams JC, Sheldon JR, Imlay HD, Dutter BF, Draelos MM, Skaar EP, Sulikowski GA (2019) Synthesis of the siderophore coelichelin and its utility as a probe in the study of bacterial metal sensing and response. Org Lett 21(3):679–682. https://doi.org/10.1021/acs.orglett.8b03857

  • Xu L, Dong Z, Fang L, Luo Y, Wei Z, Guo H, Zhang G, Gu YQ, Coleman-Derr D, Xia Q, Wang Y (2019) OrthoVenn2: a web server for whole-genome comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res 47(W1):W52–w58. https://doi.org/10.1093/nar/gkz333

  • Yaakop AS, Chan K-G, Ee R, Lim YL, Lee S-K, Manan FA, Goh KM (2016) Characterization of the mechanism of prolonged adaptation to osmotic stress of Jeotgalibacillus malaysiensis via genome and transcriptome sequencing analyses. Sci Rep 6(1):33660. https://doi.org/10.1038/srep33660

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Zerouki C, Bensalah F, Kuittinen S, Pappinen A, Turunen O (2021) Whole-genome sequencing of two Streptomyces strains isolated from the sand dunes of Sahara. BMC Genomics 22(1):578. https://doi.org/10.1186/s12864-021-07866-x

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Zin NM, Ismail A, Mark DR, Westrop G, Schniete JK, Herron PR (2021) Adaptation to endophytic lifestyle through genome reduction by Kitasatospora sp. SUK42. Front Bioeng Biotechnol 9:916. https://doi.org/10.3389/fbioe.2021.740722

Download references

Acknowledgements

The authors would like to thank the support of VAST–Culture Collection of Microorganisms, Institute of Biotechnology, Vietnam Academy of Science and Technology (www.vccm.vast.vn).

Funding

This study was financially supported by the Graduate University of Science and Technology, Vietnam Academy of Science and Technology under Grant number GUST.STS.ĐT2020-SH04.

Author information

Authors and Affiliations

Authors

Contributions

NTQ and THNV conceived of this study. NTQ, TTAN, and TLB designed and performed the experiments. NTQ, CCN and TTXL supervised and implemented the statistical analysis. NTQ and THNV wrote the manuscript. QTP improved the writing of the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Quyet-Tien Phi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

The participant has consented to the submission of this article to the journal. We confirm that the manuscript, or part of it, has neither been published nor is currently under consideration for publication. This work and the manuscript were approved by all co-authors.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

13213_2022_1684_MOESM1_ESM.docx

Additional file 1: Table S1. Screening of antioxidant activities of methanolic extracts (concentration-1.6 mg/mL) of 15 endophytic actinomycetes isolated from Bruguiera gymnorrhiza. Table S2. Similar characteristics of Streptomyces sp. VCCM 22513 and S. parvulus strain sankarensis-A10. Table S3. Genomic features of S. parvulus VCCM 22513. Table S4. Genes involved in mycothiol and ergothioneine biosynthesis identified in S. parvulus VCCM 22513. Table S5. Genetic determinants contributing to oxidative stress resistance of S. parvulus VCCM 22513. Table S6. Secondary metabolite biosynthetic gene clusters related to antioxidant activities. Table S7. Osmolyte transport and biosynthesis genes prediction in S. parvulus VCCM 22513. Table S8. Genes involved in heat and cold shock response of S. parvulus VCCM 22513. Table S9. Genes involved in degradation of aromatic compounds. Fig. S1. Phylogenetic tree based on full-length 16S rRNA gene sequences derived from strain VCCM 22513 and other closely-related Streptomyces strains. Scale bar, 0.005 substitutions per nucleotide position. Fig. S2. Effect of different H2O2 concentrations (0.6-2M H2O2) to the growth of mangrove endophytic S. parvulus VCCM 22513 and endophytic S. cavourensis strain YBQ59. Fig. S3. The growth of strain VCCM 22513 on CM basal medium un-supplemented and supplemented with gentisate (60 μM), cinnamate (80 μM), styrene (80 μM), or protocatechuate (60 μM).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Quach, N.T., Vu, T.H.N., Bui, T.L. et al. Genomic and physiological traits provide insights into ecological niche adaptations of mangrove endophytic Streptomyces parvulus VCCM 22513. Ann Microbiol 72, 27 (2022). https://doi.org/10.1186/s13213-022-01684-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13213-022-01684-6

Keywords

  • Antioxidant
  • Bruguiera gymnorrhiza
  • Environmental stressors
  • Genome mining
  • Genome-wide comparison
  • Streptomyces parvulus