Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships

Bushley, Kathryn E; Turgeon, B Gillian

doi:10.1186/1471-2148-10-26

Research article
Open access
Published: 26 January 2010

Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships

Kathryn E Bushley¹ &
B Gillian Turgeon¹

BMC Evolutionary Biology volume 10, Article number: 26 (2010) Cite this article

13k Accesses
161 Citations
1 Altmetric
Metrics details

Abstract

Background

Nonribosomal peptide synthetases (NRPSs) are multimodular enzymes, found in fungi and bacteria, which biosynthesize peptides without the aid of ribosomes. Although their metabolite products have been the subject of intense investigation due to their life-saving roles as medicinals and injurious roles as mycotoxins and virulence factors, little is known of the phylogenetic relationships of the corresponding NRPSs or whether they can be ranked into subgroups of common function. We identified genes (NPS) encoding NRPS and NRPS-like proteins in 38 fungal genomes and undertook phylogenomic analyses in order to identify fungal NRPS subfamilies, assess taxonomic distribution, evaluate levels of conservation across subfamilies, and address mechanisms of evolution of multimodular NRPSs. We also characterized relationships of fungal NRPSs, a representative sampling of bacterial NRPSs, and related adenylating enzymes, including α-aminoadipate reductases (AARs) involved in lysine biosynthesis in fungi.

Results

Phylogenomic analysis identified nine major subfamilies of fungal NRPSs which fell into two main groups: one corresponds to NPS genes encoding primarily mono/bi-modular enzymes which grouped with bacterial NRPSs and the other includes genes encoding primarily multimodular and exclusively fungal NRPSs. AARs shared a closer phylogenetic relationship to NRPSs than to other acyl-adenylating enzymes. Phylogenetic analyses and taxonomic distribution suggest that several mono/bi-modular subfamilies arose either prior to, or early in, the evolution of fungi, while two multimodular groups appear restricted to and expanded in fungi. The older mono/bi-modular subfamilies show conserved domain architectures suggestive of functional conservation, while multimodular NRPSs, particularly those unique to euascomycetes, show a diversity of architectures and of genetic mechanisms generating this diversity.

Conclusions

This work is the first to characterize subfamilies of fungal NRPSs. Our analyses suggest that mono/bi-modular NRPSs have more ancient origins and more conserved domain architectures than most multimodular NRPSs. It also demonstrates that the α-aminoadipate reductases involved in lysine biosynthesis in fungi are closely related to mono/bi-modular NRPSs. Several groups of mono/bi-modular NRPS metabolites are predicted to play more pivotal roles in cellular metabolism than products of multimodular NRPSs. In contrast, multimodular subfamilies of NRPSs are of more recent origin, are restricted to fungi, show less stable domain architectures, and biosynthesize metabolites which perform more niche-specific functions than mono/bi-modular NRPS products. The euascomycete-only NRPS subfamily, in particular, shows evidence for extensive gain and loss of domains suggestive of the contribution of domain duplication and loss in responding to niche-specific pressures.

Background

Nonribosomal peptide synthetases (NRPSs) are multimodular megasynthases which catalyze biosynthesis of small bioactive peptides (NRPs) via a thiotemplate mechanism independent of ribosomes [1–5]. NRPS encoding genes (NPSs) are plentiful in fungi and bacteria but are not known in plants or animals. The enzymes they encode biosynthesize a staggering diversity of chemical products because their substrates can include both D and L forms of the 20 amino acids used in ribosomal protein synthesis, as well as non-proteinogenic amino acids such as ornithine, imino acids, and hydroxy acids such as α-aminoadipic and α-butyric acids [1]. The natural functions of most NRPs for producing organisms are largely unknown, although recently it has become clearer that they play fundamental roles in fungal reproductive and pathogenic development, morphology, cell surface properties, stress management, and nutrient procurement [6–15] in addition to better-known roles as toxins/mycotoxins involved in plant or animal pathogenesis or as life-saving pharmaceuticals such as antibiotics, immunosuppressants, and anticancer agents.

NRPSs use a set of core domains, known as a module, to accomplish peptide synthesis. A minimal module consists of three core domains: an adenylation (A) domain which recognizes and activates the substrate via adenylation with ATP, a thiolation (T) or peptidyl carrier protein (PCP) domain which binds the activated substrate to a 4'-phosphopantetheine (PP) cofactor via a thioester bond and transfers the substrate to a condensation (C) domain which catalyzes peptide bond formation between adjacent substrates on the megasynthase complex [1]. Several specialized C-terminal domains involved in chain termination and release of the final peptide product have also been identified [16, 17]. In bacteria, chain release is most commonly effected by a thioesterase (TE) domain [18], which releases the peptide by either hydrolysis or internal cyclization [16, 17, 19]. In fungi, only a few NRPSs, such as the ACV synthetases, are known to release products via a TE domain and chain release is carried out by a variety of mechanisms, two of which predominate and occur less frequently in bacterial systems: 1) a terminal C domain, which catalyzes release by inter- or intra-molecular amide bond formation [16], and 2) a thioesterase NADP(H) dependent reductase (R) domain [20–23], which catalyzes reduction with NADPH to form an aldehyde. An additional mechanism, which has been reported only in biosynthesis of fungal ergot alkaloids, involves nonenzymatic cyclization by formation of a diketopiperazine ring [16, 24].

NRPSs may contain additional modifying domains which alter the substrate during NRPS biosynthesis: 1) an epimerization (E) domain which catalyzes epimerization of an amino acid from the L to the D configuration [25], 2) an N-methylation (M) domain (methyltransferase) which catalyzes transfer of a methyl group from an S-adenosylmethionine to an α-amino of the amino acid substrate, and 3) a specialized C domain termed a cyclization (Cyc) domain which catalyzes formation of oxazoline or thiazoline rings by internal cyclization of cysteine, serine, or threonine residues [26]. Additional tailoring enzymes which are not part of the NRPS may modify either the substrate or the final peptide product by glycosylation, hydroxylation, acylation, or halogenation [27, 28].

NRPSs may be monomodular, consisting of a single A-T-C module, or multimodular, consisting of repeated A-T-C modules. The suite of 14 NRPSs found in the genome of the Dothideomycete Cochliobolus heterostrophus is representative of the diversity of NPS genes in filamentous ascomycetes in that it contains a representative from most currently recognized groups of fungal NRPSs [10, 6], and, with the exception of duplicated copies of ChNPS12, the modular domain architectures of each encoded enzyme are distinct (Additional file 1). In addition to mono- and multi-modular NPSs, a hybrid gene (ChNPS7; PKS24) encoding an incomplete NRPS module (A-T) fused to a polyketide synthase (PKS) unit is present [10, 29]. Hybrid PKS;NRPS synthetases (e.g. ACE1, SYN2 in Magnaporthe oryzae), the reverse organization of ChNPS7;PKS24, are more common in filamentous fungi as well as in bacteria [30–34], although C. heterostrophus lacks a representative. PKSs, found in fungi, bacteria, and plants are large megasynthases, related to fatty acid synthases, that biosynthesize small molecule polyketides with as diverse natural functions as NRPS metabolites.

The evolutionary mechanisms giving rise to genes encoding enzymes with such diverse modular architectures are clearly complex. Likely mechanisms include: 1) tandem duplication and loss of individual modules or domains, 2) gene fusion/fission, and 3) recombination and/or gene conversion of individual modules or domains either within the same NPS or between different NPSs. It has been suggested that genes involved in secondary metabolite (small molecule) biosynthesis tend to be located in subtelomeric regions, a factor which may contribute to their rapid evolution by the aforementioned mechanisms [35, 36].

NPSs are generally recognized as a rapidly evolving gene class in fungi leading to few clearly identifiable orthologs between species and highly discontinuous distributions [10, 37, 38]. However, as has been observed for members of other eukaryotic gene families (e.g., major histocompatibility complex [39], immune response [40], zinc-finger [41], reproductive [42], olfactory/chemosensory [43–47], MADS-box [48], and F-box gene families [49] among others), within each family, conservation and rates of gene duplication and loss are likely to vary among subgroups of genes encoding proteins of different function. In fact, some C. heterostrophus NPSs (NPS2, NPS4, NPS6 and NPS10), are conserved or moderately conserved across euascomycote fungi [8, 10, 50] and their NRP products are involved in basic cellular functions such as growth and development, reproduction, and pathogenesis [6–8]. The majority of NPSs, however, are highly discontinuously distributed across fungal taxa and even closely related species may share only a few homologs. Some, e.g., Cochliobolus carbonum HTS1, the gene encoding the NRPS for biosynthesis of HC-toxin [51], and Alternaria alternata apple pathotype AMT, the gene encoding the NRPS for biosynthesis of AM-toxin [52], appear unique even to one race or pathotype within a single species. These lineage-specific synthetases tend to have more specialized, niche-specific functions.

Higher rates of gene duplication and loss may reflect an adaptive response to selective pressure from pathogens, interactions with other organisms, or other environmental pressures. Recent work suggests that, in fungi, genes involved in responses to stress are more likely to undergo duplication and loss than growth related genes [53]. Thus, we hypothesize that NRPSs with conserved functions involved in growth and development will show less variation in gene copy number and maintain a relatively conserved domain architecture in comparison with NRPSs with more niche-specific functions.

The multimodular structure of NRPSs and the complex mechanisms by which they evolve present challenges to phylogenetic analysis and consequently little work has been done to characterize phylogenetic relationships across this large class of megasynthases or to ask whether subclasses of common function can be identified, based on close relationships with NRPSs whose chemical products are known. In this study, we undertook phylogenomic analyses on a comprehensive dataset of fungal NRPS proteins to: 1) identify major subfamilies, 2) analyze patterns of distribution of these major subfamilies across fungal taxonomic groups, 3) understand relationships among selected bacterial NRPSs, fungal monomodular NRPS/NRPS-like proteins, fungal multimodular NRPSs, and related adenylating enzymes, including α-aminoadipate reductases involved in lysine biosynthesis in fungi, 4) consider mechanisms of evolution of multimodular NRPSs, and 5) analyze patterns of NRPS gene and A domain duplication and loss across fungi.

Results and Discussion

Identification and domain structure of candidate NRPSs

Candidate NRPSs extracted from each sequenced genome are listed in Additional file 2. Genus and species abbreviations for all organisms mentioned in this study are shown in Additional file 3. The proposed domain structure for each NRPS, based on searches with our fungal-specific HMMER models (Additional file 4) and the PFAM and Interpro databases, is shown in Additional file 2. The majority of multimodular NRPSs were composed of one or more standard NRPS modules (A-T-C) with or without modifying domains (E, M, etc), while most monomodular NRPSs lacked complete A-T-C modules and consisted of a single A domain or an A-T unit followed by a variety of C-terminal domains, several of which have not previously been identified as core NRPS domains (Additional file 2).

Phylogenomic analysis and subfamily identification

All known NRPS/NRPS-like proteins formed a monophyletic group supported by greater than 90% bootstrap support in ML analyses and greater than 50% bootstrap support in the NJ analysis (Fig. 1), separating them from most other known adenylating enzymes selected as potential outgroups, e.g., Acyl AMP ligases (AAL), CPS1 [54], Long Chain Fatty Acid ligases (LCFAL), Acetyl-CoA synthetases (ACoAS), and Ochratoxin synthetases (OCHRA)(Fig. 1, Additional file 5). The α-aminoadipate reductases (AAR), homologs of S. cerevisiae Lys2 [23, 55, 56, 56, 57], grouped within this well-supported clade of NRPS/NRPS-like proteins rather than with the other adenylating enzymes (Figs. 1, 2, Additional file 6), suggesting that AARs are more closely related to NRPSs than to other adenylating enzymes.

