# Multiple Streptomyces species with distinct secondary metabolomes have identical 16S rRNA gene sequences

Sanjay Antony-Babu, D. Stien, Véronique Eparvier, Delphine Parrot, Sophie Tomasi, Marcelino Suzuki

► To cite this version:

Sanjay Antony-Babu, D. Stien, Véronique Eparvier, Delphine Parrot, Sophie Tomasi, et al. Multiple Streptomyces species with distinct secondary metabolomes have identical 16S rRNA gene sequences. Scientific Reports, 2017, 7 (1), pp.11089. 10.1038/s41598-017-11363-1. hal-01590493

HAL Id: hal-01590493

https://hal.sorbonne-universite.fr/hal-01590493

Submitted on 19 Sep 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Multiple *Streptomyces* species with distinct secondary metabolomes have identical 16S rRNA gene sequences

Received: 14 September 2016

Accepted: 23 August 2017

Published online: 11 September 2017

Sanjay Antony-Babu<fnref n="1" />,⁴, Didier Stien<fnref n="1" />, Véronique Eparvier², Delphine Parrot³, Sophie Tomasi³ & Marcelino T. Suzuki<fnref n="1" />

Microbial diversity studies using small subunit (SSU) rRNA gene sequences continue to advance our understanding of biological and ecological systems. Although a good predictor of overall diversity, using this gene to infer the presence of a species in a sample is more controversial. Here, we present a detailed polyphasic analysis of 10 bacterial strains isolated from three coastal lichens *Lichina confinis*, *Lichina pygmaea* and *Roccella fuciformis* with SSU rRNA gene sequences identical to the type strain of *Streptomyces cyaneofuscatus*. This analysis included phenotypic, microscopic, genetic and genomic comparisons and showed that despite their identical SSU rRNA sequences the strains had markedly different properties, and could be distinguished as 5 different species. Significantly, secondary metabolites profiles from these strains were also found to be different. It is thus clear that SSU rRNA based operational taxonomy units, even at the most stringent cut-off can represent multiple bacterial species, and that at least for the case of *Streptomyces*, strain de-replication based on SSU gene sequences prior to screening for bioactive molecules can miss potentially interesting novel molecules produced by this group that is notorious for the production of drug-leads.

It is now clear that almost all eukaryotic organisms have an associated prokaryotic microflora. Thanks to advances in massively parallel sequencing technology, accessibility to high performance computing (HPC), and the development of comparative analytical software, the number of prokaryotic microflora surveys have increased in recent years<fnref n="1" />. These studies are mostly based on the analysis of small subunit (SSU) rRNA diversity² and take advantage of the correlation between the gene sequence and prokaryotic taxonomical hierarchy³. In this type of studies sequences are clustered together at a certain level of identity (97% in the vast majority of the cases) in what are termed Operational Taxonomic Units (OTUs)⁴. While such OTU-based analysis is very useful for describing major patterns of microbial diversity, they are increasingly being used as proxies to bacterial species⁵. It is also common to find ecological or physiological traits attributed to OTUs of uncultured microorganisms based on similarities of their SSU rRNA gene sequences to those of bacterial strains with known ecological or physiological traits⁶. However, it is not clear whether these OTUs are functionally coherent, that is, whether strains assigned to a single OTU represent a single species and whether strains with identical SSU rRNA gene sequences exhibit similar ecological or physiological functions.

To answer the above-mentioned questions, one approach would be to study a bacterial taxon that is widespread across different ecosystems and niches, potentially represents keystone organisms and that are metabolically active. Members of the genus *Streptomyces* represent good model organisms with such properties⁷⁻¹¹ and well-known for their ability to produce vast numbers of biologically-active secondary metabolites¹². Since

