Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

Comparative genomics of actinomycetes with a focus on natural product biosynthetic genes

James R Doroghazi1* and William W Metcalf12

Author Affiliations

1 Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Champaign, IL 61801, USA

2 Department of Microbiology, University of Illinois at Urbana-Champaign, Champaign, IL 61801, USA

For all author emails, please log on.

BMC Genomics 2013, 14:611  doi:10.1186/1471-2164-14-611

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/14/611


Received:2 April 2013
Accepted:4 September 2013
Published:11 September 2013

© 2013 Doroghazi and Metcalf; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Actinomycetes are a diverse group of medically, industrially and ecologically important bacteria, studied as much for the diseases they cause as for the cures they hold. The genomes of actinomycetes revealed that these bacteria have a large number of natural product gene clusters, although many of these are difficult to tie to products in the laboratory. Large scale comparisons of these clusters are difficult to perform due to the presence of highly similar repeated domains in the most common biosynthetic machinery: polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs).

Results

We have used comparative genomics to provide an overview of the genomic features of a set of 102 closed genomes from this important group of bacteria with a focus on natural product biosynthetic genes. We have focused on well-represented genera and determine the occurrence of gene cluster families therein. Conservation of natural product gene clusters within Mycobacterium, Streptomyces and Frankia suggest crucial roles for natural products in the biology of each genus. The abundance of natural product classes is also found to vary greatly between genera, revealing underlying patterns that are not yet understood.

Conclusions

A large-scale analysis of natural product gene clusters presents a useful foundation for hypothesis formulation that is currently underutilized in the field. Such studies will be increasingly necessary to study the diversity and ecology of natural products as the number of genome sequences available continues to grow.

Keywords:
Actinomycetes; Natural products; Genomics; Secondary metabolism

Background

The class Actinobacteria is the largest within the phylum Actinobacteria and contains many bacteria relevant to human health and industry (see [1] for review). These bacteria are Gram-positive with genomic GC content generally over 55%. Some of them, such as the Streptomyces, were originally mistaken for fungi, as evidenced by the name of the group (myces is derived from the Greek word for fungus) and were once considered relatives of fungi based on morphology and life cycle. The existence of a life cycle involving multiple, distinct stages and morphologies has also made some actinomycetes, such as “Streptomyces coelicolor” A3(2), important model systems for studying differentiation and the signaling pathways involved therein.

The class Actinobacteria, or the actinomycetes, contains both the most deadly bacterial pathogen and the organisms that are the most important for antibiotic production. Mycobacterium tuberculosis is the second leading cause of death worldwide due to an infectious agent (after HIV/AIDS [2]), while the genus Streptomyces is the source of over half of the bioactive metabolites from bacteria [3]. The genus Corynebacterium contains deadly pathogens but also includes non-pathogens that are the leading producers of L-amino acids, which represent some of the most important microbial products in terms of both volume and value [4]. Numerous other pathogens and pharmaceutical producers, as well as ecologically and industrially important taxa are also found among this important microbial group.

Actinomycetes have historically been a leading source for natural product discovery [5]. These compounds, also called secondary metabolites, have a wide range of industrial uses, including as antineoplastic, antifungal, antimicrobial, herbicidal and plant growth promoting agents. They are also important components of iron-acquisition systems and signaling molecules important for development. Production of secondary metabolites may also be important adaptations to environments such as soil, and may aid competition for resources such as plant matter. Whatever their use, the genes that are responsible for production of individual secondary metabolites are almost always located together in the genome and are referred to as biosynthetic gene clusters. The co-location and horizontal transfer of these gene clusters is fascinating in and of itself, but is also a trait that aids in discovery, characterization and comparison of the genes responsible for secondary metabolite biosynthesis (see [6,7] for an overview and discussion of evolutionary implications).

Many researchers have voiced optimism that genome mining for novel secondary metabolites will result in a renaissance of discovery and fill the innovation gap that has left the pipelines at low levels [8-10]. The main reason for this is that Streptomyces and related genera, the traditional focus of discovery, rarely express their full inventory of chemical weapons when cultivated in the lab. For example, “Streptomyces coelicolor” A3(2) was a genetic workhorse for some 40 years before having its genome sequenced and was known to make only four secondary metabolites. The genome sequence revealed an additional 18 biosynthetic gene clusters [11]. Biosynthetic gene clusters which are present but not known to produce any secondary metabolites are referred to as cryptic clusters. There have been no systematic studies to date, however, on whether a cryptic biosynthetic gene cluster in one species is also likely to be cryptic in a second species, and therefore the fraction of undiscovered secondary metabolites based solely on genetic capacity may tend to overestimate the number of pathways that are cryptic. With this in mind, being able to classify and compare biosynthetic gene clusters, and thus systematically catalog the extent of natural product diversity, is an important first step towards a full exploitation of secondary metabolites in bacteria. This is, however, a difficult bioinformatics task for the two most common classes of natural products, type I polyketide synthases (PKS), and nonribosomal peptide synthetases (NPRS), due to the multiple similar domains present in both (see [12] for a review).

Currently, there are six actinomycete genera with sufficient numbers of completed genomes to allow an in-depth analysis of secondary metabolic diversity. We compared the genomes within these six, Mycobacterium, Corynebacterium, Rhodococcus, Arthrobacter, Frankia, and Streptomyces, in detail to determine the extent to which natural product gene clusters are conserved within each genus. We also present a broad, genome-scale comparison of complete genomes across the class Actinobacteria.

Methods

All genomes were downloaded from NCBI on September 21, 2011. An attempt was made to include all species for which publicly available closed genomes were available within the order Actinomycetales as shown within NCBI taxonomy browser, although this taxonomic group has been re-ordered recently to compose the class Actinobacteria[1]. Plasmids were omitted from the analysis to prevent skewing long term evolutionary trends. Predicted proteins were used as annotated, and an all-v-all BLAST comparison was performed using BLAST v2.2.26+ [13].

Phylogeny and whole genome comparisons

