# Massive Gene Flux Drives Genome Diversity between Sympatric Streptomyces Conspecifics

Abdoul-Razak Tidjani, Jean-Noël Lorenzi, Maxime Toussaint, Erwin van Dijk, Delphine Naquin, Olivier Lespinet, Cyril Bontemps, Pierre Leblond

► To cite this version:

Abdoul-Razak Tidjani, Jean-Noël Lorenzi, Maxime Toussaint, Erwin van Dijk, Delphine Naquin, et al.. Massive Gene Flux Drives Genome Diversity between Sympatric Streptomyces Conspecifics. mBio, 2019, 10.1128/mBio.01533-19. hal-02278804

HAL Id: hal-02278804

https://hal.univ-lorraine.fr/hal-02278804

Submitted on 4 Sep 2019

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
# Massive Gene Flux Drives Genome Diversity between Sympatric *Streptomyces* Conspecifics

Abdoul-Razak Tidjani,ª Jean-Noël Lorenzi,ª,b Maxime Toussaint,ª Erwin van Dijk,ᵇ Delphine Naquin,ᵇ Olivier Lespinet,ᵇ Cyril Bontemps,ª Pierre Leblondª

**ABSTRACT** In this work, by comparing genomes of closely related individuals of *Streptomyces* isolated at a spatial microscale (millimeters or centimeters), we investigated the extent and impact of horizontal gene transfer in the diversification of a natural *Streptomyces* population. We show that despite these conspecific strains sharing a recent common ancestor, all harbored significantly different gene contents, implying massive and rapid gene flux. The accessory genome of the strains was distributed across insertion/deletion events (indels) ranging from one to several hundreds of genes. Indels were preferentially located in the arms of the linear chromosomes (ca. 12 Mb) and appeared to form recombination hot spots. Some of them harbored biosynthetic gene clusters (BGCs) whose products confer an inhibitory capacity and may constitute public goods that can favor the cohesiveness of the bacterial population. Moreover, a significant proportion of these variable genes were either plasmid borne or harbored signatures of actinomycete integrative and conjugative elements (AICEs). We propose that conjugation is the main driver for the indel flux and diversity in *Streptomyces* populations.

**IMPORTANCE** Horizontal gene transfer is a rapid and efficient way to diversify bacterial gene pools. Currently, little is known about this gene flux within natural soil populations. Using comparative genomics of *Streptomyces* strains belonging to the same species and isolated at microscale, we reveal frequent transfer of a significant fraction of the pangenome. We show that it occurs at a time scale enabling the population to diversify and to cope with its changing environment, notably, through the production of public goods.

**KEYWORDS** *Streptomyces*, conjugation, gene transfer, plasticity, population

*Streptomyces* organisms are prominent soil-dwelling bacteria found in all terrestrial ecosystems, typified by their complex differentiation cycle and their central roles in biogeochemical cycles and the homeostasis of the soil. They are known as the most prolific bacterial genus for the production of specialized metabolites (e.g., antibiotics and antifungals) and enzymes of high biotechnological and medical interest (1). In soil, these compounds are involved in interactions with the surrounding microbial communities and contribute to plant health and growth (2). Each species or strain was reported to produce one or only a few compounds, but access to the whole-genome sequences revealed a much more important biosynthetic gene reservoir (3, 4).

The genome of *Streptomyces* is remarkably large (6 to 12 Mb) and possesses one of the few linear bacterial chromosomes ever described. The genome is highly compartmentalized, with genes conserved across the species distributed at the center of the chromosome and variable genes in the terminal parts (5, 6). The biosynthetic gene clusters (BGCs) are consistently enriched in the variable part of the chromosome. The proportion of the flexible genome is impressive and is reflected at the genus level by

**Citation** Tidjani A-R, Lorenzi J-N, Toussaint M, van Dijk E, Naquin D, Lespinet O, Bontemps C, Leblond P. 2019. Massive gene flux drives genome diversity between sympatric *Streptomyces* conspecifics. mBio 10:e01533-19. https://doi.org/10.1128/mBio.01533-19.

**Editor** Julian E. Davies, University of British Columbia

**Copyright** © 2019 Tidjani et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Cyril Bontemps, cyril.bontemps@univ-lorraine.fr, or Pierre Leblond, pierre.leblond@univ-lorraine.fr.

**Received** 13 June 2019

**Accepted** 25 July 2019

**Published** 3 September 2019

<footnote>ªUniversité de Lorraine, INRA, DynAMic, Nancy, France</footnote>

<footnote>ᵇInstitute for Integrative Biology of the Cell (I2BC), CEA, CNRS, University Paris-Sud, University Paris-Saclay, Gif-sur-Yvette, France</footnote>
a wide-open pangenome, i.e., the entire gene set of all species of the genus (7–9). The *Streptomyces* genomes can be completed by large linear and circular plasmids. Linear chromosomes and plasmids share the same invertron structure that consists of long terminal inverted repeats (TIRs) ending with telomeric sequences (10–12). In addition, the *Streptomyces* genomes are recombinogenic, as revealed either by the high frequency of spontaneous rearrangements of large DNA fragments (13) or by multilocus sequence analysis-based studies at the intra- and interspecific levels showing that recombination rates exceed those observed within many bacterial species (14).

Gene flux, i.e., the gain and loss of genetic material, crucially depends first on the access to exogenous information that can be achieved thanks to horizontal gene transfer (HGT) and second on the recombination capacities of the recipient genome that promotes the integration of genetic material in the genome (15). Regarding HGT, *Streptomyces* exhibits little talent for natural competence, and transducing phages have remained elusive (16, 17) despite extensive research. Thus, conjugation remains the only gene exchange mechanism described in *Streptomyces*, and the genomes of these bacteria are known to be rich in conjugative elements (integrated in the chromosome or of plasmid origin). Two types of conjugative elements coexist: integrated and conjugative element (ICE) depending on a type IV secretion system, and actinomycete ICE (AICE) whose transfer depends on a DNA translocase (TraB) (18, 19). The latter is the most widespread in *Streptomyces* and is capable of mobilizing chromosomal DNA in *trans*, i.e., not physically linked to the element (20).

