Diversity of Myxobacteria Isolated from Indonesian Mangroves and Their Potential for New Antimicrobial Sources
Mangroves are unique environments of small forest trees in brackish water with transitional zones in the coastal area. They tolerate a wide range of salinity, oxygen, nutrient and harbor many types of microorganisms . Different microorganisms origin from this habitat have been documented with different their bioactivity and compounds . However, there is a limited information about the antimicrobial potential of myxobacteria isolated especially from Indonesian mangroves, which are known as the world’s largest mangroves . A preliminary study considered that marine environments are promising habitats for the isolation of antibiotic producing bacteria including myxobacteria. Myxobacteria are Gram-negative Deltaproteobacteria and it is one of the most fascinating groups producing natural microbial products . They produced more than a hundred new carbon skeletons and derivatives of compounds with antibacterial, antifungal, antimalarial, antioxidative, or antiviral activity [ – ]. Four genera of rare marine/saline origin are known so far and from two ( Haliangium and Enhygromyxa ) numerous (bioactive) secondary metabolites had been isolated like haliangicin, salimabromide, salimyxins, enhygrolides, and haliamide. Sharma & Thakur  noted that crude extracts with strong antimicrobial activities hint to the presence of promising antimicrobial compounds against infectious disease and treatment. Furthermore, these results were consistent with a previous study of Linares-Otoya et al.  on antimicrobial screening of predatory bacteria including myxobacteria from coastline. Therefore, it can be assumed that also saline environments like mangroves harbor an enormous potential of new antibiotic producing myxobacteria  and an exotic myxobacteria such as species Racemicystis for example from uncommon habitat should be explored. Myxobacteria possess a social lifestyle with moving to get nutrients, preying to another microorganisms, and surviving from lack environment cooperatively in predatory groups . Their communities are commonly found in various environments, even in extreme conditions [ – ]. A preliminary study considered the use of specific primers and probes to analyze myxobacterial communities in soil samples  in addition Li et al.  and Brinkhoff et al.  carried out PCR-DGGE (Denaturing Gradient Gel Electrophoresis) and quantitative real time PCR to define myxobacteria diversity. Clone bank analyses, cultivation approaches were also used to compare their diversity from different habitats . Presently, next generation sequencing (NGS) of hypervariable regions of the 16S rRNA gene, is a tool that gives profound insights into microbial communities. Jiang et al.  stated tremendous success using this method in comparing the bacterial community diversity and composition between mudflat, edge, bulk, and rhizosphere of mangrove’s samples in Hong Kong, China. Furthermore, Linares-Otoya et al.  were able to detect predatory bacteria with NGS from the Peruvian coastline stating that the microbiome present in this region is a promising source for heterotrophic bacterial strains and has potential for bioprospecting of antibiotics. Nevertheless, there are no specific information about myxobacterial communities in mangrove habitat using the NGS method’s. The use of culture-independent methods is of great interest, as it has identified a larger number of novel taxa when compared to culture-dependent methods under standard laboratory condition [ , – ]. However, the culture-dependent methods is important step to get novel antimicrobial compound, because without the successful isolated novel taxa from nature, the production of secondary metabolites never happened . The aim of this study was to evaluate the myxobacterial community composition in Indonesian mangroves by Illumina sequencing of 16S rRNA genes using specific primers, to isolate myxobacteria, and to discover their potential for antimicrobial activity. 