OrthoMCL version 2.0 with default settings was used for further analysis of BLAST results [14]. OrthoMCL similarity groups with “S. coelicolor” A3(2) genes annotated as ribosomal proteins were used for phylogenetic analysis. Only ribosomal protein genes in similarity groups containing a single gene from each species were used for this analysis. The complete list of genes used is: L1, L2, L3, L4, L5, L6, L7/L12, L9, L10, L11, L13, L14, L15, L16, L17, L18, L19, L20, L21, L22, L23, L24, L25p, L27, L29, L35, S1, S3, S5, S6, S7, S8, S9, S10, S11, S12, S13, S15, S17, S19, S20. The amino acid sequences of these genes were aligned with Clustal W 1.83 [15] and concatenated for phylogenetic analysis. The concatenated gene tree was made using FastTree 2.1.5 run with the Gamma20 model [16]. A NeighborNet network was created using the same data in the program SplitsTree 4.11.3 [17].

Groups of similar genes as output by OrthoMCL were parsed with custom Perl scripts to calculate pairwise genome similarity. Similarity was calculated as Sij/Gi, where Sij is the number of similar genes between genomes i and j, and Gi is the total number of genes in genome i. When multiple genes from the organisms being compared appeared in one similarity group, the count for number of similar genes was determined by whichever genome has fewer copies. Dividing by the total number of genes in only one genome means that there are two similarity measures presented for each pairwise comparison.

Biosynthetic gene cluster discovery and comparison

Signature enzymes for major classes of secondary metabolites were found using profile Hidden Markov Models (pHMMs) and the program HMMER [18]. The pHMMs used are a mixture of those reported by Medema et al. [19] with the same cut-offs mentioned therein for PKS I, PKS II, PKS III, NRPS, indolocarbazoles, aerobactin-like siderophores, butyrolactones, aminoglycosides, and β-lactams, including screening for fatty acid synthases that are hit by the PKS models. New pHMMs were made for discovery of terpene synthases based on the sequences published in [20], lanthipeptides based on the required cyclase domain, see [21] for review, and thiazole-oxazole modified microcins, or TOMMs based on the YcaO domain [22]. The new pHMMs and alignments are presented in a stand-alone website (see Additional file 1). Phosphonates were found using a BLAST search and screening for sequences containing the EDK-X(5)-NS motif present in all verified PepM sequences (see [23] for review). Gene clusters were defined by extending six genes to either side of a significant pHMM hit (past the specified cut-off), joining additional hits within that window into the same cluster, and re-initiating the six gene count after encountering additional hits. The six gene extension was a practical choice; when we defined gene clusters with longer extensions the comparisons included more noise (divergent genomic neighborhoods not related to biosynthetic genes), and fewer genes in each cluster resulted in too little data for comparisons. This choice was made with future automation in mind. Similar gene clusters were found using an array of tools including phylogenetic comparisons and Mauve [24] alignments after concatenation of all gene clusters in each strain into one sequence. A website showing all gene clusters are included as Additional file 1. Gene cluster diagrams also include domain annotations, but these are not manually curated and some domains are incorrectly split in half. Gene annotation and domain names are available on mouseover.

Additional file 1. A stand-alone website showing all natural product gene clusters analyzed in this study, along with separate files for conserved clusters mentioned in the text and pHMM files. Use of the HTML files requires Javascript. Homologous genes are shown in the same color. All homologous genes on a page are highlighted upon mouseover of any of them. Mouseover also produces a description containing the locus tag and annotation for each gene. Mouseover for a domain box above the gene arrows shows the domain name. Clicking on a gene arrow produces a page with the amino acid sequence and a link to BLAST the nr protein database.

Format: ZIP Size: 6.7MB Download fileOpen Data

Results and discussion

102 closed actinomycete genomes were grouped into seven broad categories according to isolation source, smear-ripened cheese being the most narrowly defined (Figure 1). The two most common isolation sources for actinomycetes are animal hosts and soil, although recently marine actinomycetes have garnered significant interest. Obligate pathogens, which by definition live in a well-defined and constant niche, tend to have undergone genome reduction, a trend not limited to actinomycetes [25]. Bacteria that dwell in soil, a very diverse and changing habitat, may benefit from a larger repertoire of genes that allows acclimation, response and adaptation to changing conditions and hence have much larger genomes.

thumbnailFigure 1. Genome size, isolation source and number of secondary metabolite gene clusters. The phylogenetic tree shown is calculated on concatenated ribosomal proteins and rooted with Bifidobacterium bifidum as an outgroup. The two bar plots are presented with species in the same order as the phylogenetic tree, representing genome size in thousands of genes on the left, colored by habitat, and number of secondary metabolite gene clusters on the right. Any combinations of cluster types found together count independently, e.g. an NRPS/PKS hybrid would be counted once as an NRPS and once as a PKS. The colors corresponding to habitat type and secondary metabolite class are shown in the key below the bar plots.

To provide context for the gene cluster comparisons, we constructed a phylogenetic tree using concatenated amino acid sequences from 41 ribosomal proteins shared by all strains (Figure 1). This is tree in good agreement with the phylogeny published by Gao and Gupta using 35 conserved genes from 98 actinobacterial genomes [26], although there are a couple of notable differences. In our tree Nakamurella multipartita DSM 44233T is found outside of the Pseudonocardiales, where it was within Pseudonocardiales based on their tree. Geodermatophilus obscurus G-20T was found to branch with Frankia, whereas their analysis suggested that it lay outside of the Frankiales. We also show that the groups they refer to as Micrococcales I and II group together, from Leifsonia xyli to Arthrobacter chlorophenolicus on our tree. Because it has already been shown that there can be extensive horizontal gene transfer within the actinomycetes [27,28], and that genome-based trees can differ from 16S and concatenated gene trees [29], we tested for recombination in the data set using the PHI test implemented in SplitsTree (p=1.0). A NeighborNet analysis was also not largely reticulate (Additional file 2), as one would expect for a data set impacted by homologous recombination. The secondary metabolite classes examined are also shown in Figure 1. While this is not an exhaustive list, it does cover all common secondary metabolites of actinomycetes. As might be expected, genome size and number of secondary metabolite biosynthetic gene clusters are positively correlated, as larger genomes can accommodate more gene clusters devoted to secondary metabolism (Figure 1 and Figure 2). This has also been noted in genomes of anaerobic microbes [30]. Interestingly, for genomes containing between 2000 and 6000 genes, pathogens tend to have a larger number of secondary metabolite biosynthetic gene clusters than free-living isolates from soil. This trend may not continue as more genomes from this order are sampled, however, as most of the pathogen genomes supporting this trend are from Mycobacterium. The same may be true with other patterns relating to isolation source.