These gene fluxes contribute to increase the flexibility of the genome and explain the diversity observed at the genus level. McDonald and Currie (21) recently showed that the *Streptomyces* genus is ancient (380 million years old) and that acquisition and retention of genes through HGT seem rare among nonclosely related lineages. Yet, nothing is known regarding the intensity and short-term impact of HGT at the population scale. Only some recent examples have shown that *Streptomyces* sharing an identical 16S rRNA gene sequences can exhibit differential phenotypes, including antimicrobial activities (22–24), suggesting rapid genome diversification.

Genome sequencing of closely related strains locally cooccurring has already revealed the existence of a large diversity in terms of gene content within different bacterial species (25–29). However, little is known regarding bacteria living in soil. Only a few examples are documented, and similarly, they seemed to indicate that variability in gene or allelic content can impact the ecology of soil bacteria. For instance, the polymorphism in quorum-sensing genes in *Bacillus subtilis* living in the same cubic centimeter impacted the ability to communicate and led to kin differentiation (30). In the same way, genomic comparison of two closely related groups of *Myxococcus* strains sampled at a centimeter scale enabled the identification of a 150-kb region involved in their sexual isolation (31). It seems thus relevant to assess genome diversification among soil inhabiting congenics, i.e., closely related strains, to unravel the genetic bases underlying the phenotypic diversity that supports niche differentiation and adaptation. The deeper the phylogenetic relationships (infraspecific level) are, the more recent the revealed events will be. Beyond analysis of 16S rRNA gene sequences, multilocus sequence analysis (MLSA), average nucleotide identity (ANI) across the genome, or phylogenomic tree reconstruction enable us to reach close phylogenetic relationships and to infer recent molecular events.

In soil ecosystems, the heterogeneity of the matrix constitutes a powerful driver of microbial growth and dynamics (32, 33). Biotic interactions in soil are most likely to occur between genetically related cells, produced by cell division, living close together on the same soil particle (34). Thus, bacteria in soil do not live as individual cells but are more likely structured in populations composed of closely related strains deriving from a recent common ancestor. The physical and phylogenetic proximities within these populations maximize the potentiality of gene exchanges. These latter reduce genetic isolation, disrupting the confinement of bacteria resulting from their reproduction scheme. However, the true extent of the gene flux that breaks clonality, and the impact of this on the function of natural bacterial populations, remains to be explored.
For that purpose, we undertook the first comparative genomics study of *Strepto-
myces* conspecifics. We favored the isolation of such closely related strains by sampling
*Streptomyces* at the microscale and sequenced and compared their genomes. Next, we
identified a high level of divergence in terms of the presence/absence of genes
between them and mapped the events along the linear chromosome. This allowed us
to identify recombination hot spots, some of which include BGCs. Next, we revealed the
high prevalence of actinomycete integrative and conjugative elements (AICEs). We then
consider the role of AICEs as motor of genome diversification and how this can later
contribute to population functioning.

RESULTS

**High levels of genetic diversity within a *Streptomyces* soil population.** To uncover microevolutionary processes, we isolated and sequenced the genomes of strains from a *Streptomyces* population inhabiting soil aggregates. For that purpose, four grains of soil, on the order of cubic millimeters in size and distant from each other from 2 cm to a maximum of 8 cm, were sampled from a clod of soil from which we isolated 129 strains. Based on identical 16S rDNA-coding gene sequences and a high degree of identity (>99.8%, nucleotide [nt]) measured with a five-gene multilocus sequence typing (MLST) scheme (see Fig. S1A in the supplemental material), we identified 32 highly related strains at the infraspecific level. The closest species to our group of strains was *Streptomyces olivochromogenes*, with in average 99.93% identity in the 16S rRNA gene and 94.5% ANIblast with the strain DSM40451. The genomes of 11 strains, representative of the different MLSTs identified in the cluster 6 of the population (Fig. S1A), were fully sequenced by a combination of Nanopore (Oxford Nanopore Technologies) and Illumina technologies (see Table S1). For each strain, a single scaffold for each replicon (i.e., chromosome and plasmids) was obtained with high-quality standards enabling full genome comparison. Phylogenomic tree construction (Fig. S1B) and calculation of ANIb percentages (ranging from 98.68% to 99.99%) between all isolates confirmed that they belonged to the same species and diverged recently from a common ancestor. The genomes ranged from 11.75 Mb to 12.44 Mb in size, positioning our strains among the biggest bacterial sequenced genomes (35). The chromosomes are linear with terminal inverted repeats (TIRs) of various sizes (from 311 kb to 587 kb) (see Table S2), as previously reported at the intra- and interspecific levels (12, 36). Following genome annotation (RAST), we determined that the pangenome of our population consisted of 13,814 genes (Fig. 1A). By extrapolating, we found that it reaches a maximum of 23,772 genes by 1,000 genomes. In contrast, for 11 strains randomly chosen in the genus, the pangenome was found to be approximately twice that of the population (23,672 versus 13,814) (Fig. 1B). Regarding the population pangenome, approximately one-third (5,036 genes) was not shared across the population and constituted the accessory genome (Fig. 1A). No two strains shared the same gene content; the closest differed by 12 predicted genes, while the two more distant isolates possessed 1,014 and 1,393 specific genes. Approximately one-quarter (27.6%; 1,392 genes) of the accessory genes were harbored on circular or linear extrachromosomal elements (98.8 kb to 394 kb), which were predicted *in silico* or experimentally observed by pulsed-field gel electrophoresis (PFGE) (see Fig. S2 and Table S2). These elements were found in single strains, and their mobile nature highlights the role of HGT as a driver of specificity, enabling the transfer of numerous genes simultaneously.

**Genome diversity is unevenly distributed in numerous indels along the linear chromosome.** Pairwise genome comparisons (n = 65) among the 11 sequenced genomes enabled us to identify disruptions in chromosomal synteny with insertion or deletion events (indels) of a minimum of one gene. The numbers of breaks differing between two strains ranged from 1 to 124 and increased linearly with the phylogenetic distance, indicating that they accumulate over time (see Fig. S3). When a break was shared by at least two strains, we considered it a single genetic event that may have been vertically inherited. This allowed consideration of 452 independent events that occurred during the population diversification (see Table S3). In extreme cases, the
<figure><img src="image_2.png" /><figcaption>**FIG 1** Pangenome analysis of the *Streptomyces* population. (A) Comparative Anvi'o genomic analysis of the 11 conspecific strains isolated in this study. The inner layers represent individual genomes organized regarding their phylogenetic relationships as indicated by the dendrogram. In the layers, dark colors indicate the presence of a gene group and light color its absence. The core (8,778 genes) and the accessory (5,036 genes) genomes are indicated in blue and pink, respectively, in the outmost layer. The blue layer represents the number of genomes among the population contributing to each gene group, and the green layer describes the gene groups in which at least one gene was functionally annotated using cluster of orthologous genes (COGs). (B) Comparison of the *Streptomyces* pangenome evolution at the genus and the population levels. The evolution of the pangenome of 59 complete *Streptomyces* genomes at the genus level is represented in blue. Its extrapolation to 1,000 genomes did not show any shift in the trend of the curve, indicating an open pangenome. The evolution of the population pangenome is represented in green. Its extrapolation rapidly reaches a plateau (see zoomed section).</figcaption></figure>