The tree topologies resulting from phylogenetic analyses of individual A domains revealed two major groups of fungal NRPSs (Fig. 1, Additional file 6). The first group (Fig. 1, light blue rectangle) consists of primarily mono- or bi-modular fungal NRPSs which group with bacterial NRPS A domains. Exceptions to the predominately mono/bi-modular fungal NRPS structures include the ACV synthetases and the clade containing A domains from the eleven modules of SimA (cyclosporin biosynthesis) [58] and from several related fungal NRPSs. The other large group contains exclusively fungal and primarily multimodular NRPSs and includes siderophore synthetases and a group we term the Euascomycete-only synthetases, as its members are restricted to euascomycetes. Both grouped together with greater than 97% bootstrap support in analyses of a reduced dataset which included selected representatives from each subfamily (Fig. 2, red arrow, Additional file 7).

Phylogenetic analyses identified nine major subfamilies of fungal NRPSs. Subfamilies were defined as the most internal branch from the root node that formed a monophyletic group which was supported by greater than >70% bootstrap support, shared identical taxon composition across all three phylogenetic methods, and contained a representative fungal NRPS. These groups were named after a representative C. heterostrophus or other fungal NRPSs of well-known function in the group (Figs. 1, 2, Additional file 6). Subfamilies include: 1) fungal PKS;NRPS hybrids, 2) ChNPS11/ETP toxin module 1 synthetases, ChNPS12-like/ETP module 2 toxin-like synthetases, 4) ChNPS10-like synthetases, 5) Cyclosporin synthetases (CYCLO), 6) α-aminoadipate reductases (AAR), 7) ACV synthetases (ACV), 8) siderophore synthetases (SID), and 9) the Euascomycete-only synthetases (EAS). Deep phylogenetic relationships among mono/bi-modular subfamilies were unresolved and lacked bootstrap support (Figs. 1, 2, Additional file 6A-C). A domains from a few ascomycete (BC1G11613.1, MGG 14967.5, MGG07803.5) and several urediniomycete (UM05245.1, Sr31423, and PGTG06519.1) proteins did not group with any of the major subfamilies and were not placed consistently in the trees when assessed by different phylogenetic methods. Homologs of bimodular A. fumigatus SidE, a putative siderophore synthetase [37], formed two clades corresponding to each module and consistently grouped with the SID subfamily but without bootstrap support in the larger phylogeny and with low bootstrap support (>50%) in the reduced phylogeny. We term this group SIDE but do not consider them as a major subfamily (Fig. 2, Additional files 6A-C, 7)

Relationships between fungal and bacterial NRPSs: horizontal transfer or vertical transmission and massive loss?

The majority of bacterial sequences (Additional file 8), identified as top hits in blast searches using a representative from each of the major fungal NRPS subfamilies to query the public databases, were eubacterial in origin and formed a monophyletic group (although lacking bootstrap support), which we term the major bacterial clade (Figs. 1, 2, gray, Additional file 6). This clade contains two fungal representatives suspected of being horizontally transmitted from bacteria to fungi. One is the fungal ChNPS7;PKS24 hybrid NRPS;PKS synthetase which is nested within this clade; previous independent analyses of both the NRPS [10] and the PKS portion of this protein [29] found the same placement (Fig. 2, Additional file 6). The other is the ACV synthetases, a group postulated to have been horizontally transferred from bacteria to fungi [59–64], which groups as sister to, or within, the major bacterial clade (Figs. 1, 2, Additional file 6). Our analysis also shows that each of the three fungal ACV synthetase A domains groups with the corresponding bacterial A domain rather than forming separate clades of fungal and bacterial A domains. These results support previous claims of horizontal transfer based on observations of closer sequence similarity than expected between these fungal and bacterial genes [61–64] (Fig. 2, Additional file 6).

In contrast, bacterial siderophore synthetases (eg. Pyoverdine (PvdD, PvdI, PvdJ, PvdL), yersiniabactin (ybtE), and Pyochelin (PchE, PchF)) group separately from fungal NRPSs (SID) that biosynthesize intracellular siderophores and fungal NRPSs (NPS6) in the EAS subfamily that biosynthesize extracellular siderophores (Fig. 2, Additional file 6). This suggests that fungal and bacterial capacities to chelate iron via small molecule siderophores have evolved independently (Fig. 2).

The remaining bacterial A domains included in this study that grouped with high bootstrap support with fungal A domains were associated with the ChNPS12-like/ETP module 2 and AAR subfamilies (Fig. 2, Additional file 6). Bacterial representatives of the former included representatives from the gammaproteobacterium (Hahella chejuensis) and two closely related species of actinobacteria (Salinospora sp.); of the latter, all are closely related species of Pseudomonads. In the two cases of proposed horizontal transfer discussed above (e.g., ChNPS7 [10, 29], and ACV [59–64] synthetases), the fungal genes are nested within a large clade of bacterial sequences. The reverse phylogenetic situation is observed for bacterial genes grouping with the AAR and ChNPS12 subfamilies as, in these cases, bacterial NRPSs, are nested within a large clade of fungal NRPSs (ChNPS12) or group as sister to fungal NRPSs (AARs). These placements suggest that either the fungal genes were transferred to bacteria or that the origin of these groups predates the divergence of eukaryotes and prokaryotes and the observed pattern reflects extensive loss or incomplete sampling from bacteria. Clearly, further sampling of bacterial sequences is needed to adequately address these hypotheses, but we favor the theory that these NRPS subfamilies may have originated prior to the divergence of prokaryotes and eukaryotes. We hypothesize that the lack of phylogenetic signal for resolving relationships among the fungal mono/bi-modular subfamilies may in part reflect an ancient and rapid radiation of these groups.

Distribution of NRPS subfamilies across fungal taxonomic groups

The distribution of fungal NRPS subfamilies across the major fungal taxonomic groups supports previous findings that NRPSs are much more abundant in Euascomycetes than in Basidiomycetes and are scarce in Chytridiomycota, Zygomycota, Schizosaccharomycota, and Hemiascomycota [10, 65, 66]. The number and distribution of NRPSs in each subfamily are shown in Table 1. EAS and PKS;NRPS subfamilies were significantly overrepresented in Euascomycete taxa when evaluated by Fisher's exact tests, while ChNPS12-like synthetases were statistically overrepresented in Basidiomycete taxa (Table 1, asterisks). The Chytridiomycota, Zygomycota, Schizosaccharomycota, and Hemiascomycota contained only a few NRPSs. All Zygomycota and Hemiascomycota lacked genes encoding NRPS-type proteins other than a single AAR. The chytrid genome contained two additional NRPS-like proteins grouping with the ChNPS12/ETP module 2 subfamily, and the two Schizosaccharomycota taxa examined contained one additional NRPS for siderophore biosynthesis (Table 1). No subfamilies were statistically overrepresented in these groups.

Table 1 Numbers of NRPS genes in each subfamily across fungal taxonomic groups^a

Full size table

Lineage specific expansions and contractions

When patterns of gene duplication and loss were analyzed for the total number of NRPSs/genome (combining all subfamilies) over the tree of fungi (Fig. 3; Additional file 9), a highly significant expansion was found on the branch leading to Euascomycetes (p = 7 × 10^-5). Significant expansions were also found within euascomycetes on the branches leading to the Aspergillus species (p = .028), to F. graminearum (p = .011) and to M. oryzae (p = .032). N. crassa showed a highly significant (p = 5 × 10^-5) contraction in total number of NRPS genes (Fig. 3), likely due to the efficiency of RIP and/or other genome defense mechanisms, that reduce the rate of fixation of duplicates [67, 68].

Our data support previous findings [66], including our own [10], that unicellular fungi have few, if any, genes for secondary metabolism (Table 1, Fig. 3). Ancestral reconstructions show that in hemiascomycete yeasts, this is due to loss of all NRPSs, except for a single AAR encoding gene, that were present in basidiomycetes and inferred to be present in the ancestor of ascomycetes (Fig. 3). However, both the fission yeast S. pombe and the unicellular basidiomycete yeast Sporobolomyces roseus contain one additional gene encoding a NRPS (a siderophore synthetase and an unknown, respectively) in additional to the single AAR encoding gene, suggesting that a unicellular habit may not preclude the existence of secondary metabolite genes such as NRPSs. Patterns of expansion and contraction also do not seem to occur preferentially in fungal pathogens versus nonpathogens. While a number of pathogenic fungi (e.g., F. graminearum, A. fumigatus, and M. oryzae) do show evidence for expansions in numbers of NRPS, we also see expansion in the nonpathogen, A. nidulans.

Subfamily distribution

AAR

A single ortholog of S. cerevisiae Lys2, the AAR involved in reduction of α-aminoadipic acid in the fungal lysine biosynthetic pathway [23, 69], was found in all fungi surveyed except the Microsporidian, Enchephalitozoon cuniculi, an intracellular parasite which has lost the majority of genes involved in amino acid biosynthesis [70] and the heterokaryotic basidiomycete Postia placenta, which contains two (Table 1)(Additional file 2).

ChNPS11/ETP/ChNPS12

In a phylogeny of a reduced set of representative A domains from each subfamily (Fig. 2), homologs of ChNPS11, ChNPS12, and the ETP toxin synthetases, GliP for Gliotoxin and SirP for Sirodesmin production, group together with strong bootstrap support (>80%), suggesting all share a common evolutionary origin. In the larger phylogeny of the complete dataset (Fig. 1, Additional file 6), they formed two separate clades each supported by >70% bootstrap support, but lacked this level of support for the entire group. The first clade (ChNPS11/ETP module 1) includes the first module of the ETP toxin synthetases and monomodular ChNPS11. The second module of the ETP toxin synthetases, however, groups within a larger clade containing the two NRPSs from the chytrid genome, several eubacterial NRPSs, and a clade containing both euascomycete and basidiomycete homologs of ChNPS12 (ChNPS12/ETP module 2). While fungal NRPSs associated with ChNPS11 and ETP toxin synthetases are found only in Euascomycota, NRPSs from both eubacteria and from the most basal fungal group, Chytridiomycota, were nested within this larger clade with high bootstrap support (>80%) (Figs. 1, 2).

ChNPS10, CYCLO, SID

Three subfamilies, monomodular ChNPS10, NRPSs grouping with SIMA (CYCLO), and NRPSs (SID) involved in intracellular (primarily) siderophore biosynthesis, contain representatives from both Basidiomycota and Euascomycota. While all euascomycetes and many basidiomycetes examined contain at least one representative from the SID subfamily (Table 1) [65], ChNPS10 and CYCLO are more discontinuously distributed and a representative is not found in all taxa (Table 1, Additional file. 6).

ACV and PKS;NRPS

PKS;NRPSs were restricted to and statistically overrepresented in euascomycetes. As has been noted previously [29, 71], all fungal PKS;NRPS hybrids fall into a single, well supported, monophyletic group, which suggests a single origin (Table 1). However, not all ascomycetes have a representative of this group and the number of corresponding genes varies widely among taxa (Table 1). C. heterostrophus, for example, lacks a representative but M. oryzae has six. While ACV synthetases are found in both bacteria and fungi, within fungi, they appear restricted to Eurotiomycete and Hypocrealean taxa. This study did not identify any additional ACV synthetases in fungi apart from the known ones in Penicillium chrysogenum, Aspergillus nidulans, and Cephalosporium acremonium (Additional files 2, 6), supporting previous conclusions that their distribution is likely the product of one or more isolated horizontal transfer events [59, 61–64].

Euascomycete only (EAS)

The EAS subfamily contains by far the greatest number of NRPSs and is both restricted to and statistically overrepresented in Euascomycetes (Table 1).

Hypothesized origins based on taxonomic distribution

Fig. 4 shows the hypothesized origins of each subfamily based on taxonomic distribution of the oldest member of each group. By this criterion, the presence of bacterial sequences grouping within the ChNPS11/ETP module 1 and ChNPS12/ETP module 2 clades suggests that the origins of these groups may predate the divergence of eubacteria and eukaryotes (Figs. 2, 4). The AAR subfamily must have arisen also either prior to or very early in the origin of the fungi as a representative is present in all fungi, including the most basal group, the Chytridiomycota (Table 1, Figs. 2, 4). Since the SID, CYCLO, and ChNPS10 subfamilies all contain representatives from both Euascomycota and Basidiomycota, these groups must have evolved prior to the divergence of the Dikarya (Fig. 4). The EAS, PKS;NRPS, and ACV synthetases contained only euascomycete representatives. Both PKS;NRPS and EAS may thus have originated in the ancestor of euascomycetes (Fig. 4). As discussed above, the grouping of fungal ACV synthetase A domains with the corresponding A domains of bacterial ACV synthetases within a large clade of bacterial sequences provides evidence for horizontal transfer and suggests that this group originated within prokaryotes (Fig. 4).