Additional file 2. A NeighborNet analysis on concatenated ribosomal proteins.

Format: PDF Size: 23KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 2. Number of gene clusters per genome compared to genome size and habitat. This compares the total number of gene clusters with the total number of genes. Each point is colored according to Figure 1. Note that the number of gene clusters for soil-isolated strains with genomes between 4–6000 genes declines before the number of secondary metabolite gene clusters from pathogens.

To examine the overall similarity of the genomes between these organisms, we performed an all-vs-all BLAST search and grouped the results into sets of homologs using OrthoMCL. Two comparisons are shown in Figure 3. Both axes are ordered in the same way, based on the ribosomal protein tree. Each pairwise comparison is a tally of the homologs shared by two genomes. If multiple homologs were listed for each organism (e.g. T. whipplei has two copies of a gene and S. bingchinggensis has four) then the smaller number was counted for that single comparison. The total number of homologs for each pair of organisms was then divided by the total number of genes. This was done such that every vertical column is divided by the corresponding strain on the top, horizontal tree. For example, Tropheryma whipplei has only 783 protein coding genes due to reductive evolution as an intracellular pathogen. Therefore, T. whipplei shares nearly all its gene set with other strains (vertical column); while containing only a fraction of the genes present in other strains (horizontal row). In contrast, S. bingchenggensis has the largest number of protein-coding genes (10,022), so the many smaller genomes contain only a small fraction of the genes held by S. bingchenggensis, and this is reflected by a dark-colored vertical column.

thumbnailFigure 3. Whole genome similarity. The order for this comparison is the same as the phylogenetic tree in Figure 1, which has been shrunk and placed upon both axes for orientation. The heatmap legend is shown on the right. All comparisons between a genome and itself occur on a line stretching from the top left to the bottom right corners. The number of similar genes between two genomes is the numerator for each comparison and the genome represented by each column is used as the denominator for each comparison. The column divided by the size of T. whipplei, the smallest genome, is marked with ^ and the largest genome, S. bingchenggensis, is indicated with *.

Overall genome similarity clearly reflects the organismal phylogeny when distinguishing genera and large branches within a genus; however, the taxonomic level of genus is not uniformly applied. For example, Salinispora, Verrucosispora, and Micromonospora strains clearly show genomic similarities on the same degree as the other genera analyzed here and, thus, could be considered a single genus. The oldest of these genera, and therefore the one with precedence in naming, is Micromonospora[31]. Verrucosispora was described as a novel genus on the basis of a lack of arabinose in whole cell sugars, the presence of 10-methyl C17:0 fatty acids, and a 16S rRNA gene sequence not previously found in the family Micromonosporaceae[32]. The genus Salinispora was differentiated from other genera based largely on 16S rRNA gene diversity, a unique combination of fatty acid type and major menaquinones, and the requirement of sea water for growth [33]. It also appears that the genus Arthrobacter, which has long been divided into two groups, should be represented by two genera and Renibacterium should also remain separate. The case for Arthrobacter groups remaining in the same genus, however, was systematically considered and the two groups were determined to be members of the same genus with two “nuclei” [34]. A broader utilization of genomic data by the taxonomic community would assist in the creation of universal criteria for both species and genera definitions [35,36]. The genomes generated for research on natural products are very useful for improving actinobacterial systematics. Because taxonomy impacts both research focus and the interpretation of results, scientists with an interest in natural products should in turn not ignore the impact their data can have on taxonomy.

The whole genome comparisons also show a noticeable, but somewhat uneven, difference between rapid and slow-growing mycobacteria. It appears that the rate of genomic change leading to the branch containing Mycobacterium leprae and the M. tuberculosis strains has affected genomic content more than the change from rapid-growing nonpathogens to the slow-growing pathogens Mycobacterium sp. JDM601, Mycobacterium avium subsp. paratuberculosis K-10, Mycobacterium avium 104, Mycobacterium ulcerans Agy99 and Mycobacterium marinum M. In other words, the switch to pathogenicity itself did not require rapid genomic change because such rapid change is isolated to the M. leprae and M. tuberculosis branch of the tree. Unlike with Mycobacterium strains, the Corynebacterium isolates do not show such a large change between pathogens and nonpathogens. This is also reflected by what is known about the evolution of pathogenicity in Corynebacterium, as many pathogenicity factors appear to be acquired through recent horizontal gene transfer [37].

Gene cluster diversity

Given the diversity of lifestyles and habitats of actinomycetes it should be expected that discrete genera use secondary metabolites differently. For many of the genera examined, the most conserved secondary metabolite clusters are siderophores, whether they are NRPS products or NRPS-independent. 41 out of 102 genomes contain at least one gene cluster for NRPS-independent siderophore biosynthesis (aerobactin-like), but 31/34 in the Corynebacterium, Mycobacterium, Nocardia group do not have this class of siderophores. The Corynebacterium, Mycobacterium, Nocardia group (from Mycobacterium africanum to Segniliparus rotundus DSM 44985 in Figure 1), all contain the gene cluster for mycolic acid, with the exception of Corynebacterium kroppenstedtii (see Additional file 1, Conserved Clusters). In general, the genera with more pathogenic members, Corynebacterium and Mycobacterium, have higher proportions of conserved secondary metabolite gene clusters than the essentially saprophytic genera Streptomyces and Rhodococcus (Figure 4). This may be due to the increased homogeneity of environments inhabited by pathogens compared to free-living bacteria. This pattern based on host-association is broken with the Frankia, however, as Frankia species have almost no overlap in their secondary metabolic capabilities. All gene cluster families (GCFs) are shown in Additional file 3, and a stand-alone website is provided in Additional file 1 that contains all gene clusters found in the complete set of genomes. All conserved clusters mentioned are also present on the website provided under the “Conserved Clusters” link.