closest strains differed by only two events, while the most distinct strains differed by 235.

Comparative genomics of the 11 strains shows that indels were scattered along the chromosome but tended to form hot spots (Fig. 2A) with locations where at least three distinct gene assemblies were encountered (Fig. 3 and 4). Remarkably, indel distribution was heterogeneous along the chromosome, with increased density in the terminal regions of the chromosome, and was exaggerated as the phylogenetic distance of strains increased (e.g., RLB1-8 versus RLB3-6). Mapping the core genome of the *Streptomyces* genus (971 genes [J.N. Lorenzi, A. Thibessard, O. Lespinet, P. Leblond, unpublished]) identified the two distal positions, which delineate the core region from chromosomal arms of approximately 2.4 Mb and 1.4 Mb which are devoid of conserved genes on the left and right arms, respectively. The indel density was approximately 5-fold higher in the arms than in the core region (Fig. 2B). Moreover, the compilation of all the pairwise comparisons shows that the density of indels increases from these positions toward the ends of the chromosome, forming a gradient of events. Although the size of the indel in the core region appears larger than in the arms, the total numbers of variable genes follow the same trend, i.e., a strong increase in the arms. The highest variability was noticed at the very ends of the chromosome, as no two strains were identical in terms of the size of their terminal inverted repeats (Table S2). As previously reported, this variability might result from the formation of indels in the TIRs followed by chromosome arm replacement (36–38).

**Conjugative elements as the motor of rapid genome diversification.** Genome mining for the presence of integrative and conjugative elements in our population allowed the identification of one ICE and 25 AICEs. These 26 elements are likely to be functional (i.e., excisable and self-transmissible), since all the key functions (integrase
<figure><img src="image_3.png" /><figcaption>**FIG 2** Distribution of insertion and deletion events along the linear chromosome of *Streptomyces*. (A) The scheme at the top represents the *Streptomyces* chromosome with a megabase scale. The position of the core genome of the *Streptomyces* genus is highlighted by a light gray</figcaption></figure>

(Continued on next page)
<figure><img src="image_4.png" /><figcaption>**FIG 3** ICE/AICE insertion hot spot. The scheme is illustrating a hot spot for conjugative and integrative element insertion. Two different elements, AICE and ICE, are inserted in the tRNA<sub>Thr</sub> genes of strains RLB1-9 and RLB1-8, respectively. The tRNA<sub>Thr</sub> insertion site remains empty in RLB3-17. Direct repeats (DR) flanking the mobile element are represented in red (43 nucleotides for RLB1-9 and 42 nucleotides for RLB1-8). Only the key genes used to identify the elements are shown. Regions in light green correspond to syntenic regions.</figcaption></figure>

and Tra translocase) and *cis* sequences (flanking *att* sites) are present in each of them. Each genome has on average five elements, which make our population particularly rich compared to two elements found on average in a previous study on *Actinobacteria* (18). All the identified elements in our population were present in indels, meaning that they were not part of the core genome, highlighting their remarkable dynamics among the population. They were mainly inserted within the central region of the chromosome (Fig. 2B) and tended to form hots pots with several distinct elements being aggregated in the same target sequence. As illustrated in Fig. 3, the same tRNA target site was occupied by two distinct elements; while an AICE was present in RLB1-9, an ICE was observed in RLB1-8. In the third strain (RLB3-17), the target site was empty, which suggests that at least two of these strains experienced different conjugative events. While conjugative transfer is acknowledged as a powerful driver of diversity in telluric bacteria, their mobility in the soil remains difficult to demonstrate (39, 40). Here, the fact that each strain of the population is unique in terms of AICE/ICE content suggests that the exchanges are intense in the soil.

**HGT fosters social cohesiveness of bacterial population.** *Streptomyces* organisms are prolific antibiotic producers, and specialized metabolite (SM) production has already been reported to be highly variable among related strains of *Streptomyces* (22–24). We tested our *Streptomyces* isolates in pairwise inhibitory and resistance bioassays and found that none was able to inhibit its conspecifics. On the other hand, some but not all produced antimicrobial activity against a *Bacillus* strain isolated from the same soil, named "killer strains." To correlate the inhibitory phenotype with the presence/absence of biosynthetic gene clusters (BGCs), we performed comparative genomic analyses using antiSMASH (41). We identified 51 unique BGCs within the population, each strain possessing around 35 clusters (not shown). Interestingly, some

**FIG 2** Legend (Continued)

frame. The terminal inverted repeats (TIRs) are shown as red arrows. The bottom portion shows illustrations of pairwise genome comparisons (among the 10 possible pairs) within the population using strain RLB1-8 as a reference. Rectangles represent the linear chromosome of each strain. The strains are ordered from top to bottom relative to their phylogenetic distance to the reference. Triangles above the chromosome represent insertions in the reference strain, while triangles below correspond to insertions in the compared strain. For the sake of clarity, only the insertions of at least 10 predicted genes are shown. The height of a triangle reflects the number of genes involved in the insertion (the scale on the right of the chromosome indicates the number of genes). The colored triangles correspond to examples of insertion hot spots (the yellow and red triangles correspond to the hot spots depicted in Fig. 3 and 4, respectively). (B) The scheme represents the distribution of all the indel events identified within the population by genome pairwise comparisons. Each histogram bar corresponds to the number of indels within a 0.5-Mb window. The proportion of ICE/AICE insertions within a window is shown in black. The dotted blue line ("smiley" curve) corresponds to the polynomial trend curve (order 2, R² value) of the indel distribution. The red curve shows the number of genes involved in indels within a window.
<figure><img src="image_5.png" /><figcaption>**FIG 4** Identification of biosynthetic gene clusters in a variability hot spot. The scheme is illustrating a hot spot of variability observed between three strains (RLB3-17, RLB1-8, and RLB1-9). In two of them, different biosynthetic gene clusters were predicted by antiSMASH: one including a tRNA-dependent cyclodipeptide synthase in RLB1-8 and one including a nonribosomal peptide synthase in RLB1-9. Regions in light green correspond to syntenic regions.</figcaption></figure>