Thus, taxonomic distributions suggest a more ancient origin of one or more of the mono/bi-modular NRPS subfamilies (ChNPS11/ETP/ChNPS12, ACV), possibly predating the divergence of eubacteria and fungi (Table 1, Fig. 4). The strongly supported co-grouping of fungal and bacterial sequences in the ChNPS11/ETP/ChNPS12 group, as in the outgroup adenylating enzymes (Fig. 2, Additional file 6A-C) suggests this is a tenable hypothesis. ACV is a special case and likely the result of horizontal transfer from bacteria to fungi. In contrast, the fungal-specific multimodular groups (SID and EAS), which group together with high bootstrap support in the reduced phylogeny (Table 1, Fig. 2, Additional file 6A-C), appear to be of more recent origin and are restricted to and highly expanded in fungi.

Mono- and bi- modular NRPS subfamilies

Unlike many of the multimodular NRPSs, most monomodular subfamilies lack a complete NRPS module (A-T-C) and consist of a single A domain or an A-T domain combination followed by a variety of C-terminal domains (Fig. 5). Many of the mono/bi-modular groups show a conserved domain architecture across all members in a subfamily, suggesting their domain architectures may be functionally constrained. Available functional data suggest that the NRP products of several of these groups may play more central roles in cellular metabolism related to responses to oxidative stress and growth and development.

Whether monomodular NRPSs may act alone or in concert with non-NRPS proteins is currently unknown. However, in bacterial systems, both single A domains as well as A-T domain units, known as initiation modules, can interact with other NRPS proteins and accomplish biosynthesis by first activating and then transfering the activated substrate either to a C domain in the same NRPS or to a C domain in a different NRPS (nonlinear biosynthesis) [5].

AARs and Lysine biosynthesis

AARs are conserved not only taxonomically but also in terms of domain structure. All have an identical structure consisting of an A-T unit followed by a thioester reductase (R) domain (IPR010080), a member of the NAD(P)-binding Rossman fold domain superfamily (SSF51735). There are two primary pathways for lysine biosynthesis, the diaminopimelic acid pathway (DAP), found predominantly in bacteria and plants, and the α-aminoadipate pathway (AAA), found primarily in fungi and a few bacteria [69]. As noted above, AARs catalyze reduction of α-aminoadipic acid in the AAA pathway [69]. The fact that AARs have a C-terminal R domain in common with several other NRPS subfamilies (PKS;NPRS, ChNPS10, EAS, discussed below) supports our conclusions based on phylogenetic relationships that AARs are more closely related to NRPSs than other adenylating enzymes (Fig. 5).

Bacterial sequences grouping with fungal AAR are comprised of a single A domain followed by an acyl-transferase domain (PFAM01757) but lack the C-terminal R domain found in fungal AARs. We conclude that they are likely not involved in lysine biosynthesis in bacteria. Although there is evidence for the existence of lysine biosynthesis through the AAA pathway in some prokaryotes [72], current data suggests that these pathways do not include a step involving reduction of α-aminoadipic acid [72]. Thus, our data support previous conclusions that AARs are fungal-specific enzymes [73–75].

PKS;NRPS

Nearly all fungal PKS;NRPS hybrids have the same domain structure (KS-AT-M-KR-ACP-C-A-T-R) (Fig. 5, Additional file 2). The terminal R domain has been reported previously in several PKS;NRPS hybrids [76–78].

ChNPS10

The ChNPS10 subfamily also has a conserved domain architecture across all genes in the subfamily, consisting of an A-T unit followed by two additional C-terminal domains. The first is a NAD(P) binding domain (IPR016040) also showing closest similarity to thioesterase reductase (R) domains and the second is a dehydrogenase domain with closest hits to ADH short chain dehydrogenases (IPR002198) (Fig. 5).

ChNPS11/ETP/ChNPS12

The large and highly diverse clade of ChNPS11/ETP/ChNPS12 homologs reveals the diversity of C-terminal domains that can follow A-T units and shows that, as for some bacterial NRPSs, fungal NRPS or NRPS-like proteins can consist of single A domains (Figs. 5, 6).

At the base of this group are monomodular ChNPS11 and module 1 of the bimodular ETP toxin synthetases, SirP and GliP, which contain complete A-T-C modules (Fig. 6). Module 2 of SirP and GliP groups at the base of the ChNPS12/ETP module 2 clade. The second module of the ETP toxin synthetases contains a complete module followed by an additional T domain (A-T-C-T) (Fig. 6). This group also contains several fungal proteins with an incomplete (MGG15248.6) or a degenerate (BC1G07441_07442.1) first module (Figs. 2, 5).

Nested within this clade is a group of bacterial NRPSs with a single A domain and two NRPS-like proteins from the chytrid B. dentrobatidis (Fig. 6). One of the chytrid NRPSs (BDEG_03514.1) has a T-C-T-A-T domain architecture followed by a domain with similarity to FSH1 (IPR005645), a serine hydrolase domain. The other chytrid protein (BDEG_08447.1) has an A-T unit followed by two additional domains. The first shows closest similarity to polynucleotidyl transferase, Ribonuclease H fold (IPR012337), a domain associated with nucleic acid binding functions and found in a variety of proteins including HIV RNase H, transposases, and exonucleases [79, 80] (Fig. 6). The second domain shows closest similarity to the membrane-associated domain LPS-induced tumor necrosis factor alpha factor (LITAF, IPR006629, PF10601), which contains a characteristic cysteine rich zinc-binding motif found also in intracellular Zn²⁺ binding proteins and animal transcription factors. The zinc and DNA-binding domains found in the chytrid NRPSs are intriguing (Fig. 5). Gliotoxin and Sirodesmin PL have been shown to inhibit viral reverse transcriptase [81] and general transcription [82], respectively. In the case of Sirodesmin PL, the addition of zinc and other IIB series metals (Hg and Cd) both decreases toxin production in Leptosphaeria maculans and also reverses the inhibition of transcription, suggesting interactions of Sirodesmin PL with either cellular zinc or zinc-containing metalloenzymes such as RNA polymerases [82, 83]. Whether these phenotypes relate to our identification of Zn-binding domains in the corresponding chytrid NRPS is unknown.

ChNPS12 (CocheC5_118012), and its paralog (CocheC5_116719) contain a single A domain followed by a domain showing closest similarity to a ferric reductase transmembrane domain (IPR013130). The closest homologs of ChNPS12 (Fig. 6, euascomycete ChNPS12 group 1) are present in both euascomycete and basidiomycete group 1 and have the same domain structure as the C. heterostrophus NPS12 proteins (Fig. 6). Sister to all group 1 NPS12-like proteins is a group of proteins consisting of standalone A domains (Fig. 6, basidiomycete NPS12 group 2). These were found only in the brown-rot heterokaryotic fungus, P. placenta, which carries eight closely related copies.

The monomodular bacterial NRPSs nested within the ChNPS12/ETP module 2 subfamily also consist of a standalone A domains. As noted earlier, for many bacterial NRPS systems (e.g., VibE, MxcE, and YbtE), single A domains may be involved in NRPS biosynthesis by activating and transferring the activated substrate to a different NRPS [5]. Only one example of this type of synthesis has been reported for fungi (e.g., C. purpurea ergot alkaloid biosynthesis) [5, 84], but our identification of these single fungal A domains grouping with other known NPRSs (e.g., ETP toxins) (Figs. 2, 6, Additional file 6) suggests that this mechanism could be more common in fungi than previously appreciated.

The diversity of domain structures found within the ChNPS11/ETP/ChNPS12 group leads us to hypothesize that there may be several distinct functional groups within this clade.

Multimodular NRPS subfamilies

The majority of multimodular NRPSs are found in the SID and EAS subfamilies. These subfamilies group together with high bootstrap support (>97%) in analyses of the reduced dataset (Fig. 2). Analyses that included a larger number of bacterial sequences (KE Bushley and BG Turgeon, unpublished) support our phylogenetic and distribution data that the SID and EAS subfamilies are restricted to fungi. As noted above, two subfamilies containing genes encoding multimodular NRPSs, the CYCLO and ACV synthetases, group with the primarily mono/bi-modular suite of NRPSs. (Table 1, Fig. 2). SID synthetases show a relatively conserved domain architecture, are present in the majority of euascomycetes sampled, and are thought to have evolved by module duplication and selective loss of A domains or complete modules, as described in detail in Bushley et al. [65].

Diversity within the EAS subfamily

The EAS subfamily, in additional to containing the vast majority of fungal NRPSs, also shows the greatest diversity of both domain architecture and function (Figs. 2, 7, Additional file 2). It includes proteins that are both structurally and functionally conserved (e.g. homologs of ChNPS6 which biosynthesize extracellular siderophores), as well as those that are highly lineage specific (e.g. HTS1 [51] and AMT [52] synthetases for host selective toxins, Tex1 [85] and other peptaibol synthetases, and ergot alkaloid synthetases). The highly diverse domain architectures and discontinuous distribution of corresponding A domains make the identification of orthologs across species extremely challenging.

ChNPS6/PerA

Perhaps the only group for which orthologs can be clearly identified are homologs of the most conserved NRPS in the EAS clade, ChNPS6, which biosynthesizes an extracellular iron scavenging siderophore that serves as a virulence factor for several fungi and is also involved in combating oxidative stress [10, 6] (Figs. 2, 6, 7). Although ChNPS6 appears to have undergone a gene duplication event, it is single copy in all species examined except Trichoderma reesii (Fig. 7), which contains two paralagous copies. All ChNPS6 homologs have a highly conserved domain structure consisting of a single A-T-C module followed by a module with a degenerate A domain (dA-T-C) [10]. Sister to the ChNPS6 group is a clade containing both ChNPS8 and an Epichloe festucae NRPS, PerA; the latter NRPS mediates symbiotic interactions of E. festuca with its grass host by producing an NRP insect deterrent, peramine [86] (Fig. 7, arrow).

Ergot alkaloid synthetases

NRPSs synthesizing ergot alkaloids consistently grouped sister to the ChNPS6 and ChNPS8/PerA clade but without bootstrap support (Figs. 2, 7). These synthetases were found only in animal pathogens in the Eurotiales and grass endophytes such as C. purpurea (Figs. 2, 7). Given that grass endophytes such as C. purpurea are thought to have an animal pathogenic ancestor [87] and that their ergot alkaloid NRP products have toxic effects on livestock and other animals [88–91], we hypothesize that NRPSs synthesizing ergot alkaloids originally evolved to function in animal pathogenesis.

Peptaibol synthetases

Peptaibol synthetases, which were restricted to the Hypocrealean taxa examined in this study (Trichoderma/Hypocrea), also formed a well supported group. However, as discussed below, several modules of each peptaibol synthetase group outside of the main clade (Table 1, Figs. 2, 7)

Dothideomycete host-selective toxin synthetases

A domains of the A. alternata apple pathotype-specific AMT synthetase which produces the host-selective toxin, AM toxin, grouped consistently with modules 1 and 3 of ChNPS1 and ChNPS3 (discussed below). Modules of tetramodular C. carbonum HTS1, responsible for biosynthesis of another host selective toxin, the cyclic tetrapeptide, HC-toxin, grouped in disparate locations in the EAS clade such that clear homologs of HTS1 A domains were not recognizable in any of the species in our dataset (Figs. 2, 7).

ChNPS4