thumbnailFigure 4. Gene cluster conservation by class and genus. Histograms showing the conservation of lanthipeptide, PKS, NRPS and NRPS-independent siderophores are shown for each genus. For example, in Arthrobacter there are three PKS gene clusters that are unique unto themselves and one type of PKS gene cluster that shows up four times. To emphasize the abundance of some classes in certain genera, the number 10 is highlighted in red on the y-axis when present.

Additional file 3. List of genes grouped together within the genera of interest, gene range is separated by commas and gene groups are separated by semicolons.

Format: DOCX Size: 16KB Download fileOpen Data

One use for GCFs is the potential for cluster boundary delineation. Over evolutionary time natural product gene clusters will change their location on genomes and phylogenetic trees through horizontal gene transfer and genome rearrangements [6,7]. This mobility changes the surrounding genes, and if the GCF is found in enough genomic backgrounds, then the genes surrounding the cluster will change. The drop in gene content similarity is used to determine gene cluster boundaries shown in Figure 5. Knowing the genes involved in biosynthesis is essential for synthetic biologists and geneticists attempting to refactor pathways or to attempt heterologous expression of natural products in a new host.

thumbnailFigure 5. Conserved gene clusters in Streptomyces and Frankia. Four biosynthetic gene clusters were found to be well conserved in Frankia and three in Streptomyces. (A) The conservation of gene content was used to find putative gene cluster boundaries. The first and last conserved gene numbers are shown, along with the cluster name from our analysis. The standard six gene window extending from a pHMM hit are shown, with the exception of FRAAL_16 (Frankia alni ACN14a cluster 16) where additional genes are shown downstream of the original cluster due to continuing gene cluster content similarity. SCO_17 (“Streptomyces coelicolor” A3(2) cluster 17) occurs in a large, well conserved region, so no cluster boundaries were determined. (B) shows the sequence logo of all putative lanthipeptide precursors in Frankia clusters similar to FRAAL_9. (C) and (D) The occurrence of the well conserved gene clusters in Frankia and Streptomyces are shown on subtrees of the ribosomal protein phylogeny from Figure 1, along with the closest outgroup taxa. FRAAL_16 was found in all strains shown in (C) and is therefore not shown. Open circles show the start of branches lacking the relevant gene cluster. Open squares show the start of branches containing clusters with the same basic machinery but large scale changes compared to the other genomes. These branches are also shown with dashed lines. The number of lanthipeptide clusters similar to FRAAL_9 present in each Frankia genome is shown in red above the terminal branch on the tree.

Another use for GCFs is in correlating with molecular families through MS analyses. The basis for this work is that similar gene clusters should produce similar natural products [38]. The gene cluster families presented here can be correlated with the presence of such similar products, or molecular families, to uncover novel associations and find new natural products that would otherwise remain hidden in the analysis of a single sample.

Mycobacterium

Within Mycobacterium, many of the PKS gene clusters are well conserved in large phylogenetic groups, Figure 6, which are largely accounted for by differences in the complicated cell wall of the mycobacteria. For example, the gene cluster for the production of mycolic acid is shared by all strains, whereas the genes for production of phthiocerol are only present in slow-growing, pathogenic strains. In contrast, the NRPS clusters, with one exception, are either unique or shared with only a single close relative. The single exception is the gene cluster for mycobactin synthesis, a characterized siderophore, which is found in all strains except M. leprae. Two scotochromogenic strains, Mycobacterium gilvum and Mycobacterium sp. Spyr1 (which is proposed as synonymous with M. gilvum[39]) share a lycopene cyclase not found in the other strains that is possibly the source of their coloration (Mflv_0944-0956, Mspyr1_50120-50240).

thumbnailFigure 6. Distribution of abundant gene clusters in Mycobacterium. Six common secondary metabolite gene clusters with possibly informative phylogenetic distributions are shown from Mycobacterium. The tree is a subtree of that shown in Figure 1. Strains shown in blue contain the gene cluster indicated and strains shown in black do not. All strains above JDM601 are nonpathogenic or rarely opportunistic pathogens; branches leading to pathogenic strains are colored red. While none of these gene clusters have identical distribution patterns, the data suggest that some are important enough to be conserved for pathogenesis and others may be more useful in soil or aquatic habitats. Genome IDs used in this figure are: MAB, Mycobacterium abscessus ATCC 19977T; MSMEG, Mycobacterium smegmatis str. MC2 155; Mflv, Mycobacterium gilvum PYR-GCK; Mspyr1, Mycobacterium sp. Spyr1; Mjls, Mycobacterium sp. JLS; Mkms, Mycobacterium sp. KMS; Mmcs, Mycobacterium sp. MCS; JDM601, Mycobacterium sp. JDM601; MAP, Mycobacterium avium subsp. Paratuberculosis K-10; MAV, Mycobacterium avium 104; MUL, Mycobacterium ulcerans Agy99; MMAR, Mycobacterium marinum M; MLBr, Mycobacterium leprae Br4923; MCAN, Mycobacterium canettii CIPT 140010059; Rv, Mycobacterium tuberculosis H37RvT; Mb, Mycobacterium bovis AF2122/97; and MAF, Mycobacterium africanum GM041182.

