- Research article
- Open Access
Genome based analysis of type-I polyketide synthase and nonribosomal peptide synthetase gene clusters in seven strains of five representative Nocardia species
BMC Genomics volume 15, Article number: 323 (2014)
Actinobacteria of the genus Nocardia usually live in soil or water and play saprophytic roles, but they also opportunistically infect the respiratory system, skin, and other organs of humans and animals. Primarily because of the clinical importance of the strains, some Nocardia genomes have been sequenced, and genome sequences have accumulated. Genome sizes of Nocardia strains are similar to those of Streptomyces strains, the producers of most antibiotics. In the present work, we compared secondary metabolite biosynthesis gene clusters of type-I polyketide synthase (PKS-I) and nonribosomal peptide synthetase (NRPS) among genomes of representative Nocardia species/strains based on domain organization and amino acid sequence homology.
Draft genome sequences of Nocardia asteroides NBRC 15531T, Nocardia otitidiscaviarum IFM 11049, Nocardia brasiliensis NBRC 14402T, and N. brasiliensis IFM 10847 were read and compared with published complete genome sequences of Nocardia farcinica IFM 10152, Nocardia cyriacigeorgica GUH-2, and N. brasiliensis HUJEG-1. Genome sizes are as follows: N. farcinica, 6.0 Mb; N. cyriacigeorgica, 6.2 Mb; N. asteroides, 7.0 Mb; N. otitidiscaviarum, 7.8 Mb; and N. brasiliensis, 8.9 - 9.4 Mb. Predicted numbers of PKS-I, NRPS, and PKS-I/NRPS hybrid clusters ranged between 4–11, 7–13, and 1–6, respectively, depending on strains, and tended to increase with increasing genome size. Domain and module structures of representative or unique clusters are discussed in the text.
We conclude the following: 1) genomes of Nocardia strains carry as many PKS-I and NRPS gene clusters as those of Streptomyces strains, 2) the number of PKS-I and NRPS gene clusters in Nocardia strains varies substantially depending on species, and N. brasiliensis strains carry the largest numbers of clusters among the species studied, 3) the seven Nocardia strains studied in the present work have seven common PKS-I and/or NRPS clusters, some of whose products are yet to be studied, and 4) different N. brasiliensis strains have some different gene clusters of PKS-I/NRPS, although the rest of the clusters are common within the N. brasiliensis strains. Genome sequencing suggested that Nocardia strains are highly promising resources in the search of novel secondary metabolites.
Actinomycetous strains of the genus Nocardia usually live in soil or water and play saprophytic roles in the environment, but also are opportunistic human pathogens, infecting the respiratory tract, skin, brain, and other organs of both immunocompromised and immunocompetent patients. To date, more than 80 species have been established in the genus Nocardia, and approximately one-third to one-half of the species have been reported as human pathogens [1–3]. Because of their medical importance, Nocardia strains have accumulated in microbial collections as a resource for clinical and scientific studies in the last few decades (e.g., [4–7]).
Although Nocardia strains belong to the Order Actinomycetales together with Streptomyces strains, the latter being known as a rich resource for discovery of secondary metabolites, few studies have been focused on secondary metabolites and their synthetic genes in Nocardia strains.
Type I polyketide synthase (PKS-I) and nonribosomal peptide synthetase (NRPS) gene clusters are two of the major secondary metabolite-producing clusters in bacteria and are involved in the biosynthesis of polyketide chains and nonribosomal peptides, respectively. It has been found that these clusters produce several medically and industrially important compounds, such as pathogenic factors, avermectin, erythromycin, and vancomycin.
In the present paper, we searched for PKS-I and NRPS genes in the genomes of representative Nocardia strains and analyzed their sequence similarities and differences in domain/module structures. While we were sequencing and analyzing Nocardia draft genomes, two new Nocardia genomes of N. cyriacigeorgica GUH-2 [8, 9] and N. brasiliensis HUJEG-1  were published. We included them in the present analysis together with N. farcinica genome, which our group has published previously .
N. otitidiscaviarum IFM 11049 and N. brasiliensis IFM 10847 were from the IFM culture collections of MMRC, Chiba University, Japan . N. asteroides NBRC 15531T and N. brasiliensis NBRC 14402T were from the NBRC culture collection . Cells were cultured in brain heart infusion liquid culture medium (Difco) in the conventional manner.
Acquisition of whole-genome sequences
Genomic DNA of N. otitidiscaviarum IFM 11049, N. brasiliensis (IFM 10847, NBRC 14402T), and N. asteroides NBRC 15531T was prepared as described previously . Genome sequences were read by the pyrosequencing method using genome sequencer GS FLX Instruments and GS FLX Titanium Kits (Roche Applied Science, Japan). The read redundancy for the four draft genomes ranged between 55 and 104. We assembled the sequence reads of N. otitidiscaviarum IFM 11049, N. brasiliensis IFM 10847, N. brasiliensis NBRC 14402T, and N. asteroides NBRC 15531T, and obtained 65, 223, 115, and 39 contigs, which were longer than 500 bp. The estimated genome sizes of N. otitidiscaviarum IFM 11049, N. brasiliensis IFM 10847, N. brasiliensis NBRC 14402T, and N. asteroides NBRC 15531T were 7.9 Mb, 9.2 Mb, 8.9 Mb, and 7.0 Mb, respectively. The draft genome sequences of N. otitidiscaviarum IFM 11049, N. brasiliensis (IFM 10847, NBRC 14402T), and N. asteroides NBRC 15531T are available at GenBank/EMBL/DDBJ under the accession numbers BATZ01000001–BATZ01000065, BAUA01000001–BAUA01000223, BAFT01000001–BAFT01000128, and BAFO01000001–BAFO01000049, respectively. The complete genome sequences of N. cyriacigeorgica GUH-2, N. brasiliensis HUJEG-1 (=ATCC 700358), and N. farcinica IFM 10152 were downloaded from DDBJ , with accession numbers FO082843, CP0033876, and AP006618, respectively.
Analysis of PKS-I and NRPS gene clusters
The assembled contig sequences were submitted to the auto-annotation pipeline MiGAP [15, 16] at DDBJ as described previously . Assigned ORFs were further searched for signature domains of PKS-I and NRPS genes using the InterPro domain database [18, 19]. ORFs having ketosynthase (KS) domain (IPR014030, IPR014031, IPR020841) or condensation (C) domain (IPR001242) were identified, and their adjacent genes were further analyzed as PKS-I and NRPS gene candidates. Module organizations were determined manually based on search results using InterPro database, results using PKS/NRPS analysis website , and signature sequences deduced using MOTIF search . We also used antiSMASH [22, 23], a website for antibiotics and secondary metabolite analysis, for finding orthologous clusters and predicting substrates for adenylation domains. PKS-I and NRPS gene clusters of N. farcinica IFM 10152, N. cyriacigeorgica GUH-2, and N. brasiliensis HUJEG-1 were also identified using the N. farcinica genomic database [11, 24]. We assumed that two or more PKS-I and/or NRPS genes that were adjacent to each other constitute one cluster for secondary metabolite production (See Additional file 1: Table S1, for details and exceptions). We also assumed that one multi-domain PKS-I or NRPS gene that was not accompanied by adjacent PKS-I/NRPS genes constitute one independent cluster. However, genes having only a single PKS-I or NRPS domain were excluded from the present analysis because we considered them atypical, and focused on multi-domain clusters. The contig sequences containing PKS-I and NRPS gene clusters are available at GenBank/EMBL/DDBJ under the following accession numbers: [AB700569 - AB700587] (N. otitidiscaviarum IFM 11049), [AB701575 - AB701605] (N. brasiliensis IFM 10847), [AB701607 – AB701636] (N. brasiliensis NBRC 14402T), and [AB685274], [AB700124 - AB700133], [AB700557 - AB700568] (N. asteroides NBRC 15531T).
Search for orthologous gene clusters among species and strains
BLASTP search was performed using the NCBI Protein BLAST program against the non-redundant protein sequence database [25, 26]. We considered Nocardia genes homologous to other genes when they have more than 70% sequence similarity in BLASTP search, and also when their domain organizations have high similarity. We also compared clusters with domain organizations that only partially match each other, as described in the text.
Results and discussion
The two leftmost columns in Table 1 list Nocardia strains studied in the present paper and their exact (complete genome) or estimated (draft) genome sizes. The genome sizes ranged between 6.0 and 9.4 Mb, similar to those of representative Streptomyces strains (5.0 - 11.9 Mb), the most abundant sources of secondary metabolites [27–29]. The fourth column indicates that the strains are from different clinical origins.
Figure 1 illustrates phylogenetic positions of the five Nocardia species (seven strains) studied in the present paper among 78 other established Nocardia strains. It also includes Streptomyces coelicolor and Mycobacterium tuberculosis for comparison. Four out of the five species are located in different clades of the 16S rRNA phylogenetic tree, indicating that the present analysis is based on information from a wide range of Nocardia species. N. asteroides and N. cyriacigeorgica are in the same clade. We also included three strains from N. brasiliensis to elucidate intra-species variations.
PKS-I, NRPS, and PKS-I/NRPS hybrid gene clusters from the Nocardia strains were predicted as described in Methods. Numbers of the three different types of clusters and the total number of clusters in each strain are listed in the four rightmost columns in Table 1. Among the seven strains, the numbers of PKS-I, NRPS, PKS-I/NRPS hybrid clusters, and their total number increased proportionally to the genome size, except for N. otitidiscaviarum and N. asteroides (Table 1) as reported in other genera . N. farcinica had the least while N. brasiliensis HUJEG-1 had the highest number of the gene clusters. The total number of clusters within the three N. brasiliensis genomes differed (27 to 30), suggesting that different strains of the same species potentially produce their own unique products (see below).
We also counted the numbers of type-II PKS, type-III PKS and terpene synthesis clusters in each genome. The numbers ranged between 0 and 3, except in N. otitidiscaviarum and N. brasiliensis strains, which have five and eight clusters for terpene synthesis per genome, respectively. In the present paper, however, we focused on PKS-I and NRPS secondary metabolite clusters because their products usually have larger molecular weights with more complex chemical structures than the others and have unique pharmacological activities.
Figure 2 shows all the clusters found in each genome. Presumptive orthologous clusters, as defined in Methods, are aligned in the same row of the table. The rightmost column shows secondary metabolites referred from the database, e.g., , and also those inferred using the tools described in Methods.
Clusters common among the seven strains
Figure 2 suggests that seven presumable products (lines #1, #2, #4, #5, #25, #27, and #35) are common among the seven strains belonging to the five species. It is noteworthy that clusters #1, #2, #4, and #5 reside close to the original points of replication (ori.) in the three species whose completed genome sequences are known, in accordance with a report showing that conserved genes reside in the internal core region of actinomycete genomes .
Mycolic acid (pks13)
We predicted that the products of the PKS-I genes in line #1 were mycolic acids, cell wall components in members belonging to Corynebacterineae, because these PKS-Is showed the same domain organization as those of pks13 in Mycobacterium tuberculosis for the synthesis of mycolic acids , and also showed over 80% sequence similarities to PKS-Is of N. farcinica annotated for mycolic acid synthesis .
NRPS genes of line #2 were predicted to be for poly-lysine synthesis because their module organizations are the same as that of poly-lysine synthetase (Pls) in Streptomyces albulus. The corresponding gene, nfa3790, is also annotated as a Pls homolog in the N. farcinica database . The sequence identity between Pls in S. albulus and nfa3790 in N. farcinica was 55%.
Ser/The-rich nonribosomal peptides
NRPS gene clusters in line #4 were present in all strains examined, but only a partial sequence was found in N. otitidiscaviarum (Additional file 2: Figure S1A). Zoropogui et al.,  has suggested that #4 in N. cyriacigeorgica was for synthesis of 2-amino-9,10 epoxi-8-oxodecanoic acid, a component of HC-toxin, which is produced by a plant pathogen causing corn leaf spots (reviewed in ). However, we suggest two other possibilities based on the domain organization of the cluster. One is that the intact #4 cluster is for synthesis of a serine/threonine-rich peptide composed of 11 amino acids. The second possibility is that the same sequence consists of two different NRPS clusters, since two thioesterase domains are present within the sequence, and accordingly, produces two peptide chains. In the latter case, the products of cluster #4 in N. asteroides would be two peptides: one composed of four amino acids (NCAST_11_00880 & NCAST_11_00870) and the other composed of seven amino acids (NCAST_11_00860 & NCAST_11_00850). Likewise, in N. brasiliensis, the products of cluster #4 would also be two peptides: one composed of six amino acids (O3I_004080 & O3I_004085) and the other composed of five amino acids (O3I_004090) (Figure 2, #4, rightmost column).
A possible reason why N. otitidiscaviarum, unlike the other species, has only a partial #4 cluster is that this species is phylogenetically distant from the other strains, as shown in Figure 1. The relationships among the products of these gene clusters, the phylogenetic positions of the strains, and their pathogenicity to plants and animals are interesting issues to clarify.
NRPS genes in line #35 were also present in all the strains examined, although those of N. otitidiscaviarum and N. brasiliensis IFM 10847 were partial compared with those of the other strains (Additional file 2: Figure S1B). The genes in five strains, except for N. otitidiscaviarum and N. brasiliensis IFM 10847, have 12 modules, and many of their adenylation domains were predicted to select Ser as the substrate. Hence, we assumed the products would be Ser-rich 12 aa peptides. Besides the clusters common among the seven strains in line #35, similar NRPS clusters were also present in the adjacent clusters of line #36, but only in four species excluding N. brasiliensis strains. Interestingly, the predicted products are Ser-rich 13 aa peptides.
Although the conservation of Ser- (and Thr-) rich large peptides synthesized by clusters of #4, #35, and #36 in the Nocardia strains suggests that they have important roles, physiological roles of the products remain to be investigated.
Nocobactin (nbt)-like siderophore
The PKS-I/NRPS hybrid cluster #5 in N. farcinica produces nocobactin, a siderophore and a pathogenic factor, as proven by Hoshino et al., . Figure 3 compares clusters in line #5, which are candidates of nocobactin-like siderophore-producing genes. N. asteroides has a full set of genes required for siderophore synthesis (NCAST_11_00460 through NCAST_11_00420 and NCAST_13_1740 in #5). NCAST_13_1740 contains an nbtF-like gene , but is found separated from the rest of nbt genes in contig 11 by more than 180 Mb based on the contig sequences. In N. brasiliensis, the corresponding cluster structure is different from that in N. farcinica in terms of module number and domain organization, suggesting that a functionally similar but structurally different molecule is synthesized in N. brasiliensis. In N. otitidiscaviarum, cluster #5 lacks genes corresponding to nbtA-C and nbtE, which are required for nocobactin synthesis in N. farcinica, suggesting that cluster #5 in N. otitidiscaviarum may not be able to produce a nocobactin-like siderophore. Interestingly, however, nbtD- and nbtF-like genes, which both have 75% similarities to nbtD and nbtF genes in N. farcinica, respectively, are found in the middle of two different contigs in the N. otitidiscaviarum genome (contigs 3 and 41, respectively; listed in cluster #5 in Figure 2; see also Figure 3), suggesting gene loss/gain and recombination during evolution within the genus Nocardia. Such gene loss/gain is not only observed in the nbt-like gene cluster #5, but also found in nfa7170-7200 homologs (#4), nfa50330 homologs (cluster #35), and nfa50630-50620 homologs (#36), as shown in Additional file 2 (Figure S1). It is further noteworthy that N. otitidiscaviarum has another candidate gene cluster for siderophore synthesis, which is shown in line #40 of Figure 2. The domain and module organization of this cluster are shown in Figure 3F, which shows the differences from the one for nocobactin synthesis.
Other common gene clusters
Two NRPS genes in lines #25 and #27 (Figure 2) were common in all the strains sequenced, having four and three modules, respectively. However, the amino-acid composition of the products could not be predicted in silico using antiSMASH. Chemical structures and physiological roles of the products remain to be investigated.
Clusters missing in a few strains but found in others
PKS-I clusters of line #13 (Figure 2) were present in the six strains except N. otitidiscaviarum. Other PKS-I clusters in line #30 were present in the four strains but not in N. asteroides and in the two N. brasiliensis strains (NBRC 14402T, IFM 10847). Nfa30250 (Figure 2 #13) has been predicted to be involved in mycocerosic acid synthesis in the N. farcinica genome project . On the other hand, O3I_032485 in N. brasiliensis HUJEG-1 (Figure 2 #30) has been putatively annotated as mycocerosate synthase . Both PKS-Is in N. farcinica (#13) and in N. brasiliensis (#30) showed approximately 47% amino acid similarities to Mycobacterium tuberculosis PKS-I [GenBank/EMBL/DDBJ accession number: CCP46654] for the synthesis of mycocerosic acid, a pathogenic factor [40, 41]. All the PKS-Is listed in line #13 and #30 showed sequence similarities ranging between 62 and 75%, having almost the same protein length (approximately 2200 aa) and identical domain organization (KS/AT/DH/KR/ER/ACP) (Additional file 1: Table S1).
It is possible that in N. otitidiscaviarum, cluster #30 is a substitute for cluster #13. Interestingly, N. farcinica and N. cyriacigeorgica possess both #13 and #30, which could be related to the strong pathogenicity of the two strains.
Polyunsaturated fatty acid (pfaA)
Orthologous genes of cluster #11 are found in N. asteroides, N. otitidiscaviarum, and N. brasiliensis IFM 10487 and HUJEG-1, but not in N. farcinica, N. cyriacigeorgica, and N. brasiliensis NBRC 14402 (Figure 2). The PKS-I sequence identities among N. asteroides, N. otitidiscaviarum, and N. brasiliensis ranged from 66 to 78%. The modular structures of PKS-Is in line #11 (KS/AT/ACP/ACP/KR) are unusual because they have two tandem ACP domains, and one KR domain is located after ACP (Additional file 3: Figure S2). These features have been known in polyunsaturated fatty acid (PUFA) synthase, PfaA, of marine bacteria, such as those belonging to the genus Shewanella[42, 43]. The module organizations are similar between #11 and PfaA (Figure S2), and their amino-acid sequence similarities are over 50%. Hence, we predict that the products of PKS-I in line #11 are polyunsaturated fatty acids, as already reported in N. brasiliensis HUJEG-1 genome . However, reports of production of polyunsaturated fatty acids have been limited only in some psychrophilic, piezophilic, or halophilic bacteria in prokaryotes [44, 45] and have not been found in Nocardia strains. Thus, future chemical and synthetic analysis are required to explore the potential of Nocardia strains as industrial producers of these pharmaceutically and nutraceutically valuable compounds.
Only a few species-specific clusters are found in N. farcinica (cluster #29), N. cyriacigeorgica (#7, #8, #24, #34), and N. otitidiscaviarum (#40, #41, #42), suggesting that the evolution of their polyketides and nonribosomal metabolites are not as dynamic as in other species (Figure 2). On the other hand, N. asteroides has nine species-specific clusters (PKS-I, #43 - #46; NRPS, #28, #47 - #49; PKS-I/NRPS hybrid, #50). Among them, we selected #46 as a representative unique cluster structure in N. asteroides (Figure 4A). The cluster consists of four adjacent genes, namely, NCAST33_02270, _02280, _02290, and _02300, each of which has the top BLAST-hit to a gene sequence in different Streptomyces strains (see Additional file 1: Table S1, for more details). The cluster consisted of twelve modules and was larger than any other cluster in the Nocardia strains studied in this paper. The orders of domains within the modules follow the accepted theory for PKS-I gene cluster structures, i.e., having repeats of “KS-AT-optional domains-ACP” [46, 47]. The same module organization, however, could not be found in available public databases, and the cluster has an amino-acid sequence similarity of less than 53% to any other known PKS-I clusters. This suggests that the gene cluster may be involved in production of a novel secondary metabolite. We propose a chemical structure of the polyketide chain synthesized by this cluster as shown in Figure 4B based on the PKS-I assembly line rule  as follows. The presence of a ketosynthaseQ (KSQ) domain at the N terminus of NCAST_33_02270 and a thioesterase (Te) domain at the C terminus of NCAST_33_02300 indicates that NCAST_33_02270 and NCAST_33_02300 contain the modules that initiate and terminate PKS-I assembly line, respectively. Among the eleven AT domains of module 1 (m1) to 11 (m11), ten AT domains, except that in m8, had the HAFHS signature amino-acid sequence, which is specific for malonyl-CoA in substrate recognition [48, 49]. The substrate of the AT domain in m8, which has IASHS amino-acid sequence, and the starter molecules loaded on LM could not be predicted using bioinformatic approach. Hence, this polyketide backbone is predicted to be Cx-C2-C2-C2-C2-C2-C2-C2-Cy-C2-C2-C2, where C2 is a unit derived from malonyl-CoA, and Cx and Cy are carbon backbones derived from presently unknown substrates. An inactive DH (dh)-KR pair, which is responsible for formation of a hydroxyl group, is present as an optional domain in m1, m3, and m4 modules, while the optional domain present in m2 is a DH-ER-KR trio, which completely reduces the ketone residue formed by the m2 module. Optional domains of m5 through m11 are DH-KR pairs, which form double bonds from the ketone residues produced by these modules. The resulting molecule has a polyketide backbone consisting of more than 23 carbons.
The N. brasiliensis strains had the largest number of species-specific clusters among the strains studied in the present work: five to six PKS-I clusters (#15, #16, #19, #22, #26, #51, #52), six NRPS clusters (#3, #6, #17, #18, #21, #32), and two to four PKS-I/NRPS hybrid clusters (#9, #10, #12, #23), depending on the strains. Among them, three PKS-I clusters (#22, #51, #52) consisted of a single module, suggesting that their final or intermediate products may be small in accordance with the assembly line rule [46, 47], unless the modules are used iteratively, as has been reported in actinomycetes (e.g., [50, 51]). Interestingly, PKS-Is in #22 shows 61% sequence similarity to PksE of Streptomyces griseus [GenBank/EMBL/DDBJ accession number: AAO25858], whose product is an unusual polyketide compound including a 9-membered enediyne core .
The rest of the clusters with multiple modules may possibly produce large species-specific products. In particular, nonribosomal peptides produced by clusters #18 and #21 are, respectively, predicted to consist of nine and six (to eight) amino acids.
Intra-species variations of clusters
Strain-specific as well as species-specific PKS-I clusters were found in N. brasiliensis strains (#19, #30, #52), indicating different strains of the same species may potentially produce different products. Because there are large numbers of Nocardia strains stored in several bio-resource centers (e.g., [5–7]), and are rapidly accumulating , these Nocardia strains constitute a highly promising future resource for exploring secondary metabolites.
Other unique examples of clusters
Figure 5A shows a module structure of PKS-I/NRPS hybrid cluster #20 (Figure 2) in N. otitidiscaviarum. The left half (approximately 8,300 aa) of the cluster has high similarity to N. brasiliensis hybrid cluster #20 (75 - 82% similarities), but the right half (7,746 aa of NOTIT_47_07220 plus 9,157 aa of NOTIT_47_07270) has top hits to two proteins in Rhodococcus opacus (7,746 aa) [GenBank/EMBL/DDBJ accession number: YP_002777453] (Figure 5D) and Gordonia aichiensis (9,517 aa) [GenBank/EMBL/DDBJ accession number: WP_005170336], with 59% (4,667/7,781) and 54% (4,543/8,322) amino acid similarities, respectively. The right-half consisting of two NRPSs is widely conserved in Rhodococcus spp. including R. jostii as shown in Figure 5E.
Because the gene cluster of N. otitidiscaviarum #20 contains one PKS-I module and 19 NRPS modules, the product has been tentatively predicted to include one polyketide chain and 19 amino acids (Additional file 4: Figure S4). It should also be mentioned, however, that the two genes in cluster #20, NOTIT_47_07270 and NOTIT_47_07220, each contain thioesterase domains (Te) at their C-terminal ends, suggesting another possibility that the product may contain 12 amino acids instead of 19.
All cluster structures we analyzed in the present work are listed in (Additional file 1: Table S1).
We conclude the following: 1) genomes of Nocardia strains carry as many PKS-I and NRPS clusters as Streptomyces strains, 2) the number of PKS-I and NRPS gene clusters in Nocardia strains varies substantially depending on species, and N. brasiliensis strains carry the largest number of clusters among the species studied, 3) the seven Nocardia strains studied in the present work have six common PKS-I/NRPS clusters, some of whose products are yet to be studied, and 4) different N. brasiliensis strains have a few different clusters for secondary metabolite synthesis. Also, the following are suggested: 1) there is no clear relation between genome size and pathogenicity in Nocardia strains, e.g. N. farcinica and N. brasiliensis are both prevalent pathogens, but their genome sizes are 6.0 Mb and 9.4 Mb, respectively, the minimum and maximum among the strains studied, and 2) some genes (e.g. cluster #17 in Figure 2) are likely to have been horizontally transferred from (or to) other actinomycetous strains, such as Rhodococcus spp.
To summarize, in this study, we compared complete and draft genome sequences of seven strains from five representative Nocardia species. The sequences we obtained provided useful information for inferring numbers and molecular structures of secondary metabolites potentially produced by Nocardia strains. Genome sequencing revealed the possibility that Nocardia strains are as attractive resources as Streptomyces strains, the largest resource of natural compounds, in the search for new useful secondary metabolites.
HK: Senior Chief, Biological Resource Center, National Institute of Technology and Evaluation (NBRC).
NI: Chief, NBRC.
AH: Senior Chief, NBRC.
AT-N: Postdoctoral Fellow, Medical Mycology Research Center (MMRC), Chiba University.
TM: Research Associate, MMRC, Chiba University.
KS: Senior Director, NBRC.
NF: Senior Director, NBRC.
TG: Professor, MMRC, Chiba University.
Type-I polyketide synthase
Nonribosomal peptide synthetase
Biological Resource Center, National Institute of Technology and Evaluation
Medical Mycology Research Center
Acyl carrier protein
Beaman BL, Beaman L: Nocardia species: host-parasite relationships. Clin Microbiol Rev. 1994, 7 (2): 213-264.
Kageyama A, Yazawa K, Ishikawa J, Hotta K, Nishimura K, Mikami Y: Nocardial infections in Japan from 1992 to 2001, including the first report of infection by Nocardia transvalensis. Eur J Epidemiol. 2004, 19 (4): 383-389.
Brown-Elliott BA, Brown JM, Conville PS, Wallace RJ: Clinical and Laboratory Features of the Nocardia spp. Based on Current Molecular Taxonomy. Clin Microbiol Rev. 2006, 19 (2): 259-282. 10.1128/CMR.19.2.259-282.2006.
Medical Mycology Research Center: Japan: Chiba University, http://www.pf.chiba-u.ac.jp/eng/index.html,
NBRC culture collection: http://www.nbrc.nite.go.jp/e/index.html,
American Type Culture Collection, http://www.atcc.org/,
Leibniz Institute DSMZ-German Collection of Microorganisms and cell culture: https://www.dsmz.de/,
Zoropogui A, Pujic P, Normand P, Barbe V, Beaman B, Beaman L, Boiron P, Colinon C, Deredjian A, Graindorge A, Mangenot S, Nazaret S, Neto M, Petit S, Roche D, Vallenet D, Rodriguez-Nava V, Richard Y, Cournoyer B, Blaha D: Genome sequence of the human- and animal-pathogenic strain Nocardia cyriacigeorgica GUH-2. J Bacteriol. 2012, 194 (8): 2098-2099. 10.1128/JB.00161-12.
Zoropogui A, Pujic P, Normand P, Barbe V, Belli P, Graindorge A, Roche D, Vallenet D, Mangenot S, Boiron P, Rodriguez-Nava V, Ribun S, Richard Y, Cournoyer B, Blaha D: The Nocardia cyriacigeorgica GUH-2 genome shows ongoing adaptation of an environmental Actinobacteria to a pathogen’s lifestyle. BMC Genomics. 2013, 14 (1): 286-10.1186/1471-2164-14-286.
Vera-Cabrera L, Ortiz-Lopez R, Elizondo-Gonzalez R, Perez-Maya AA, Ocampo-Candiani J: Complete genome sequence of Nocardia brasiliensis HUJEG-1. J Bacteriol. 2012, 194 (10): 2761-2762. 10.1128/JB.00210-12.
Ishikawa J, Yamashita A, Mikami Y, Hoshino Y, Kurita H, Hotta K, Shiba T, Hattori M: The complete genomic sequence of Nocardia farcinica IFM 10152. Proc Natl Acad Sci U S A. 2004, 101 (41): 14925-14930. 10.1073/pnas.0406410101.
IFM culture collections of MMRC: Japan: Chiba University, http://www.pf.chiba-u.ac.jp/eng/bioresoures/index.html,
Ichikawa N, Oguchi A, Ikeda H, Ishikawa J, Kitani S, Watanabe Y, Nakamura S, Katano Y, Kishi E, Sasagawa M, Ankai A, Fukui S, Hashimoto Y, Kamata S, Otoguro M, Tanikawa S, Nihira T, Horinouchi S, Ohnishi Y, Hayakawa M, Kuzuyama T, Arisawa A, Nomoto F, Miura H, Takahashi Y, Fujita N: Genome sequence of Kitasatospora setae NBRC 14216 T: an evolutionary snapshot of the family Streptomycetaceae. DNA Res. 2010, 17 (6): 393-406. 10.1093/dnares/dsq026.
DNA Data Bank of Japan: http://www.ddbj.nig.ac.jp/index-e.html,
Microbial Genome Annotation Pipeline: http://www.migap.org/index.php/en/aboutpipeline,
Sugawara H, Ohyama A, Mori H, Kurokawa K: Microbial Genome Annotation Pipeline (MiGAP) for diverse users [abstract]. The 20th International Conference on Genome Informatics (GIW2009). 2009, 20 (S001): 1-2.
Komaki H, Ichikawa N, Oguchi A, Hanamaki T, Fujita N: Genome-wide survey of polyketide synthase and nonribosomal peptide synthetase gene clusters in Streptomyces turgidiscabies NBRC 16081. J Gen Appl Microbiol. 2012, 58 (5): 363-372. 10.2323/jgam.58.363.
InterPro domain database: http://www.ebi.ac.uk/interpro/,
Mulder NJ, Kersey P, Pruess M, Apweiler R: In silico characterization of proteins: UniProt, InterPro and Integr8. Mol Biotechnol. 2008, 38 (2): 165-177. 10.1007/s12033-007-9003-x.
PKS/NRPS analysis Web-site: http://nrps.igs.umaryland.edu/nrps,
KEGG MOTIF Search: http://www.genome.jp/tools/motif/,
Antibiotics and Secondary Metabolite Analysis Shell (antiSMASH): http://www.secondarymetabolites.org/,
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 (Web Server issue): W339-W346.
Nocardia farcinica genome project page. http://nocardia.nih.go.jp/,
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1016/S0022-2836(05)80360-2.
NCBI/BLAST Home: http://blast.ncbi.nlm.nih.gov/,
Nett M, Ikeda H, Moore BS: Genomic basis for natural product biosynthetic diversity in the actinomycetes. Nat Prod Rep. 2009, 26 (11): 1362-1384. 10.1039/b817069j.
Berdy J: Bioactive microbial metabolites. J Antibiot (Tokyo). 2005, 58 (1): 1-26. 10.1038/ja.2005.1.
Berdy J: Thoughts and facts about antibiotics: where we are now and where we are heading. J Antibiot (Tokyo). 2012, 65 (8): 385-395. 10.1038/ja.2012.27.
Eppinger H: Über eine neue pathogene Cladothrixund eine durch sie hervorgerufene Pseudotuberculosis (Cladothrichica). Beitrage zur pathologischen Anatomie. 1891, 9: 287-328.
Snijders EP: Cavia-scheefkopperij, een nocardiose. Geneeskundig Tijdschrift voor Nederlandsch-Indie. 1924, 64: 85-87.
Lindenberg A: Un nouveau mycétome. Archives deParasitologie. 1909, 13: 265-282.
Gordon RE, Mihm JM: A Comparison of Nocardia asteroides and Nocardia brasiliensis. J Gen Microbial. 1959, 20: 129-135. 10.1099/00221287-20-1-129.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Jenke-Kodama H, Sandmann A, Müller R, Dittmann E: Evolutionary implications of bacterial polyketide synthases. Mol Biol Evol. 2005, 22 (10): 2027-2039. 10.1093/molbev/msi193.
Portevin D, De Sousa-D'Auria C, Houssin C, Grimaldi C, Chami M, Daffé M, Guilhot C: A polyketide synthase catalyzes the last condensation step of mycolic acid biosynthesis in mycobacteria and related organisms. Proc Natl Acad Sci U S A. 2004, 101 (1): 314-31937. 10.1073/pnas.0305439101.
Yamanaka K, Murayama C, Takagi H, Hamano Y: Epsilon-poly-L-lysine dispersity is controlled by a highly unusual nonribosomal peptide synthetase. Nat Chem Biol. 2008, 4 (12): 766-772. 10.1038/nchembio.125.
Walton JD: HC-toxin. Phytochem. 2006, 67 (14): 1406-1413. 10.1016/j.phytochem.2006.05.033.
Hoshino Y, Chiba K, Ishino K, Fukai T, Igarashi Y, Yazawa K, Mikami Y, Ishikawa J: Identification of nocobactin NA biosynthetic gene clusters in Nocardia farcinica. J Bacteriol. 2011, 193 (2): 441-448. 10.1128/JB.00897-10.
Mathur M, Kolattukudy PE: Molecular cloning and sequencing of the gene for mycocerosic acid synthase, a novel fatty acid elongating multifunctional enzyme, from Mycobacterium tuberculosis var. bovis Bacillus Calmette-Guerin. J Biol Chem. 1992, 267 (27): 19388-19395.
Brennan PJ, Nikaido H: The envelope of mycobacteria. Annu Rev Biochem. 1995, 64: 29-63. 10.1146/annurev.bi.64.070195.000333.
Kaulmann U, Hertweck C: Biosynthesis of polyunsaturated fatty acids by polyketide synthases. Angew Chem Int Ed Engl. 2002, 41 (11): 1866-1869. 10.1002/1521-3773(20020603)41:11<1866::AID-ANIE1866>3.0.CO;2-3.
Jiang H, Zirkle R, Metz JG, Braun L, Richter L, Van Lanen SG, Shen B: The role of tandem acyl carrier protein domains in polyunsaturated fatty acid biosynthesis. J Am Chem Soc. 2008, 130 (20): 6336-6337. 10.1021/ja801911t.
Wallis JG, Watts JL, Browse J: Polyunsaturated fatty acid synthesis: what will they think of next?. Trends Biochem Sci. 2002, 27 (9): 467-10.1016/S0968-0004(02)02168-0.
Okuyama H, Orikasa Y, Nishida T: Significance of antioxidative functions of eicosapentaenoic and docosahexaenoic acids in marine microorganisms. Appl Environ Microbiol. 2008, 74 (3): 570-574. 10.1128/AEM.02256-07.
Fischbach MA, Walsh CT: Assembly-line enzymology for polyketide and nonribosomal Peptide antibiotics: logic, machinery, and mechanisms. Chem Rev. 2006, 106 (8): 3468-3496. 10.1021/cr0503097.
Weissmam KJ: Introduction to polyketide biosynthesis. Methods Enzymol. 2009, 459: 3-16.
Del Vecchio F, Petkovic H, Kendrew SG, Low L, Wilkinson B, Lill R, Cortes J, Rudd BA, Staunton J, Leadlay PF: Active-site residue, domain and module swaps in modular polyketide synthases. J Ind Microbiol Biotechnol. 2003, 30: 489-494.
Kakavas SJ, Katz L, Stassi D: Identification and characterization of the niddamycin polyketide synthase genes from Streptomyces caelestis. J Bacteriol. 1997, 179: 7515-7522.
Gaisser S, Trefzer A, Stockert S, Kirschning A, Bechthold A: Cloning of an avilamycin biosynthetic gene cluster from Streptomyces viridochromogenes Tu57. J Bacteriol. 1997, 178 (20): 6271-7278.
Van Lanen SG, Oh TJ, Liu W, Wendt-Pienkowski E, Shen B: Characterization of the maduropeptin biosynthetic gene cluster from Actinomadura madurae ATCC 39144 supporting a unifying paradigm for enediyne. Biosynthesis. 2007, 129 (43): 13082-13094.
Zazopoulos E, Huang K, Staffa A, Liu W, Bachmann BO, Nonaka K, Ahlert J, Thorson JS, Shen B, Farnet CM: A genomics-guided approach for discovering and expressing cryptic metabolic pathways. Nat Biotechnol. 2003, 21 (2): 187-190. 10.1038/nbt784.
A genome database of microorganisms sequenced at NITE: http://www.bio.nite.go.jp/dogan/project/view/OPACUS,
McLeod MP, Warren RL, Hsiao WW, Araki N, Myhre M, Fernandes C, Miyazawa D, Wong W, Lillquist AL, Wang D, Dosanjh M, Hara H, Petrescu A, Morin RD, Yang G, Stott JM, Schein JE, Shin H, Smailus D, Siddiqui AS, Marra MA, Jones SJ, Holt R, Brinkman FS, Miyauchi K, Fukuda M, Davies JE, Mohn WW, Eltis LD: 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. 10.1073/pnas.0607048103.
This work was supported by research grants of the Ministry of Education, Culture, Sports, Science, & Technology in Japan [grants #21406003], by the National BioResource Project (http://www.nbrp.jp/), and by the Cooperative Research Grant of NEKKEN, 2012 to TG. We thank to Mr. Syuji Yamazaki and Dr. Atsushi Yamazoe for launching the massive genome sequencing project of the genus Nocardia and for providing us with the sequences.
The GenBank/EMBL/DDBJ accession numbers for the sequences containing the type-I PKS and NRPS gene clusters are: AB685274, AB700124-AB700133, AB700557-AB700568 (Nocardia asteroides NBRC 15531T), AB700569-AB700587 (Nocardia otitidiscaviarum IFM 11049), AB701575-AB701605 (Nocardia brasiliensis IFM 10847), and AB701607-701636 (Nocardia brasiliensis NBRC 14402T).
The authors declare that they have no competing interests.
HK and NI: analyzed gene clusters and predicted their functions; AH, AT-N, and TM: performed sequencing; KS and NF: organized the sequencing project and edited the manuscript; TG: wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: ORFs and module/domain structures of PKS-I, NRPS, and PKS-I/NRPS hybrid gene clusters in genomes. Data for N. otitidiscaviarum IFM 11049, N. asteroides NBRC 15531 T, N. brasiliensis NBRC 14402 T and IFM 10847 are shown. (XLSX 33 KB)
Additional file 2: Figure S1: Representative NRPS gene clusters in N. farcinica and their homologs in other strains. A. N. asteroides has a cluster with an overall similarity to nfa7170-7200; but the third ORF, NCAST_11_00860, is similar to the first ORF nfa7170, rather than the third ORF nfa7190 of the corresponding position. N. brasiliensis NBRC 14402T and IFM 10847 lack ORFs corresponding to nfa7180. B. N. brasiliensis IFM 10847 has only partial sequences of N. farcinica nfa50330-homologous gene, while the homolog in N. otitidiscaviarum is not only partial but also distantly located in the genome. C. N. asteroides possesses an nfa50630 homolog, but lacks an nfa50620 homolog. N. brasiliensis strains have no homologs. (PPTX 105 KB)
Additional file 3: Figure S2: Comparison of putative polyunsaturated fatty acid synthase (PfaA) genes between the genus Nocardia (Figure 2 #11) and the genus Shewanella.(PPTX 76 KB)
Additional file 4: Figure S3: Predicted chemical structure of the product from PKS-I/NRPS hybrid gene cluster #20 in N. otitidiscaviarum.(PPTX 107 KB)
About this article
Cite this article
Komaki, H., Ichikawa, N., Hosoyama, A. et al. Genome based analysis of type-I polyketide synthase and nonribosomal peptide synthetase gene clusters in seven strains of five representative Nocardia species. BMC Genomics 15, 323 (2014). https://doi.org/10.1186/1471-2164-15-323
- Nocardia asteroides
- Nocardia otitidiscaviarum
- Nocardia brasiliensis
- Nocardia farcinica
- Nocardia cyriacigeorgica
- Genome sequence
- Type-I polyketide synthase
- Nonribosomal peptide synthetase