Both HTS1 and ChNPS4 A domain relationships exemplify the challenges of identifying orthologs within the rapidly evolving EAS clade. ChNPS4 has been shown to play a role in C. heterostrophus conidial cell surface hydrophobicity and a homolog in the related Dothideomycete Alternaria brassicicola plays a role in conidial wall development and integrity [9]. Each of the four A domains from each module of tetramodular ChNPS4 groups with strong support with the corresponding A domains of tetramodular AbNPS1 in the closely related Dothideomycete, A. brassicae. These A domains group within a larger clade containing Metarhizium anisopliae NRPS PesA although without bootstrap support (Figs. 2, 7, Additional file 6). However, A domains from NRPSs found in other euascomycetes that group with each of the ChNPS4 modules contain from two to six modules. While some of these A domains are clearly related to those of ChNPS4, module duplication and loss obscure the history of this group.

Evolutionary mechanisms giving rise to multimodular NRPSs

The greater diversity of domain architectures seen in multimodular NRPSs is likely due to the multiplicity of evolutionary mechanisms which may generate the corresponding multimodular genes. The EAS subfamily, in particular, contains NRPSs varying from monomodular proteins involved in ergot alkaloid biosynthesis (PS2 and PS4) and ChNPS6 (which has one complete and one degenerate A domain) to the eighteen module TEX1 synthetase responsible for peptaibol biosynthesis in Trichoderma virens (Hypocrea virens) [85] (Figs. 2, 7, Additional file 6). Several subgroups within the EAS illustrate some of the mechanisms by which the diverse domain architectures of multimodular NRPSs may arise.

Tandem Duplication

Cyclosporin synthetase (SimA) is a clear example of tandem duplication of modules of an NRPS in a single species (Tolypocladium inflatum). All eleven A domains from this protein group together as a single well supported monophyletic group (Fig. 2) which also includes certain A domains from other fungal NRPSs, such as ChNPS1 module 2 and ChNPS3 modules 2 and 4.

Peptaibol synthetases illustrate a more complex process of tandem duplication of modules of an NRPS. Peptaibol synthetases are highly lineage specific and found only within the Hypocreales to date. Using H. virens TEX1 as a point of reference, we found that all modules of TEX1 group together in three separate, well-supported clades with modules of two peptaibol synthetases (Trire2_23171 and Trire2_123786) in the related species, Trichoderma reesii (Figs. 2, 7). TEX1 module 13 falls outside of the other two TEX1 clades (Figs. 2, 7, 8), The nearly one-to-one relationship between modules of TEX1 and modules of T. reesii Trire2_23171 suggests that tandem duplication of modules giving rise to these orthologous genes must have occurred prior to divergence of these two species (Fig. 8). However, at least one additional internal duplication has occurred since divergence from an ancestral species (e.g., note the relationship between T. reesii Trire2_23171 modules 18 and 19) (Fig. 8). The relationship of these two peptaibol synthetases with the T. reesii 14 module peptaibol synthetases, Trire2_123786 is less straightforward. However, we note that certain A domains from Trire2_123786 modules 2, 6, and 11 form widowed branches at the base of clades which contain A domains of at least two, and more often, all three peptaibol synthetases (Fig. 8, stippled boxes). We hypothesize that these may be ancestral domains. Previous studies suggest that like T. reesii, T. virens also harbors additional NRPSs involved in peptaibol biosynthesis [92].

Recombination

Two NRPSs found in C. heterostrophus, ChNPS1 and ChNPS3, demonstrate the potential role of recombination and modular rearrangement in the generation of multimodular NRPSs. Modules 1 and 3 of both ChNPS1 and ChNPS3 group within the EAS subfamily with AMT synthetase, a lineage specific NRPS found only in a single strain of related A. alternata [52] (Figs. 2, 7, Additional file 6A-C). Module 2 of ChNPS1 and modules 2 and 4 of ChNPS3, however, group with the CYCLO synthetases among the mono/bi-modular NRPS subfamilies (Fig. 2, Additional file 6). The phylogenetically unlinked locations of ChNPS1 and ChNPS3 modules in the larger phylogeny suggests that a recombination event must have given rise to the extant genes in C. heterostrophus (Fig. 9). A domains of several other euascomycete NRPSs, for example, bimodular Fusarium equiseti Enniatin synthetase (FeESYN1) and trimodular M. oryzae, MGG00022, also show recombinant structures. Module 1 A domains of both proteins group in the EAS clade with the C. heterostrophus pseudogene ChNPS13, but without bootstrap support (Fig. 2), at positions distinct from modules 1 and 3 of ChNPS1 and ChNPS3. The C-terminal A domain of ESYN1 and the A domains of the final two modules of MGG00022 group in the CYCLO clade (Figs. 2, 9), like module 2 of ChNPS1 and modules 2 and 4 of ChNPS3. Thus, homologs of modules of ChNPS1 and ChNPS3 appear in different combinations in other fungi and demonstrate that recombination plays an important role in the evolution of multimodular NRPSs.

Stability of NRPS gene copy number and domain architectures across subfamilies

Many multigene families experience gene duplication and loss and evolve by a birth-death process [93–96]. Variation in gene copy number resulting from gene duplication and loss is thought to be influenced by both functional and dosage requirements as well as random processes such as genomic drift [43, 44, 97, 98]. Recent studies suggest that functionally conserved genes, such as those involved in growth and development or other basic cellular processes, tend to experience both less variation in copy number [53] and more stable domain organizations [49] than genes involved in environmental and stress responses [53, 99].

For multimodular genes such as NRPSs, duplication and loss or birth-death evolution [93–95] can occur at two hierarchical levels: 1) at the level of the whole gene, and 2) at the level of domains within a gene (intragenic). In the latter case, genes encoding NRPSs whose products are involved in more conserved functions, such as the AARs, would be expected to have more stable domain architectures than those encoding proteins with niche-specific functions. The latter may experience less functional constraint allowing for flexible gain and loss of domains leading to diversity of domain structures. Because NRPS A domains are involved in substrate selection [100, 101], their loss or gain could result in a rapid change in the chemical product of an NRPS.

The range of variation in copy number of NRPS-encoding genes and in number of A domains/NRPS for each subfamily is shown for Euascomycete taxa only in Fig. 10. Variation in gene copy number is the highest for the EAS subfamily but both the PKS;NRPS and ChNPS12 subfamilies also show substantial variation (Fig. 10A). The EAS subfamily also shows by far the greatest variation in number of A domains/NRPS, followed by CYCLO and SID subfamilies, suggestive of less stable domain architectures and higher rates of intragenic domain duplication for these three groups. All of the remaining mono/bi-modular subfamilies show remarkably conserved domain architectures (Fig. 5, 10B), supporting available functional data which suggests these groups may have more central conserved roles in metabolism.

When we compared gene and domain duplication and loss in different subfamilies across euascomycetes, no particular subfamily showed significant evidence for nonrandom expansion or contraction of number of genes. When patterns of the total number of A domains per subfamily were analyzed, the EAS subfamily was the only group which showed highly significant (P < .00001) deviation from a random birth-death process (data not shown). These results support other observations that gain and loss of domains is an important evolutionary force within the EAS subfamily and may represent an adaptive response to niche-specific environmental pressures.

Chain termination mechanisms

Our survey revealed that fungal NRPSs have a variety of C-terminal domains involved in chain termination. The most common for multimodular NRPSs is a C domain while for monomodular NRPSs it is an R domain (Additional file 2). R domains have previously been identified and shown to play a role in peptide release in fungal AARs [23, 56], a number of fungal PKS;NRPSs [76–78], and in a minority of bacterial NRPSs including SafA and MxcG [20, 21] and the PKS;NRPS hybrid, myxalamid [22]. Some multimodular NRPSs, however, also have a terminal R domain suggesting this may be a common release mechanism for fungal NRPSs. (Additional file 2). Two different release mechanisms have been identified for R type domains in fungal NRPSs, indicating the possibility of R domains subtypes. In fungal AAR's, the R domain reduces the enzyme bound α-aminoadipic acid [23]. The C-terminal R domain in the fungal PKS;NRPSs for Equisetin biosynthesis (EqiS), however, catalyzes a Dieckmann condensation reaction, thus performing a function similar to bacterial TE domains [76]. Some mono- and some multi-modular NRPSs terminate in T domains (Additional file 2) although these have not been implicated previously in chain release.

As noted previously, bacterial NRPSs generally have a TE domain at the C-terminal end for peptide release but TE domains have been found only in a few fungal NRPSs, notably the ACV synthetases [16]. We identified several other fungal NRPSs (AN2621.4, FGSG_11989.3, and Phchr1_2706), grouping with modules of Cyclosporin synthetase, which also contain a C-terminal TE domain (Additional file 2). However, our data suggest that TE domains are indeed rare in fungal NRPSs providing further support for the claims of horizontal transfer from bacteria to fungi of genes encoding ACV synthetases, and possibly these other fungal genes with TE domains.

Conclusions

Phylogenomic analysis identified nine major subfamilies of fungal NRPSs which fall into two main groups: 1) a group of primarily mono/bi-modular enzymes (ChNPS10, AAR, ChNPS12, ChNPS11/ETP, PKS:NRPS, and CYCLO subfamilies) that group with bacterial NRPSs, and 2) a group of primarily multimodular proteins (EAS, SID) which appear both restricted to and highly expanded within fungi. Analyses demonstrate that α-aminoadipate reductases are more closely related to NRPSs than to other adenylating enzymes and provide further support for previous claims of horizontal transfer of certain NRPSs from bacteria to fungi. In addition, phylogenomic relationships among subfamilies, taxonomic distributions, structural conservation of domain architecture, and known data on function suggest that several of the mono/bi-modular groups are older in origin and play more central roles in cellular metabolism. The highly expanded group of fungal multimodular NRPSs, particularly the EAS subfamily, have less conserved domain architectures due to domain/module duplication and loss, and tend to perform more niche-specific functions, typically considered the realm of "secondary" metabolites.

Methods

Identification of putative NRPSs in fungal genomes

A set of fungal NRPSs with known chemical products was extracted from the NCBI database (Additional file 10), aligned using MUSCLE [102] with the 13 NRPSs identified previously in the Dothideomycete, C. heterostrophus C4 strain [10], and used to construct an initial HMMER model of fungal NRPS A domains using HMMER 2.0 http://hmmer.janelia.org (Additional file 11). This model was tested for specificity and ability to identify NRPSs proteins in fungal genomes for which NRPSs have been well characterized (e.g., C. heterostrophus and Gibberella zeae/Fusarium graminearum) and was found to correctly identify all known NRPSs in the genomes of these species as top hits. Protein datasets of a taxonomically representative sample of fungal genomes (Additional file 12) were downloaded and searched using both a local and global version of the fungal NRPS HMMER model. Proteins that were hit by our A domain model with an e-value less than 1 were considered possible NRPSs. A similar search strategy was employed on the nucleotide genome sequences using GENEWISE [103] and the same HMMER model to identify candidates that might have been missed or mis-annotated by automated gene calling programs. This approach did not identify any additional genes but did identify missed domains and also revealed a number of split gene annotations in the automated protein calls which we have reannotated. These included BC1G09040_09041.1, BC1G07441_07442.1, and FGSG11659.3 and FGSG11630.3 which we conclude represents a single gene corresponding to the MIPS and version 2 broad annotation (FG_00042.1), (Additional file 2).