Sample Collection The mangrove samples were collected from Taman Muara Tawar, Bekasi (MT) (6°088772′N 107°002393′E), Taman Muara Angke, Jakarta (MA) (6°105321′N 106°735578′E), and Taman Mangrove Api-Api, Yogyakarta (MK) (7°894662′S 110°02554′E), at Java, Indonesia in November 2018 (Fig. S1). Total of six samples including sediment and leaf flakes were taken from each location. Approximately 50 g of the upper sediment with leaf, which fell on it, were placed in a sterile zip lock plastic bag. All samples were air dried at 30 °C for minimalizing contamination from fungi. Isolation of Myxobacteria The strains were isolated as previously described by Mohr et al.  using water agar with autoclaved Escherichia coli strain K12 and Stan 21 agar with filter paper . Swarming colonies or fruiting bodies were observed under dissecting microscope (Olympus SZX10) every 5–15 days and transferred to new water agar plates with Escherichia coli strain K12 and finally to VY/2 agar plates . The pure cultures from VY/2 agar plates were transferred into 20 ml CY/H liquid medium  and 1.5 ml of well-grown cultures were directly conserved at −80 °C. Identification of Pure Cultures by 16S rRNA Gene Sequences The bacterial DNA was extracted using Invitek Spin Plant Mini Kit (Invisorb) following the manufacturer’s instruction. One microliter of template DNA was directly amplified using F27/R1525 as described . The PCR products were checked by 0.8% agarose gel electrophoresis at 70 V for 40 min and purified by NucleoSpin Gel and PCR Clean up Kit (Macherey–Nagel, Düren, Germany). The forward and reverse sequences of 16S rRNA gene fragments were assembled with the BioEdit program  and closely related type strains were identified using the NCBI 16S rRNA gene database. Sequences were deposited at NCBI database under accession number MW199130, MW182265-MW182288. Microbiome Analysis The total DNA was extracted from 250 mg of sample using the MOBio PowerSoil® Kit following the manufacturer’s instructions. The quality of DNA was measured in Nano-photometer IMPLAN. Amplification was performed using PrimerSTAR HS DNA Polymerase (Takara, Otsu, Shigu, Japan) following the manufacturer’s instructions. Amplification for analyzing the microbial community composition was described by Rath et al.  with modifications. Forward primers W2 and W5  specifically targeting Cystobacterineae and Sorangiineae / Nannocystineae (Table S1) of the Myxococcales were separately used in conjunction with reverse primer R1525  in a first PCR reaction. Samples were pre-denaturated at 96 °C for 3 min following 25 cycles of denaturation at 94 °C for 1 min, annealing at 56.6 °C for 1 min, extension at 72 °C for 2 min. The PCR product was checked on 2.0% agarose gel electrophoresis and purified using NucleoSpin Gel PCR Clean up Kit (Macherey–Nagel, Düren, Germany). One microliter of purified PCR product was used as template in a second PCR with primers 807F and 1050R containing part of the sequencing primer sites as short overhangs (Table S1) for 20 cycles to enrich for target sequences. A third amplification step of 10 cycles added the two indices and Illumina adapters to amplicons . The second and third PCR reaction following a first step PCR reaction. Obtained products were pooled in equimolar ratios and sequenced on an Illumina Miseq (2 × 300 bases, San Diego, USA). The bioinformatic processing was performed as previously described . Raw reads were merged with the Ribosomal Database Project (RDP) assembler . Sequences were aligned within MOTHUR  (gotoh algorithm using the SILVA database  and subjected to preclustering (diffs = 2) yielding so-called phylotypes that were filtered for an average abundance of ≥ 0.001% and a sequence length ≥ 250 bp before analysis. The potential duplications of same Myxococcales sequences between both sets of primers pair were counted and were included to each primer pair analysis. Phylotypes were assigned to a taxonomic affiliation based on the naive Bayesian classification with a pseudo-bootstrap threshold of 80% . All sequences not matching to the order Myxococcales were deleted before further analysis. The relative abundance of genera was plotted using pivot table on Microsoft Excel. Sequences were deposited at NCBI database under accession number SRX9502207-SRX9502210 (Project PRJNA678217). Preparation for Crude Extracts Seed cultures from the 25 isolates were prepared from myxobacterial swarming colonies after seven and before 15 days of incubation (depending on the myxobacteria growth rate) on VY/2 agar medium  by inoculation into 20 ml CY/H liquid medium . The cultures were shaken at 160–180 rpm. After seven days, 10% of the seed cultures were transferred into 100 ml screening VY/2 liquid medium , which contained 2% XAD-absorber resin and incubated at room temperature for 7 up to 15 days. After 2% XAD-absorber was collected, eluate molecules bound on XAD was released using 70 ml using 70 ml acetone (J. T. Baker) and the organic phase was evaporated (Heidolph, Laborota 4003 control) at a temperature of 38–40 °C. The residue was mixed with 1 ml methanol (J. T. baker) and stored in a freezer at −20 °C. Screening for Antimicrobial Activity All 25 crude extracts from 25 strains were evaluated by a test panel using the following pathogenic microorganisms: Escherichia coli WT-BW 25113, Escherichia coli JW0451-02, Acinetobacteria baumanii DSM 30008, Citrobacter freundii DSM 30039, Pseudomonas aeruginosa Pa14, Staphylococcus aureus Newmann, Mycobacterium smegmatis ATCC 700084, Candida albicans DSM 1665, Wickerhamomyces anomalus DSM 6766, and Mucor hiemalis DSM 2656. The antimicrobial activity was assessed by a growth inhibition test using a serial dilution of each crude extract against different pathogenic microorganisms in a 96-well plate . The highest dilution, at which growth inhibition was observed, was noted. The antimicrobial activity was visualized quantitatively by heat map with Heatmapper software . Identification of Active Compound The crude extract of Racemicystis sp. strain 503MSO was fractionated on 96-well plate to analyze the peak-activity correlation. Separation was conducted on an Agilent 1260 high performance liquid chromatography (HPLC) (Agilent Technology USA) using XBridge C 18 column (Waters, Milford, Massachusetts, USA) (100 × 2.1 mm). The solvents consist of 5% acetonitrile (MeCN) in water (H 2 O) added with 5 mmol ammonium acetate (NH4OAc) and 0.04 ml/l acetic acid (CH3COOH) as solvent A and solvent B consist of 95% MeCN, 5 mmol (NH4OAc,) and 0.04 ml/l CH3COOH. Gradient system starting from 10% B to 100% B in 30 min and maintaining at 100% for 10 min, followed by post-run in 10 min for column re-conditioning. Flow rate was 0.3 ml/min and column temperature maintained at 40 °C. The fractions was collected every 30 s . In paralel, the crude extract was analyzed by using diode array detector—high-resolution electrospray ionization mass-spectrometry (DAD-UV-HRESIMS) (MaXis ESI TOF, Bruker Daltonik GmbH, Bremen, Germany) using bridge ethylene hybrid (BEH) C 18 column (Waters ACQUITY, Milford, MA, USA) (1.7 µm 2.1 × 50 mm). A linear gradient from 95% H 2 O and 5% MeCN to 5%H 2 O and 95% MeCN suplemented with 0.1% formic acid was used . The eluate on 96-well plate was dried under streaming nitrogen (N 2 ) gas and the respective bacteria was incubated directly to the plate. The clear zone on the well plate indicating the active compounds as well as the retention time in the respective mass chromatogram. The data were analyzed following Krug and Müller  and Hoffmann et al.  using Data Analysis 4.2 B383 (Bruker Daltonics) and identified with the in-house compounds library (Myxobase). Statistical Analysis The results were analyzed using The IBM SPSS Statistics version 26.0 . It was expressed as mean ± SEM (Standard Error of Mean). The differences between mean were statistically analyzed using one-way analysis of ANOVA with test of homogeneity of variances Bonferroni. P < 0.05 considered to statistical significance. 

The mangrove samples were collected from Taman Muara Tawar, Bekasi (MT) (6°088772′N 107°002393′E), Taman Muara Angke, Jakarta (MA) (6°105321′N 106°735578′E), and Taman Mangrove Api-Api, Yogyakarta (MK) (7°894662′S 110°02554′E), at Java, Indonesia in November 2018 (Fig. S1). Total of six samples including sediment and leaf flakes were taken from each location. Approximately 50 g of the upper sediment with leaf, which fell on it, were placed in a sterile zip lock plastic bag. All samples were air dried at 30 °C for minimalizing contamination from fungi. 

The strains were isolated as previously described by Mohr et al.  using water agar with autoclaved Escherichia coli strain K12 and Stan 21 agar with filter paper . Swarming colonies or fruiting bodies were observed under dissecting microscope (Olympus SZX10) every 5–15 days and transferred to new water agar plates with Escherichia coli strain K12 and finally to VY/2 agar plates . The pure cultures from VY/2 agar plates were transferred into 20 ml CY/H liquid medium  and 1.5 ml of well-grown cultures were directly conserved at −80 °C. 

The bacterial DNA was extracted using Invitek Spin Plant Mini Kit (Invisorb) following the manufacturer’s instruction. One microliter of template DNA was directly amplified using F27/R1525 as described . The PCR products were checked by 0.8% agarose gel electrophoresis at 70 V for 40 min and purified by NucleoSpin Gel and PCR Clean up Kit (Macherey–Nagel, Düren, Germany). The forward and reverse sequences of 16S rRNA gene fragments were assembled with the BioEdit program  and closely related type strains were identified using the NCBI 16S rRNA gene database. Sequences were deposited at NCBI database under accession number MW199130, MW182265-MW182288. 

The total DNA was extracted from 250 mg of sample using the MOBio PowerSoil® Kit following the manufacturer’s instructions. The quality of DNA was measured in Nano-photometer IMPLAN. Amplification was performed using PrimerSTAR HS DNA Polymerase (Takara, Otsu, Shigu, Japan) following the manufacturer’s instructions. Amplification for analyzing the microbial community composition was described by Rath et al.  with modifications. Forward primers W2 and W5  specifically targeting Cystobacterineae and Sorangiineae / Nannocystineae (Table S1) of the Myxococcales were separately used in conjunction with reverse primer R1525  in a first PCR reaction. Samples were pre-denaturated at 96 °C for 3 min following 25 cycles of denaturation at 94 °C for 1 min, annealing at 56.6 °C for 1 min, extension at 72 °C for 2 min. The PCR product was checked on 2.0% agarose gel electrophoresis and purified using NucleoSpin Gel PCR Clean up Kit (Macherey–Nagel, Düren, Germany). One microliter of purified PCR product was used as template in a second PCR with primers 807F and 1050R containing part of the sequencing primer sites as short overhangs (Table S1) for 20 cycles to enrich for target sequences. A third amplification step of 10 cycles added the two indices and Illumina adapters to amplicons . The second and third PCR reaction following a first step PCR reaction. Obtained products were pooled in equimolar ratios and sequenced on an Illumina Miseq (2 × 300 bases, San Diego, USA). The bioinformatic processing was performed as previously described . Raw reads were merged with the Ribosomal Database Project (RDP) assembler . Sequences were aligned within MOTHUR  (gotoh algorithm using the SILVA database  and subjected to preclustering (diffs = 2) yielding so-called phylotypes that were filtered for an average abundance of ≥ 0.001% and a sequence length ≥ 250 bp before analysis. The potential duplications of same Myxococcales sequences between both sets of primers pair were counted and were included to each primer pair analysis. Phylotypes were assigned to a taxonomic affiliation based on the naive Bayesian classification with a pseudo-bootstrap threshold of 80% . All sequences not matching to the order Myxococcales were deleted before further analysis. The relative abundance of genera was plotted using pivot table on Microsoft Excel. Sequences were deposited at NCBI database under accession number SRX9502207-SRX9502210 (Project PRJNA678217). 

Seed cultures from the 25 isolates were prepared from myxobacterial swarming colonies after seven and before 15 days of incubation (depending on the myxobacteria growth rate) on VY/2 agar medium  by inoculation into 20 ml CY/H liquid medium . The cultures were shaken at 160–180 rpm. After seven days, 10% of the seed cultures were transferred into 100 ml screening VY/2 liquid medium , which contained 2% XAD-absorber resin and incubated at room temperature for 7 up to 15 days. After 2% XAD-absorber was collected, eluate molecules bound on XAD was released using 70 ml using 70 ml acetone (J. T. Baker) and the organic phase was evaporated (Heidolph, Laborota 4003 control) at a temperature of 38–40 °C. The residue was mixed with 1 ml methanol (J. T. baker) and stored in a freezer at −20 °C. 

All 25 crude extracts from 25 strains were evaluated by a test panel using the following pathogenic microorganisms: Escherichia coli WT-BW 25113, Escherichia coli JW0451-02, Acinetobacteria baumanii DSM 30008, Citrobacter freundii DSM 30039, Pseudomonas aeruginosa Pa14, Staphylococcus aureus Newmann, Mycobacterium smegmatis ATCC 700084, Candida albicans DSM 1665, Wickerhamomyces anomalus DSM 6766, and Mucor hiemalis DSM 2656. The antimicrobial activity was assessed by a growth inhibition test using a serial dilution of each crude extract against different pathogenic microorganisms in a 96-well plate . The highest dilution, at which growth inhibition was observed, was noted. The antimicrobial activity was visualized quantitatively by heat map with Heatmapper software . 

The crude extract of Racemicystis sp. strain 503MSO was fractionated on 96-well plate to analyze the peak-activity correlation. Separation was conducted on an Agilent 1260 high performance liquid chromatography (HPLC) (Agilent Technology USA) using XBridge C 18 column (Waters, Milford, Massachusetts, USA) (100 × 2.1 mm). The solvents consist of 5% acetonitrile (MeCN) in water (H 2 O) added with 5 mmol ammonium acetate (NH4OAc) and 0.04 ml/l acetic acid (CH3COOH) as solvent A and solvent B consist of 95% MeCN, 5 mmol (NH4OAc,) and 0.04 ml/l CH3COOH. Gradient system starting from 10% B to 100% B in 30 min and maintaining at 100% for 10 min, followed by post-run in 10 min for column re-conditioning. Flow rate was 0.3 ml/min and column temperature maintained at 40 °C. The fractions was collected every 30 s . In paralel, the crude extract was analyzed by using diode array detector—high-resolution electrospray ionization mass-spectrometry (DAD-UV-HRESIMS) (MaXis ESI TOF, Bruker Daltonik GmbH, Bremen, Germany) using bridge ethylene hybrid (BEH) C 18 column (Waters ACQUITY, Milford, MA, USA) (1.7 µm 2.1 × 50 mm). A linear gradient from 95% H 2 O and 5% MeCN to 5%H 2 O and 95% MeCN suplemented with 0.1% formic acid was used . The eluate on 96-well plate was dried under streaming nitrogen (N 2 ) gas and the respective bacteria was incubated directly to the plate. The clear zone on the well plate indicating the active compounds as well as the retention time in the respective mass chromatogram. The data were analyzed following Krug and Müller  and Hoffmann et al.  using Data Analysis 4.2 B383 (Bruker Daltonics) and identified with the in-house compounds library (Myxobase). 

The results were analyzed using The IBM SPSS Statistics version 26.0 . It was expressed as mean ± SEM (Standard Error of Mean). The differences between mean were statistically analyzed using one-way analysis of ANOVA with test of homogeneity of variances Bonferroni. P < 0.05 considered to statistical significance. 

Mangrove Myxobacteria Isolates Seventy strains were successfully isolated from three different sources of Indonesian mangroves. All of the strains showed fruiting bodies formation and swarming on a surface agar medium. Therefore, based on different morphologies and 16S rRNA gene sequences analysis, the number of myxobacteria replicates from same source of organisms and location were reduced and 25 of 70 isolates were obtained. Twenty-five isolates, 16 from MA sampling site, 3 from MT sampling site and 6 from MK sampling site were selected for further analysis. Based on 16S rRNA gene sequence analysis, members of three suborders of Myxococcales were identified (see Table ). Some isolates showed less than 98.60% similarity to close related type strains. Therefore, full genome sequencing is needed for their further characterization and whether they comprise novel species. As main characteristics of myxobacteria as previously described by Reichenbach , some of the swarming and fruiting bodies of myxobacterial isolates from mangroves can be seen in Fig. . Myxococcus sp. strains 431MSO and 451MSO have spherical fruiting bodies with yellow or oranges red colors on VY/2 agar medium. Coralococcus sp. strain 412MSO move with swarm colonies and forms fruiting bodies with coralloid-branched shapes. Chondromyces sp. strain 151MSO builds tree shaped fruiting bodies on Stan21  agar medium and Archangium sp. strain 455MSO makes swarm colonies with branched radial veins on VY/2 agar medium. Racemycistis sp. strain 503MSO has swarming area like the genus Sorangium on VY/2 agar medium. Myxobacterial Diversity Myxobacterial diversity was evaluated after amplification using two different primer pairs W2/R1525 and W5/R1525, targeting the suborders Cystobacterineae and Sorangiineae/Nannocystineae , respectively [ , , ]. A total of 20,057 myxobacterial sequences (1761 ± 1002 Standard Error of Mean/SEM per sample) from 12 samples, where each sample was amplified twice with the two different primer pairs, were obtained. The success of the W2/R1525 primer to amplify members of the order Myxococcales order varied and between 0.56% and 70.42% of the obtained reads belonged to this order (see Table ). With the W5/R1525 primer pair, 7.27–27.03% of amplified sequences were belonging to Myxococcales , indicating a low abundance of those bacteria in that sample. In contrast to the W5/R1525 primer pair, 70.42% of Myxococcales sequences with the W2/R1525 primer pair in one Yogyakarta mangrove sample indicate high abundance of Myxococcales . Moreover, sample MK from Yogyakarta contained more diverse myxobacteria than other and sample MT from Bekasi showed less diversity according to the results. Thirteen major genera were identified from three sampling sites analyzed in this study (Fig. ). The W2/R1525 primer pair could also amplify sequences indicating the presence of Cystobacter , Myxococcus, Stigmatella, Archangium and Anaeromyxobacter of the Cystobacterineae. However, also Haliangium sequences was amplified by this primer. The W5/R1525 primer pair amplified sequence indicating the presence of members of the genera Haliangium , Kloferia and Nannocystis of the Nannocystineae and Chondromyces , Labilithrix , Phaselicystis, Polyangium and Sandaracinus of the Sorangiineae. In addition, some sequences of Myxococcus, Stigmatella and Cystobacter of the Cystobacterineae were observed as being amplified by this primer. Overall, the primers showed good specificity for their targets. Clear differences in myxobacterial diversity could be observed in the sampling sites. Stigmatella spp. exclusively existed in MK samples, Yogyakarta mangroves whereas Nannocystis spp . , and Labilithrix spp . were present only in Jakarta and Yogyakarta mangroves. Only Myxococcus spp . were observed in all six samples, indicating that it is a common myxobacterium in mangroves samples (Table ). Antimicrobial Activity of Myxobacteria Twenty-five crude extracts from the 25 isolates were tested for their antimicrobial activity against ten pathogenic microorganisms, including five isolates Gram-negative bacteria, two-isolates Gram-positive bacteria and three fungi isolates (Table ). Thirteen crude extract showed activity against at least one of the human pathogenic microorganisms tested. Two of them showed activity against both Gram-negative strains Escherichia coli and Citrobacter freundii and nine were active against Gram-positive Staphylococcus aureus. Furthermore, three of the extracts showed activity against the yeasts Candida albicans and Wickerhamomyces anomalus. The crude extract of Racemicystis sp. strain 503MSO was fractionated in order to identify the responsible active compounds. We considered and selected Racemicystis sp. strain 503MSO out of 13 active strains because of ≤ 98.60% 16S rRNA sequence similarity to the type strain and because of missing information about compounds from genus Racemicystis in the Myxobase database. The crude extract of Racemicystis sp. strain 503MSO was fractionated in order to identify the responsible active compounds. The compound identification was done by comparing detected mass of the parent ions of the active fractions with in-house Myxobase database. The Myxobase is a database to support research with myxobacteria, which are increasingly recognized as producers of secondary metabolites. Within Myxobase, the information of bioactivity, retention time, UV spectrum, molecular mass and elemental formula of the molecule responsible for the active peak and HPLC chromatogram were provided. A summary of the active compounds from strain 503MSO was provided in Figure S1. Two active fractions from Racemicystis sp . strain 503MSO were identified comprising compounds with masses of m/z 375.2531 and 604.3857 [M + H] + , respectively after the high-resolution mass-spectrometry (HR-MS) analysis of these active fractions. The peak chromatogram which indicated the presence of compound from the crude extract was mentioned in Fig. S1. 

Seventy strains were successfully isolated from three different sources of Indonesian mangroves. All of the strains showed fruiting bodies formation and swarming on a surface agar medium. Therefore, based on different morphologies and 16S rRNA gene sequences analysis, the number of myxobacteria replicates from same source of organisms and location were reduced and 25 of 70 isolates were obtained. Twenty-five isolates, 16 from MA sampling site, 3 from MT sampling site and 6 from MK sampling site were selected for further analysis. Based on 16S rRNA gene sequence analysis, members of three suborders of Myxococcales were identified (see Table ). Some isolates showed less than 98.60% similarity to close related type strains. Therefore, full genome sequencing is needed for their further characterization and whether they comprise novel species. As main characteristics of myxobacteria as previously described by Reichenbach , some of the swarming and fruiting bodies of myxobacterial isolates from mangroves can be seen in Fig. . Myxococcus sp. strains 431MSO and 451MSO have spherical fruiting bodies with yellow or oranges red colors on VY/2 agar medium. Coralococcus sp. strain 412MSO move with swarm colonies and forms fruiting bodies with coralloid-branched shapes. Chondromyces sp. strain 151MSO builds tree shaped fruiting bodies on Stan21  agar medium and Archangium sp. strain 455MSO makes swarm colonies with branched radial veins on VY/2 agar medium. Racemycistis sp. strain 503MSO has swarming area like the genus Sorangium on VY/2 agar medium. 

Myxobacterial diversity was evaluated after amplification using two different primer pairs W2/R1525 and W5/R1525, targeting the suborders Cystobacterineae and Sorangiineae/Nannocystineae , respectively [ , , ]. A total of 20,057 myxobacterial sequences (1761 ± 1002 Standard Error of Mean/SEM per sample) from 12 samples, where each sample was amplified twice with the two different primer pairs, were obtained. The success of the W2/R1525 primer to amplify members of the order Myxococcales order varied and between 0.56% and 70.42% of the obtained reads belonged to this order (see Table ). With the W5/R1525 primer pair, 7.27–27.03% of amplified sequences were belonging to Myxococcales , indicating a low abundance of those bacteria in that sample. In contrast to the W5/R1525 primer pair, 70.42% of Myxococcales sequences with the W2/R1525 primer pair in one Yogyakarta mangrove sample indicate high abundance of Myxococcales . Moreover, sample MK from Yogyakarta contained more diverse myxobacteria than other and sample MT from Bekasi showed less diversity according to the results. Thirteen major genera were identified from three sampling sites analyzed in this study (Fig. ). The W2/R1525 primer pair could also amplify sequences indicating the presence of Cystobacter , Myxococcus, Stigmatella, Archangium and Anaeromyxobacter of the Cystobacterineae. However, also Haliangium sequences was amplified by this primer. The W5/R1525 primer pair amplified sequence indicating the presence of members of the genera Haliangium , Kloferia and Nannocystis of the Nannocystineae and Chondromyces , Labilithrix , Phaselicystis, Polyangium and Sandaracinus of the Sorangiineae. In addition, some sequences of Myxococcus, Stigmatella and Cystobacter of the Cystobacterineae were observed as being amplified by this primer. Overall, the primers showed good specificity for their targets. Clear differences in myxobacterial diversity could be observed in the sampling sites. Stigmatella spp. exclusively existed in MK samples, Yogyakarta mangroves whereas Nannocystis spp . , and Labilithrix spp . were present only in Jakarta and Yogyakarta mangroves. Only Myxococcus spp . were observed in all six samples, indicating that it is a common myxobacterium in mangroves samples (Table ). 

Twenty-five crude extracts from the 25 isolates were tested for their antimicrobial activity against ten pathogenic microorganisms, including five isolates Gram-negative bacteria, two-isolates Gram-positive bacteria and three fungi isolates (Table ). Thirteen crude extract showed activity against at least one of the human pathogenic microorganisms tested. Two of them showed activity against both Gram-negative strains Escherichia coli and Citrobacter freundii and nine were active against Gram-positive Staphylococcus aureus. Furthermore, three of the extracts showed activity against the yeasts Candida albicans and Wickerhamomyces anomalus. The crude extract of Racemicystis sp. strain 503MSO was fractionated in order to identify the responsible active compounds. We considered and selected Racemicystis sp. strain 503MSO out of 13 active strains because of ≤ 98.60% 16S rRNA sequence similarity to the type strain and because of missing information about compounds from genus Racemicystis in the Myxobase database. The crude extract of Racemicystis sp. strain 503MSO was fractionated in order to identify the responsible active compounds. The compound identification was done by comparing detected mass of the parent ions of the active fractions with in-house Myxobase database. The Myxobase is a database to support research with myxobacteria, which are increasingly recognized as producers of secondary metabolites. Within Myxobase, the information of bioactivity, retention time, UV spectrum, molecular mass and elemental formula of the molecule responsible for the active peak and HPLC chromatogram were provided. A summary of the active compounds from strain 503MSO was provided in Figure S1. Two active fractions from Racemicystis sp . strain 503MSO were identified comprising compounds with masses of m/z 375.2531 and 604.3857 [M + H] + , respectively after the high-resolution mass-spectrometry (HR-MS) analysis of these active fractions. The peak chromatogram which indicated the presence of compound from the crude extract was mentioned in Fig. S1. 

Due to the increasing need for new antibiotics, the present study focused on the evaluation of the myxobacterial community as well as the isolation of novel members of the order Myxococcales and their potential for antimicrobial activity. Recent isolations of members of the myxobacteria have frequently resulted in the identification of new secondary metabolites with antibiotic activities. Untapped habitats such as various marine ecosystem have been proven to be a valuable resources for new microorganisms [ , , ]. Indonesian mangroves explored in the present study were unexploited marine habitats and have a huge potential for the isolation of new uncommon species of myxobacteria. Twenty-five strains were isolated from three different sites of Indonesian mangroves by conventional methods for isolation. Myxococcus strains were presented in all six samples and it reflected both myxobacteria isolates versus myxobacterial communities that were amplified based on 16S rRNA gene sequences. According to the fruiting body formation of Myxococcus , they were easy to be purified and cultivated under laboratory conditions . Previously, this genus was used as model for wide-ranging studies of biotechnology [ , – ]. The analysis of myxobacterial communities by culture-independent methods revealed more abundances than myxobacterial isolates by culture-dependent method. For instance, the culture-independent community analyses revealed the genus Haliangium to be one of the most dominant and common myxobacteria in MA and MT samples, whereas it was not found by the culture-dependent method in our study. The faster growing fungal or bacterial contaminants exhaust nutrients before slow growing myxobacteria could be observed and further, the insufficient growth conditions in laboratory missing requirements such as changing tidal-caused salinity might be reasons that rare myxobacterial species were not isolated . Thirteen out of 25 crude extracts from myxobacterial isolates yielded different activities against pathogenic microorganisms. These predatory bacteria including myxobacteria from coastline may be promising sources for novel antibiotics . In this study, we only focused on the active crude extract of Racemicystis sp . strain 503MSO, an isolate from Yogyakarta mangroves. Racemicystis sp . strain 503MSO was selected for further analysis because there was no information about potential bioactive compounds from the genus Racemicystis in the Myxobase database yet. Therefore, it has a higher possibility to get new compounds from this genus. The active fraction of Racemicystis sp. strain 503MSO with m/z 375.2531 [M + H] + and 604.3847 [M + H] + were detected in this study and may represent so far, an undescribed compound. Therefore, Racemicystis sp. strain 503MSO should be further explored including substantial work for optimizing isolation, structure elucidation and biological activity of the pure substance. To the best of our knowledge no bioactive compounds from the same genus and other sources were reported in the scientific publication to date. Indonesian mangroves were inhabited by complex myxobacterial communities as shown by amplification using primer targeting Cystobacterinae (W2/R1525) and Sorangiineae/Nannocystineae (W5/R1525) [ , , ]. Previously, Linares-Otoya  observed that the use of universal primers (F515/R806) yielded only a neglectable amount of myxobacterial sequences (equivalent to 0.0002%) as also indicated by our work. The rapid sequencing-based analysis of myxobacterial communities has advantages compared to previous studies comprising DGGE (Denaturing Gel Gradient Electrophoresis) , hybridization analysis of a 16S rRNA gene library  or library analysis  and by use of two PCR reactions an overview of all three suborders, i.e., Cystobacterineae , Sorangiineae and Nannocystineae were obtained . However, with the current method, a high amount of non- Myxococcales genera were amplified and target sequences often comprised < 10% of sequences. Evidently, optimization of primers is still necessary. Moreover, a species level information is only possible in exceptional cases using the V5V6 hypervariable region as performed here. Overall, mangrove habitats are a rich source for antimicrobial compounds produced by myxobacteria, highlighting the necessity to explore their diversity. 

This study confirmed that neglected areas such as mangroves are promising habitats for the isolation of novel myxobacteria strains and potential sources for unknown secondary metabolites producers. Based on comparison of its monoisotopic mass and retention
time of compound described in Myxobase database, we suggest further studies for isolating and upscaling active compound of Racemicystis sp. strain 503MSO and applying novel models of biological test as well. 

Below is the link to the electronic supplementary material. Supplementary file1 (TIFF 111 KB) Supplementary file2 (TIFF 1804 KB) Supplementary file3 (DOCX 14 KB) 