<footnote n="1">¹Sorbonne Universités, UPMC Univ Paris 06, CNRS, Laboratoire de Biodiversité et Biotechnologies Microbiennes (LBBM), Observatoire Océanologique, F-66650 Banyuls/Mer, France CNRS, USR 3579, LBBM, Observatoire Océanologique, 66650, Banyuls-sur-Mer, France. ²CNRS, Institut de Chimie des Substances Naturelles, 91198, Gif-sur-Yvette cedex, France. ³UMR CNRS 6226, Institut des Sciences Chimiques de Rennes, Equipe CORINT "Chimie Organique et Interface", UFR Sciences Pharmaceutiques et Biologiques, Univ. Rennes 1, Université Bretagne Loire, 2 Avenue du Pr. Léon Bernard, F-35043, Rennes, France. ⁴Present address: School of Biological Sciences, University of Nebraska-Lincoln, Lincoln, Nebraska, 68588, United States of America. Correspondence and requests for materials should be addressed to M.T.S. (email: suzuki@obs-banyuls.fr)</footnote>
<table><caption>**Table 1.** Source of the 10 selected strains from Parrot *et al.* (2015) used in this study, and their phenotypic and genomic characteristics. <sup>a</sup>MA- marine agar, AIA – actinomycetes isolation agar. <sup>b</sup>H- homogenate, W- Wash.</caption><thead><tr><th rowspan="2">Isolate code</th><th rowspan="2">Lichen source</th><th rowspan="2">Medium of isolation<sup>a</sup></th><th rowspan="2">Inoculum type<sup>b</sup></th><th colspan="3">Halotolerance - NaCl %</th><th colspan="3">Growth at different pH</th><th rowspan="2">Color group number</th><th rowspan="2">GC %</th></tr><tr><th>5%</th><th>10%</th><th>20%</th><th>5</th><th>7</th><th>9</th></tr></thead><tbody><tr><td>MOLA1420</td><td rowspan="3">Lichina pygmaea</td><td>MA</td><td>H</td><td></td><td>++</td><td></td><td></td><td>++</td><td>+</td><td>1</td><td>64%</td></tr><tr><td>MOLA1421</td><td>MA</td><td>H</td><td></td><td>++</td><td></td><td></td><td>++</td><td>+</td><td>1</td><td>64%</td></tr><tr><td>MOLA1427</td><td>AIA</td><td>H</td><td></td><td>++</td><td></td><td></td><td>++</td><td></td><td>2</td><td>64%</td></tr><tr><td>MOLA1488</td><td rowspan="4">Lichina confinis</td><td>AIA</td><td>W</td><td>+</td><td>++</td><td></td><td></td><td>++</td><td>+</td><td>3</td><td>64%</td></tr><tr><td>MOLA1493</td><td>AIA</td><td>W</td><td>+</td><td>++</td><td></td><td></td><td>++</td><td>++</td><td>4</td><td>66%</td></tr><tr><td>MOLA1578</td><td>AIA</td><td>H</td><td>++</td><td>++</td><td></td><td></td><td>++</td><td>++</td><td>5</td><td>66%</td></tr><tr><td>MOLA1596</td><td>AIA</td><td>W</td><td>++</td><td>++</td><td></td><td></td><td>++</td><td>++</td><td>5</td><td>66%</td></tr><tr><td>MOLA1615</td><td rowspan="3">Roccella fuciformis</td><td>AIA</td><td>W</td><td>++</td><td>++</td><td></td><td></td><td>++</td><td>++</td><td>5</td><td>66%</td></tr><tr><td>MOLA1611</td><td>AIA</td><td>W</td><td>++</td><td>++</td><td>+</td><td>+</td><td>++</td><td>++</td><td>6</td><td>64%</td></tr><tr><td>MOLA1612</td><td>AIA</td><td>W</td><td>++</td><td>++</td><td>+</td><td>+</td><td>++</td><td>++</td><td>6</td><td>64%</td></tr></tbody></table>

these molecules play central roles in microbial competition and inter- and intra- species signaling<sup>13</sup> they represent a phenotypic characteristic with niche determining potential<sup>14</sup>, and due to their bioactivity have extensively been developed as drugs<sup>12</sup>. To add further significance to *Streptomyces* as model for the current study, systematics of these bacteria have had a “turbulent” past<sup>15</sup> where close to 3000 named-species existed the 1970s albeit, almost solely on patent literature<sup>16</sup>. After several efforts to amend their taxonomy<sup>17</sup>, around 600 species of *Streptomyces* are currently officially recognized<sup>18</sup> and those are classified in 130 supra-specific clades<sup>17, 19, 20</sup>, regardless of the diversity observed among strains with similar SSU rRNA gene sequences<sup>20-24</sup>.

In the study presented here, the phenotypes of 10 *Streptomyces* strains, isolated as part of an earlier microbial diversity study<sup>25</sup> from three different lichens (*Lichina confinis*, *Lichina pygmeae* and *Roccella fuciformis*), using different treatments, and that had 100% similar SSU rRNA gene sequences were compared. The strains were subject of morphological, phenotypical, molecular taxonomic and metabolic analyses to test whether stains with identical SSU rRNA gene sequences represented the same species, and presented similar functional traits.

## Results

**Stress resistance phenotypic assays.** Stress resistance was evaluated from as increase in dry weight after growth in International Streptomyces Project media 2 (ISP2) that had been either amended with NaCl (5, 10 and 20%) or adjusted to the different pH levels (5.0, 7.0 & 9.0; Table 1). The individual growth curves are presented in supplementary Figure S1.

**Halotolerance.** All 10 tested strains grew well in 10% NaCl and the strains could be clustered into 4 groups based on the growth kinetics at different concentrations of NaCl (Table 1). Three strains, MOLA1420, MOLA1421 and MOLA1427, all isolated from *Lichina pygmaea*, only grew with 10% NaCl. Interestingly, although MOLA1420 and MOLA1421 were isolated on marine agar they were not particularly more halotolerant when compared to some strains isolated using actinomycetes isolation agar. MOLA1488 and MOLA1493, both strains that originated from *Lichina confinis*, showed moderate growth at 5% salt (in addition to optimal growth at 10%). Five strains isolated from *Roccella fuciformis* could be separated into two groups based on their growth at different salt concentrations and all five strains grew well in 5% and 10% salts. MOLA1578, MOLA1596 and MOLA1615 could not grow at 20% that differentiated them from the MOLA1611 and MOLA1612, which could.

