Skip to main content

Genus level analysis of PKS-NRPS and NRPS-PKS hybrids reveals their origin in Aspergilli

Abstract

Background

Filamentous fungi produce a vast amount of bioactive secondary metabolites (SMs) synthesized by e.g. hybrid polyketide synthase-nonribosomal peptide synthetase enzymes (PKS-NRPS; NRPS-PKS). While their domain structure suggests a common ancestor with other SM proteins, their evolutionary origin and dynamics in fungi are still unclear. Recent rational engineering approaches highlighted the possibility to reassemble hybrids into chimeras — suggesting molecular recombination as diversifying mechanism.

Results

Phylogenetic analysis of hybrids in 37 species – spanning 9 sections of Aspergillus and Penicillium chrysogenum – let us describe their dynamics throughout the genus Aspergillus. The tree topology indicates that three groups of PKS-NRPS as well as one group of NRPS-PKS hybrids developed independently from each other. Comparison to other SM genes lead to the conclusion that hybrids in Aspergilli have several PKS ancestors; in contrast, hybrids are monophyletic when compared to available NRPS genes — with the exception of a small group of NRPSs. Our analysis also revealed that certain NRPS-likes are derived from NRPSs, suggesting that the NRPS/NRPS-like relationship is dynamic and proteins can diverge from one function to another. An extended phylogenetic analysis including bacterial and fungal taxa revealed multiple ancestors of hybrids. Homologous hybrids are present in all sections which suggests frequent horizontal gene transfer between genera and a finite number of hybrids in fungi.

Conclusion

Phylogenetic distances between hybrids provide us with evidence for their evolution: Large inter-group distances indicate multiple independent events leading to the generation of hybrids, while short intra-group distances of hybrids from different taxonomic sections indicate frequent horizontal gene transfer. Our results are further supported by adding bacterial and fungal genera. Presence of related hybrid genes in all Ascomycetes suggests a frequent horizontal gene transfer between genera and a finite diversity of hybrids — also explaining their scarcity. The provided insights into relations of hybrids and other SM genes will serve in rational design of new hybrid enzymes.

Background

Secondary metabolites (SMs), non-growth associated compounds, have been subject to research efforts due to their wide range of bioactivities. Polyketides like sterigmatocystin and aflatoxin, two potent mycotoxins [1], cause food spoilage; while others like the cholesterol lowering lovastatins can be used as medical drugs [2]. Many SMs are promising leads for anti-cancer drugs as e.g. the non-ribosomal peptide malformins [3]. The enzyme classes producing these distinct compounds — polyketide synthases (PKSs) and non-ribosomal peptide synthetases (NRPSs) — are also seen in PKS-NRPS or NRPS-PKS hybrids, seemingly chimeric genes creating a chimeric compound. The products of hybrids are often bioactive, e.g. the mycotoxins cyclopiazonic acid, pyranonigrin, and cytochalasin. [46].

The evolutionary events leading to new enzymes and hence compounds have been described in detail for PKSs and NRPSs. PKSs diversify by exchange of initiation modules for modification of primer units, module duplication, horizontal gene transfer [7, 8]. Other studies suggest a burst of PKS duplications in the early Pezizomycotina, a predecessor of mainly Ascomycota, [9] as major driver for PKS diversity. For NRPSs, studies suggest duplication and loss of NRPSs, horizontal gene transfer (HGT) from bacteria to fungi, and gain and loss of domains as driver for diversity [10, 11].

In contrast, hybrids have been neglected by phylogenetic studies, although their combination of NRPS and PKS domains suggests an interesting evolutionary history. The existing studies have focused on known hybrids [12]. Lawrence et al. [13] have shown that one Cochliobolus heterostrophus NRPS-PKS hybrid gene originates from Burkholderiales; which they suggest to be acquired by HGT in the early evolution of the Pezizomycotina. A model for the avirulence factor ACE1 gene [14] by Khaldi et al. [15] shows gene duplication, loss, and horizontal gene transfer — a common event between fungi [16] — as the driver in diversification of ACE1 hybrids. ACE1 duplicated in ancestors of Eurotiomycetes and Dothideomycetes giving rise to an ACE1-like hybrid, then during diversification of species. Within Aspergillus species, only A. clavatus preserved the ACE1-like hybrid and in addition received ACE1 from Magnaporthe grisea through HGT (ccsA).

Using genome data from 38 strains of the SM-rich Aspergillus genus and Penicillium chrysogenum we describe the phylogenetic dynamics of PKS-NRPS and NRPS-PKS hybrids and relate them to PKSs and the structurally similar PKS-likes, as well as NRPSs and the structurally similar NRPS-likes (often the likes show shorter domain arrangement). The genus Aspergillus can be divided in taxonomic sections, i.e. groups of fungi with similar morphological and metabolic characteristics, as e.g. the black Aspergilli of section Nigri. Due to their morphology, the species of this section were divided in biseriates and uniseriates, which reflects their metabolic capabilites [17] and genetic diversity [18, 19].

Furthermore, we identify origins of hybrids in bacteria and fungal genera. Understanding hybrid evolution and diversity will provide insights into molecular evolution and put rational engineering of these proteins within our grasp.

Results

Genus wide analysis identifies independent groups of hybrids

Recent work has highlighted the dynamics of SM genes in fungi and their diversifying mechanisms [711]. In this study, we are describing the diversity of rare PKS-NRPS and NRPS-PKS hybrids and compare them to related classes like PKSs and NRPSs. Due to the similarity of these enzymes, we expect that a fusion of NRPSs and PKSs could have occurred during early fungal evolution.

In order to investigate this, we created a Maximum Likelihood phylogeny (ML) of hybrid proteins from a selection of Aspergilli of section Nigri (including biseriates and uniseriates), Circumdati, Candidi, Flavi, Fumigati, Ochracerosii, Terrei, and P. chrysogenum to cover several Eurotiomycetes (Fig. 1).

Fig. 1
figure1