of these BGCs were present in a variability hot spot (Fig. 4), with three strains exhibiting different gene content at these locations. The strain RLB3-17 contained no BGC, while RLB1-8 harbored a cyclodipeptide synthase BGC and RLB1-9 a nonribosomal peptide synthetase (NRPS) BGC. The disruption of a key gene belonging to the NRPS BGC (Fig. 4) only present in the *Streptomyces* killer strains (including RLB1-9) abolished the inhibitory activity against *Bacillus* (see Fig. S4), proving that this BGC is responsible for the antimicrobial activity of the killer strains under these conditions.

## DISCUSSION

In this work, we isolated *Streptomyces* conspecifics at a soil microscale, and we revealed by genome comparison the extent of the massive gene flux that occurs in this natural population. This is reflected by a large pangenome, with almost clonal strains harboring strain-specific gene content. The size of the pangenome is influenced by three parameters: (i) the capacity to capture new DNA sequences, (ii) ecological selection of adaptive gene sets, and (iii) the evolutionary time enabling the accumulation of genome divergence (42). The fact that the population pangenome appears reduced in comparison with that of the genus (Fig. 1B) might reflect a constrained flux of genetic information within the microhabitat that is dependent of the diversity of the incoming information mostly arising from sister strains. However, the genetic diversity remained high in our population, implying an intense rate of gain and loss of genome regions over a short evolutionary time.

In *Streptomyces*, we have shown recently that repair of chromosomal double-strand breaks (DSB) experimentally induced in the terminal regions triggers the formation of chromosomal rearrangements (43). The formation of short insertions was also associated with DSB repair depending of nonhomologous end joining (NHEJ). Although the underlying mechanisms are mostly unknown, the plasticity of the subtelomeric regions appears as a hallmark of chromosome linearity in many organisms, from the limited range of bacterial taxa harboring linear chromosomes (44) to eukaryotes (45), including humans (46). In *Streptomyces*, tolerance to DNA rearrangements in the terminal regions could find its origins in the ancestral chromosomal linearization event. It is assumed that the ancestral actinomycete chromosome (i.e., harboring the ancestral core genes) was circular and linearized by recombination with an incoming linear replicon, leading to the addition of chromosomal arms harboring mostly contingency genes (47). However, if the ancestral gain of chromosomal arms can explain the tolerance to DNA rearrangements, it does not shed light on the evolutionary mechanisms that led to the contemporary variability of the *Streptomyces* genome (8, 9, 48). With respect to this, previous work has shown that chromosomal regions separating the core region from chromosomal arms were riddled with short indels (1 to 10 genes) forming a gradient
toward the terminal ends (5). We speculate that the indels, observed at the population level, when cumulated over evolutionary times, gave rise to the genetic compartmentalization observed within the genus (5).

Considering the high prevalence and plasticity of AICEs, we further propose that conjugative mechanisms are responsible for the massive gene flux observed in indels. To mobilize DNA, the translocase TraB binds to specific 8-mer repeats (named *clt* for *cis* acting locus of transfer) present on the AICE and assumed to be also present on the chromosome (*clt* like). Thus, beyond promoting self-transfer of the elements, TraB is able to transfer chromosomal markers in *trans* (20). This property was characterized and exploited in early studies on *Streptomyces* genetics (49). While conjugative transfer of chromosomal DNA was proved to be efficient under laboratory conditions (49, 50), the extent (frequency and size of transferred DNA stretches) remains to be elucidated. The indels ranging from 1 to 241 genes give an indication of the size of the chromosomal DNA stretches that may be transferred through conjugational processes. This is reminiscent of the distributive conjugal transfer in mycobacteria that enables the transfer of chromosomal fragments creating chimeric chromosomes in the recipient cells (51, 52). It is tempting to speculate that several DNA fragments could be transferred at once and generate multiple and diverse recombination events in the recipient *Streptomyces* hyphae.

At the ecological level, genome plasticity increases the functional potential at the strain but also the population level. For instance, it can enable the specific expression of specialized metabolites that can constitute public goods. Bacterial populations (i.e., sympatric closely related strains) can act as social units, meaning that their competitive interactions are more intense toward individuals not belonging to the population rather than between members of the population. Variability of SM production and resistance has already been reported as a driver of ecological cohesive units in *Vibrio* (53) and *Streptomyces* (23). Accordingly, differential SM production and resistance can be used as a proxy to understand the ecological dynamics of the bacterial population (53–55). Here, we confirm that the intense gene flux that *Streptomyces* conspecifics experienced promotes this variability, as previously proposed by Vetsigian et al. (23).

The production of SM within a bacterial population may underpin two principal ecological outcomes, bearing in mind that actual interactions in the soil may be more complex or different than under laboratory conditions. The production of antimicrobial activity and its respective resistance may benefit the sole carrier of the antibiotic production gene and mediate competition with conspecifics as well as with other species. Alternatively, an antibiotic produced by only some individuals of the population where the rest are resistant supports the hypothesis that antibiotics can constitute public goods, benefiting nonproducing but resistant conspecifics. Since none of our isolates was able to inhibit its conspecifics, it supports the hypothesis that they formed a cooperative population. The acquisition of the capacity to produce antimicrobial activity together with the capacity to resist its toxicity can result from a single HGT event thanks to the clustering of biosynthetic and resistance genes in the same genetic region, which is a typical trait of BGCs in actinomycetes. Therefore, the emergence of individuals able to overcome their conspecifics by producing a lethal activity is easily achieved by HGT. However, as we show in this work, our *Streptomyces* population experienced a massive and intense gene flux, which, in addition to favoring the acquisition of BGCs by individuals can also promote a rapid dissemination of all or parts of BGCs to the rest of the population. The replacement of a complete BGC obviously diversifies the antimicrobial arsenal, while its loss accompanied or not by that of the resistant determinant can result in a high intrapopulational SM variability. Further, recombination within or between BGCs may be enhanced thanks to their highly redundant structure (i.e., modular organization of NRPS and polyketide synthase [PKS]) and foster the emergence of new BGCs and the production of new antimicrobial
compounds. The intense turnover of the BGC and resistance gene sets would constantly modulate the interactions within the population and, through this dynamic state, favor the social cohesiveness of the population.