**Optimal growth pH.** All test strains grew well in pH 7.0. Based on the growth kinetics at different pHs, the strains could be clustered into 4 groups. Unlike with their salt tolerance, pH growth kinetics had little correlation with the lichens of origin. Particularly two groups could be defined in isolates from different lichens: MOLA1420 and MOLA1421 that were isolated from *Lichina pygmaea* and MOLA1488 from *L. confinis* showed trace growth at pH 9.0 along with optimal growth at pH 7.0; MOLA1493 from *L. confinis* and MOLA1578, MOLA1596 and MOLA1615 from *R. fuciformis* showed good growth in both pH 7.0 and 9.0, but no growth at 5.0. MOLA1611 and MOLA1612 from *R. fuciformis* had the widest growth range with trace growth at pH 5.0 and good growth at pH 7.0 and 9.0. Such wide pH range is typical of alkaliphilic members of the genus *Streptomyces*<sup>26</sup>. MOLA1427 from *L. pygmaea* was the only strain that grew solely at pH 7.0 and was the most stringent in terms of pH requirement.

**Color groups and morphology.** Color group properties are listed in supplementary Table S1 and the numbers are correlated with the isolates in Table 1. A total of 6 color groups were found, with 3 color groups represented by multiple strains and 3 represented by a single strain. Figure 1 shows the colony characteristics and Figure S2 spore chain morphology respectively, that were associated to the 6 different color groups. All strains isolated from *Lichina* spp. could produce melanin on ISP 6 while the 5 strains from *Roccella fuciformis* could not. Remarkably, and in partial agreement with the halotolerance experiments, each color group included strains isolated from a single lichen species and isolates from each of the three lichens belonged to two different color groups. Except for MOLA1488, all strains produced spores on ISP 3 media. All strains from *Lichina pygmaea* in color groups 1 and 2 showed the same spore chain arrangement while the two isolates from *L. confinis* showed
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Colony morphology on oatmeal agar (ISP 3) plates of each of the color group recorded from the 10 test strains. (A) Edge of colony showing aerial spore mass & diffusible pigments, and (B) Colony reverse.</figcaption></figure>

very different colony characteristics on ISP 3 with one of them, MOLA1488 not producing any spores whereas MOLA1496 showed spores held in clusters. The five tested isolates from *Roccella fuciformis* belonged to color groups 5 and 6, and had spores arranged in chains. Members of color group 5 have spirally arranged spores in contrast to the straight chains observed in group 6.

**Protein coding genes.** The tree obtained from the sequences of five concatenated housekeeping gene sequences is presented in Fig. 2. Using a 100% nucleotide identity as the clustering cut-off, the ten test strains were grouped into 7 types that clustered similarly when the different genes were analyzed individually. Dissimilarities between the sequences are shown in supplementary Figure S3 along with phylogenetic trees constructed using individual genes. Three pairs of strains, MOLA1420 - MOLA1421, MOLA1578 - MOLA1596, and MOLA1611 - MOLA1612 had identical sequences for all genes sequenced in this study. Isolate MOLA1427 from *Lichina pyg-maea* branched with two other isolates from the same lichen (MOLA1420 and MOLA1421) in all phylogenetic reconstructions. The five isolates from *Roccella fuciformis* always clustered together based on all 5 genes, with MOLA1615 mostly grouping with the pair MOLA1578-MOLA1596. As a general feature MOLA1488 and MOLA1493 tended to be more distant to the remainder of the strains. Finally, no differences in topology were observed when using DNA or their translated sequences (data now shown).
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Tree showing diversity between the strains with identical SSU rRNA gene sequences. Branch colors represent the lichen sources as follows: blue *R. fuciformis*, green *L. confinis*, red *L. pygmaea*. Phylogenetic tree inferred on concatenated sequences of the 5 protein coding sequences (atpD, gyrB, recA, rpoB and trpD) using neighbor-joining analyses. Confidence in tree topology was calculated by bootstrap test amongst 1000 replicates. There were a total of 3844 positions in the final dataset.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Upper triangular matrix of mean Δ*T*m values. Δ*T*m is the difference in melting temperature of the homologous genomic DNA and hybrid DNA preparation corresponding to their 50% decrease in fluorescence in melting curves. The values in the brackets are standard error of the mean between 3 replicates. Hybrid pairs showing less than 5 degree difference are considered as same genomic species.</figcaption></figure>

**Genomic similarities.** GC content in the genomes was estimated using fluorometrically measured DNA melting curves. The GC% values of each strain are listed in Table 1. Six strains had 64% GC while the other 4 had 66%, and no systematic correlation with species, color groups or other phenotypic characteristics could be established. Based on a criteria recommended by Stackebrandt et al., (2002)²⁷ of 5 °C difference cut-off in melting temperature between homologous and heterologous DNA-DNA hybrids the strains could be classified into 5 different genomic species (Fig. 3). The three isolates from *Lichina pygmaea* belongs to same species, whereas the two strains of *Lichina confinis* represent two additional species and finally the 5 strains from *Roccella fuciformis* represent two additional separate species.

**Metabolic profiling.** Metabolic profiles were distinct among the different strains as indicated in Table S2 and Figures S4–S5. Differences between the metabolic profiles of the 10 strains are represented in Figure S4 and S5. As expected from the remainder of the analyses, MOLA1596 and MOLA 1611 were most closely related with 6 common metabolites. No obvious relationship was found between the other strains. It should be noted that MOLA1488 was again the most distinct of all, with a very high relative proportion of specific secondary metabolites not present in any other strain. High-resolution mass spectrometry allowed the putative identification of 38 natural products, while 14 metabolites were possibly novel. Among the latter, remarkably, 11 metabolites were unique to a single strain.

**Discussion**