For each fungal genome, A domains from all candidate NRPSs were aligned, using MUSCLE [102], with A domains from the 12 NRPSs previously identified from C. heterostrophus [10] (Additional file 1) and with A domains from related adenylating enzymes in the AMP-binding family (PFAM PF00501) [e.g., acyl CoA ligases (ACoAL), acetyl CoA synthetases (ACoAS), acyl AMP ligases (AAL), homologs of C. heterostrophus CPS1 (CPS1) [54], long chain fatty acid ligases (LCFAL), and homologs of Ochratoxin synthetase (OCHRA) [104] (Additional file 5). An initial phylogenetic analysis was conducted using the WAG+G model in PhyML to define a set of candidate NRPS proteins for each genome. Proteins from each genome grouping within a monophyletic group containing A domains of the known C. heterostrophus NRPS proteins and separated from the outgroup proteins with consistently high bootstrap support (>90), were retained in the dataset as candidate NRPSs or NRPS-like proteins. We chose to use individual A domains, rather than to include only proteins containing a complete A-T-C module as has been used in previous studies [105] because the latter would miss several putative NRPS or NRPS-like proteins (e.g. C. heterostrophus NPS10 and NPS12 [10]) that lack a complete A-T-C module. In addition, freestanding A domains in bacterial NRPSs have been shown to catalyze NRPS biosynthesis by activating and transferring substrates in trans to separate NRPSs [5] and the evolutionary relationship between monomodular NRPS-like proteins and multimodular NRPSs was also of interest.

Annotation of domain architectures

All candidate proteins were annotated with our initial fungal NRPS A model and the PFAM models for C (PF00668) and T (PF00550) domains. Using the domains identified in the dataset from this search, a refined set of fungal specific NRPS HMMER models was built for the A (FungalNPSAMP.hmm), C (FungalNPSCON.hmm), and T (FungalNPSTHIOL.hmm) domains (Additional file 4). These models more accurately identified C and T domains in NRPSs with known/manually curated annotations than the generic PFAM models and were thus used to annotate A-T-C domain structures of all candidate fungal NRPSs. In addition, all candidate proteins were used as queries against the PFAM and INTERPRO domain databases to identify additional non-canonical NRPS domains present in these proteins. A complete domain architecture was compiled for each protein by merging these two approaches (Additional file 2).

Phylogenomic analyses

Representatives of both fungal and bacterial adenylating enzymes used as outgroups (Additional file 5) in identification of putative NRPSs were also used as outgroups in phylogenomic analyses. While all AARs grouped as putative NRPSs, to reduce the size of the dataset, only a taxonomically representative sample of the fungal AARs were included in the full phylogenetic analyses. Fungal A domains from NRPSs with known function and/or chemical products present in GenBank were also included (Additional file 10). To select a diverse group of bacterial proteins, a representative A domain of each subfamily of fungal NRPSs was used to query the nr protein database at NCBI and the top 5 bacterial protein hits for each, as well as a number of bacterial proteins with known chemical products, were selected (Additional file 8). The complete set of A domains were extracted from these 58 bacterial proteins for a total of 99 A domains.

All candidate NRPS and outgroup A domains were aligned with MUSCLE [102]. Portions of ambiguous alignment were first adjusted manually and then masked to remove columns in the alignment with > 30% gaps prior to phylogenetic analysis (Additional file 13). A few candidate A domains were partial (BC1G15479, FG11319, AN8504, and Pa3740) and were removed from the final analysis because they did not align well with other NRPSs. ProtTest [106] was used to identify an appropriate protein substitution matrix as it has been shown that spurious choice of a matrix can lead to inaccurate phylogenies [107]. The RtREV+G+F model had the best likelihood score for all criteria (AIC and BIC) except for AIC-1 with sample size corrected for the number of sites in the alignment, which identified WAG+G as the best model. Three methods were used for phylogeny construction: 1) Maximum likelihood (ML) using RaxML [108] with the RtREV+G+F substitution model, 2) ML using PhyML with the WAG+G model [109], and 3) Neighbor joining (NJ) using NEIGHBOR in PHYLIP [110] and a distance matrix created in TREEPUZZLE [111] with the WAG+G substitution model. We used a Gamma distribution with four rate categories to model rate variation in all analyses. Bootstrapping was performed to assess the robustness of the phylogeny. Bootstrap datasets of 500 replicates for ML analysis and 200 replicates for the NJ analyses were created using SEQBOOT in PHYLIP and analyzed by the respective methods.

Because bootstrap support has been observed to decline in larger datasets [112–114], we also performed analyses on a subset of the data containing representatives from each of the major subfamilies identified. This dataset was aligned separately with MUSCLE and also masked with slightly less stringent conditions to remove columns containing greater than 50% gaps (Additional file 14). Phylogenetic analyses were performed on this dataset using the same methods described above.

Alignments have been deposited in TREEBASE (http://www.treebase.org/treebase/index.html, Study accession number = S2573 Matrix accession number = M4916).

Subfamily identification and modelling

Fungal NRPS subfamilies were characterized as monophyletic groups defined by the most internal branch from the root above a bootstrap cutoff level (we chose 70%) [115, 116] that also shared identical taxon composition across all three phylogenetic methods and had fungal NRPS representation (Additional file 6). The SID group was a single exception in that in the full phylogenies (Fig. 1, Additional file 6) maximum likelihood methods supported this clade with 68% and 74% bootstrap support while NJ did not provide support above 50% (Fig. 1, Additional file 6). This clade is, however, supported by >80% bootstrap support in all phylogenetic methods in analysis of the reduced dataset (Fig. 2, Additional file 7).

Distribution of NRPS subfamilies across fungal taxonomic groups

To address patterns of distribution of NRPSs across fungal taxonomic groups, we tallied NRPS counts in Chytridiomycota, Zygomycota, Basidiomycota, Schizosaccharomycota, Hemiascomycota, and Euascomycota. Fisher's exact tests were used to test for associations between taxonomic groups and the proportion of genes in each NRPS subfamily.

Lineage specific expansions and variation in birth-death rates

We calculated and graphed the average and range of the number of genes encoding NRPSs in each subfamily per euascomycete genome and the number of A domains per NRPS for each subfamily to assess broad patterns of variation in numbers of genes and numbers of A domains/gene across subfamilies (Fig. 10)

We used the method of Hahn et al. [117, 118], which applies a stochastic birth and death process along a phylogeny to test for statistically significant lineage specific expansions and contractions of 1) number of NRPS genes and 2) numbers of NRPS A domains/subfamily. For these analyses, we created an ultrametric species tree with the PL method in r8s [119] using the phylogeny of the concatenated protein dataset of Fitzpatrick et al. [120] (Additional file 9).

We performed two separate analyses using CAFÉ [118] to look at patterns of gene and A domain expansions. The first analysis looked at patterns of the total number of NRPSs (e.g. all subfamilies combined) to look for broad patterns of expansions and contractions across the full tree of fungi (excluding B. dendrobatidis). The second analysis analyzed duplications and losses in each subfamily separately and was restricted to the euascomycete taxa because the birth-death model assumes that at least one gene of each subfamily is present in the common ancestor of all taxa. The ACV synthetase subfamily was excluded because parsimony inferred that this family had zero genes at the root. For all analyses, we used 1000 re-samplings and significant deviations from a random birth-death model were determined by viterbi p-values below .05.

All additional file figure legends and notes are in Additional file 15.

References