Hybrid dynamics throughout Aspergilli. (A) Maximum Likelihood (ML) phylogeny of PKS-NRPS and NRPS-PKS hybrid proteins was created on aligned and trimmed protein sequences. Branches shown in grey (A-E) are shown in Additional file 9: Figure S1, Additional file 10: Figure S2, Additional file 11: Figure S3, Additional file 12: Figure S4, Additional file 13: Figure S5. Sections and species groups are indicated by tip color; the orientation of hybrids N-type (NRPS-PKS) and P-type (PKS-NRPS) is indicated by tip shape. Percentage of 1000 times bootstrap values (red) are indicated for major branches, the remaining values are shown in Additional file 9: Figure S1, Additional file 10: Figure S2, Additional file 11: Figure S3, Additional file 12: Figure S4, Additional file 13: Figure S5. The phylogenetic tree is available as Additional file 1. (B) Unrooted view of the phylogenetic tree shown in A. Tip colors indicate sections. Branches highlighted in A are marked in B. (C) Classification of nearest neighbors of ML phylogeny. A matrix of tip distances was extracted from the tree and nearest neighbors classified according to their section, thus the barplot shows the origin of hybrid genes by a species. The background color indicates the section of Aspergillus species while the bar color indicates the section of the nearest neighbor hybrid ortholog

In this phylogeny, NRPS-PKS and PKS-NRPS hybrids form several distinct groups, with the PKS-NRPS orientation being more abundant than the NRPS-PKS orientation (Fig. 1, Additional file 1). The analysis indicates related compounds and which hybrids are conserved throughout Aspergilli in different sections.

A group of hybrids containing the cytochalasin producing hybrid A. clavatus 6366 (Additional file 9: Figure S1) shows a large phylogenetic distance to other hybrids, indicating that ACE1-like hybrids are rare in Aspergilli — sustaining the hypothesis by Khaldi et al. [15]. Since A. sclerotioniger 605326 is the nearest neighbour of A. clavatus 6366 we predict this hybrid to produce sclerotionigrin [20], a cytochalasan.

The synteny plot (Fig. 2) further sustains that A. clavatus 6366 originates from another genus, since the synteny to related cytochalasan producing hybrids from Aspergillus species (e.g. A. sclerotioniger 605326) is low.

Fig. 2
figure2

Synteny plots of hybrid gene clusters. A Synteny plot of cytochalasin the cytochalasin hybrid from A. clavatus and related hybrids from the phylogeny (for phylogeny see Additional file 1). Numbers of species names indicate protein id of the colorized hybrid gene. Color indicates section. A. sclerotioniger is known to produce sclerotionigrin — a cytochalasan. The synteny plot shows that the sequence of A. clavatus 6366 differs from the other hybrids. B Synteny plot of the cyclopiazonic acid hybrid from A. flavus. The hybrid gene shows high conservation between Aspergilli from sections Flavi, Nigri and Circumdati (for phylogeny see Additional file 12: Figure S4)

Another group of hybrids is conserved in biseriate Nigri species and A. homomorphus, A. clavatusA. campestris and A. ochraceoroseus (Additional file 10: Figure S2). The tree topology indicates this as common hybrid duplication in Aspergilli through its conservation in many species. The short phylogenetic distance from A. ibericus 400692 and A. sclerotiicarbonarius 380544 to A. campestris 310784 and A. ochraceoroseus 492959 hybrids is surprising, as these species are from different sections.

In one case hybrid 370420 of A. homomorphus — a member of uniseriate species — is forming a subgroup with hybrids conserved in biseriates of section Nigri (Additional file 10: Figure S2). Hence the hybrid has been retained by A. homomorphus or gained by horizontal gene transfer (HGT) from another species.

Our results show similar cases to the ACE1 scenario suggested by Khaldi et al. [15]: A. ibericus hybrids 443386 and 469268 and A. steynii hybrids 454498 and 477231 are duplications with larger phylogenetic distances which suggest duplication, loss, and HGT to happen frequently (Additional file 10: Figure S2).

Hybrids with known compounds show that substrate specificity is unrelated to phylogenetic proximity. Hybrids producing pseurotin A and isoflavipucine (using different substrates) are located in sisterclades (Additional file 11: Figure S3). The broad substrate acceptance of isoflavipucine – shown to create 63 diverse compounds [21] – supports a common origin. Biosynthetically related hybrids like chaetoglobosin and cytochalasin hybrids seem to have evolved in parallel as indicated by the phylogenetic distance in the tree (Additional file 12: Figure S4).

Sparse hybrids indicate HGT

Following are four subgroups, two consisting of PKS-NRPS, including pyranonigrin related hybrids, and two groups of NRPS-PKS orientation (Additional file 13: Figure S5). Pyranonigrin-related hybrids are, with the exception of hybrid A. steynii 463238 unique for section Nigri (at least in the scope of our dataset).

Notably, NRPS-PKS hybrids (Additional file 13: Figure S5 indicated by tip labels) are rare among the analyzed species and are only present in a few species: the biseriates of section Nigri, A. indologenus, A. steynii, A. campestris, and Penicillium chrysogenum (Additional file 13: Figure S5). Their absence in other species and their scarcity point towards recent acquisition by HGT of all NRPS-PKS hybrids.

The phylogeny indicates two major groups of NRPS-PKS hybrids and two hybrids as outgroups (A. campestris 323099 and A. steynii 418130). While one major group is biseriate specific, the other group consists of P. chrysogenum 85311, hybrids from biseriates, and A. indologenus 482416. A. luchuensis and A. piperis are the only species that carry hybrids from both major groups of NRPS-PKS hybrids, pointing towards a HGT before their speciation, or retention of a hybrid. The position of the P. chrysogenum 85311 in the phylogeny points towards HGT as well.

Genus wide analysis provides evidence for HGT

Hybrid diversity in fungi is mostly driven by evolution followed by purifying selection [15] and HGT [13, 15]. With the ML phylogeny established, it was an obvious step to extend our analysis for detection of potential HGT.

If all hybrids were inherited vertically and variation in hybrid content caused by purifying selection, we would expect branches of the phylogenetic tree to only contain hybrids from Aspergillus species of the same section and show longer phylogenetic distances due to accumulation of mutations. Hence, identifying homologs of hybrids with short phylogenetic distance from different sections of Aspergilli indicates HGT.