The study presented here examined the diversity between 10 isolates that were affiliated to the genus *Streptomyces* and exhibited identical SSU rRNA gene sequences. Phenotypic characteristics such as color grouping, halotolerance, pH growth range and protein coding gene sequences showed that the strains could be consistently segregated into 6 groups. Two of these groups, however, represented the same microscopic morphology and genomic species (by percentage genome similarities). This is, to best of our knowledge, the first report of bacterial strains that have 100% identical SSU rRNA sequences and that are yet diverse enough to be classified as 5 different species. Phylogenetic and metabolic diversity between *Streptomyces* strains with highly similar SSU rRNA gene sequences has been studied previously²⁴, ²⁸, ²⁹ and such metabolic diversity within phylogenetically close strains have been also demonstrated in other actinomycetes such as *Salinispora*³⁰. The metabolic and genomic variations between strains with highly similar SSU rRNA gene sequences appear to be particularly high in the bacterial family Actinomycetaceae as noted in this study, and also in strains of family Vibrionaceae³¹. Additionally, globally widespread planktonic cyanobacteria *Synechococcus* and *Prochlorococcus* that belong to different genera, bear
completely different photosynthetic apparatus despite the high SSU rRNA similarity across the two (>98%)³². Nevertheless, none of these reports examined strains with identical SSU rRNA sequences, nor did they find that these bacterial strains belonged to different species.

Consequences of under-appreciating diversity between strains with identical or nearly identical SSU rRNAs genes are high. This is further exacerbated by the fact that in most studies using new generation sequencing, only a partial sequence of the gene is obtained, and sequences are clustered at operational taxonomy units at 97% sequence similarity. Here we would highlight two major repercussions that directly affect linked to human welfare. Firstly, the test organisms in this study belong to the family Actinomycetaceae, that represents the main natural source of bioactive natural products³³. These molecules make up the largest source of pharmacological compounds (such as antibiotics) that are currently used³⁴ and thus have been the main subjects of natural products drug discovery screenings. In these discovery campaigns, it is common to de-replicate the strains based on SSU rRNA gene sequences. In the likelihood that our results with S. cyanofuscatus is also true to other actinomycetes, de-replication based on SSU rRNA sequences would result in overlooking possible novel chemical structures with drug potential, as it was the case for 14 metabolites from 6 analyzed metabolomes. In the case of Vibrionaceae, the SSU sequences of several species that are pathogenic to humans and to marine invertebrates of commercial importance, are in some case nearly identical to those of non-pathogenic species (i.e. V. parahaemolyticus and V. natriegens³⁵), which severely hampers the use of SSU rRNA and its coding gene for environmental monitoring.

Finally, it was very interesting to note that many of the phenotypical characteristics of the strains, protein coding sequences and the secondary metabolic patterns were most similar between strains isolated from the same lichen species. This is in line with earlier studies that have shown that phenotypical (and by extension and most likely genomic) similarity might be higher in sympatry²⁴, ³⁶, and they could indicate niche-specific adaptation²⁴, ³⁷, ³⁸. Alternatively (or in addition) this phenotypical and genomic similarity between sympatric strains might be due to an increased opportunity for horizontal gene exchange, a phenomenon that is well known to occur in *Streptomyces*³⁹.

## Conclusion

A polyphasic approach was used to assess phenotypic, morphology and genomic relationship between 10 strains with identical SSU rRNA gene sequences. We demonstrate that despite this identity, the strains group as 5 different species. The results hold significance to the increased use of SSU rRNA gene sequences as proxies of prokaryotic diversity. We also show that the biosynthetic properties varied between the 5 species.

## Methods

**Source of isolates.** A total of ten bacterial strains with identical SSU rRNA gene sequence were used in this study. The strains were isolated from 3 different coastal lichens; *Lichina pygmaea*, *Lichina confinis* and *Roccella fuciformis* (Table 1). Details of isolation protocols are described by Parrot and collaborators²⁵. Briefly, the lichen samples were aseptically crushed, washed with sterile seawater and ground to form a homogenate. The wash water and the homogenate were used as different inoculant on actinomycetes isolation agar (AIA, Difco™, France) and marine agar (MA, Difco™, France), both media supplemented with nalidixic acid (Sigma-Aldrich, France) and cycloheximide (Sigma-Aldrich, France). The isolates were deposited in the MOLA culture collection (WDCM911, Collection of Microbial Observatory Laboratoire Arago, France). All 10 isolates of interest have identical SSU rRNA gene sequences to that of *Streptomyces cyaneofuscatus* type strain JCM 4364¹.

**Preparation of spore suspension.** The strains were cultured on International Streptomyces Project medium 3 (ISP 3; oat meal agar⁴⁰, ⁴¹) for 21 to 30 days. The spores were collected by scrapping over top of the colonies and suspended in sterile phosphate buffered saline solution. The preparation was passed through sterile cotton wool in a syringe to remove the hyphal matrix. Numbers of viable spores were estimated number of colonies formed when plated on nutrient agar media. The spore suspensions were maintained at -20°C for not more than 2 weeks. Before use the suspensions were thawed at 4°C.

**Halotolerance and pH growth range.** Halotolerance of the isolates were estimated by plotting growth curves at differing pH and salt concentrations. International Streptomyces Project media 2 (ISP 2) broths were prepared with three NaCl concentrations (5%, 10% & 20%) and three different pH values (5.0, 7.0 & 9.0). The pH tests were carried out in broths that had 10% NaCl and the halotolerance tests were done at in media at pH 7.0.