Finking R, Marahiel MA: Biosynthesis of nonribosomal peptides. Annual Review of Microbiology. 2004, 58: 453-488. 10.1146/annurev.micro.58.030603.123615.
CAS PubMed Google Scholar
Sieber SA, Marahiel MA: Learning from nature's drug factories: Nonribosomal synthesis of macrocyclic peptides. Journal of Bacteriology. 2003, 185 (24): 7036-7043. 10.1128/JB.185.24.7036-7043.2003.
PubMed Central CAS PubMed Google Scholar
Grunewald J, Marahiel MA: Chemoenzymatic and template-directed synthesis of bioactive macrocyclic peptides. Microbiology and Molecular Biology Reviews. 2006, 70 (1): 121-146. 10.1128/MMBR.70.1.121-146.2006.
PubMed Central PubMed Google Scholar
Stein T, Vater J, Kruft V, Otto A, Wittmann-Liebold B, Franke P, Panico M, McDowell R, Morris HR: The multiple carrier model of nonribosomal peptide biosynthesis at modular multienzymatic templates. Journal of Biological Chemistry. 1996, 271 (26): 15428-15435. 10.1074/jbc.271.26.15428.
CAS PubMed Google Scholar
Mootz HD, Schwarzer D, Marahiel MA: Ways of assembling complex natural products on modular nonribosomal peptide synthetases. ChemBioChem. 2002, 3: 490-504. 10.1002/1439-7633(20020603)3:6<490::AID-CBIC490>3.0.CO;2-N.
CAS PubMed Google Scholar
Oide S, Moeder W, Krasnoff S, Gibson D, Haas H, Yoshioka K, Turgeon BG: NPS6, encoding a nonribosomal peptide synthetase involved in siderophore-mediated iron metabolism, is a conserved virulence determinant of plant pathogenic ascomycetes. Plant Cell. 2006, 18 (10): 2836-2853. 10.1105/tpc.106.045633.
PubMed Central CAS PubMed Google Scholar
Oide S, Krasnoff SB, Gibson DM, Turgeon BG: Intracellular siderophores are essential for ascomycete sexual development in heterothallic Cochliobolus heterostrophus and homothallic Gibberella zeae. Eukaryotic Cell. 2007, 6 (8): 1339-1353. 10.1128/EC.00111-07.
PubMed Central CAS PubMed Google Scholar
Oide S: Functional characterization of nonribosomal peptide synthetases in the filamentous ascomycete phytopathogen Cochliobolus heterostrophus. PhD. 2007, Ithaca, NY: Cornell University
Google Scholar
Kim KH, Cho Y, La Rota M, Cramer RA, Lawrence CB: Functional analysis of the Alternaria brassicicola non-ribosomal peptide synthetase gene AbNPS2 reveals a role in conidial cell wall construction. Molecular Plant Pathology. 2007, 8 (1): 23-39. 10.1111/j.1364-3703.2006.00366.x.
CAS PubMed Google Scholar
Lee BN, Kroken S, Chou DYT, Robbertse B, Yoder OC, Turgeon BG: Functional analysis of all nonribosomal peptide synthetases in Cochliobolus heterostrophus reveals a factor, NPS6, involved in virulence and resistance to oxidative stress. Eukaryotic Cell. 2005, 4 (3): 545-555. 10.1128/EC.4.3.545-555.2005.
PubMed Central CAS PubMed Google Scholar
Hahn J, Dubnau D: Growth stage signal transduction and the requirements for Srfa induction in development of competence. Journal of Bacteriology. 1991, 173 (22): 7275-7282.
PubMed Central CAS PubMed Google Scholar
Schaeffer P: Sporulation and the production of antibiotics, exoenzymes, and exotonins. Bacteriological Reviews. 1969, 33 (1): 48-71.
PubMed Central CAS PubMed Google Scholar
Horinouchi S, Beppu T: Autoregulatory factors of secondary metabolism and morphogenesis in actinomycetes. Critical Reviews in Biotechnology. 1990, 10 (3): 191-204. 10.3109/07388559009038207.
CAS PubMed Google Scholar
Marahiel MA, Stachelhaus T, Mootz HD: Modular peptide synthetases involved in nonribosomal peptide synthesis. Chemical Reviews. 1997, 97 (7): 2651-2673. 10.1021/cr960029e.
CAS PubMed Google Scholar
Challis GL, Ravel J, Townsend CA: Predictive, structure-based model of amino acid recognition by nonribosomal peptide synthetase adenylation domains. Chemistry and Biology (London). 2000, 7 (3): 211-224.
CAS Google Scholar
Keating TA, Ehmann DE, Kohli RM, Marshall CG, Trauger JW, Walsh CT: Chain termination steps in nonribosomal peptide synthetase assembly lines: Directed acyl-S-enzyme breakdown in antibiotic and siderophore biosynthesis. ChemBioChem. 2001, 2 (2): 99-107. 10.1002/1439-7633(20010202)2:2<99::AID-CBIC99>3.0.CO;2-3.
CAS PubMed Google Scholar
Keating TA, Walsh CT: Initiation, elongation, and termination strategies in polyketide and polypeptide antibiotic biosynthesis. Current Opinion in Chemical Biology. 1999, 3 (5): 598-606. 10.1016/S1367-5931(99)00015-0.
CAS PubMed Google Scholar
Schneider A, Marahiel MA: Genetic evidence for a role of thioesterase domains, integrated in or associated with peptide synthetases, in non-ribosomal peptide biosynthesis in Bacillus subtilis. Archives of Microbiology. 1998, 169 (5): 404-410. 10.1007/s002030050590.
CAS PubMed Google Scholar
Kohli RM, Trauger JW, Schwarzer D, Marahiel MA, Walsh CT: Generality of peptide cyclization catalyzed by isolated thioesterase domains of nonribosomal peptide synthetases. Biochemistry. 2001, 40 (24): 7099-7108. 10.1021/bi010036j.
CAS PubMed Google Scholar
Pospiech A, Bietenhader J, Schupp T: Two multifunctional peptide synthetases and an O-methyltransferase are involved in the biosynthesis of the DNA-binding antibiotic and antitumour agent saframycin Mx1 from Myxococcus xanthus. Microbiology. 1996, 142: 741-746. 10.1099/00221287-142-4-741.
CAS PubMed Google Scholar
Silakowski B, Kunze B, Nordsiek G, Blocker H, Hofle G, Muller R: The myxochelin iron transport regulon of the myxobacterium Stigmatella aurantiaca Sg a15. European Journal of Biochemistry. 2000, 267 (21): 6476-6485. 10.1046/j.1432-1327.2000.01740.x.
CAS PubMed Google Scholar
Silakowski B, Nordsiek G, Kunze B, Blocker H, Muller R: Novel features in a combined polyketide synthase/non-ribosomal peptide synthetase: the myxalamid biosynthetic gene cluster of the myxobacterium Stigmatella aurantiaca Sga15. Chemistry & Biology. 2001, 8 (1): 59-69. 10.1016/S1074-5521(00)00056-9.
CAS Google Scholar
Ehmann DE, Gehring AM, Walsh CT: Lysine biosynthesis in Saccharomyces cerevisiae: Mechanism of alpha-aminoadipate reductase (Lys2) involves posttranslational phosphopantetheinylation by Lys5. Biochemistry. 1999, 38 (19): 6171-6177. 10.1021/bi9829940.
CAS PubMed Google Scholar
Walzel B, Riederer B, Keller U: Mechanism of alkaloid cyclopeptide synthesis in the ergot fungus Claviceps purpurea. Chemistry & Biology. 1997, 4 (3): 223-230. 10.1016/S1074-5521(97)90292-1.
CAS Google Scholar
Pfeifer E, Pavelavrancic M, Vondohren H, Kleinkauf H: Characterization of tyrocidine synthetase 1 (TY1): Requirement of posttranslational modification for peptide biosynthesis. Biochemistry. 1995, 34 (22): 7450-7459. 10.1021/bi00022a019.
CAS PubMed Google Scholar
Rausch C, Hoof I, Weber T, Wohlleben W, Huson DH: Phylogenetic analysis of condensation domains in NRPS sheds light on their functional evolution. BMC Evolutionary Biology. 2007, 7: 78-10.1186/1471-2148-7-78.
PubMed Central PubMed Google Scholar
Walsh CT, Chen HW, Keating TA, Hubbard BK, Losey HC, Luo LS, Marshall CG, Miller DA, Patel HM: Tailoring enzymes that modify nonribosomal peptides during and after chain elongation on NRPS assembly lines. Current Opinion in Chemical Biology. 2001, 5 (5): 525-534. 10.1016/S1367-5931(00)00235-0.
CAS PubMed Google Scholar
Samel SA, Marahiel MA, Essen LO: How to tailor non-ribosomal peptide products - new clues about the structures and mechanisms of modifying enzymes. Molecular Biosystems. 2008, 4 (5): 387-393. 10.1039/b717538h.
CAS PubMed Google Scholar
Kroken S, Glass NL, Taylor JW, Yoder OC, Turgeon BG: Phylogenomic analysis of type I polyketide synthase genes in pathogenic and saprobic ascomycetes. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (26): 15670-15675. 10.1073/pnas.2532165100.
PubMed Central CAS PubMed Google Scholar
Collemare J, Pianfetti M, Houlle AE, Morin D, Camborde L, Gagey MJ, Barbisan C, Fudal I, Lebrun MH, Boehnert HU: Magnaporthe grisea avirulence gene ACE1 belongs to an infection-specific gene cluster involved in secondary metabolism. New Phytologist. 2008, 179 (1): 196-208. 10.1111/j.1469-8137.2008.02459.x.
CAS PubMed Google Scholar
Maiya S, Grundmann A, Li X, Li SM, Turner G: Identification of a hybrid PKS/NRPS required for pseurotin A biosynthesis in the human pathogen Aspergillus fumigatus. ChemBioChem. 2007, 8 (14): 1736-1743. 10.1002/cbic.200700202.
CAS PubMed Google Scholar
Bergmann S, Schumann J, Scherlach K, Lange C, Brakhage AA, Hertweck C: Genomics-driven discovery of PKS-NRPS hybrid metabolites from Aspergillus nidulans. Nature Chemical Biology. 2007, 3 (4): 213-217. 10.1038/nchembio869.
CAS PubMed Google Scholar
Brendel N, Partida-Martinez LP, Scherlach K, Hertweck C: A cryptic PKS-NRPS gene locus in the plant commensal Pseudomonas fluorescens Pf-5 codes for the biosynthesis of an antimitotic rhizoxin complex. Organic & Biomolecular Chemistry. 2007, 5 (14): 2211-2213. 10.1039/b707762a.
CAS Google Scholar
Rees DO, Bushby N, Cox RJ, Harding JR, Simpson TJ, Willis CL: Synthesis of [1,2-C-13(2), N-15]-L-homoserine and its incorporation by the PKS-NRPS system of Fusarium moniliforme into the mycotoxin fusarin C. ChemBioChem. 2007, 8 (1): 46-50. 10.1002/cbic.200600404.
CAS PubMed Google Scholar
Fedorova ND, Khaldi N, Joardar VS, Maiti R, Amedeo P, Anderson MJ, Crabtree J, Silva JC, Badger JH, Albarraq A, et al: Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. PLoS Genetics. 2008, 4 (4): e1000046-10.1371/journal.pgen.1000046.
PubMed Central PubMed Google Scholar
Chang PK, Horn BW, Dorner JW: Sequence breakpoints in the aflatoxin biosynthesis gene cluster and flanking regions in nonaflatoxigenic Aspergillus flavus isolates. Fungal Genetics and Biology. 2005, 42 (11): 914-923. 10.1016/j.fgb.2005.07.004.
CAS PubMed Google Scholar
Cramer RA, Stajich JE, Yamanaka Y, Dietrich FS, Steinbach William JS, Perfect JR: Phylogenomic analysis of non-ribosomal peptide synthetases in the genus Aspergillus. Gene. 2006, 383 (15): 24-32. 10.1016/j.gene.2006.07.008.
CAS PubMed Google Scholar
Johnson R, Voisey C, Johnson L, Pratt J, Fleetwood D, Khan A, Bryan G: Distribution of NRPS gene families within the Neotyphodium/Epichloe complex. Fungal Genetics and Biology. 2007, 44 (11): 1180-1190. 10.1016/j.fgb.2007.04.009.
CAS PubMed Google Scholar
Nei M, Hughes AL: Balanced polymorphism and evolution by the birth-and-death processin the MHC loci. 11th Histocompatibility Workshop and Conference: 1992. 1992, Oxford Univ. Press
Google Scholar
de Bono B, Madera M, Chothia C: V-H gene segments in the mouse and human genomes. Journal of Molecular Biology. 2004, 342 (1): 131-143. 10.1016/j.jmb.2004.06.055.
CAS PubMed Google Scholar
Hamilton AT, Huntley S, Tran-Gyamfi M, Baggott DM, Gordon L, Stubbs L: Evolutionary expansion and divergence in the ZNF91 subfamily of primate-specific zinc finger genes. Genome Research. 2006, 16 (5): 584-594. 10.1101/gr.4843906.
PubMed Central CAS PubMed Google Scholar
Tian X, Pascal G, Fouchecourt S, Pontarotti P, Monget P: Gene birth, death, and divergence: The different scenarios of reproduction-related gene evolution. Biology of Reproduction. 2009, 80 (4): 616-621. 10.1095/biolreprod.108.073684.
CAS PubMed Google Scholar
Nozawa M, Nei M: Evolutionary dynamics of olfactory receptor genes in Drosophila species. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (17): 7122-7127. 10.1073/pnas.0702133104.
PubMed Central CAS PubMed Google Scholar
Nozawa M, Kawahara Y, Nei M: Genomic drift and copy number variation of sensory receptor genes in humans. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (51): 20421-20426. 10.1073/pnas.0709956104.
PubMed Central CAS PubMed Google Scholar
Niimura Y: Evolutionary dynamics of olfactory receptor genes in mammals. Genes & Genetic Systems. 2007, 82 (6): 503-503.
Google Scholar
Niimura Y, Nei M: Evolutionary dynamics of olfactory and other chemosensory receptor genes in vertebrates. Journal of Human Genetics. 2006, 51 (6): 505-517. 10.1007/s10038-006-0391-8.
PubMed Central CAS PubMed Google Scholar
Niimura Y, Nei M: Comparative evolutionary analysis of olfactory receptor gene clusters between humans and mice. Gene. 2005, 346: 13-21. 10.1016/j.gene.2004.09.025.
CAS PubMed Google Scholar
Nam J, Kim J, Lee S, An GH, Ma H, Nei MS: Type I MADS-box genes have experienced faster birth-and-death evolution than type II MADS-box genes in angiosperms. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101 (7): 1910-1915. 10.1073/pnas.0308430100.
PubMed Central CAS PubMed Google Scholar
Xu GX, Ma H, Nei M, Kong HZ: Evolution of F-box genes in plants: Different modes of sequence divergence and their relationships with functional diversification. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (3): 835-840. 10.1073/pnas.0812043106.
PubMed Central CAS PubMed Google Scholar
Turgeon BG, Oide S, Bushley K: Creating and screening Cochliobolus heterostrophus non-ribosomal peptide synthetase mutants. Mycological Research. 2008, 112: 200-206. 10.1016/j.mycres.2007.10.012.
CAS PubMed Google Scholar
Scottcraig JS, Panaccione DG, Pocard JA, Walton JD: The cyclic peptide synthetase catalyzing HC-toxin production in the filamentous fungus Cochliobolus carbonum is encoded by a 15.7-kilobase open reading frame. Journal of Biological Chemistry. 1992, 267 (36): 26044-26049.
CAS Google Scholar
Johnson RD, Johnson L, Itoh Y, Kodama M, Otani H, Kahmoto K: Cloning and characterization of a cyclic peptide synthetase gene from Alternaria alternata apple pathotype whose product is involved in AM-toxin synthesis and pathogenicity. Molecular Plant-Microbe Interactions. 2000, 13 (7): 742-753. 10.1094/MPMI.2000.13.7.742.
CAS PubMed Google Scholar
Wapinski I, Pfeffer A, Friedman N, Regev A: Natural history and evolutionary principles of gene duplication in fungi. Nature. 2007, 449 (7158): 54-U36. 10.1038/nature06107.
CAS PubMed Google Scholar
Lu SW, Kroken S, Lee BN, Robbertse B, Churchill ACL, Yoder OC, Turgeon BG: A novel class of gene controlling virulence in plant pathogenic ascomycete fungi. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (10): 5980-5985. 10.1073/pnas.0931375100.
PubMed Central CAS PubMed Google Scholar
Suvarna K, Seah L, Bhattacherjee V, Bhattacharjee JK: Molecular analysis of the LYS2 gene of Candida albicans: homology to peptide antibiotic synthetases and the regulation of the alpha-aminoadipate reductase. Current Genetics. 1998, 33 (4): 268-275. 10.1007/s002940050336.
CAS PubMed Google Scholar
Eibel H, Philippsen P: Identification of the Cloned Saccharomyces cerevisiae LYS2 Gene by an Integrative Transformation Approach. Molecular & General Genetics. 1983, 191 (1): 66-73. 10.1007/BF00330891.
CAS Google Scholar
Sinha AK, Bhattach JK: Lysine biosynthesis in Saccharomyces - Conversion of alpha-aminoadipate into alpha-aminoadipic delta-semialdehyde. Biochemical Journal. 1971, 125 (3): 743-
PubMed Central CAS PubMed Google Scholar
Weber G, Schorgendorfer K, Schneiderscherzer E, Leitner E: the peptide synthetase catalyzing Cyclosporine production in Tolypocladium niveum is encoded by a giant 45.8-kilobase open reading frame. Current Genetics. 1994, 26 (2): 120-125. 10.1007/BF00313798.
CAS PubMed Google Scholar
Aharonowitz Y, Cohen G, Martin JF: Penicillin and Cephalosporin biosynthetic genes - structure, organization, regulation, and evolution. Annual Review of Microbiology. 1992, 46: 461-495. 10.1146/annurev.mi.46.100192.002333.
CAS PubMed Google Scholar
Brakhage AA, Al-Abdallah Q, Tuncher A, Sprote P: Evolution of beta-lactam biosynthesis genes and recruitment of trans-acting factors. Phytochemistry. 2005, 66 (11): 1200-1210. 10.1016/j.phytochem.2005.02.030.
CAS PubMed Google Scholar
Liras P, Martin JF: Gene clusters for beta-lactam antibiotics and control of their expression: why have clusters evolved, and from where did they originate?. International Microbiology. 2006, 9 (1): 9-19.
CAS PubMed Google Scholar
Buades C, Moya A: Phylogenetic analysis of the isopenicillin-N-synthetase horizontal gene transfer. Journal of Molecular Evolution. 1996, 42 (5): 537-542. 10.1007/BF02352283.
CAS PubMed Google Scholar
Landan G, Cohen G, Aharonowitz Y, Shuali Y, Graur D, Shiffman D: Evolution of Isopenicillin-N synthase genes may have involved horizontal gene-transfer. Molecular Biology and Evolution. 1990, 7 (5): 399-406.
CAS PubMed Google Scholar
Penalva MA, Moya A, Dopazo J, Ramon D: Sequences of Isopenicillin-N synthetase genes suggest horizontal gene-transfer from prokaryotes to eukaryotes. Proceedings of the Royal Society of London Series B-Biological Sciences. 1990, 241 (1302): 164-169. 10.1098/rspb.1990.0081.
CAS Google Scholar
Bushley KE, Ripoll DR, Turgeon BG: Module evolution and substrate specificity of fungal nonribosomal peptide synthetases involved in siderophore biosynthesis. BMC Evolutionary Biology. 2008, 8: 328-10.1186/1471-2148-8-328.
PubMed Central PubMed Google Scholar
von Dohren H: Biochemistry and general genetics of nonribosomal pepticle synthetases in fungi. Molecular Biotechnology of Fungal Beta-Lactam Antibiotics and Related Peptide Synthetases. 2004, 88: 217-264.
Google Scholar
Selker EU: Genome defense and DNA methylation in Neurospora. Cold Spring Harbor Symposia on Quantitative Biology. 2004, 69: 119-124. 10.1101/sqb.2004.69.119.
CAS PubMed Google Scholar
Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, et al: The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003, 422 (6934): 859-868. 10.1038/nature01554.
CAS PubMed Google Scholar
Velasco AM, Leguina JI, Lazcano A: Molecular evolution of the lysine biosynthetic pathways. Journal of Molecular Evolution. 2002, 55 (4): 445-459. 10.1007/s00239-002-2340-2.
CAS PubMed Google Scholar
Katinka MD, Duprat S, Cornillot E, Metenier G, Thomarat F, Prensier G, Barbe V, Peyretaillade E, Brottier P, Wincker P, et al: Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi. Nature. 2001, 414 (6862): 450-453. 10.1038/35106579.
CAS PubMed Google Scholar
Bohnert HU, Fudal I, Dioh W, Tharreau D, Notteghem JL, Lebrun MH: A putative polyketide synthase peptide synthetase from Magnaporthe grisea signals pathogen attack to resistant rice. Plant Cell. 2004, 16 (9): 2499-2513. 10.1105/tpc.104.022715.
PubMed Central PubMed Google Scholar
Nishida H, Nishiyama M, Kobashi N, Kosuge T, Hoshino T, Yamane H: A prokaryotic gene cluster involved in synthesis of lysine through the amino adipate pathway: a key to the evolution of amino acid biosynthesis. Genome Research. 1999, 9: 1175-1183. 10.1101/gr.9.12.1175.
CAS PubMed Google Scholar
Xu HY, Andi B, Qian JH, West AH, Cook PF: The alpha-aminoadipate pathway for lysine biosynthesis in fungi. Cell Biochemistry and Biophysics. 2006, 46 (1): 43-64. 10.1385/CBB:46:1:43.
CAS PubMed Google Scholar
Zabriskie TM, Jackson MD: Lysine biosynthesis and metabolism in fungi. Natural Product Reports. 2000, 17 (1): 85-97. 10.1039/a801345d.
CAS PubMed Google Scholar
Kwang-Deuk A, Nishida H, Yoshiharu M, Yokota A: Aminoadipate reductase gene: a new fungal-specific gene for comparative evolutionary analyses. BMC Evolutionary Biology. 2002, 2: 6-10.1186/1471-2148-2-6.
Google Scholar
Sims JW, Schmidt EW: Thioesterase-like role for fungal PKS-NRPS hybrid reductive domains. Journal of the American Chemical Society. 2008, 130 (33): 11149-11155. 10.1021/ja803078z.
CAS PubMed Google Scholar
Sims JW, Fillmore JP, Warner DD, Schmidt EW: Equisetin biosynthesis in Fusarium heterosporum. Chemical Communications. 2005, 186-188. 10.1039/b413523g. 2
Song ZS, Cox RJ, Lazarus CM, Simpson TJ: Fusarin C biosynthesis in Fusarium moniliforme and Fusarium venenatum. ChemBioChem. 2004, 5 (9): 1196-1203. 10.1002/cbic.200400138.
CAS PubMed Google Scholar
Chattopadhyay D, Finzel BC, Munson SH, Evans DB, Sharma SK, Strakalaitis NA, Brunner DP, Eckenrode FM, Dauter Z, Betzel C, et al: Crystallographic analyses of an active HIV-1 ribonuclease H domain show structural features that distinguish it from the inactive form. Acta Crystallographica Section D-Biological Crystallography. 1993, 49: 423-427. 10.1107/S0907444993002409.
CAS Google Scholar
Rice P, Craigie R, Davies DR: Retroviral integrases and their cousins. Current Opinion in Structural Biology. 1996, 6 (1): 76-83. 10.1016/S0959-440X(96)80098-4.
CAS PubMed Google Scholar
Declercq E, Billiau A, Ottenheijm HCJ, Herscheid JDM: Anti-reverse transcriptase activity of Gliotoxin analogs. Biochemical Pharmacology. 1978, 27 (5): 635-639. 10.1016/0006-2952(78)90497-5.
CAS Google Scholar
Rouxel T, Chupeau Y, Fritz R, Kollmann A, Bousquet JF: Biological effects of Sirodesmin-Pl, a phytotoxin produced by Leptosphaeria maculans. Plant Science. 1988, 57 (1): 45-53. 10.1016/0168-9452(88)90140-9.
CAS Google Scholar
Rouxel T, Kollmann A, Bousquet JF: Zinc suppresses sirodesmin PL toxicity and protects Brassica napus plants against the blackleg disease caused by Leptosphaeria maculans. Plant Science. 1990, 68 (1): 77-86. 10.1016/0168-9452(90)90155-H.
CAS Google Scholar
Hoffmeister D, Keller NP: Natural products of filamentous fungi: enzymes, genes, and their regulation. Natural Product Reports. 2007, 24 (2): 393-416. 10.1039/b603084j.
CAS PubMed Google Scholar
Wiest A, Grzegorski D, Xu BW, Goulard C, Rebuffat S, Ebbole DJ, Bodo B, Kenerley C: Identification of peptaibols from Trichoderma virens and cloning of a peptaibol synthetase. Journal of Biological Chemistry. 2002, 277 (23): 20862-20868. 10.1074/jbc.M201654200.
CAS PubMed Google Scholar
Tanaka A, Tapper BA, Popay A, Parker EJ, Scott B: A symbiosis expressed non-ribosomal peptide synthetase from a mutualistic fungal endophyte of perennial ryegrass confers protection to the symbiotum from insect herbivory. Molecular Microbiology. 2005, 57 (4): 1036-1050. 10.1111/j.1365-2958.2005.04747.x.
CAS PubMed Google Scholar
Spatafora JW, Sung GH, Sung JM, Hywel-Jones NL, White JF: Phylogenetic evidence for an animal pathogen origin of ergot and the grass endophytes. Molecular Ecology. 2007, 16 (8): 1701-1711. 10.1111/j.1365-294X.2007.03225.x.
CAS PubMed Google Scholar
Clay K, Cheplick GP: Effect of ergot alkaloids from fungal endophyte-infected grasses on fall armyworm (Spodoptera frugiperda). Journal of Chemical Ecology. 1989, 15 (1): 169-182. 10.1007/BF02027781.
CAS PubMed Google Scholar
Fiserova A, Pospisil M: Role of ergot alkaloids in the immune system. Ergot-The Genus Claviceps. Edited by: Kren V, Cvak L. 1999, Amsterdam: Harwood, 451-467.
Google Scholar
Panaccione DG, Cipoletti JR, Sedlock AB, Blemings KP, Schardl CL, Machado C, Seidel GE: Effects of ergot alkaloids on food preference and satiety in rabbits, as assessed with gene-knockout endophytes in perennial ryegrass (Lolium perenne). Journal of Agricultural and Food Chemistry. 2006, 54 (13): 4582-4587. 10.1021/jf060626u.
CAS PubMed Google Scholar
Cross D: Ergot Alkaloid Toxicity. Clavicipitalean Fungi: Evolutionary Biology, Chemistry, Biocontrol, and Cultural Impacts. Edited by: White JF Jr, Bacon CW, Hywel-Jones NL, Spatafora JW. 2003, New York, New York: Marcel Dekker, Inc, 475-494.
Google Scholar
Wei XY, Yang FQ, Straney DC: Multiple non-ribosomal peptide synthetase genes determine peptaibol synthesis in Trichoderma virens. Canadian Journal of Microbiology. 2005, 51 (5): 423-429. 10.1139/w05-006.
CAS PubMed Google Scholar
Hughes AL, Nei M: Evolution of the major histocompatibility complex - independent origin of nonclassical Class I genes in different groups of mammals. Molecular Biology and Evolution. 1989, 6 (6): 559-579.
CAS PubMed Google Scholar
Nei M, Gu X, Sitnikova T: Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94 (15): 7799-7806. 10.1073/pnas.94.15.7799.
PubMed Central CAS PubMed Google Scholar
Ota T, Nei M: Divergent evolution and evolution by the birth-and-death process in the Immunoglobulin V-H gene family. Molecular Biology and Evolution. 1994, 11 (3): 469-482.
CAS PubMed Google Scholar
Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annual Review of Genetics. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.
PubMed Central CAS PubMed Google Scholar
Nei M: The new mutation theory of phenotypic evolution. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (30): 12235-12242. 10.1073/pnas.0703349104.
PubMed Central CAS PubMed Google Scholar
Nei M, Niimura Y, Nozawa M: The evolution of animal chemosensory receptor gene repertoires: roles of chance and necessity. Nature Reviews Genetics. 2008, 9 (12): 951-963. 10.1038/nrg2480.
CAS PubMed Google Scholar
Korbel JO, Kim PM, Chen X, Urban AE, Weissman S, Snyder M, Gerstein MB: The current excitement about copy-number variation: how it relates to gene duplications and protein families. Current Opinion in Structural Biology. 2008, 18 (3): 366-374. 10.1016/j.sbi.2008.02.005.
PubMed Central CAS PubMed Google Scholar
Lautru S, Challis GL: Substrate recognition by nonribosomal peptide synthetase multi-enzymes. Microbiology (Reading). 2004, 150 (Part 6): 1629-1636. 10.1099/mic.0.26837-0.
CAS Google Scholar
Stachelhaus T, Mootz HD, Marahiel MA: The specificity-conferring code of adenylation domains in nonribosomal peptide synthetases. Chemistry and Biology (London). 1999, 6 (8): 493-505.
CAS Google Scholar
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
PubMed Central CAS PubMed Google Scholar
Birney E, Clamp M, Durbin R: GeneWise and Genomewise. Genome Research. 2004, 14: 988-995. 10.1101/gr.1865504.
PubMed Central CAS PubMed Google Scholar
Karolewiez A, Geisen R: Cloning a part of the ochratoxin A biosynthetic gene cluster of Penicillium nordicum and characterization of the ochratoxin polyketide synthase gene. Systematic and Applied Microbiology. 2005, 28 (7): 588-595. 10.1016/j.syapm.2005.03.008.
CAS PubMed Google Scholar
Rausch C, Weber T, Kohlbacher O, Wohlleben W, Huson DH: Specificity predictions of adenylation domains in nonribosomal peptide synthetases (NRPS) using transductive support vector machines (TSVMs). Nucleic Acids Research. 2005, 33 (18): 5799-5808. 10.1093/nar/gki885.
PubMed Central CAS PubMed Google Scholar
Abascal F, Zardoya R, Posada D: ProtTest: Selection of best-fit models of protein evolution. Bioinformatics. 2005, 21 (9): 2104-2105. 10.1093/bioinformatics/bti263.
CAS PubMed Google Scholar
Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McInerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evolutionary Biology. 2006, 6: 29-10.1186/1471-2148-6-29.
PubMed Central PubMed Google Scholar
Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML web-servers. Systematic Biology. 2008, 75 (5): 758-771. 10.1080/10635150802429642.
Google Scholar
Guindon SGO: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
PubMed Google Scholar
Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.6. 2005, Distributed by the author. Department of Genome Sciences, University of Washington, Seattle
Google Scholar
Schmidt HA, Strimmer K, Vingron M, Haeseler AV: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.
CAS PubMed Google Scholar
Bremmer B, Jansen R, Oxelman B, Backlund M, Lantz H, KJ K: More characters or more taxa for a robust phylogeny: case study from the coffee family (Rubiaceae). Systematic Biology. 1999, 48: 413-435. 10.1080/106351599260085.
Google Scholar
Mitchell A, Mitter C, Regier JC: More taxa or more characters revisited: Combining data from nuclear protein-encoding genes for phylogenetic analyses of Noctuoidea (Insecta: Lepidoptera). Systematic Biology. 2000, 49 (2): 202-224. 10.1093/sysbio/49.2.202.
CAS PubMed Google Scholar
Sanderson MJ, Wojciechowski MF: Improved bootstrap confidence limits in large-scale phylogenies, with an example from Neo astragalus (leguminosae). Systematic Biology. 2000, 49 (4): 671-685. 10.1080/106351500750049761.
CAS PubMed Google Scholar
Auwerx J, Baulieu E, Beato M, Becker-Andre M, Burbach PH, Camerino G, Chambon P, Cooney A, Dejean A, Dreyer C, et al: A unified nomenclature system for the nuclear receptor superfamily. Cell. 1999, 97 (2): 161-163. 10.1016/S0092-8674(00)80726-6.
Google Scholar
Hodge T, Cope M: A myosin family tree. Journal of Cell Science. 2000, 113 (19): 3353-
CAS PubMed Google Scholar
Hahn MW, De Bie T, Stajich JE, Nguyen C, Cristianini N: Estimating the tempo and mode of gene family evolution from comparative genomic data. Genome Research. 2005, 15 (8): 1153-1160. 10.1101/gr.3567505.
PubMed Central CAS PubMed Google Scholar
De Bie T, Cristianini N, Demuth JP, Hahn MW: CAFE: a computational tool for the study of gene family evolution. Bioinformatics. 2006, 22 (10): 1269-1271. 10.1093/bioinformatics/btl097.
CAS PubMed Google Scholar
Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19 (2): 301-302. 10.1093/bioinformatics/19.2.301.
CAS PubMed Google Scholar
Fitzpatrick DA, Logue ME, Stajich JE, Butler G: A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis. BMC Evolutionary Biology. 2006, 6: 99-10.1186/1471-2148-6-99.
PubMed Central PubMed Google Scholar
Page RDM: TreeView: An application to display phylogenetic trees on personal computers. Computer Applications in the Biosciences. 1996, 12 (4): 357-358.
CAS PubMed Google Scholar
Jeanmougin F, Thompson JD, Gouy M, Higgins DG, Gibson TJ: Multiple sequence alignment with Clustal X. Trends in Biochemical Sciences. 1998, 23 (10): 403-405. 10.1016/S0968-0004(98)01285-7.
CAS PubMed Google Scholar

Download references

Acknowledgements

BGT acknowledges, with gratitude, the US Department of Energy Joint Genome Institute (JGI) for their fungal genome program, in particular, for their support in generating the sequence of race O, strain C5, of Cochliobolus heterostrophus http://genome.jgi-psf.org/CocheC5_1/CocheC5_1.home.html. BGT and KEB are especially grateful to D. Schneider (USDA ARS) for helpful discussions and for providing computer resources for phylogenetic and other analyses. KEB would like to thank J. Doyle, S. Kroken, and A. Siepel for discussions of phylogenetic analyses, J. Stajich and M. Hahn for discussions of ultrametric tree construction and birth-death analyses, and the Cornell Computational Biology Service Unit facility and staff, A. Siepel, K. Nixon, and the CIPRES project http://www.phylo.org/sub_sections/portal/ for computer resources for phylogenetic and other computational analyses. BGT acknowledges the support of the Division of Molecular and Cellular Biosciences, National Science Foundation, the USDA Cooperative State Research Education and Extension Service, National Research Initiative and the BARD foundation.

Author information

Authors and Affiliations

Department of Plant Pathology & Plant-Microbe Biology, Cornell University, 334 Plant Science Bldg., Ithaca, NY, 14853, USA
Kathryn E Bushley & B Gillian Turgeon

Authors

Kathryn E Bushley
View author publications
You can also search for this author in PubMed Google Scholar
B Gillian Turgeon
View author publications
You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to B Gillian Turgeon.

Additional information

Authors' contributions

BGT and KEB conceived of the study, and participated in its design, coordination and data interpretation. KEB carried out the protein and domain identifications, performed alignments, phylogenetic analyses, and drafted the manuscript. BGT advised KEB on manuscript content and in critical revisions. KEB and BGT jointly wrote the final versions of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12862_2009_1242_MOESM1_ESM.PDF

Additional file 1: Diagram of Cochliobolus heterostrophus NRPSs and their domain structure. 30 individual AMP domains are indicated. See Additional file 15 for detailed description. (PDF 11 KB)

12862_2009_1242_MOESM2_ESM.PDF

Additional file 2: Annotation of all proteins used in the study. Identification, accession numbers, genomic locations, and domain architectures of NRPSs identified in 38 fungal genomes. (PDF 216 KB)

12862_2009_1242_MOESM3_ESM.PDF

Additional file 3: Genus and species abbreviations. Genus and species abbreviations for all taxa used in the study. (PDF 89 KB)

12862_2009_1242_MOESM4_ESM.ZIP

Additional file 4: Profile HMMs for fungal-specific NRPS A, T, and C domains. Zipped text file (file name extension .hmm to be used with the program package HMMER (http://hmmer.janelia.org). See Additional file 15 for detailed description. (ZIP 62 KB)

12862_2009_1242_MOESM5_ESM.PDF

Additional file 5: Fungal and bacterial AMP-binding protein outgroups. Selection of fungal and bacterial AMP-binding protein used as outgroups in phylogenetic analyses. (PDF 128 KB)

12862_2009_1242_MOESM6_ESM.PDF

Additional file 6: Phylogenies, full dataset. A. NJ, B. ML (PhyML), and C. ML (RAxML) phylogenies of the full AMP dataset. See Additional file 15 for detailed description. (PDF 167 KB)

12862_2009_1242_MOESM7_ESM.PDF

Additional file 7: Phylogenies, reduced dataset. Topologies of NJ, ML (PhyML), and ML (RAxML) phylogenetic analyses of the reduced NRPS AMP domain dataset containing selected representatives of each subfamily. See Additional file 15 for detailed description. (PDF 224 KB)

Additional file 8: Bacterial outgroup proteins. Bacterial proteins used as outgroups. (PDF 145 KB)

12862_2009_1242_MOESM9_ESM.PDF

Additional file 9: Newick phylogenetic tree. Text file containing newick phylogenetic tree for opening in tree visualization programs such as Treeview [121]. See Additional file 15 for detailed description. (PDF 81 KB)

12862_2009_1242_MOESM10_ESM.PDF

Additional file 10: Known fungal NRPSs used for constructing the initial HMMER model. Known fungal NRPS AMP domains used for constructing the initial HMMER model. (PDF 152 KB)

12862_2009_1242_MOESM11_ESM.HMM

Additional file 11: HMMER AMP domain models used as the initial model for NRPS identification. HMMER AMP domain models used as the initial model for NRPS identification. Zipped text files (file name extension .hmm to be used with the program package HMMER http://hmmer.janelia.org. (HMM 189 KB)

12862_2009_1242_MOESM12_ESM.PDF

Additional file 12: Fungal protein datasets used in phylogenomic analyses. Fungal protein (AMP domain) datasets used in phylogenomic analyses. (PDF 142 KB)

12862_2009_1242_MOESM13_ESM.ZIP

Additional file 13: MUSCLE alignment for complete AMP domain dataset. MUSCLE alignment of 558 fungal and bacterial AMP domains used in phylogenetic analyses of the complete dataset. Zipped text file containing alignment in fasta format for visualization in sequence alignment editor such as ClustalX [122]. See Additional file 15 for detailed description. (ZIP 123 KB)

12862_2009_1242_MOESM14_ESM.ZIP

Additional file 14: MUSCLE alignment for reduced AMP domain dataset. MUSCLE alignment of the reduced dataset of fungal and bacterial AMP domains containing selected representatives of each major fungal subfamily and bacterial clades. Zipped text file containing alignment in fasta format for visualization in sequence alignment editor such as ClustalX [122]. See Additional file 15 for detailed description. (ZIP 70 KB)

12862_2009_1242_MOESM15_ESM.PDF

Additional file 15: Additional file figure legends and notes. Selected additional file figure legends and notes. (PDF 57 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Authors’ original file for figure 9

Authors’ original file for figure 10

Rights and permissions

This article is published under license to 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.

Reprints and permissions

About this article

Cite this article

Bushley, K.E., Turgeon, B.G. Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships. BMC Evol Biol 10, 26 (2010). https://doi.org/10.1186/1471-2148-10-26

Download citation

Received: 06 August 2009
Accepted: 26 January 2010
Published: 26 January 2010
DOI: https://doi.org/10.1186/1471-2148-10-26

Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships

Abstract

Background

Results

Conclusions

Background

Results and Discussion

Identification and domain structure of candidate NRPSs

Phylogenomic analysis and subfamily identification

Relationships between fungal and bacterial NRPSs: horizontal transfer or vertical transmission and massive loss?

Distribution of NRPS subfamilies across fungal taxonomic groups

Lineage specific expansions and contractions

Subfamily distribution

AAR

ChNPS11/ETP/ChNPS12

ChNPS10, CYCLO, SID

ACV and PKS;NRPS

Euascomycete only (EAS)

Hypothesized origins based on taxonomic distribution

Mono- and bi- modular NRPS subfamilies

AARs and Lysine biosynthesis

PKS;NRPS

ChNPS10

ChNPS11/ETP/ChNPS12

Multimodular NRPS subfamilies

Diversity within the EAS subfamily

ChNPS6/PerA

Ergot alkaloid synthetases

Peptaibol synthetases

Dothideomycete host-selective toxin synthetases

ChNPS4

Evolutionary mechanisms giving rise to multimodular NRPSs

Tandem Duplication

Recombination

Stability of NRPS gene copy number and domain architectures across subfamilies

Chain termination mechanisms

Conclusions

Methods

Identification of putative NRPSs in fungal genomes

Annotation of domain architectures

Phylogenomic analyses

Subfamily identification and modelling

Distribution of NRPS subfamilies across fungal taxonomic groups

Lineage specific expansions and variation in birth-death rates

References

Acknowledgements

Author information

Authors and Affiliations

Corresponding author

Additional information

Authors' contributions

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

About this article

Cite this article

Share this article

Keywords

BMC Ecology and Evolution

Contact us