To find the best homologs of hybrids for each species, we extracted distances of hybrids from the ML phylogeny and classified them according to origin (Fig. 1). This analysis works best with the hybrid-rich section Nigri, as this group includes many closely related species, which reduce bias, but also the other Aspergillus species and P. chrysogenum provide insights into hybrid dynamics.

Our analysis reveals that biseriates of section Nigri contain mostly conserved groups of hybrids, but can contain some hybrids derived from other sections (Fig. 1 panel C). A. sclerotiicarbonarius contains the hybrid 361763 related to the aspyridone hybrid from A. nidulans (see also Additional file 11: Figure S3) and A. heteromorphus contains a majority of hybrid homologs from uniseriate Aspergillus species. A. ibericus contains one third of hybrid homologs from other sections. A. sclerotioniger 605326 is homolog to hybrid 376297 from uniseriate A. saccharolyticus.

Uniseriates of section Nigri only contain a low number of hybrid genes, some of them showing orthology to hybrids from other sections (Fig. 1 panel C). A. saccharolyticus, a uniseriate, contains one hybrid from biseriate species (376297) and one hybrid (388526) which shows high conservation to a hybrid from A. steynii and A. oryzae. The latter is responsible for cyclopiazonic acid (CPA) synthesis, a mycotoxin [22].

In this analysis, we included only few non-Nigri species with often only one representative per section. Hence, for these cases, Fig. 1 will show that all hybrids in these species have homologs in species of other sections than their own. This is of course biased due to the selection of species, and further genome sequencing in the future will do much to deconvolute this. The analysis does however still give a good indication of the origin of hybrids. E.g. the isoflavipucine hybrid (325) from A. terreus and a hybrid (260046) from A. campestris show high conservation and thus a short phylogenetic distance which indicates HGT between these species rather than hybrid conservation. Although not derived by HGT, but still worth mentioning, are hybrids from A. steynii. Its hybrids are representative of almost every subgroup of hybrids in the dataset, showing a high diversity in this species.

A. steynii and section Nigri species contain a large number of diverse hybrids which are related to most subgroups of the dataset. If new lineages of hybrids would frequently emerge throughout sections we would expect more section specific hybrids and A. steynii as well as Nigri species would cover less of the hybrid groups. This suggests that the evolutionary events leading to hybrid generation happened before species diversification in the genus Aspergillus. Another observation is that closely related hybrids are present in many phylogenetically distant sections. This points to diversification of hybrids occuring through recombination events after HGT of hybrids.

Additionally, we expect that NRPS-PKS hybrids were either derived by joining of independent NRPS and PKS genes or acquired independently from another source, since they show large phylogenetic distance to PKS-NRPS hybrids. Since the phylogenetic distance could be biased by the amount of structurally similar PKS-NRPS hybrids, we created further comparisons on basis of single domains.

PKS analysis shows common ancestors for PKSs and hybrids

Since intrinsically, hybrids did show large phylogenetic distances (Fig. 1, Additional files 2 and 3), we hypothesized their origin from related SM genes. Previous studies prove hybrid parts as exchangeable [23, 24], hence, we proposed that other SM genes could join together in filamentous fungi to form a hybrid. In order to study this, we created a ML phylogeny of 1369 ketosynthase (KS) domains of PKS-like, PKS and hybrids to elucidate their phylogenetic relations (Fig. 3).

Fig. 3
figure3

Phylogeny of PKS, PKS-like and hybrid proteins. The maximum likelihood phylogeny was created from KS domains of PKS, PKS-like and hybrid proteins. Tip color shows SM gene type (red: hybrid, green: PKS, blue: PKS-like). Hybrids linked to compounds are labelled with the compound name. Grey highlighted groups A-C are shown in Additional file 14: Figure S6, Additional file 15: Figure S7, Additional file 16: Figure S8. Phylogenetic tree available as Additional file 2

The tree topology shows multiple groups consisting of PKS only, mixed PKS and hybrids and PKS-likes. PKS-likes form two unrelated groups, suggesting that they are largely unrelated to PKSs. Hybrids are separated into two groups. NRPS-PKS hybrids are located as sister clade to 6-methylsalicylic acid (6-MSA) PKS related genes — including PKSs for synthesis of yanuthones, terreic acid and patulin (Fig. 3 branch A, Additional file 14: Figure S6). PKS-NRPSs are clustering together with other PKSs that frequently break into hybrid clades and separate known examples from each other (Fig. 3 branch B, Additional file 15: Figure S7). Thus, we suggest these PKSs and PKS-NRPS hybrids to have common ancestors in fungi. PKSs linked to citreoviridin and pyripyropene are lcoated in a sister clade to hybrids (Fig. 3 branch B, Additional file 16: Figure S8). The pyripyropene PKS has an adjacent adenylation domain in its cluster, thus these PKSs could be the ideal precursor for the molecular evolution of hybrids.

In summary, hybrids do not form a monophyletic clade inside the ML phylogeny — rather, clades contain mixes of PKSs and hybrids. Hence we can hypothesize that PKSs and hybrids had common ancestors — distinct ones for NRPS-PKS and PKS-NRPS genes. Additionally, the analysis shows that NRPS-PKS and PKS-NRPS hybrids are unrelated as indicated by the phylogeny of hybrids (Fig. 1).

Phylogeny of NRPSs and hybrids reveals monophyletic clade of hybrids

Hybrids incorporate amino acids (e.g. tyrosine in case of the cytotoxic aspyridone or L-phenylalanine in case of cytochalasins) into compounds in a manner similar to NRPS and NRPS-likes. Thus we sought to investigate the phylogenetic relationship of these proteins.