Spore suspensions (~10² viable spores) of the strain were inoculated into 9 tubes of each media and incubated in an orbital shaker at 28 °C at 120 rpm. For each strain, 3 tubes of every media were removed for analyses after 7, 14 and 21 days. Each tube represented a replicate (among triplicates) for each strain per condition. The broth was individually filtered through pre-weighed Whatman® cellulose filter paper (Grade 1). The harvested biomass oven dried at 55 °C for 24 h. Total dry weights of each of the triplicates per condition for each strain were recorded and plotted against time to evaluate tolerance to pH and salt.

**Color grouping.** The test strains were grouped according to their pigment production characteristics as recommended by Williams *et al.*, 1969. A computer-assisted numerical method recommended by Antony-Babu and colleagues⁴² was used since this method could be more objective than the protocol first proposed in 1969²³. The 10 test isolates were inoculated as spore suspension (~50 spores) onto ISP 3 and International *Streptomyces* Project media 6 (ISP 6, Difco, France⁴³) as triplicate streak cultures. ISP 3 plates were incubated 21 days and ISP 6 for 4 days, both at 28 °C. At the end of 4 days ISP 6 media cultures were checked for melanin pigment production. The ISP 3 plates were examined by eye for aerial spore mass color, substrate mycelium pigmentation and the color of any soluble pigments after 7, 14 and 21 days. The isolates were grouped based on similarities in the colony color properties.
**Morphological examination.** Spore chain arrangements of the strains were examined microscopically. To encourage spore formation, the strains were inoculated (~50 spores) onto ISP 3 media. Sterile coverslips were inserted into the agar surface in a 45° angle, with the bottom face leaning over the inoculum. The set-up was incubated at 28 °C for 21 days. After the incubation period, the coverslips were carefully removed and stained with methyl violet stain. The spore chain morphology was examined under 100 X magnification using a light microscope.

**Genomic DNA isolation.** In order to extract genomic DNA, the test strains were cultured on media described by Antony-Babu and Goodfellow²⁶. This media allows growth of streptomycetes on agar surfaces that are easily detached from the substrate mycelium. The aerial biomass was harvested after 5 days using sterile tooth-picks without touching the agar. The cells were first briefly rinsed with 5 mM EDTA (pH 8.0) and then with 99% ethanol. The pelleted biomass was resuspended in 10% SDS and mixed thoroughly. The suspension was heat-lysed at 80 °C for 20 minutes. The lysate was extracted using a conventional phenol: chloroform: isoamylalcohol protocoll, followed by chloroform: isoamyl alcohol extraction and ethanol precipitation. The dried pellets were suspended in molecular grade DNase and RNase free water. The genomic DNA extracts were checked for quality & quantity on agarose gel and NanoVue spectrophotometer (GE Healthcare) and stored at -20 °C until further use.

**Protein coding gene sequencing.** Genomic DNA samples were used as templates for PCR amplification of the protein-coding housekeeping genes: *atpD*, *gyrB*, *recA*, *rpoB* and *trpB* using PCR primers previously described by Guo and colleagues²¹. The original conditions described by the authors were modified to accommodate the use of a fast PCR reaction mix KAPA2G (Clinisciences, Nanterre France), where the denaturation, annealing and amplification steps were carried out for 15, 15 and 30 seconds respectively. Amplicons were purified using Agencourt® AMPure® XP Kit (Beckman Coulter, Villepinte, France) and used in sequencing reactions were carried out with Big Dye Terminator (V3.1; Life Technologies) using internal sequencing primers²¹ as well as the PCR primers used to synthetize the amplicons. Cycle sequencing products were cleaned using Agencourt® CleanSeq® Kit (Beckman Coulter) and sequenced using ABI 3130xl genetic Analyser Sequencer (Life Technologies). Partial sequences were assembled into contigs using MEGA 5.1 and for each of the housekeeping gene, the sequences available on the Genbank database that were less than 2000 bp and attributed to genus *Streptomyces* were retrieved. Contigs were aligned along with the database sequences using CLUSTAL X⁴⁴. Neighbor-joining phylogenetic trees were built using distance matrices obtained using the Kimura-2 parameter model and their robustness evaluated from 1000 trees built using bootstrap resampling.

**Genomic GC%.** G + C content of the strains were obtained fluorimetrically based on the procedure described by De Ley and Tijtgat (1970)⁴⁵. The original method was downscaled to 20 µl with 0.5 µg of each genomic DNA suspended in a 0.1 ml optical clear tubes along with 0.1x SSC buffer (pH 8) and 1x SYBR Gold. A thermal ramp of 25°C to 100°C with a 1 °C rise per minute was performed using a StepOnePlus™ Real-Time PCR system (Applied Biosystems). The fluorescence at 300 nm was measured at each step of the thermal ramp and the T~m~ of each strain calculated from the minimum value of the slope of the tangent to the melting curve of fluorescence versus temperature. The %G + C of each strain was calculated using the equation: %G + C = (1.99T~m~) - 71.08.