## MATERIALS AND METHODS

**Strain isolation and maintenance.** Four grains of soil on the order of cubic millimeters in size and distant from each other from 2 cm to a maximum of 8 cm were sampled from a clod of soil collected in the Montiers-sur-Saulx forest in France (GPS coordinates: 48°32′37.248″N, 5°18′21.946″E) and stored for 48 h at 4°C before processing. For bacterial isolation, each soil aggregate was dissolved 1:100 (wt/vol) in sterile water by vortexing for 15 min. Serial dilutions from 10⁻¹ to 10⁻³ were spread on *Streptomyces* isolation medium (SIM) (56), and bacterial colonies showing a typical *Streptomyces* phenotype were randomly picked after 7 days of incubation at 30°C. After three consecutive subcultures on mannitol soy agar plates (57), strains were stored as spore suspensions in 20% glycerol (57).

**Genome sequencing, annotation, and DNA manipulation.** Total genomic DNA was isolated as previously described by Kieser et al. (57). Each genome was sequenced by one-dimension MinION (Oxford Nanopore Technologies, UK) (mean coverage of 108×) and further corrected with MiSeq sequencing (Illumina, CA, USA) (mean coverage 50×), enabling the acquisition of each chromosome in one scaffold and the identification of extrachromosomal elements when present. Sequencing and assembling were performed via the I2BC NGS platform (France). Genome sequencing data and genome accession numbers are listed in Table S2 in the supplemental material. Coding sequence prediction and annotation were performed using the NCBI Prokaryotic Genome Annotation Pipeline (58). Genes involved in SM biosynthesis were predicted with antiSMASH (41). To disrupt a gene cluster putatively involved in the biosynthesis of a NRPS in strain RLB1-9, a 897-bp sequence internal to the peptide synthetase gene (open reading frame [ORF]-07845) was amplified and cloned into the suicide vector pIJ8668 to give pIJ8668-intORF-07845 (59). Intergeneric *Escherichia coli/Streptomyces* conjugation (60) followed by selection for apramycin resistance (pIJ8668) allowed us to select NRPS-encoding gene disruption by homologous recombination. Independent mutant clones (i.e., from distinct transformation assays) were tested for their antimicrobial activities. Pulsed-field gel electrophoresis (PFGE) was performed on whole genomic DNA in agarose plugs as described previously (61).

**Phylogenetic analyses.** 16S rRNA genes and MLST genes (*atpD*, *gyrB*, *recA*, *rpoB*, and *trpB*) were amplified using primers and PCR conditions described by Rintala et al. (62) and Guo et al. (63), respectively. Neighbor-joining phylogenetic trees were built using single or concatenated MLST gene (3,178 positions) sequences with a Kimura's two-parameter distance correction and 100 bootstrap replicates (MEGA7 software [64]). Phylogenomic analyses were performed with the nucleotidic sequences of the core genome identified by a Best BLAST Hit analysis (65), subsequently aligned with Muscle (v3.8.31 [66]), and corrected with Gblocks (67). From the 8,778 genes of the population core genome, we selected the 5,213 genes (5,149,602 nucleotide positions) shared with *Streptomyces avermitilis*, the reference strain used as a root in the phylogenetic reconstruction. Maximum likelihood phylogenomic tree (GTR, 100 bootstrap replicates) was built with RaxML (68). Python scripts (unpublished) were used to calculate the average nucleotide identity (ANIb) (69) and to identify and to extrapolate pangenomes. Pangenome visualization was performed using Anvi'o (70), and extrapolation graphs were computed with the UpsetR package in script R. Functional AICES and ICEs were identified first by searching signature genes (replication, transfer, excision, and integration) using BLASTP (65) as described by Bordeleau et al. (18) and then by identifying direct repeats flanking the element. Breaks of synteny were identified using an in-lab developed software (unpublished and available upon request).

**Antibacterial activity assays.** Drops of *Streptomyces* spore suspensions (5 µl of spore stock, 10⁷ spores per ml) were spotted and incubated for 5 days at 30°C on GA medium plates (54). Plates were overlaid with a Luria-Bertani top (0.4% agar), inoculated with *Bacillus* (optical density at 600 nm [OD₆₀₀] of 0.1), and incubated at 4°C for 2 h and then at 30°C overnight. Growth inhibition of *Bacillus* was assessed by the observation of a zone of clearance. To test *Streptomyces-Streptomyces* inhibitory capacities, spores of a first strain were spotted on mannitol soy medium plates for a 3-day incubation at 30°C. Then, a similar deposit of spores of the second strain (indicator) was adjacency spotted (5 mm) and incubated for 5 days at 30°C. The absence of phenotypic variation in comparison with control strains inoculated alone was considered a neutral interaction, and growth or sporulation alterations were recorded as antagonism.

**Accession number(s).** The genome accession numbers for the strains in this study are CP041650 to CP041654 and CP041601 to CP041613.

## SUPPLEMENTAL MATERIAL

Supplemental material for this article may be found at https://doi.org/10.1128/mBio .01533-19.

*FIG S1*, PDF file, 0.6 MB.

*FIG S2*, PDF file, 0.4 MB.

*FIG S3*, PDF file, 0.4 MB.

*FIG S4*, PDF file, 0.4 MB.
TABLE S1, PDF file, 0.4 MB.

TABLE S2, PDF file, 0.5 MB.

TABLE S3, PDF file, 0.8 MB.

## ACKNOWLEDGMENTS

This work was supported by a grant overseen by the French National Research Agency (ANR) as part of the Investissements d'Avenir program (ANR-11-LABX-0002-01, Lab of Excellence ARBRE).

We thank the research system SOERE-OPE (Système d'Observation et d'Experimentation au long terme pour la Recherche en Environnement, Observatoire Pérenne de l'Environnement) for access to the sampling site.

## REFERENCES

1. Barka EA, Vatsa P, Sanchez L, Gaveau-Vaillant N, Jacquard C, Meier-Kolthoff JP, Klenk H-P, Clément C, Ouhdouch Y, van Wezel GP. 2016. Taxonomy, physiology, and natural products of *Actinobacteria*. Microbiol Mol Biol Rev 80:1–43. https://doi.org/10.1128/MMBR.00019-15.

2. Olanrewaju OS, Babalola OO. 2019. *Streptomyces*: implications and interactions in plant growth promotion. Appl Microbiol Biotechnol 103: 1179–1188. https://doi.org/10.1007/s00253-018-09577-y.

3. Aigle B, Lautru S, Spiteller D, Dickschat JS, Challis GL, Leblond P, Pernodet J-L. 2014. Genome mining of *Streptomyces ambofaciens*. J Ind Microbiol Biotechnol 41:251–263. https://doi.org/10.1007/s10295-013-1379-y.

4. Baltz RH. 2017. Gifted microbes for genome mining and natural product discovery. J Ind Microbiol Biotechnol 44:573–588. https://doi.org/10.1007/s10295-016-1815-x.

5. Choulet F, Aigle B, Gallois A, Mangenot S, Gerbaud C, Truong C, Francou F-X, Fourrier C, Guérineau M, Decaris B, Barbe V, Pernodet J-L, Leblond P. 2006. Evolution of the terminal regions of the *Streptomyces* linear chromosome. Mol Biol Evol 23:2361–2369. https://doi.org/10.1093/molbev/msl108.

6. Redenbach M, Kieser HM, Denapaite D, Eichner A, Cullum J, Kinashi H, Hopwood DA. 1996. A set of ordered cosmids and a detailed genetic and physical map for the 8 Mb *Streptomyces coelicolor* A3(2) chromosome. Mol Microbiol 21:77–96. https://doi.org/10.1046/j.1365-2958.1996.6191336.x.

7. 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.

8. Kim J-N, Kim Y, Jeong Y, Roe J-H, Kim B-G, Cho B-K. 2015. Comparative genomics reveals the core and accessory genomes of *Streptomyces* species. J Microbiol Biotechnol 25:1599–1605. https://doi.org/10.4014/jmb.1504.04008.

9. Zhou Z, Gu J, Li Y-Q, Wang Y. 2012. Genome plasticity and systems evolution in *Streptomyces*. BMC Bioinformatics 13 Suppl 10:S8. https://doi.org/10.1186/1471-2105-13-S10-S8.

10. Lin YS, Kieser HM, Hopwood DA, Chen CW. 1993. The chromosomal DNA of *Streptomyces lividans* 66 is linear. Mol Microbiol 10:923–933. https://doi.org/10.1111/j.1365-2958.1993.tb00964.x.

11. Huang CH, Lin YS, Yang YL, Huang SW, Chen CW. 1998. The telomeres of *Streptomyces* chromosomes contain conserved palindromic sequences with potential to form complex secondary structures. Mol Microbiol 28:905–916. https://doi.org/10.1046/j.1365-2958.1998.00856.x.

12. Chen CW, Huang C-H, Lee H-H, Tsai H-H, Kirby R. 2002. Once the circle has been broken: dynamics and evolution of *Streptomyces* chromosomes. Trends Genet 18:522–529. https://doi.org/10.1016/S0168-9525(02)02752-X.

13. Thibessard A, Leblond P. 2014. Subtelomere plasticity in the bacterium *Streptomyces*, p 243–258. In Louis EJ, Becker MM (ed), Subtelomeres. Springer, Berlin, Germany.

14. Doroghazi JR, Buckley DH. 2010. Widespread homologous recombination within and between *Streptomyces* species. ISME J 4:1136–1143. https://doi.org/10.1038/ismej.2010.45.

15. Ochman H, Lawrence JG, Groisman EA. 2000. Lateral gene transfer and the nature of bacterial innovation. Nature 405:299–304. https://doi.org/10.1038/35012500.

16. Stuttard C. 1982. Temperate phages of *Streptomyces venezuelae*: lysogeny and host specificity shown by phages SV1 and SV2. Microbiology 128:115–121. https://doi.org/10.1099/00221287-128-1-115.

17. Burke J, Schneider D, Westpheling J. 2001. Generalized transduction in *Streptomyces coelicolor*. Proc Natl Acad Sci U S A 98:6289–6294. https://doi.org/10.1073/pnas.101589398.

18. Bordeleau E, Ghinet MG, Burrus V. 2012. Diversity of integrating conjugative elements in actinobacteria: coexistence of two mechanistically different DNA-translocation systems. Mob Genet Elements 2:119–124. https://doi.org/10.4161/mge.20498.

19. Thoma L, Muth G. 2016. Conjugative DNA-transfer in *Streptomyces*, a mycelial organism. Plasmid 87–88:1–9. https://doi.org/10.1016/j.plasmid.2016.09.004.

20. Pettis GS. 2018. Spreading the news about the novel conjugation mechanism in *Streptomyces* bacteria. Environ Microbiol Rep 10:503–510. https://doi.org/10.1111/1758-2229.12659.

21. McDonald BR, Currie CR. 2017. Lateral gene transfer dynamics in the ancient bacterial genus *Streptomyces*. mBio 8:e00644-17. https://doi.org/10.1128/mBio.00644-17.

22. Antony-Babu S, Stien D, Eparvier V, Parrot D, Tomasi S, Suzuki MT. 2017. Multiple *Streptomyces* species with distinct secondary metabolomes have identical 16S rRNA gene sequences. Sci Rep 7:11089. https://doi.org/10.1038/s41598-017-11363-1.

23. Vetsigian K, Jajoo R, Kishony R. 2011. Structure and evolution of *Streptomyces* interaction networks in soil and in silico. PLoS Biol 9:e1001184. https://doi.org/10.1371/journal.pbio.1001184.

24. Vicente CM, Thibessard A, Lorenzi J-N, Benhadj M, Hôtel L, Gacemi-Kirane D, Lespinet O, Leblond P, Aigle B. 2018. Comparative genomics among closely related *Streptomyces* strains revealed specialized metabolite biosynthetic gene cluster diversity. Antibiotics (Basel) 7:e86. https://doi.org/10.3390/antibiotics7040086.

25. Harris HMB, Bourin MJB, Claesson MJ, O'Toole PW. 2017. Phylogenomics and comparative genomics of *Lactobacillus salivarius*, a mammalian gut commensal. Microb Genom 3:e000115. https://doi.org/10.1099/mgen.0000115.

26. Park CJ, Andam CP. 2019. Within-species genomic variation and variable patterns of recombination in the tetracycline producer *Streptomyces rimosus*. Front Microbiol 10:552. https://doi.org/10.3389/fmicb.2019.00552.

27. Croucher NJ, Coupland PG, Stevenson AE, Callendrello A, Bentley SD, Hanage WP. 2014. Diversification of bacterial genome content through distinct mechanisms over different timescales. Nat Commun 5:5471. https://doi.org/10.1038/ncomms6471.

28. Levade I, Terrat Y, Leducq J-B, Weil AA, Mayo-Smith LM, Chowdhury F, Khan AI, Boncy J, Buteau J, Ivers LC, Ryan ET, Charles RC, Calderwood SB, Qadri F, Harris JB, LaRocque RC, Shapiro BJ. 2017. *Vibrio cholerae* genomic diversity within and between patients. Microb Genom 3:e000142. https://doi.org/10.1099/mgen.0000142.

29. Bruns H, Crüsemann M, Letzel A-C, Alanjary M, McInerney JO, Jensen PR, Schulz S, Moore BS, Ziemert N. 2018. Function-related replacement of bacterial siderophore pathways. ISME J 12:320–329. https://doi.org/10.1038/ismej.2017.137.

30. Stefanic P, Kraigher B, Lyons NA, Kolter R, Mandic-Mulec I. 2015. Kin discrimination between sympatric *Bacillus subtilis* isolates. Proc Natl Acad Sci U S A 112:14042–14047. https://doi.org/10.1073/pnas.1512671112.
31. Wielgoss S, Didelot X, Chaudhuri RR, Liu X, Weedall GD, Velicer GJ, Vos M. 2016. A barrier to homologous recombination between sympatric strains of the cooperative soil bacterium *Myxococcus xanthus*. ISME J 10:2468-2477. https://doi.org/10.1038/ismej.2016.34.

32. Vos M, Wolf AB, Jennings SJ, Kowalchuk GA. 2013. Micro-scale determinants of bacterial diversity in soil. FEMS Microbiol Rev 37:936-954. https://doi.org/10.1111/1574-6976.12023.

33. Wolf AB, Vos M, de Boer W, Kowalchuk GA. 2013. Impact of matric potential and pore size distribution on growth dynamics of filamentous and non-filamentous soil bacteria. PLoS One 8:e83661. https://doi.org/10.1371/journal.pone.0083661.

34. Raynaud X, Nunan N. 2014. Spatial ecology of bacteria at the microscale in soil. PLoS One 9:e87217. https://doi.org/10.1371/journal.pone.0087217.

35. Sayers EW, Agarwala R, Bolton EE, Brister JR, Canese K, Clark K, Connor R, Fiorini N, Funk K, Hefferon T, Holmes JB, Kim S, Kimchi A, Kitts PA, Lathrop S, Lu Z, Madden TL, Marchler-Bauer A, Phan L, Schneider VA, Schoch CL, Pruitt KD, Ostell J. 2019. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res 47:D23-D28. https://doi.org/10.1093/nar/gky1069.

36. Wenner T, Roth V, Fischer G, Fourrier C, Aigle B, Decaris B, Leblond P. 2003. End-to-end fusion of linear deleted chromosomes initiates a cycle of genome instability in *Streptomyces ambofaciens*. Mol Microbiol 50: 411-425. https://doi.org/10.1046/j.1365-2958.2003.03698.x.

37. Fischer G, Wenner T, Decaris B, Leblond P. 1998. Chromosomal arm replacement generates a high level of intraspecific polymorphism in the terminal inverted repeats of the linear chromosomal DNA of *Streptomyces ambofaciens*. Proc Natl Acad Sci U S A 95:14296-14301. https://doi.org/10.1073/pnas.95.24.14296.

38. Uchida T, Miyawaki M, Kinashi H. 2003. Chromosomal arm replacement in *Streptomyces griseus*. J Bacteriol 185:1120-1124. https://doi.org/10.1128/jb.185.3.1120-1124.2003.

39. Clerc S, Simonet P. 1996. Efficiency of the transfer of a pSAM2-derivative plasmid between two strains of *Streptomyces lividans* in conditions ranging from agar slants to non-sterile soil microcosms. FEMS Microbiol Ecol 21:157-165. https://doi.org/10.1111/j.1574-6941.1996.tb00343.x.

40. Sullivan JT, Ronson CW. 1998. Evolution of rhizobia by acquisition of a 500-kb symbiosis island that integrates into a phe-tRNA gene. Proc Natl Acad Sci U S A 95:5145-5149. https://doi.org/10.1073/pnas.95.9.5145.

41. Blin K, Wolf T, Chevrette MG, Lu X, Schwalen CJ, Kautsar SA, Suarez Duran HG, de los Santos ELC, Kim HU, Nave M, Dickschat JS, Mitchell DA, Shelest E, Breitling R, Takano E, Lee SY, Weber T, Medema MH. 2017. antiSMASH 4.0-improvements in chemistry prediction and gene cluster boundary identification. Nucleic Acids Res 45:W36-W41. https://doi.org/10.1093/nar/gkx319.

42. McInerney JO, McNally A, O'Connell MJ. 2017. Why prokaryotes have pangenomes. Nat Microbiol 2:17040. https://doi.org/10.1038/nmicrobiol.2017.40.

43. Hoff G, Bertrand C, Piotrowski E, Thibessard A, Leblond P. 2018. Genome plasticity is governed by double strand break DNA repair in *Streptomyces*. Sci Rep 8:5272. https://doi.org/10.1038/s41598-018-23622-w.

44. Casjens S, Palmer N, van Vugt R, Huang WM, Stevenson B, Rosa P, Lathigra R, Sutton G, Peterson J, Dodson RJ, Haft D, Hickey E, Gwinn M, White O, Fraser CM. 2002. A bacterial genome in flux: the twelve linear and nine circular extrachromosomal DNAs in an infectious isolate of the Lyme disease spirochete *Borrelia burgdorferi*. Mol Microbiol 35:490-516. https://doi.org/10.1046/j.1365-2958.2000.01698.x.

45. Ricchetti M, Dujon B, Fairhead C. 2003. Distance from the chromosome end determines the efficiency of double strand break repair in subtelomers of haploid yeast. J Mol Biol 328:847-862. https://doi.org/10.1016/S0022-2836(03)00315-2.

46. Linardopoulou EV, Williams EM, Fan Y, Friedman C, Young JM, Trask BJ. 2005. Human subtelomers are hot spots of interchromosomal recombination and segmental duplication. Nature 437:94-100. https://doi.org/10.1038/nature04029.

47. Bentley SD, Chater KF, Cerdeño-Tárraga A-M, Challis GL, Thomson NR, James KD, Harris DE, Quail MA, Kieser H, Harper D, Bateman A, Brown S, Chandra G, Chen CW, Collins M, Cronin A, Fraser A, Goble A, Hidalgo J, Hornsby T, Howarth S, Huang C-H, Kieser T, Larke L, Murphy L, Oliver K, O'Neil S, Rabbinowitsch E, Rajandream M-A, Rutherford K, Rutter S, Seeger K, Saunders D, Sharp S, Squares R, Squares S, Taylor K, Warren T, Wietzorrek A, Woodward J, Barrell BG, Parkhill J, Hopwood DA. 2002.

Complete genome sequence of the model actinomycete *Streptomyces coelicolor* A3(2). Nature 417:141-147. https://doi.org/10.1038/417141a.

48. Hsiao N-H, Kirby R. 2008. Comparative genomics of *Streptomyces avermitilis*, *Streptomyces cattleya*, *Streptomyces maritimus* and *Kita-satospora aureofaciens* using a *Streptomyces coelicolor* microarray system. Antonie Van Leeuwenhoek 93:1-25. https://doi.org/10.1007/s10482-007-9175-1.

49. Hopwood DA. 2006. Soil to genomics: the *Streptomyces* chromosome. Annu Rev Genet 40:1-23. https://doi.org/10.1146/annurev.genet.40.110405.090639.

50. Hopwood DA, Kieser T, Wright HM, Bibb MJ. 1983. Plasmids, recombination and chromosome mapping in *Streptomyces lividans* 66. J Gen Microbiol 129:2257-2269. https://doi.org/10.1099/00221287-129-7-2257.

51. Gray TA, Krywy JA, Harold J, Palumbo MJ, Derbyshire KM. 2013. Distributive conjugal transfer in mycobacteria generates progeny with meiotic-like genome-wide mosaicism, allowing mapping of a mating identity locus. PLoS Biol 11:e1001602. https://doi.org/10.1371/journal.pbio.1001602.

52. Derbyshire KM, Gray TA. 2014. Distributive conjugal transfer: new insights into horizontal gene transfer and genetic exchange in mycobacteria. Microbiol Spectr 2:MGM2-0022-2013. https://doi.org/10.1128/microbiolspec.MGM2-0022-2013.

53. Cordero OX, Wildschutte H, Kirkup B, Proehl S, Ngo L, Hussain F, Le Roux F, Mincer T, Polz MF. 2012. Ecological populations of bacteria act as socially cohesive units of antibiotic production and resistance. Science 337:1228-1231. https://doi.org/10.1126/science.1219385.

54. Abrudan MI, Smakman F, Grimbergen AJ, Westhoff S, Miller EL, van Wezel GP, Rozen DE. 2015. Socially mediated induction and suppression of antibiosis during bacterial coexistence. Proc Natl Acad Sci U S A 112:11054-11059. https://doi.org/10.1073/pnas.1504076112.

55. Butaité E, Baumgartner M, Wyder S, Kümmerli R. 2017. Siderophore cheating and cheating resistance shape competition for iron in soil and freshwater *Pseudomonas* communities. Nat Commun 8:414. https://doi.org/10.1038/s41467-017-00509-4.

56. D'Costa VM, McGrann KM, Hughes DW, Wright GD. 2006. Sampling the antibiotic resistome. Science 311:374-377. https://doi.org/10.1126/science.1120800.

57. Kieser T, Bibb MJ, Buttner MJ, Chater KF, Hopwood DA. 2000. Practical Streptomyces genomics. John Innes Foundation, Norwich, United Kingdom.

58. 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:6614-6624. https://doi.org/10.1093/nar/gkw569.

59. Sun J, Kelemen GH, Fernández-Abalos JM, Bibb MJ. 1999. Green fluorescent protein as a reporter for spatial and temporal gene expression in *Streptomyces coelicolor* A3(2). Microbiology 145:2221-2227. https://doi.org/10.1099/00221287-145-9-2221.

60. Kim M-K, Ha H-S, Choi S-U. 2008. Conjugal transfer using the bacteriophage phiC31 att/int system and properties of the attB site in *Streptomyces ambofaciens*. Biotechnol Lett 30:695-699. https://doi.org/10.1007/s10529-007-9586-0.

61. Leblond P, Francou FX, Simonet J-M, Decaris B. 1990. Pulsed-field gel electrophoresis analysis of the genome of *Streptomyces ambofaciens* strains. FEMS Microbiol Lett 72:79-88. https://doi.org/10.1111/j.1574-6968.1990.tb03866.x.

62. Rintala H, Nevalainen A, Rönkä E, Suutari M. 2001. PCR primers targeting the 16S rRNA gene for the specific detection of streptomyces. Mol Cell Probes 15:337-347. https://doi.org/10.1006/mcpr.2001.0379.

63. Guo Y, Zheng W, Rong X, Huang Y. 2008. 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-159. https://doi.org/10.1099/ijss.0.65224-0.

64. Kumar S, Stecher G, Tamura K. 2016. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol 33: 1870-1874. https://doi.org/10.1093/molbev/msw054.

65. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. 1990. Basic local alignment search tool. J Mol Biol 215:403-410. https://doi.org/10.1016/S0022-2836(05)80360-2.

66. Edgar RC. 2004. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5:113. https://doi.org/10.1186/1471-2105-5-113.

67. Talavera G, Castresana J. 2007. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequences.
quence alignments. Syst Biol 56:564–577. https://doi.org/10.1080/10635150701472164.

68. Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30:1312–1313. https://doi.org/10.1093/bioinformatics/btu033.

69. Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. 2007. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol 57:81–91. https://doi.org/10.1099/ijs.0.64483-0.

70. Eren AM, Esen ÖC, Quince C, Vineis JH, Morrison HG, Sogin ML, Delmont TO. 2015. Anvi’o: an advanced analysis and visualization platform for ‘omics data. PeerJ 3:e1319. https://doi.org/10.7717/peerj.1319.