We created a ML tree of 2428 adenylation domains from NRPS, NRPS-like and hybrid proteins which interestingly led to mostly monophyletic groups (Fig. 4, Additional file 4 and 5): NRPS-likes form two groups which are monophyletic. Other groups comprise NRPSs which appear to have a common ancestor, with few NRPS-likes forming a sister clade. This indicates that NRPS-likes developed from NRPSs in certain cases. Hybrids form a monophyletic group, they are however located in a sister clade with a group of NRPS and NRPS-likes conserved in uniseriate Nigri species (Additional file 17: Figure S9), NRPS homologs are also found in A. heteromorphus and A. ellipticus). These proteins could possibly have a common ancestor. Overall, NRPS and hybrid evolution seems to be largely independent. Thus domains seem to be specific for either NRPS or hybrid proteins.

Fig. 4
figure4

ML phylogeny of NRPS, NRPS-like and hybrid A domains. Tip colors indicate SM protein type; Tip labels show associated compound (if applicable). The phylogenetic tree shows that hybrids are monophyletic when compared to NRPSs and NRPS-likes. Phylogenetic tree available as additional file 3. Branch A shown in Additional file 17: Figure S10

Extended analysis of hybrids shows two events leading to hybrid evolution

We used protein blast on the NCBI non-redundant protein database to find homologs of Aspergillus hybrid genes. Adenylation domains from 288 best hits were extracted and added to the Aspergillus hybrid adenylation domain set. Subsequent alignment and ML analysis generated the phylogeny in Fig. 5.

Fig. 5
figure5

ML phylogeny of A domains from PKS-NRPS and NRPS-PKS hybrid genes and their best blast hits from the NCBI non-redundant protein database. Protein sequences were aligned with clustal omega, trimmed and filtered using trimal, and a ML phylogeny generated using iqtree. Compounds of hybrids, if present in MIBiG are shown above tips. Thanamycin is shown multiple times because each domain of the lipopeptide synthetase is annotated. Tip colors indicate whether the hybrid domain is derived from the genera Aspergillus, Penicillium or a different taxonomical class. Phylogenetic tree available as additional file 6

Best blast hits were mostly originating from Ascomycete classes: Dothideomycetes, Eurotiomycetes, Leotiomycetes, Sordariomycetes, Xylonomycetes, one Orbilliomycete, and one Exobasidiomycete hybrid were included. Hybrids from bacterial classes include Proteobacteria, Terrabacteria and Planctomycetes.

We found fungal sequences distributed throughout the tree, and although many ascomycete taxa are included, the tree topology indicates that hybrids are conserved throughout these taxa. Certainly our blast search might bias the tree topology. Nonetheless, if PKSs and NRPSs would recombine frequently in fungi, we would expect more intermediates. There are some NRPSs, mostly from Sordariomycetes and Dothideomycetes, which are related to PKS-NRPS hybrids. These could be remnants of ancestral NRPSs which have been donors for hybrids in fungi.

The majority of the tree consists of PKS-NRPS hybrids, while NRPS-PKS hybrids from fungi and bacterial NRPSs and hybrids (from Terrabacteria and Proteobacteria), are co-clustering in one location (Fig. 5 branch A, Additional file 18: Figure S10, Additional files 6, 7, and 8). Inside the cluster, we can identify the thanamycin hybrid gene from Pseudomonas sp. SHC52, a lipopeptide. What’s more, we can identify hybrids KPC78190.1, APD71785.1, WP_023586037.1 from Streptomyces sp. in a sister clade to hybrids from multiple fungal genera. This indicates that lipopeptides and hybrids from Streptomyces could be horizontally transferred to fungi — giving rise to NRPS-PKS hybrids in fungi.

The phylogeny also shows related hybrids of different sections in the same branch, as in the case for genes of aspyridone and fumosorinone — two compounds similar in structure [25]. This supports that the structural diversity of hybrids throughout Ascomycetes might be limited. PKS-NRPS homologs of Aspergilli are co-clustering with many hybrids from Sordariomycetes and Eurotiomycetes. Thus, the distances of pyranonigrin associated hybrids observed earlier (Fig. 1) can now be explained with the added dataset (Fig. 5 branch B, Additional file 19: Figure S11). The A. ellipticus hybrid (460246) is clustering closer together with Sordariomycetes; the same for A. steynii which carries an Eurotiomycete -related hybrid. Recurrence of the same genera emphasizes that hybrid diversity might be limited in fungi, which is why there are usually so few.

Discussion

Analysis of the SM protein repertoire of a genus-wide dataset led us to discover dynamics between NRPSs, NRPS-likes, PKSs, PKS-likes and NRPS-PKS as well as PKS-NRPS hybrids. While previous studies included hybrids related to known examples [12], focused on NRPS domains [11], or were considering single hybrids for analysis [13, 15], we combined analysis of A and KS-domains through all hybrids in a genus wide species set and related them to other SM enzymes to focus on their evolutionary history.

Gallo et al. [12] also showed that the hybrid and NRPS adenylation domains are monophyletic — which our analysis supports. Yet, we could identify a sister clade of NRPS and NRPS-likes which groups together with hybrids. These genes are valuable leads in the investigation of SM gene evolution in the genus Aspergillus. Additionally, one common ancestor of PKSs and hybrid PKS-NRPSs had been hypothesized [12]. Our results indicate multiple PKSs which could be the ancestor of hybrids. We can confirm that despite structural (or biosynthetic) similarity of cytochalasins and chaetoglobosisn, the genes for their synthetases have a distinct phylogenetic history — ccsA, the cytochalasin hybrid, is more closely related to the equisetin hybrid than to the chaetoglobosin hybrid.

Previous studies indicated that NRPS-PKSs are of bacterial origin, while PKS-NRPSs, due to their abundance in fungi, are of fungal origin [13]. Our analysis shows that only the A and KS domains of NRPS-PKS hybrids have similar phylogenetic histories in fungi (as indicated in [13]). The extended comparison to hits from the NCBI non-redundant protein database in this study revealed that the collective of NRPS-PKSs in Ascomycetes is related to bacterial hybrids and lipopeptide synthetases, suggesting a bacterial origin. PKS-NRPS hybrids on the other hand, show similarity to different fungal PKSs in the ML phylogeny. Hence, we suspect that A and KS domains have a different phylogenetic history in this type of hybrids. Their monophyly in A domain comparisons suggests that one NRPS ancestor was able to recombine with different PKSs. Hence our analysis indicates that multiple events rather than one event gave rise to hybrid evolution.

According to previous studies, the ccsA ancestor was duplicated during Ascomycete evolution, lost in Aspergilli, and reacquired by HGT in A. clavatus from Magnaporthe grisea [15]. In our analysis of SM genes throughout 38 Aspergillus strains — with addition of best hits from the NCBI non-redundant protein database — we can further provide evidence for HGT of hybrids through Ascomycetes. Biseriate species A. sclerotioniger, A. heteromorphus as well as uniseriate species A. saccharolyticus, all of which belong to section Nigri, contain conserved ccsA homologs (Fig. 1). A. sclerotioniger is producer of sclerotionigrin, a SM related to cytochalasins from A. clavatus.

Duplication, loss and reacquisition of related hybrids appears to be common since A. ibericus and A. steynii contain duplications with larger phylogenetic distances (Aspibe1_443386 and Aspibe1_469268; Aspste1_454498 and Aspste1_477231). This duplication, loss, and reacquisition pattern is similar to the one proposed for ace1. The ace1 gene is duplicated during Sordariomycete evolution to ace1 and syn2. In an HGT event, syn2 is transferred from Magnaporthe grisea to A. clavatus resulting in ccsA — an organism which lost the original ace1 ancestor [15]. Additionally, we find several cases where phylogenetic distances to hybrids of other sections are shorter than to hybrids of their own section.

For example, A. heteromorphus, although biseriate, contains mostly hybrid homologs from uniseriate species and A. saccharolyticus contains a homolog to a hybrid from A. flavus. This is further supported by our extended analysis using additional hybrids from other fungi and bacteria. A. ellipticus and A. steynii show homologs to Sordariomycetes and Eurotiomycetes, respectively — indicating a loss and reaquisition of a pyranonigrin-related hybrid. The phylogenetic relationship of hybrids points to few events which yielded novel hybrids, the phylogenetic distances however, exclude a loss only diversification. We suggest that hybrids are frequently transferred, hence the short distances between hybrids of distantly related species and the finite number of clades/their low amount.

Aspergilli show vast amounts of SM genes; hybrids however, are the minority (Table 1). These low counts, their dynamics and their unrelated groups (Fig. 1) pose the question whether independent events lead to their evolution and whether their diversity is finite.

Table 1 Number of hybrids and total number of predicted SM proteins in Aspergillus species and Penicillium chrysogenum

Hybrids were subject to many engineering efforts since their structure suggested inpendendent units, which could be recombined. The genes tenS and dmbS, part of the hybrids producing the similar compounds tenellin and desmethylbassianin, respectively, could be fused together and resulted in the production of tenellin A [23]. Nielsen et al. [24] created a fusion of CcsA and Syn2 resulting in niduporthin, a novel chimeric compound. Despite these successful recombinations, creating a chimera of hybrids remained challenging as other studies showed a limitation to recombination efforts [26]. Our study will facilitate further efforts to engineer hybrids since we could show their dynamics inside the genus Aspergillus and relate it to other fungal genera. Especially NRPSs and PKSs, which have been identified in this study as hybrid related, might be amenable to modification.

Conclusion

In this study, we were able to show that several hybrid homologs are shared throughout different sections which points to HGT and frequent gene gain and loss events. What’s more, several PKSs were found together with hybrids in clades, while NRPSs and hybrids are mostly monophyletic. A more extensive search against the non-redundant protein database and subsequent phylogenetic analysis shows that hybrids with different orientation have a separate evolutionary origin, they are frequently subject to HGT between taxa and their diversity in Ascomycetes appears to be finite, as shown by the clusters of the phylogenetic tree of hybrids. All in all, this study contributes to our understanding of the events which led to hybrid evolution and provide cases for recombination experiments.

Methods

Data collection

Protein sequences and SMURF annotations for Aspergillus and Penicillium species were downloaded from JGI (http://genome.jgi.doe.gov/).

Genetic dereplication

Secondary metabolite proteins were annotated with known compounds through BLAST against examples from the Minimum information on biosynthetic gene clusters (MIBiG) database [27].

Identifying best hits in the NCBI non-redundant protein database

Adenylation domains of hybrids were blasted against the NCBI non-redundant protein database (downloaded from http://ftp://ftp.ncbi.nlm.nih.gov/blast/db/) using protein basic local alignment search tool (BLAST) [28]. The best 15 hits which suffice a query coverage cutoff of 80% were retained. Taxonomic labels for best hits were added using ETE 3 toolkit [29] (Additional file 5).

Protein domains, alignment and maximum likelihood analysis

InterproScan5 (version 5.22-61.0) [30] was used to identify domains on protein sequences. Domains of the same type which were under 350 amino acids long and less than 100 amino acids apart were merged before proceeding. Sequences were handled using Biopython [31]. Resulting domain sequences were aligned using Clustal Omega [32] and cut using trimal [33]. In case of the full hybrid protein tree, full protein sequences were aligned and trimmed. IQ-tree [34] was used on aligned sequences using Model Finder Plus [35] and 1000 times ultra fast bootstrap [36].

Finding best homologs in trees

A distance matrix was extracted from the Maximum Likelihood tree using the cophenetic function of the ape package [37] in R [38]. Best homologs were plotted using ggplot2 [39].

Visualization

Phylogenetic trees were visualized using ggtree [40] in combination with the ape package [37] for unrooted trees and ggstance [41]. Gene clusters were plotted using Easyfig [42].

References

  1. 1

    Bennett JW, Klich M, Mycotoxins M. Mycotoxins. Clin Microbiol Rev. 2003; 16(3):497–516. https://doi.org/10.1128/CMR.16.3.497.

  2. 2

    Rubinstein A, Lurie Y, Groskop L, Weintrob M. Cholesterol-Lowering Eeffects of a 10 mg Daily Dose of Lovastatin in Patients with Initial Total Cholesterol Levels 200 to 240 mg/dl (5.18 to 6.21 mmol/liter). Am J Cardiol. 1991; 68(11):1123–6.

    CAS  Article  Google Scholar 

  3. 3

    Wang J, Jiang Z, Lam W, Gullen EA, Yu Z, Wei Y, Wang L, Zeiss C, Beck A, Cheng EC, Wu C, Cheng YC, Zhang Y. Study of malformin C, a fungal source cyclic pentapeptide, as an anti-cancer drug. PLoS ONE. 2015; 10(11):1–19. https://doi.org/10.1371/journal.pone.0140069.

    Article  Google Scholar 

  4. 4

    Awakawa T, Yang XL, Wakimoto T, Abe I. Pyranonigrin E: A PKS-NRPS hybrid metabolite from Aspergillus niger identified by genome mining. ChemBioChem. 2013; 14(16):2095–9. https://doi.org/10.1002/cbic.201300430.

    CAS  Article  Google Scholar 

  5. 5

    Tokuoka M, Seshime Y, Fujii I, Kitamoto K, Takahashi T, Koyama Y. Identification of a novel polyketide synthase-nonribosomal peptide synthetase (PKS-NRPS) gene required for the biosynthesis of cyclopiazonic acid in Aspergillus oryzae. Fungal Genet Biol. 2008; 45(12):1608–15. https://doi.org/10.1016/j.fgb.2008.09.006.

    CAS  Article  Google Scholar 

  6. 6

    Casella JF, Flanagan MD, Lin S. Cytochalasin D inhibits actin polymerization and induces depolymerization of actin filaments formed during platelet shape change. Nature. 1981; 293(5830):302–5. https://doi.org/10.1038/293302a0.

    CAS  Article  Google Scholar 

  7. 7

    Ridley CP, Lee HY, Khosla C. Evolution of polyketide synthases in bacteria. Proc Nat Acad Sci USA. 2008; 105(12):4595–600. https://doi.org/10.1073/pnas.0710107105.

    CAS  Article  Google Scholar 

  8. 8

    Throckmorton K, Wiemann P, Keller NP, Vol. 7. Evolution of Chemical Diversity in a Group of Non-reduced Polyketide Gene Clusters: Using Phylogenetics to Inform the Search for Novel Fungal Natural Products; 2015, pp. 3572–607. https://doi.org/10.3390/toxins7093572.

    CAS  Article  Google Scholar 

  9. 9

    Koczyk G, Dawidziuk A, Popiel D. The distant siblings - A phylogenomic roadmap illuminates the origins of extant diversity in fungal aromatic polyketide biosynthesis. Genome Biol Evol. 2015; 7(11):3132–54. https://doi.org/10.1093/gbe/evv204.

    CAS  Article  Google Scholar 

  10. 10

    Bushley KE, Ripoll DR, Turgeon BG. Module evolution and substrate specificity of fungal nonribosomal peptide synthetases involved in siderophore biosynthesis. BMC Evol Biol. 2008; 8(1):328. https://doi.org/10.1186/1471-2148-8-328.

    Article  Google Scholar 

  11. 11

    Bushley KE, Turgeon BG. Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships. BMC Evol Biol. 2010; 10:26. https://doi.org/10.1186/1471-2148-10-26.

    Article  Google Scholar 

  12. 12

    Gallo A, Ferrara M, Perrone G. Phylogenetic study of polyketide synthases and nonribosomal peptide synthetases involved in the biosynthesis of mycotoxins. Toxins. 2013; 5(4):717–742. https://doi.org/10.3390/toxins5040717.

    CAS  Article  Google Scholar 

  13. 13

    Lawrence DP, Kroken S, Pryor BM, Arnold AE. Interkingdom gene transfer of a hybrid NPS/PKS from bacteria to filamentous Ascomycota. PloS one. 2011; 6(11):28231. https://doi.org/10.1371/journal.pone.0028231.

    CAS  Article  Google Scholar 

  14. 14

    Böhnert HU, Fudal I, Dioh W, Tharreau D, Notteghem J-L, Lebrun M-H. A putative polyketide synthase/peptide synthetase from Magnaporthe grisea signals pathogen attack to resistant rice. Plant Cell. 2004; 16(9):2499–513. https://doi.org/10.1105/tpc.104.022715.

    Article  Google Scholar 

  15. 15

    Khaldi N, Collemare J, Lebrun M-H, Wolfe KH. Evidence for horizontal transfer of a secondary metabolite gene cluster between fungi. Genome Biol. 2008; 9(1):18. https://doi.org/10.1186/gb-2008-9-1-r18.

    Article  Google Scholar 

  16. 16

    Fitzpatrick Da. Horizontal gene transfer in fungi. FEMS Microbiol Lett. 2012; 329(1):1–8. https://doi.org/10.1111/j.1574-6968.2011.02465.x.

    Article  Google Scholar 

  17. 17

    Lamboni Y, Nielsen KF, Linnemann AR, Gezgin YK, Hell K, Nout MJR, Smid EJ, Tamo M, Van Boekel MAJS, Hoof JB, Frisvad JC. Diversity in secondary metabolites including mycotoxins from strains of aspergillus section nigri isolated from raw cashew nuts from benin, west africa. PLoS ONE. 2016; 11(10):1–14. https://doi.org/10.1371/journal.pone.0164310.

    Article  Google Scholar 

  18. 18

    Theobald S, Vesth TC, Rendsvig JK, Nielsen KF, Riley R, Abreu LMD, Salamov A, Frisvad JC, Larsen TO, Andersen MR, Hoof JB. Uncovering secondary metabolite evolution and biosynthesis using gene cluster networks and genetic dereplication. Sci Rep. 2018; 8(17957):1–12. https://doi.org/10.1038/s41598-018-36561-3.

  19. 19

    Vesth TC, Nybo JL, Theobald S, Frisvad JC, Larsen TO, Nielsen KF, Hoof JB, Brandl J, Salamov A, Riley R, Gladden JM, Phatale P, Nielsen MT, Lyhne EK, Kogle ME, Strasser K, McDonnell E, Barry K, Clum A, Chen C, LaButti K, Haridas S, Nolan M, Sandor L, Kuo A, Lipzen A, Hainaut M, Drula E, Tsang A, Magnuson JK, Henrissat B, Wiebenga A, Simmons BA, Mäkelä MR, de Vries RP, Grigoriev IV, Mortensen UH, Baker SE, Andersen MR. Investigation of inter- and intraspecies variation through genome sequencing of Aspergillus section Nigri. Nature Genetics. 2018. https://doi.org/10.1038/s41588-018-0246-1.

    CAS  Article  Google Scholar 

  20. 20

    Petersen LM, Bladt TT, Dürr C, Seiffert M, Frisvad JC, Gotfredsen CH, Larsen TO. Isolation, structural analyses and biological activity assays against chronic lymphocytic leukemia of two novel cytochalasins - Sclerotionigrin A and B. Molecules. 2014; 19(7):9786–97. https://doi.org/10.3390/molecules19079786.

    Article  Google Scholar 

  21. 21

    Qiao K, Zhou H, Xu W, Zhang W, Garg N, Tang Y. A fungal nonribosomal peptide synthetase module that can synthesize thiopyrazines. Organ Lett. 2011; 13(7):1758–61. https://doi.org/10.1021/ol200288w.

    CAS  Article  Google Scholar 

  22. 22

    Sorenson WG, Tucker JD, Simpson JP. Mutagenicity of the tetramic mycotoxin cyclopiazonic acid. Appl Environ Microbiol. 1984; 47(6):1355–57.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. 23

    Heneghan MN, Yakasai AA, Williams K, Kadir KA, Wasil Z, Bakeer W, Fisch KM, Bailey AM, Simpson TJ, Cox RJ, Lazarus CM. The programming role of trans-acting enoyl reductases during the biosynthesis of highly reduced fungal polyketides. Chem Sci. 2011; 2(5):972. https://doi.org/10.1039/c1sc00023c.

    CAS  Article  Google Scholar 

  24. 24

    Nielsen ML, Isbrandt T, Petersen LM, Mortensen UH, Andersen MR, Hoof JB, Larsen TO. Linker flexibility facilitates module exchange in fungal hybrid PKS-NRPS engineering. PLoS ONE. 2016; 11(8):1–18. https://doi.org/10.1371/journal.pone.0161199.

    Article  Google Scholar 

  25. 25

    Liu L, Zhang J, Chen C, Teng J, Wang C, Luo D. Structure and biosynthesis of fumosorinone, a new protein tyrosine phosphatase 1B inhibitor firstly isolated from the entomogenous fungus Isaria fumosorosea. Fungal Genet Biol. 2015; 81:191–200. https://doi.org/10.1016/j.fgb.2015.03.009.

    CAS  Article  Google Scholar 

  26. 26

    Boettger D, Bergmann H, Kuehn B, Shelest E, Hertweck C. Evolutionary Imprint of Catalytic Domains in Fungal PKS-NRPS Hybrids. ChemBioChem. 2012; 13(16):2363–73. https://doi.org/10.1002/cbic.201200449.

    CAS  Article  Google Scholar 

  27. 27

    Medema MH, Kottmann R, Yilmaz P, Cummings M, Biggins JB, Blin K, De Bruijn I, Chooi YH, Claesen J, Coates RC, Cruz-Morales P, Duddela S, Düsterhus S, Edwards DJ, Fewer DP, Garg N, Geiger C, Gomez-Escribano JP, Greule A, Hadjithomas M, Haines AS, Helfrich EJN, Hillwig ML, Ishida K, Jones AC, Jones CS, Jungmann K, Kegler C, Kim HU, Kötter P, et al. Minimum Information about a Biosynthetic Gene cluster. Nat Chem Biol. 2015; 11(9):625–31. https://doi.org/10.1038/nchembio.1890.

    CAS  Article  Google Scholar 

  28. 28

    Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10(1):421. https://doi.org/10.1186/1471-2105-10-421.

    Article  Google Scholar 

  29. 29

    Huerta-Cepas J, Serra F, Bork P. ETE 3: Reconstruction, Analysis, and Visualization of Phylogenomic Data. Mole Biol Evol. 2016; 33(6):1635–8. https://doi.org/10.1093/molbev/msw046.

    CAS  Article  Google Scholar 

  30. 30

    Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, Pesseat S, Quinn AF, Sangrador-Vegas A, Scheremetjew M, Yong SY, Lopez R, Hunter S. InterProScan 5: Genome-scale protein function classification. Bioinformatics. 2014; 30(9):1236–40. https://doi.org/10.1093/bioinformatics/btu031.

    CAS  Article  Google Scholar 

  31. 31

    Cock PJA, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, Friedberg I, Hamelryck T, Kauff F, Wilczynski B, De Hoon MJL. Biopython: Freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009; 25(11):1422–23. https://doi.org/10.1093/bioinformatics/btp163.

    CAS  Article  Google Scholar 

  32. 32

    Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thompson JD, Higgins DG. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mole Syst Biol. 2011; 7(1):539. https://doi.org/10.1038/msb.2011.75.

    Article  Google Scholar 

  33. 33

    Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009; 25(15):1972–73. https://doi.org/10.1093/bioinformatics/btp348.

    CAS  Article  Google Scholar 

  34. 34

    Nguyen LT, Schmidt HA, Von Haeseler A, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mole Biol Evol. 2015; 32(1):268–74. https://doi.org/10.1093/molbev/msu300.

    Article  Google Scholar 

  35. 35

    Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates. Nature Methods. 2017; 14(6):587–91. https://doi.org/10.1038/nmeth.4285.

    CAS  Article  Google Scholar 

  36. 36

    Minh BQ, Nguyen MAT, Von Haeseler A. Ultrafast approximation for phylogenetic bootstrap. Mole Biol Evol. 2013; 30(5):1188–95. https://doi.org/10.1093/molbev/mst024.

    CAS  Article  Google Scholar 

  37. 37

    Paradis E, Claude J, Strimmer K. APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics (Oxford, England). 2004; 20(2):289–90.

    CAS  Article  Google Scholar 

  38. 38

    R Core Team. R: A Language and Environment for Statistical Computing. 2017. https://www.r-project.org/.

  39. 39

    Wickham H. ggplot2: Elegant Graphics for Data Analysis. 2009. http://ggplot2.org.

  40. 40

    Yu G, Smith D, Zhu H, Guan Y, Lam TT-Y. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017; 8:28–36. https://doi.org/10.1111/2041-210X.12628.

    Article  Google Scholar 

  41. 41

    Henry L, Wickham H, Chang W. ggstance: Horizontal ’ggplot2’ Components. 2016.

  42. 42

    Sullivan MJ, Petty NK, Beatson SA. Easyfig: A genome comparison visualizer. Bioinformatics. 2011; 27(7):1009–10. https://doi.org/10.1093/bioinformatics/btr039.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Thomas O. Larsen for discussions on the topic.

Funding

ST, TCV, and MRA acknowledge funding from The Villum Foundation, grant VKR023437.

Author information

Affiliations

Authors

Contributions

ST performed data retrieval, processing, performed ML analyses and wrote the manuscript. ST, TV and MR analyzed results. MR and TV advised ST on manuscript content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mikael R. Andersen.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1

ML phylogeny of NRPS-PKS and PKS-NRPS hybrids in Aspergilli.

Additional file 2

ML phylogeny of KS-domains of PKS, PKS-like and hybrid proteins in Aspergilli.

Additional file 3

List of collections of identical sequences from Additional file 2.

Additional file 4

ML phylogeny of A domains of NRPS, NRPS-like and hybrid proteins in Aspergilli.

Additional file 5

Data frame for A-domain tree in additional file 4.

Additional file 6

ML phylogeny of Aspergillus A-domains of hybrid proteins and hybrid proteins from best hit analysis against the non-redundant protein database.

Additional file 7

Data frame for A-domain tree in additional file 6.

Additional file 8

NCBI identifier to taxonomical class translation table.

Additional file 9: Figure S1

Branch A from hybrid maximum likelihood phylogeny (Fig. 1). Sections and species groups indicated by tip color; Orientation of hybrids N-type (NRPS-PKS) and P-type (PKS-NRPS) indicated by tip shape. Tip labels constist of jgi organism name, jgi protein id and associated compound (if applicable). Percentage values of 1000 times bootstrap are indicated in red next to the node.

Additional file 10: Figure S2

Branch B from hybrid maximum likelihood phylogeny (Fig. 1). Sections and species groups indicated by tip color; Orientation of hybrids N-type (NRPS-PKS) and P-type (PKS-NRPS) indicated by tip shape. Tip labels constist of jgi organism name, jgi protein id and associated compound (if applicable). Percentage values of 1000 times bootstrap are indicated in red next to the node.

Additional file 11: Figure S3

Branch C from hybrid maximum likelihood phylogeny (Fig. 1). Sections and species groups indicated by tip color; Orientation of hybrids N-type (NRPS-PKS) and P-type (PKS-NRPS) indicated by tip shape. Tip labels constist of jgi organism name, jgi protein id and associated compound (if applicable). Percentage values of 1000 times bootstrap values are indicated in red next to the node.

Additional file 12: Figure S4

Branch D from hybrid maximum likelihood phylogeny (Fig. 1). Sections and species groups indicated by tip color; Orientation of hybrids N-type (NRPS-PKS) and P-type (PKS-NRPS) indicated by tip shape. Tip labels constist of jgi organism name, jgi protein id and associated compound (if applicable). Percentage values of 1000 times bootstrap are indicated in red next to the node.

Additional file 13: Figure S5

Branch E from hybrid maximum likelihood phylogeny (Fig. 1). Sections and species groups indicated by tip color; Orientation of hybrids N-type (NRPS-PKS) and P-type (PKS-NRPS) indicated by tip shape. Tip labels constist of jgi organism name, jgi protein id and associated compound (if applicable). Percentage values of 1000 times bootstrap are indicated in red next to the node.

Additional file 14: Figure S6

Branch A from phylogeny of PKS, PKS-like and hybrid proteins (Fig. 3). Percentage values of 1000 times bootstrap below 100 are shown in red.

Additional file 15: Figure S7

Branch B from phylogeny of PKS, PKS-like and hybrid proteins (Fig. 3). Percentage values of 1000 times bootstrap below 100 are shown in red.

Additional file 16: Figure S8

Branch C from phylogeny of PKS, PKS-like and hybrid proteins (Fig. 3). Percentage values of 1000 times bootstrap below 100 are shown in red.

Additional file 17: Figure S9

Branch A from phylogeny of NRPS, NRPS-like and hybrid proteins (Fig. 4). Tip color indicates section/subgroup; tip shape indicates SM protein type; tip label shows jgi protein id and associated compound (if applicable). A group of NRPS and NRPS-likes from uniseriate Nigri species forms a sister clade to the monophyletic hybrids. Percentage values of 1000 times bootstrap below 100 are shown in red.

Additional file 18: Figure S10

ML phylogeny of fungal and bacterial hybrids. Subtree extracted from Fig. 5. Tip labels show species name and NCBI identifier/ JGI organism and protein identifier. Tip color indicates genus or class; Tip shape indicates SM protein type. Percentage values of 1000 times bootstrap below 100 are shown in red. Hybrids from Streptomyces form a sister clade to fungal hybrids NRPS-PKS hybrids, indicating this class from bacterial origin.

Additional file 19: Figure S11

ML phylogeny of pyranonigrin associated hybrids. Subtree extracted from Fig. 6. Tip labels show species name and NCBI identifier/ jgi organism and protein identifier. Tip color indicates genus or class; Tip shape indicates SM protein type. Percentage values of 1000 times bootstrap below 100 are shown in red. Additional tip label shows associated compound (if applicable).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Theobald, S., Vesth, T.C. & Andersen, M.R. Genus level analysis of PKS-NRPS and NRPS-PKS hybrids reveals their origin in Aspergilli. BMC Genomics 20, 847 (2019). https://doi.org/10.1186/s12864-019-6114-2

Download citation

Keywords

  • Aspergillus
  • PKS-NRPS hybrids
  • Secondary metabolites
  • Gene clusters