**Genomic similarities.** The optimum temperature of renaturation (T~or~) was calculated using the GC% obtained earlier, by applying the formula T~or~ = 0.51(%G + C) + 47. As with the GC content protocol, the reaction volumes were downscaled to 20 µl. Genomic DNA hybrids were generated by mixing 0.5 µg each of genomic DNA from the reference and the test strain and suspended with 0.1x SSC buffer (pH 8) in a 0.1 ml optical clear tubes (for a total volume of 18 µl). The homologous DNA control and hybrid DNA samples were denatured and re-annealed using in a StepOnePlus™ Real-Time PCR system (Applied Biosystems) using the following conditions: 99 °C for 10 minutes, T~or~ for 8 hours followed by progressive steps in which the temperature was dropped by 10 °C then held for 1 hour until 25 °C was achieved. 2 µl SYBR Gold was added to this reaction mix to achieve a final 0.1 X concentration of the dye and 1 X concentration of SSC in the final mix. The thermal denaturation was achieved as follows: 25 °C for 15 minutes, followed by a thermal ramp from 25–100 °C with a 1 °C rise per minute, fluorescence measurements were taken at each step of the ramp. The T~m~ of the homologous and hybrid DNA preparations were calculated by measuring the temperature corresponding to a 50% decrease in fluorescence in the melting curve of fluorescence versus temperature. ΔT~m~ is the difference between the temperatures of homologous and heterologous genomic DNA mixes; differences of 5 °C or more are considered to show that the tested strains belong to different genomic species⁴⁶, ⁴⁷.

**Metabolic profiling.** *Extraction of metabolites.* Representative strains from the 6 color groups (MOLA1420, MOLA1427, MOLA1488, MOLA1493, MOLA1596 & MOLA1611) were used in this analysis. Inocula (50 µl) from PBS spore suspensions were spread over trypic soy agar plates and incubated at 28 °C for 15 days. Entire plate contents were shredded into small pieces (~0.5 cm³) using a sterile scalpel. The agar blocks were suspended in 15 ml of HPLC grade ethyl acetate (Carlo Erba, France). The preparation was shaken overnight in the dark at 5 °C on an orbital shaker. The liquid phase was removed, dehydrated with addition of sodium sulfate (~1.0 g) and filtered through cotton wool. The extracts were evaporated in a HT4X rotatory vacuum concentrator (Genevac, Ipswich, UK) system over night at 5 °C. Ethyl acetate extracts were treated prior to analysis using solid phase extraction (Waters Sep-Pak C18, Vac 1 cm³, 100 mg). The stationary SPE phase was activated with 2 mL of acetonitrile. The cartridge was loaded with <5.0 mg of each extract, which was solubilized in 500 µL of acetonitrile. The elution was performed with 1 mL acetonitrile to eliminate the most lipophilic metabolites. The eluates were dried in a Genevac system and were dissolved in acetonitrile to afford yield 1.5 mg/mL solutions for profiling analyses.
*Profiling.* The composition of the extracts was analyzed with a Waters Acquity Ultra High Pressure Liquid Chromatography system coupled to a Waters Micromass LCT Premier Time-of-Flight mass spectrometer (Milford, MA, USA) equipped with an electrospray interface. The UHPLC system was equipped with a Waters Acquity UPLC BEH C18 column (150 mm × 2.1 mm, 1.7 µm). The oven temperature was maintained at 40 °C and the solvent flow was 0.6 mL/min. The mobile phases were H₂O (A) and acetonitrile (B), each containing 0.1% formic acid (v/v), with the following gradient: 90:10 (0–1 min), 10 to 100% B (1–20 min), and 100% B for 10 min. Analysis of each extract (1.0 µL injected) was performed in both positive (ESI⁺) and negative ionization (ESI⁻) modes in the 80–1500 Da range. The other ESI conditions were set as follows: capillary voltage 2800 V, cone voltage 40 V, source temperature 120 °C, desolvation temperature 330 °C, cone gas flow 20 L/h, desolvation gas flow 600 L/h, and MCP (micro-channel plate) detector voltage 2400 V. Profiles with ESI⁺ were best for automated analysis.

The UHPLC-MS chromatograms were converted into the cdf format using the Databridge™ application provided with MassLynx™ (Waters, Guyancourt, France). Peak selection and peak alignment were performed with MzMine 2.10⁴⁸. The extract from cultured medium was included in the analysis to exclude its peaks from the profiles. The removal of peaks was performed manually, and in some cases, peaks that were too low to be detected were added manually if present because the objective of the profiling analysis was to detect biosynthetic pathways no matter their level of expression. Subsequently aligned and cleaned data was exported to Microsoft® Excel (2007) yielding a 6 (N) × 52 (X) data matrix. Peak height information was removed and replaced by 1 if compound was present and 0 if absent. Chemometric data analyses were done using Addinsoft® XLSTAT version 2013.4.08. Hierarchical clustering analysis (HCA) was performed to evaluate the Euclidean distance between all profiles.

Putative identification of the metabolites produced from six *Streptomyces cyaneofuscatus* strains was performed. Presumptive molecular formula and the difference between the exact mass and the experimental data (Δppm) were determined using MassLynx™. The identification of compounds was proposed using three databases (Antibase2012, Dictionary of Natural Products and Reaxys) focusing on natural compounds already isolated from microorganisms, fungi, sponges and some few algae based on the origin of these six bacterial strains.

**Data availability statement.** The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

**References**

1. Gibbons, S. M. & Gilbert, J. A. Microbial diversity-exploration of natural ecosystems and microbiomes. *Curr. Opin. Genet. Dev.* **35**, 66–72 (2015).

2. Yarza, P. et al. Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. *Nat. Rev. Microbiol.* **12**, 635–645 (2014).

3. Kim, M., Oh, H.-S., Park, S.-C. & Chun, J. Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. *Int. J. Syst. Evol. Microbiol.* **64**, 346–51 (2014).

4. Schloss, P. D. & Westcott, S. L. Assessing and improving methods used in operational taxonomic unit-based approaches for 16S rRNA gene sequence analysis. *Appl. Environ. Microbiol.* **77**, 3219–26 (2011).

5. Rosselló-Mora, R. & López-López, A. In *Accessing uncultivated microorganisms* 117–130, doi:10.1128/9781555815509.ch7 (American Society of Microbiology, 2008).

6. Langille, M. G. I. et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. *Nat. Biotechnol.* **31**, 814–821 (2013).

7. Schütze, E. et al. Growth of streptomyces in soil and their impact on bioremediation. *J. Hazard. Mater.* **267**, 128–35 (2014).

8. Schrey, S. D. & Tarkka, M. T. Friends and foes: streptomyces as modulators of plant disease and symbiosis. *Antonie Van Leeuwenhoek* **94**, 11–9 (2008).

9. Bontemps, C. et al. Taxonomic and functional diversity of *Streptomyces* in a forest soil. *FEMS Microbiol. Lett.* **342**, 157–67 (2013).

10. Seipke, R. F., Kaltenpoth, M. & Hutchings, M. I. *Streptomyces* as symbionts: an emerging and widespread theme? *FEMS Microbiol. Rev.* **36**, 862–76 (2012).

11. Cockell, C. S., Kelly, L. C. & Marteinsson, V. Actinobacteria — An ancient phylum active in colcanic rock weathering. *Geomicrobiol. J.* **30**, 706–720 (2013).

12. Newman, D. J. & Cragg, G. M. Natural products as sources of new drugs over the 30 years from 1981 to 2010. *J. Nat. Prod.* **75**, 311–35 (2012).

13. Davies, J. & Ryan, K. S. Introducing the parvome: bioactive compounds in the microbial world. *ACS Chem. Biol.* **7**, 252–9 (2012).

14. Andersson, D. I. & Hughes, D. Microbiological effects of sublthal levels of antibiotics. *Nat. Rev. Microbiol.* **12**, 465–78 (2014).

15. Burkholder, P. R. & Sun, S. H. Criteria of speciation in the genus *Streptomyces*. *Ann. N. Y. Acad. Sci.* **60**, 102–123 (1954).

16. Trejo, W. H. An evaluation of some concepts and criteria used in the speciation of streptomycetes. *Trans. N. Y. Acad. Sci.* **32**, 989–997 (1970).

17. Labeda, D. P. et al. Phylogenetic study of the species within the family *Streptomycetaceae*. *Antonie Van Leeuwenhoek* **101**, 73–104 (2012).

18. List of prokaryotic names with standing in nomenclature. at http://www.bacterio.net/streptomyces.html.

19. Goodfellow, M., Kumar, Y., Labeda, D. P. & Sembiring, I. The *Streptomyces violaceusniger* clade: a home for Streptomycetes with rugose ornamented spores. *Antonie Van Leeuwenhoek* **92**, 173–99 (2007).

20. Liu, Z. *et al.* Classification of *Streptomyces griseus* (Krainsky 1914) Waksman and Henrici 1948 and related species and the transfer of *Microstreptospora cinerea* to the genus *Streptomyces* as *Streptomyces yanii* sp. nov. *Int. J. Syst. Evol. Microbiol.* **55**, 1605–10 (2005).

21. Guo, Y., Zheng, W., Rong, X. & Huang, Y. A multilocus phylogeny of the *Streptomyces griseus* 16S rRNA gene clade: use of multilocus sequence analysis for streptomycete systematics. *Int. J. Syst. Evol. Microbiol.* **58**, 149–59 (2008).

22. Rong, X. & Huang, Y. In *Methods in Microbiology* **41**, 221–251 (Elsevier, 2014).

23. Kämpfer, P., Glaeser, S. P., Parkes, L., van Geertje, K. & Dyson, P. In *The Prokaryotes* (eds Rosenberg, E., DeLong, E. F., Lory, S., Stackebrandt, E. & Thompson, F.) 889–1010, doi:10.1007/978-3-642-30138-4 (Springer Berlin Heidelberg, 2014).

24. Antony-Babu, S., Stach, J. E. M. & Goodfellow, M. Genetic and phenotypic evidence for *Streptomyces griseus* ecovars isolated from a beach and dune sand system. *Antonie Van Leeuwenhoek* **94**, 63–74 (2008).

25. Parrot, D. et al. Littoral lichens as a novel source of potentially bioactive Actinobacteria. *Sci. Rep* (2015).

26. Antony-Babu, S. & Goodfellow, M. Biosystematics of alkaliphilic streptomycetes isolated from seven locations across a beach and dune sand system. *Antonie Van Leeuwenhoek* **94**, 581–91 (2008).
27. Stackebrandt, E. et al. Report of the ad hoc committee for the re-evaluation of the species definition in bacteriology. *Int. J. Syst. Evol. Microbiol.* **52**, 1043–7 (2002).