Mycobacterium marinum is a very unique genome with regards to natural products compared to other Mycobacterium genomes. It has seven NRPS clusters, two PKS clusters, and three hybrid PKS-NRPS clusters not found in other mycobacterial genomes completed to date. This is especially surprising given the very close relationship between M. marinum and M. ulcerans, which have an average nucleotide identity of >98% [40]. Stinear et al. has shown that these clusters are not found on a single genomic island, and some of them may represent recent duplication events followed by divergence [41]. The evolution of natural product gene clusters in this group has already been mapped out in detail, including a new genome sequence for M. liflandii not included in the present study [42].

Corynebacterium

Corynebacterium is not known for its ability to produce natural products of the kind investigated here, and their genomes have not held many surprises in these regards. The most conserved cluster is that for mycolic acid as discussed above. Unlike most bacteria in the Corynebacterium-Nocardia-Mycobacterium group examined here, three pathogenic strains, Corynebacterium resistens DSM 45100T, Corynebacterium ulcerans BR-AD22 and Corynebacterium diphtheriae NCTC13129, share an aerobactin-like non-NRPS siderophore gene cluster. The ratio of isoprenoid and terpenoid biosynthesis gene clusters to PKS and NRPS clusters is high in corynebacteria compared to other genera, but this may be due simply to low overall numbers. The importance of these compounds at least to some of these strains is highlighted by the presence of the discrete mevalonate and non-mevalonate pathways for isoprene biosynthesis in Corynebacterium kroppenstedtii DSM 44385T and Corynebacterium variabile DSM 44702 [43]. Interestingly, the two mevalonate pathways seem to have reached Corynebacterium via different horizontal gene transfer routes, as they are only 54% similar to each other and more closely related to genes outside of the genus. The presence of two mevalonate pathways of different origins in Actinobacteria has been reported before, and these pathways are not unique to Corynebacterium among Actinobacteria[44].

Arthrobacter

The secondary metabolites in the Arthrobacter genomes examined here reveal little more than the divergence of Renibacterium salmoninarum ATCC 33209T from both Group I and II arthrobacteria. Overall, these strains have very few secondary metabolite gene clusters. One NRPS independent, aerobactin-like siderophore cluster is shared among all strains except Renibacterium, and a type III PKS is shared by all Group I strains. Arthrobacter arilaitensis RE117T and Arthrobacter aurescens TC1 also share a phytoene synthase gene cluster. The rest of the biosynthetic gene clusters present in this genus are unique to one strain.

Rhodococcus

The extent of secondary metabolite gene clusters revealed by Rhodococcus genome sequences was initially a surprise because no rhodococcal secondary metabolites were previously known [45]. In comparison with other actinomycete genomes, the Rhodococcus strains examined here have a skewed ratio of NRPS to PKS gene clusters. The average ratio of NRPS to PKS gene clusters for the entire data set is 0.45, but among rhodococcal genomes this ratio jumps to 2.8. In these four genomes there are only two PKS clusters that are found in only one strain, but each genome has at least four NRPS clusters that are not shared with any of the others. Despite the abundance of NRPS clusters, there are no conserved NRPS gene clusters; however, there are two conserved PKS clusters, one conserved phytoene synthase, which condenses two geranylgeranyl pyrophosphates to phytoene, one conserved lycopene cyclase, which cyclizes the ends of lycopene to the rings found in β-carotene, and a conserved butyrolactone biosynthetic gene cluster. The presence of a conserved butyrolactone biosynthetic gene cluster may indicate that a conserved cell-cell signaling pathway is important for the rhodococcal life cycle [46]. Rhodococcus strains are capable of differentiation and growth as either rods, cocci or hyphal filaments [47], but development has not been as well studied in this genus as in Streptomyces. The two strains from soil have larger genomes and more secondary metabolite biosynthetic gene clusters than Rhodococcus erythropolis PR4, a species isolated from a depth of 1,000 m in the Pacific Ocean south of Okinawa island, Japan, and Rhodococcus equi 103S, an equine pathogen.

Streptomyces

Based on solely genomic data, Streptomyces are the logical choice to mine for secondary metabolites. They have consistently high numbers of secondary metabolite biosynthetic gene clusters and a large variety of classes. Of course, streptomycetes have been the most heavily sampled historically, making rediscovery more likely when sampling from this genus. The eight genomes examined in this data set show a large diversity of gene clusters for secondary metabolism with little overlap between strains. The most common classes are PKS and NRPS, followed by terpenoids, aerobactin-like non-NRPS siderophores and lanthipeptides. All genomes contain the genes for butyrolactone biosynthesis, and in all but Streptomyces griseus at least one afsA, the central butyrolactone biosynthetic gene, homolog per genome is accompanied by a tetR family regulator immediately 5’ to afsA and in the opposite orientation (see Additional file 1, under Conserved Clusters). All eight genomes contain a non-NRPS aerobactin-like siderophore gene cluster similar to rhizobactin that is not currently tied to a product (SCO_17 in Figure 5). This gene cluster appears to be present in Catenulispora acidiphila as well, but significant changes to the gene cluster occurred between C. acidiphila and the most recent common ancestor of Streptomyces. All but Streptomyces sp. SirexAA-E contain the genes for the biosynthesis of the aerobactin-like siderophore desferrioxamine (nocardamine, SCO_10 in Figure 5). All streptomycetes, with the exception of S. griseus, contain the spore pigment type II PKS gene cluster. S. griseus contains a different spore pigment, produced instead by a type III PKS [48]. Interestingly, the lanthipeptide SapB, which was found to be required for aerial mycelia formation on rich media in “S. coelicolor” A3(2) and S. griseus[49], is only present in half of the strains.

Given the number of NRPS and PKS gene clusters in this genus, the amount of overlap with these clusters between genomes is very low. Unlike the abundance of NRPS clusters in Rhodococcus or PKS clusters in Frankia (discussed below), the ratio of NRPS to PKS clusters is also not heavily skewed in either direction and varies throughout the genus. While there has already been a significant amount of discovery of nonribosomal peptides and polyketides from Streptomyces, only a handful of terpenoids have been discovered from streptomycetes (see [20] for a review). Nevertheless, the number of terpene synthases present in these eight genomes comes close to those for PKS and NRPS biosynthesis, suggesting that a large diversity of terpenoids remain to be discovered in members of this genus.