28. Rong, X., Liu, N., Ruan, J. & Huang, Y. Multilocus sequence analysis of *Streptomyces griseus* isolates delineating intraspecific diversity in terms of both taxonomy and biosynthetic potential. *Antonie Van Leeuwenhoek* **98**, 237–48 (2010).

29. Seipke, R. F. Strain-Level diversity of secondary metabolism in *Streptomyces albus*. *PLoS One* **10**, e0116457 (2015).

30. Jensen, P. R., Williams, P. G., Oh, D.-C., Zeigler, L. & Fenical, W. Species-specific secondary metabolite production in marine actinomycetes of the genus Salinispora. *Appl. Environ. Microbiol.* **73**, 1146–52 (2007).

31. Sawabe, T. et al. Updating the Vibrio clades defined by multilocus sequence phylogeny: proposal of eight new clades, and the description of Vibrio tritonius sp. nov. *Front. Microbiol.* **4**, 414 (2013).

32. Chisholm, S. W., Moore, L. R. & Rocap, G. Physiology and molecular phylogeny of coexisting *Prochlorococcus* ecotypes. *Nature* **393**, 464–467 (1998).

33. Goodfellow, M. & Fiedler, H.-P. A guide to successful bioprospecting: informed by actinobacterial systematics. *Antonie Van Leeuwenhoek* **98**, 119–42 (2010).

34. Baltz, R. H. Renaissance in antibacterial discovery from actinomycetes. *Curr. Opin. Pharmacol.* **8**, 557–563 (2008).

35. Darshanee Ruwandeepika, H. A. et al. Pathogenesis, virulence factors and virulence regulation of vibrios belonging to the Harveyi clade. *Rev. Aquac.* **4**, 59–74 (2012).

36. Hedlund, B. P. & Staley, J. T. Phylogeny of the genus *Simonsiella* and other members of the *Neisseriaceae*. *Int. J. Syst. Evol. Microbiol.* **52**, 1377–82 (2002).

37. Schlatter, D. C., DavelosBaines, A. L., Xiao, K. & Kinkel, L. L. Resource use of soilborne *Streptomyces* varies with location, phylogeny, and nitrogen amendment. *Microb. Ecol.* **66**, 961–71 (2013).

38. Schlatter, D. et al. Resource amendments influence density and competitive phenotypes of *Streptomyces* in soil. *Microb. Ecol.* **57**, 413–20 (2009).

39. Doroghazi, J. R. & Buckley, D. H. Widespread homologous recombination within and between *Streptomyces* species. *ISME J.* **4**, 1136–43 (2010).

40. Küster, E. Outline of a comparative study of criteria used in characterisation of the actinomycetes. *Int. Bull. Bacteriol. Nomencl. Taxon.* **9**, 97–104 (1959).

41. Williams, S. T., Davies, F. L. & H. D. M. In *The soil ecosystem* (ed. Sheals, J. G.) 107–117 (The Systematics Association, 1969).

42. Antony-Babu, S., Stach, J. E. M. & Goodfellow, M. Computer-assisted numerical analysis of colour-group data for dereplication of streptomyctes for bioprospecting and ecological purposes. *Antonie Van Leeuwenhoek* **97**, 231–9 (2010).

43. Shirling, E. B. & Gottlieb, D. Methods for characterization of *Streptomyces* species. *Int. J. Syst. Bacteriol.* **16**, 313–340 (1966).

44. Larkin, M. A. et al. Clustal W and Clustal X version 2.0. *Bioinformatics* **23**, 2947–8 (2007).

45. Ley, J. D., Cattoir, H. & Reynaerts, A. The quantitative measurement of DNA hybridization from renaturation rates. *Eur. J. Biochem.* **12**, 133–142 (1970).

46. Wayne, L. G. et al. Report of the ad hoc committee on reconciliation of approaches to bacterial systematics. *Int. J. Syst. Bacteriol.* **37**, 463–464 (1987).

47. Gonzalez, D. J. et al. Observing the invisible through imaging mass spectrometry, a window into the metabolic exchange patterns of microbes. *J. Proteomics* **75**, 5069–76 (2012).

48. Pluskal, T. et al. MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. *BMC Bioinformatics* **11**, 395 (2010).

## Acknowledgements

This work was partly supported by EMR, a partnership between the UPMC, the Laboratoires Pierre Fabre and CNRS and by the project MALICA (10-INBS-02-01) of the Agence Nationale de la Recherche, France and the INSA Rennes for the ministerial grant for the PhD to Delphine Parrot. We thank Alice Rodrigues and Clémence Rohée for fruitful discussions and technical advice. We thank the BIO2MAR platform (http://bio2mar OBS-banyuls.fr) for providing technical help and support and access to instrumentation.

## Author Contributions

The study was designed by S.A.-B., M.T.S. and D.S. Experimental work with the morphological, physiological and genetic characterization of the *S. cyaneofuscatus* strains was carried out by S.A.-B. Experimental work related to the metabolome characterization of the strains was carried out by S.A.-B., D.S. and V.E. Analysis of morphological, physiological and genetic dataset was done by S.A.-B. with help from M.T.S. Analysis of metabolome data was done by D.S., V.E., D.P. and S.T. The manuscript was written by S.A.-B., D.S. and M.T.S. and revised by all authors.

## Additional Information

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-11363-1

**Competing Interests:** The authors declare that they have no competing interests.

**Accession codes:** KX582197-KX582246.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017