Frankia

Frankia strains have a large number of secondary metabolite biosynthetic gene clusters, the vast majority of which are PKS clusters not shared with other strains. There are only four unique NRPS clusters within the genus, three of which occur only once and one that is shared by two strains. There are also two hybrid NRPS/PKS clusters, both unique. Out of the PKS clusters all but three sets of clusters are unique to one strain. Of the shared PKS clusters, one is a type II PKS shared by Frankia sp. CcI3 and Frankia sp. EuI1c, and one is a type II PKS conserved by all strains. The other cluster is a type I PKS that is conserved in all strains and duplicated in Frankia sp. EuI1c and Frankia alni ACN14a. There is only one type of lanthipeptide cluster found within the genus, but it is found either twice or four times in all genomes except FsymDg (Frankia symbiont of Datisca glomerata, Figure 5B). The sequence logo for the putative precursor peptides from these twelve lanthipeptide gene clusters show two conserved cysteine residues and a conserved threonine, along with a conserved LD motif that may be related to cleavage of the leader peptide. The conservation of cysteines, threonines and serines is biologically significant in lanthipeptides, as these residues are involved in lanthionine formation and cyclization that is central to lanthipeptide function (see [21] for review).

Other genera

The marine actinomycetes in the genus Salinispora have been a recent focus of natural products research because they have been historically understudied and because they possess large numbers of secondary metabolite gene clusters [50]. Moreover, they have the genetic capacity to produce a diverse array of natural product classes, Figure 1. Of the twelve classes examined in this study, Salinispora tropica and Salinispora arenicola have gene clusters that involve seven and nine classes, respectively. Thus, of the complete genomes examined here, S. arenicola has the highest diversity of secondary metabolite classes.

The genomes of Amycolatopsis mediterranei U32 and S699 (AMED and RAM, respectively), Actinosynnema mirum 101TT (Amir), Pseudonocardia dioxanivorans CB1190T (Psed) and Saccharopolyspora erythraea NRRL 2338T (SACE) also show a large number and diversity of secondary metabolite biosynthetic gene clusters. These strains were already known to produce rifamycin (AMED and RAM), nocardicin (Amir), and erythromycin (SACE). Amycolatopsis and Saccharopolyspora in particular are heavily researched, industrially important strains. Saccharomonospora viridis DSM 43017T, a pathogen that falls within the order Pseudonocardiales, has a smaller genome compared to its closest relatives in this analysis, a common theme among pathogens, and a corresponding large decrease in secondary metabolite biosynthetic gene clusters. The order Streptosporangiales also has significant potential for secondary metabolite production based on genome mining, although this is highly variable dependent on the genus examined.

Conclusions

We have concerned ourselves here with the study of natural product genetic diversity throughout the actinomycetes because the resultant patterns and observations add depth and breadth to our understanding of their molecular biology and ecology. The work presented in this manuscript is our first step towards a systematic framework for studying natural products, a difficult bioinformatic task especially for PKS and NRPS systems. We have found patterns showing that some genera have higher prevalence of NRPS or PKS natural products compared to other genera. We have used multiple types of comparisons to group every gene cluster in each genus well-represented by complete genomes. Such gene cluster families are essential for determining cluster boundaries and as part of integrated data sets for novel natural product discovery. These groupings found conservation of the spore pigment and desferrioxamine class of siderophores in Streptomyces, along with mycolic acid, mycobactin and phthiocerol in Mycobacterium. When applied to less well-studied genera, analysis of conservation within phylogenetic groups is a first-step tool to form hypotheses about pathways that may be of similar importance. Our focus on the genomes available from Frankia has allowed us to generate hypotheses about the importance of several natural product gene cluster families that may relate to core aspects of the evolution and biology of Frankia. We also show that some mycobacterial natural product gene clusters with uncharacterized products are preferentially conserved on one of the other side of the fast or slow growing split that divides the genus. All conserved clusters are shown together on a stand-alone website, as well as the complete collection of all gene clusters found in these genomes. Our broad overview of actinomycete genomic diversity also reinforces the view that several genera within the Actinobacteria may be in need of new descriptions that take genomic diversity into account. It is our hope that this work will provide valuable leads in the field about yet unforeseen aspects of actinomycete biology and ecology.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JRD designed and performed the research and wrote the manuscript. WWM guided the research and edited the manuscript. Both authors read and approved the final manuscript.

Acknowledgments

This work was supported in part by the National Institutes of Health (GM PO1 GM077596). JRD was funded by an Institute for Genomic Biology Postdoctoral Fellowship. We would like to thank Kou-San Ju and Joel Cioni for helpful discussions, and three anonymous reviewers for their comments.

References

  1. Ludwig W, Euzéby J, Schumann P, Busse H, Trujillo M, Kämpfer P, Whitman W: Road map of the phylum Actinobacteria. In Bergey’s manual® of systematic bacteriology. 5th edition. Edited by Whitman WB, Goodfellow M, Kämpfer P, Busse H-J, Trujillo ME, Ludwig W, Suzuki K-i, Parte A. New York City, New York, USA: Springer; 2012:1-28. PubMed Abstract | Publisher Full Text OpenURL

  2. Floyd K, Dias HM, Falzon D, Fitzpatrick C, Glaziou P, Hiatt T, Lienhardt C, Nguyen L, Sismanidis C, Timimi H, Uplekar M, van Gemert W, Zignol M, Raviglione M: Global tuberculosis report 2012. Geneva, Switzerland: World Health Organization; 2012. OpenURL

  3. Bérdy J: Bioactive microbial metabolites.

    J Antibiot (Tokyo) 2005, 58(1):1-26. Publisher Full Text OpenURL

  4. Ikeda M, et al.: Amino Acid Production Processes. In Microbial Production of l-Amino Acids. 79th edition. Edited by Faurie R, Thommel J, Bathe B, Debabov VG, Huebner S, Ikeda M, Kimura E, Marx A, Möckel B, Mueller U. New York City, New York, USA: Springer Berlin Heidelberg; 2003:1-35. OpenURL

  5. Bérdy J: Thoughts and facts about antibiotics: Where we are now and where we are heading.

    J Antibiot 2012, 65(8):385-395. PubMed Abstract | Publisher Full Text OpenURL

  6. Osbourn A: Secondary metabolic gene clusters: evolutionary toolkits for chemical innovation.

    Trends Genet 2010, 26(10):449-457. PubMed Abstract | Publisher Full Text OpenURL

  7. Fischbach MA, Walsh CT, Clardy J: The evolution of gene collectives: How natural selection drives chemical innovation.

    Proc Natl Acad Sci USA 2008, 105(12):4601-4608. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Baltz RH: Renaissance in antibacterial discovery from actinomycetes.

    Curr Opin Pharmacol 2008, 8(5):557-563. PubMed Abstract | Publisher Full Text OpenURL

  9. Fischbach MA, Walsh CT: Antibiotics for emerging pathogens.

    Science 2009, 325(5944):1089-1093. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Challis GL: Mining microbial genomes for new natural products and biosynthetic pathways.

    Microbiology 2008, 154(6):1555-1569. PubMed Abstract | Publisher Full Text OpenURL

  11. Bentley SD, Chater KF, Cerdeño-Tárraga AM, Challis GL, Thomson NR, James KD, Harris DE, Quail MA, Kieser H, Harper D, et al.: Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2).

    Nature 2002, 417(6885):141-147. PubMed Abstract | Publisher Full Text OpenURL

  12. Fischbach M, Walsh C: Assembly-Line Enzymology for Polyketide and Nonribosomal Peptide Antibiotics: Logic, Machinery, and Mechanisms.

    Chem Rev 2006, 106(8):3468-3564. PubMed Abstract | Publisher Full Text OpenURL

  13. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215(3):403-410. PubMed Abstract OpenURL

  14. Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes.

    Genome Res 2003, 13(9):2178-2189. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al.: Clustal W and clustal X version 2.0.

    Bioinformatics 2007, 23(21):2947-2948. PubMed Abstract | Publisher Full Text OpenURL

  16. Price MN, Dehal PS, Arkin AP: FastTree 2–approximately maximum-likelihood trees for large alignments.

    PLoS ONE 2010, 5(3):e9490. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Huson D, Bryant D: Application of phylogenetic networks in evolutionary studies.

    Mol Biol Evol 2006, 23(2):254-321. PubMed Abstract | Publisher Full Text OpenURL

  18. Eddy SR: Profile hidden Markov models.

    Bioinformatics 1998, 14(9):755-763. PubMed Abstract | Publisher Full Text OpenURL

  19. Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, Weber T, Takano E, Breitling R: antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences.

    Nucleic Acids Res 2011, 39(suppl 2):W339-W346. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Cane D, Ikeda H: Exploration and mining of the bacterial terpenome.

    Acc Chem Res 2012, 45(3):463-472. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Arnison PG, Bibb MJ, Bierbaum G, Bowers AA, Bugni TS, Bulaj G, Camarero JA, Campopiano DJ, Challis GL, Clardy J: Ribosomally synthesized and post-translationally modified peptide natural products: overview and recommendations for a universal nomenclature.

    Nat Prod Rep 2013, 30(1):108-160. PubMed Abstract | Publisher Full Text OpenURL

  22. Dunbar KL, Melby JO, Mitchell DA: YcaO domains use ATP to activate amide backbones during peptide cyclodehydrations.

    Nat Chem Biol 2012, 8(6):569-575. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Metcalf WW, van der Donk WA: Biosynthesis of Phosphonic and Phosphinic Acid Natural Products.

    Annu Rev Biochem 2009, 78:65. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Darling ACE, Mau B, Blattner FR, Perna NT: Mauve: Multiple alignment of conserved genomic sequence with rearrangements.

    Genome Res 2004, 14(7):1394-1403. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Moran NA: Microbial minimalism: genome reduction in bacterial pathogens.

    Cell 2002, 108(5):583-586. PubMed Abstract | Publisher Full Text OpenURL

  26. Gao B, Gupta RS: Phylogenetic framework and molecular signatures for the main clades of the phylum Actinobacteria.

    Microbiol Mol Biol Rev 2012, 76(1):66-112. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Doroghazi JR, Buckley DH: Widespread homologous recombination within and between Streptomyces species.

    ISME J 2010, 4(9):1136-1143. PubMed Abstract | Publisher Full Text OpenURL

  28. Smith SE, Showers-Corneli P, Dardenne CN, Harpending HH, Martin DP, Beiko RG: Comparative Genomic and Phylogenetic Approaches to Characterize the Role of Genetic Recombination in Mycobacterial Evolution.

    PLoS ONE 2012, 7(11):e50070. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Alam MT, Merlo ME, Takano E, Breitling R: Genome-based phylogenetic analysis of Streptomyces and its relatives.

    Mol Phylogenet Evol 2010, 54(3):763-772. PubMed Abstract | Publisher Full Text OpenURL

  30. Letzel A-C, Pidot SJ, Hertweck C: A genomic approach to the cryptic secondary metabolome of the anaerobic world.

    Nat Prod Rep 2013, 30(3):392-428. PubMed Abstract | Publisher Full Text OpenURL

  31. Skerman VBD, McGowan V, Sneath PHA: Approved Lists of Bacterial Names.

    Int J Syst Bacteriol 1980, 30(1):225-420. Publisher Full Text OpenURL

  32. Rheims H, SchumannC P, Rohde M, Stackebrandt E: Verrucosispora gifhornensis gen. nov., sp. nov., a new member of the actinobacterial family Micromonosporaceae.

    Int J Syst Bacteriol 1998, 48(4):1119-1127. PubMed Abstract | Publisher Full Text OpenURL

  33. Maldonado LA, Fenical W, Jensen PR, Kauffman CA, Mincer TJ, Ward AC, Bull AT, Goodfellow M: Salinispora arenicola gen. nov., sp. nov. and Salinispora tropica sp. nov., obligate marine actinomycetes belonging to the family Micromonosporaceae.

    Int J Syst Evol Microbiol 2005, 55(5):1759-1766. PubMed Abstract | Publisher Full Text OpenURL

  34. Stackebrandt E, Fowler VJ, Fiedler F, Seiler H: Taxonomic Studies on Arthrobacter nicotianae and Related Taxa: Description of Arthrobacter uratoxydans sp. nov. and Arthrobacter sulfureus sp. nov. and Reclassification of Brevibacterium protophormiae as Arthrobacter protophormiae comb. nov.

    Syst Appl Microbiol 1983, 4(4):470-486. PubMed Abstract | Publisher Full Text OpenURL

  35. Whitman WB: Intent of the nomenclatural Code and recommendations about naming new species based on genomic sequences.

    The Bulletin of BISMiS 2011, 2(2):135-139. OpenURL

  36. Sutcliffe IC, Trujillo ME, Goodfellow M: A call to arms for systematists: revitalising the purpose and practises underpinning the description of novel microbial taxa.

    Antonie van Leeuwenhoek 2012, 101(1):13-20. PubMed Abstract | Publisher Full Text OpenURL

  37. Cerdeño-Tárraga A, Efstratiou A, Dover L, Holden M, Pallen M, Bentley S, Besra G, Churcher C, James K, De Zoysa A, et al.: The complete genome sequence and analysis of Corynebacterium diphtheriae NCTC13129.

    Nucleic Acids Res 2003, 31(22):6516-6523. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Nguyen DD, Wu CH, Moree WJ, Lamsa A, Medema MH, Zhao X, Gavilan RG, Aparicio M, Atencio L, Jackson C, et al.: MS/MS networking guided analysis of molecule and gene cluster families.

    Proc Natl Acad Sci USA 2013.

    published ahead of print June 24, 2013

    OpenURL

  39. Kallimanis A, Karabika E, Mavromatis K, Lapidus A, LaButti KM, Liolios K, Ivanova N, Goodwin L, Woyke T, Velentzas AD, et al.: Complete genome sequence of Mycobacterium sp. strain (Spyr1) and reclassification to Mycobacterium gilvum Spyr1.

    Stand Genomic Sci 2011, 5(1):144-153. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Stinear TP, Seemann T, Pidot S, Frigui W, Reysset G, Garnier T, Meurice G, Simon D, Bouchier C, Ma L, et al.: Reductive evolution and niche adaptation inferred from the genome of Mycobacterium ulcerans, the causative agent of Buruli ulcer.

    Genome Res 2007, 17(2):192-200. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Stinear TP, Seemann T, Harrison PF, Jenkin GA, Davies JK, Johnson PDR, Abdellah Z, Arrowsmith C, Chillingworth T, Churcher C, et al.: Insights from the complete genome sequence of Mycobacterium marinum on the evolution of Mycobacterium tuberculosis.

    Genome Res 2008, 18(5):729-741. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Tobias NJ, Doig KD, Medema MH, Chen H, Haring V, Moore R, Seemann T, Stinear TP: Complete Genome Sequence of the Frog Pathogen Mycobacterium ulcerans Ecovar Liflandii.

    J Bacteriol 2013, 195(3):556-564. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Hamano Y, Dairi T, Yamamoto M, Kawasaki T, Kaneda K, Kuzuyama T, Itoh N, Seto H: Cloning of a gene cluster encoding enzymes responsible for the mevalonate pathway from a terpenoid-antibiotic-producing Streptomyces strain.

    Biosci Biotechnol Biochem 2001, 65(7):1627-1635. PubMed Abstract | Publisher Full Text OpenURL

  44. Lombard J, Moreira D: Origins and early evolution of the mevalonate pathway of isoprenoid biosynthesis in the three domains of life.

    Mol Biol Evol 2011, 28(1):87-99. PubMed Abstract | Publisher Full Text OpenURL

  45. McLeod MP, Warren RL, Hsiao WWL, Araki N, Myhre M, Fernandes C, Miyazawa D, Wong W, Lillquist AL, Wang D, et al.: The complete genome of Rhodococcus sp. RHA1 provides insights into a catabolic powerhouse.

    Proc Natl Acad Sci U S A 2006, 103(42):15582-15587. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Horinouchi S: A microbial hormone, A-factor, as a master switch for morphological differentiation and secondary metabolism in Streptomyces griseus.

    Front Biosci 2002, 7:2045-2057. Publisher Full Text OpenURL

  47. Jones AL, Goodfellow M: Genus IV. Rhodococcus (Zopf 1891) emend. Goodfellow, Alderson and Chun 1998a. In Bergey’s manual® of systematic bacteriology. 5th edition. Edited by Whitman WB, Goodfellow M, Kämpfer P, Busse H-J, Trujillo ME, Ludwig W, Suzuki K-i, Parte A. New York City, New York, USA: Springer; 2012:437-464. PubMed Abstract | Publisher Full Text OpenURL

  48. Ohnishi Y, Ishikawa J, Hara H, Suzuki H, Ikenoya M, Ikeda H, Yamashita A, Hattori M, Horinouchi S: Genome Sequence of the Streptomycin-Producing Microorganism Streptomyces griseus IFO 13350.

    J Bacteriol 2008, 190(11):4050-4060. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Kodani S, Hudson M, Durrant M, Buttner M, Nodwell J, Willey J: The SapB morphogen is a lantibiotic-like peptide derived from the product of the developmental gene ramS in Streptomyces coelicolor.

    Proc Natl Acad Sci USA 2004, 101(31):11448-11453. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Udwary DW, Zeigler L, Asolkar RN, Singan V, Lapidus A, Fenical W, Jensen PR, Moore BS: Genome sequencing reveals complex secondary metabolome in the marine actinomycete Salinispora tropica.

    Proc Natl Acad Sci USA 2007, 104(25):10376-10381. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL