# Molecular and functional evolution of the fungal diterpene synthase genes

Marc Fischer, Camille Rustenhloz, Véronique Leh-Louis, Guy Perrière

► To cite this version:

Marc Fischer, Camille Rustenhloz, Véronique Leh-Louis, Guy Perrière. Molecular and functional evolution of the fungal diterpene synthase genes. BMC Microbiology, 2015, 15 (1), pp.221. 10.1186/s12866-015-0564-8. http://hal-01223553

HAL Id: hal-01223553

https://hal.science/hal-01223553

Submitted on 7 Dec 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Molecular and functional evolution of the fungal diterpene synthase genes

Marc JC Fischer¹*, Camille Rustenhloz¹, Véronique Leh-Louis² and Guy Perrière³

**Abstract**

**Background:** Terpenes represent one of the largest and most diversified families of natural compounds and are used in numerous industrial applications. Terpene synthase (TPS) genes originated in bacteria as diterpene synthase (di-TPS) genes. They are also found in plant and fungal genomes. The recent availability of a large number of fungal genomes represents an opportunity to investigate how genes involved in diterpene synthesis were acquired by fungi, and to assess the consequences of this process on the fungal metabolism.

**Results:** In order to investigate the origin of fungal di-TPS, we implemented a search for potential fungal di-TPS genes and identified their presence in several unrelated Ascomycota and Basidiomycota species. The fungal di-TPS phylogenetic tree is function-related but is not associated with the phylogeny based on housekeeping genes. The lack of agreement between fungal and di-TPS-based phylogenies suggests the presence of Horizontal Gene Transfer (HGTs) events. Further evidence for HGT was provided by conservation of synteny of di-TPS and neighbouring genes in distantly related fungi.

**Conclusions:** The results obtained here suggest that fungal di-TPSs originated from an ancient HGT event of a single di-TPS gene from a plant to a fungus in Ascomycota. In fungi, these di-TPSs allowed for the formation of clusters consisting in di-TPS, GGPPS and P450 genes to create functional clusters that were transferred between fungal species, producing diterpenes acting as hormones or toxins, thus affecting fungal development and pathogenicity.

**Keywords:** Diterpene synthase, Horizontal gene transfer, Fungi, Geranylgeranyl diphosphate synthase

**Background**

Horizontal Gene Transfer (HGT) is an exchange of genetic material between different strains or species ultimately integrated into the genome. The existence of HGTs has been perceived by the general public as a major health issue, with the recent appearance and spread of bacteria that are multiresistant to antibiotics. The rapid increase in available genomic sequences, combined with the development of analytical techniques, has enabled the detection of numerous events of recent or ancient HGTs between prokaryotes and eukaryotes, as well as within each group of organisms [1]. Following a recent review of the presence of HGT in fungi, it was found that these events are relatively common [2, 3].

The majority of these transfer events involve bacterial [4–6] or fungal donors [7–9]. The transfer of genes between plants and fungi is also possible [10–12], and the transfer of a nitrate assimilation cluster from the *Oomycota* to an ancestral *Dikarya* species was also recently documented [13].

Plants, fungi and bacteria are able to synthesize diterpenes. Sequence analyses, suggest that the di-TPSs from those organisms share a common evolutionary origin resulting from an ancient HGT [14, 15]. Two types of bacterial di-TPSs have been described. Type I di-TPSs initiate a cyclization reaction via heterolytic cleavage of the polyprenyl diphosphate, whereas type II di-TPSs initiate this reaction via protonation of a double bond or an epoxide ring. It has been suggested that both di-TPSs were horizontally transferred from soil bacteria to plants and that an early fusion of bacterial type I and II enzymes led to an increase in size of the encoded protein, from ~500 in bacteria to ~900 aa in plants [16, 17]. However, it cannot not be completely ruled out that a yet undescribed

<footnote>* Correspondence: marc.fischer@unistra.fr
¹Université de Strasbourg, INRA, Inst Natl Recherche Agron, Métab Second Vigne, Unit Mixte Recherche Santé Vigne & Qual Vins, 28 rue de Herrlisheim, F-68021 Colmar, France
Full list of author information is available at the end of the article</footnote>
bacterium, already endowed with a fused type I and II di-TPS, may have existed and may have been responsible for this HGT event. TPS genes probably played an important role in the development of land plants, as they made it possible for the gibberellin phytohormones to stimulate growth and trigger developmental processes such as seed germination, flowering, fruit formation and plant senescence [18–20]. In addition, these genes also played a role in the synthesis of molecules involved in plant defence [21]. The gibberellin phytohormone pathway was presumably not fully present in early plants. The moss Physcomitrella patens does not produce gibberellins. Instead, it utilizes the diterpene metabolite *ent*-kaurenoic acid, a gibberellin precursor, as its endogenous developmental regulator [19]. The lycophyte Selaginella moellondorffii also does not produce gibberellins. However, it encodes the gibberellin signaling pathway in its genome and uses it for spore and sexual organ development [18].

Plant diterpenoids form a very large family, with more than 3,000 distinct structures described thus far [22]. Furthermore, diterpenes belong to the very large isoprenoid family, which includes more than 40,000 described molecules [23]. Isoprenoids are classified by the number of carbons in the skeletal structure: hemi (C5), mono- (C10), sesqui- (C15), di- (C20) and tri- (C30) terpenoids are the most basic, and latex (C > 1,000) the largest representatives. Angiosperms usually have several genes encoding mono-, sesqui- and di-terpenoid synthases. The genome of the embryophyte *P. patens*, the most basal land plant species for which a genome sequence is available, together with data gathered on gymnosperms and angiosperms, have suggested part of the evolutionary history of this family of enzymes [22, 24–26]. Current evidence suggests that early land plants had a single di-TPS gene but no monoterpane or sesquiterpane synthase genes. The ancestral di-TPS gene likely underwent multiple duplication events, as revealed by the presence of only one di-TPS gene in the *P. patens* genome (a bifunctional di-TPS copalyl diphosphate synthase/kaurene synthase, [25]), as opposed to the presence of 66 TPS-like genes in the genome of the later diverging *S. moellondorffii* [27]. Most di-TPSs are approximately 210 aa longer than mono- or sesqui-TPSs, and plant mono- and sesqui-TPSs probably evolved from di-TPSs via the loss of an internal sequence [16]. Most mono- and di-terpenoid synthases also contain an N-terminal plastid transit peptide, whereas this sequence is usually absent from sesqui-TPSs.

It is also noteworthy that gibberellins are not limited to the plant and bacterial kingdoms. Gibberellin was first discovered in the fungus *Fusarium fujikuroi* (=Gibberella fujikuroi), and was shown to play an important role in the symptomology of bakanae, an *F. fujikuroi*-incited disease of rice. Gibberellins are highly diverse, with 136 analogs described from plants, fungi and

bacteria (http://www.plant-hormones.info/gibberellin_-nomenclature.htm). Fungal di-TPS enzymes are highly similar in size to plant enzymes, and the combination of biochemistry studies with molecular genetic analyses has facilitated the comparison of plant and fungal biochemical pathways leading to the formation of gibberellins in plants and fungi. The first step in gibberellin biosynthesis in both plants and fungi is catalyzed by a di-TPS, whereas later steps in the biosynthesis differ (post *ent*-7-hydroxy kaurenoic acid). The data indicate that the plant and fungal di-TPSs catalysing the initial steps in gibberellin biosynthesis are similar and therefore share a common evolutionary origin, whereas genes encoding enzymes for the late stage in the pathway most likely have different evolutionary origins in fungi and plants and were not obtained by HGT from plants or bacteria [28, 29]. The recent availability of fungal genomes in large numbers has provided us with the opportunity to investigate the multiple transfers between fungi of genes involved in diterpene synthesis and to discuss the consequences of transfers on fungal metabolism. In the present study, we sought to determine whether these transfers involved clusters of genes acting in the same functional pathway or whether only the di-TPS was transferred. We also found putative di-TPSs in Basidiomycota and sought to investigate the mechanism allowing these fungi to acquire their di-TPS.

## Results

### Distribution of di-TPS genes in fungal genomes

By querying various genome databases covering ca. 1000 fungal species and following a manual curation step of the annotations, we found 57 putative di-TPS sequences. In most cases, the gene annotation predicted/proposed by the fungal genome project was used. Otherwise, the annotation was completed by searching for exon and intron positions. BLAST searches across fungal species generated eight positive hits for putative di-TPSs that were not shortlisted for further analyses, the inferred gene being either too small (*Aulographum hederae* 408580 and 3987706, *Aspergillus kawachii* KAW_01797, *A. niger* XP_00139461, *Neofusicoccum parvum* EOD52872, *Phaeosphaeria nodorum* XP_001792462, *Punctularia strigosozonata* (12-scaffold join 470791–470964; 471015–472545; 472609–472664; 472694–472709), *Rhytidhysteron rufulum* (NODE_501 join 2571–1717; 1618–677) or too large (*Fusarium oxysporum* EGU79219) to correspond to a classical di-TPS (ca 950aa). It is also noteworthy that the sesquiterpenoid trichothe-cene synthase-like genes, found in several unrelated fungi [30], have weak similarities with di-TPS, especially in their 5' section. These sequences are very short (around 300aa) and were therefore not considered in our analyses. The list of the 57 fungal di-TPS genes examined in this study is provided in Additional file 1.
The 48 fungal species containing di-TPSs genes are not evenly distributed over the fungal phylogeny. In paerticular, none of the species in early diverging lineages of fungi contain a di-TPS (Fig. 1). Di-TPS genes are found exclusively in the *Dikarya* subkingdom (Ascomycota and Basidiomycota)é Furthermore, they are more common across Ascomycota than Basidiomycota: 30 of the 112 Ascomycota genera examined include at least one species with a di-TPS gene in its genome. In contrast, only four of 59 genera of Basidiomycota examined had di-TPS genes. Ascomycota species having di-TPS genes are not uniformly distributed in the Ascomycota tree, since none of the 31 Saccharomycetes genera gave a positive hit for di-TPS. The large number of Saccharomycetes genomes available means that it is unlikely that di-TPS genes are widespread in these organisms. It is not yet known whether di-TPSs are absent from the Pezizomycetes and Orbilomycetes species, because genome sequences for only two genera are available. In view of the small number of genomes available for Pezizomycetes and Orbilomycetes, it is

possible that species belonging to these groups may also contain di-TPS genes, which would imply the presence of these genes in all groups from Pezizomycotina. In the case of Basidiomycota, we found only five species (from four genera out of 60) having di-TPS genes, all of which belong to the Agaricomycotina subphylum, comprising four different families and genera. We were able to obtain full-length sequences for *Punctularia strigosozonata* (Corticiales), *Serpula lacrymans* (Boletales) and *Gymnopus luxurians* (Agaricales). In the case of *Moniliophthora perniciosa* and *M. rereri* (Agaricales), the sequences were not included in our analyses since these genome sequence tabases (1.9X and 15.0X respectively) have contigs that are too small to contain full-length di-TPS sequences. The *Moniliophthora* di-TPSs are nevertheless expected to be functional, as these fungi are known to produce plant-like hormones, including gibberellins and auxins [31]. According to current knowledge, Tremellomycetes, Ustilaginomycotina and Pucciniomycotina contain only a small number of species for which genomic information is

<figure><img src="image_3.png" /></figure>
available, and it is thus possible that some species be-
longing to this group also contain di-TPS genes, which
would mean that they are present in Basidiomycota.

The phylogenetic analysis of the 57 selected di-TPS
genes (Fig. 2) shows a monophyletic clade of all di-TPS
sequences, suggesting that they may have diverged from
a common ancestral sequence. Interestingly, some as-
pects of the topology of the resulting gene clade di-TPS
tree are not correlated with the species phylogeny
(Fig. 1). The di-TPS phylogeny may to a certain extent
be related to the lifestyle of the species. A large group of
27 sequences concentrates 19 of the 22 fungi living in
close relationship with plants, among which 17 are plant
parasites (Fig. 2, green box). This group does not contain
di-TPS sequences from saprophytic fungi. The 11 se-
quences from animal pathogenic fungi are classed into
two groups, one of which is scattered across saprophytic
species (six sequences). The other group forms a unique
clade in the plant parasite and plant saprophyte group.
This clade of five sequences includes the eurotiomycetous
species *T. equinum*, *T. verrucosum*, *T. rubrum*, *T. tonsur-
ans* and *A. benhamiae* sequences (Fig. 2, pink box) species
belonging to *Arthrodermataceae*, suggesting a common
ancestral origin (a plant parasitic donor) of di-TPS genes
for these species.

Catalytic sites, copy number variation, exon-intron patterns,
and functional clusters

Di-TPSs were checked for potential functional catalytic
sites. Recent phylogenies from plant TPSs concentrated at
the catalytic sites have proposed three major classes of
plant di-TPSs (Type I, II, I/II), which can be distinguished
by the presence of a certain number of typical features
[32]. Type I di-TPSs contain a DDXXD motif and NSE/
DTE signature, type II di-TPSs contain a DXDD motif,
and bifunctional type I/II di-TPSs contain both DXDD
and DDXXD motifs and a NSE/DTE signature. Such
bifunctional type I/II di-TPSs are only found in primitive
plants, such as *S. moellondorffii* (h subfamily) as well as in
some gymnosperms such as *Abies grandis* (d subfamily)
[33]. A multiple sequence alignment (Table 1) was per-
formed on the fungal sequences for the Mg²⁺-PPi binding,
H⁺-initiated cyclization and Mg²⁺-ionization/cyclization
sites already described in plants [16]. Typical bifunctional
type I/II signatures were usually found in fungal di-TPS,
which is in agreement with biochemical data that found
them to be cytosolic bifunctional enzymes [34]. Mg²⁺-PPi
binding regions contain the expected D/E residues, except
for A. oryzae_a (no D/E residue), C. purpurea_a and B.
fuckeliana_a (no alignment in this region). DXDD H⁺-ini-
tiated cyclization sites failed to be conserved in only 3 out
of 47 cases (P. anserina_a, R. rufulum_a, Z. passerinii_a).
In fungi, the Mg²⁺-ionization/cyclization DDXXD site
tends to be replaced by a DEXXE site. This DEXXE site

was not found in 11 cases (A. nidulans_a, G. lozoyensis_a,
L. palustris_a, M. oryzae_b, P. marneffei_a, R. rufulum_a,
R. rufulum_b, T. stipitatus_a, T. stipitatus_b, X. parieti-
na_a, Z. passerinii_a). The results found on fungal di-TPSs
show that in most cases, all three Mg²⁺-PPi binding, H
⁺-initiated cyclization and Mg²⁺-ionization/cyclization
sites are conserved, whereas plant di-TPSs having all three
sites are uncommon [16].

The number of fungal di-TPS genes per genome
ranges from one copy (39 species out of 48) in most
cases, to two copies (for *C. purpurea*) or paralogs (for *A.
sojae*, *M. oryzae*, *P. amygdali*, *R. rufulum*, *T. stipitatus*),
or three copies (for *G. luxurians*) or paralogs (for *A.
nidulans*). When the phylogenetic position of each of
these sequences is examined, two types of clustering are
found. 1) in the case of *C. purpurea*, both copies are
clustered together, suggesting that they have resulted
from a recent duplication event that took place specific-
ally in the genome of this species. In the case of *G. lux-
urians*, the three copies are clustered together. However,
the fact that only five Basidiomycota sequences are avail-
able does not show that they resulted from a duplication
event that took place specifically in the genome of this
species; 2) in all other cases, both paralogs are scattered
throughout the tree, showing that it is not straightfor-
ward to determine the origin of the additional sequence.
However, inter-species clustering of the di-TPS genes
can be observed for *A. nidulans* and *A. sojae*, since the
di-TPS genes A. nidulans_a and A. sojae_a belong to the
same clade, whereas the other two sequences, A. nidu-
lans_b and A. sojae_b, are both located in a different
clade.

The fungal di-TPSs also exhibit considerable variation
in their intron/exon patterns (Additional file 2). Some
species (*A. sojae*, *A. niger*) have intronless di-TPS but
also contain di-TPS genes with introns. There appears to
be no unique intron/exon pattern within each of the
Ascomycota, Eurotiomycetes, Dothideomycetes and Sor-
dariomycetes groups. On the other hand, the two Basidi-
iomycota *P. strigosozonata* and *S. lacrymans* di-TPSs
exhibit a rather similar intron/exon pattern.

The general genomic context in the chromosomic re-
gion flanking the fungal di-TPS genes was investigated
by searching for potential GGPPS or cytochrome P450
(P450) genes that could create a functional cluster
(Additional file 3). Table 2 indicates that 19 out of 56
regions contain a putative GGPPS. The frequency of
GGPPS genes per genome is not uniform in the three
Ascomycota groups represented by more than one spe-
cies: the Eurotiomycetes (4/18) have a lower frequency of
GGPPSs than the Dothideomycetes (4/11) and Sordario-
mycetes (8/19). The two Basidiomycota *P. strigosozonata*
and *S. lacrymans* di-TPSs-surrounding regions contain a
putative GGPPS, whereas no putative GGPPS were found
<figure><img src="image_4.png" /><figcaption>**Fig. 2** (See legend on next page.)</figcaption></figure>
(See figure on previous page.)

**Fig. 2** Phylogenetic tree generated from inferred amino acid sequences of fungal and plant di-TPSs. The fungal and plant di-TPS used are listed in Additional file 1. For species possessing multiple di-TPS genes, “_a”,”_b” and “_c” are added to the species designation. The colours used for the species refer to their phylogenetic positions and the Oblongs refer to their life style, as described in the legend. The tree is rooted at the level of the plant clade and branch values correspond to aLRT results (only values ≥ 0.95 are shown). The proposed HGT cases are highlighted by yellow squares

for the three *G. luxurians* di-TPSs-surrounding regions. Table 2 also indicates that P450 are more frequent and numerous than GGPPSs in Ascomycota: only five regions have no P450 (none of these 5 regions without P450 contains a GGPPS), 11 regions have one P450, 17 regions have two P450, 12 regions have three P450, seven regions have four P450, and one region has five P450. The two Basidiomycota *P. strigosozonata* and *S. lacrymans* di-TPS surrounding regions contain six and five putative P450, respectively, whereas no putative P450 were found for the three *G. luxurians* di-TPSs surrounding regions. Although some data related to the di-TPS or its immediate surrounding was available, no information could be found concerning the genomic context of L. palustris_a and C. paspali_a. The four positive hits resulting from a BLAST search across fungal species that were too small to be shortlisted did not contain GPPS or P450 genes next to the putative di-TPS (*A. kawachii* KAW_01797, *A. niger* XP_00139461, *P. nodorum* XP_001792462, *P. strigosozonata* 12-scaffold, *R. rufulum* NODE_501). The *F. oxysporum* EGU79219, considered too large to correspond to a classical di-TPS, forms a putative cluster with one GGPPS and 4 putative P450 genes.

## Horizontal gene transfers

Potential di-TPS HGTs between fungi were assessed by visual inspection of the phylogeny and further investigation of the surrounding regions, in order to highlight potential similarities in gene organization between species that are not closely related (Fig. 3). A potential HGT was then confirmed by best hit BLAST searches performed on the genes involved in a potential gene cluster (Additional files 4, 5, 6 and 7). The following set of di-TPSs having similar genes in not closely related organisms was found: A. niger_a and T. reesei_a, C. globosum_a and B. fuckelia-na_a, R. rufulum_b and G. lozoyensis_a, S. lacrymans_a and P. strigosozonata_a. Another potential HGT shown in Fig. 2 (P. teres_a and P. amygdali_b, aLRT = 1.0) is not shown in Fig. 3, since no surrounding genes other than GGPPS and P450 were involved in the potential cluster. In the case of the putative *G. lozoyensis* / *R. rufulum* HGT (aLRT = 0.99), it is noteworthy that neither di-TPSs R. rufulum_b, nor G. lozoyensis_a contain the Mg²⁺-ionization/cyclization site encountered in almost all fungal di-TPS. A. niger_a, T. reesei_a, C. globosum_a and A. nidulans_a have a rather similar gene organization, which is highlighted in Fig. 3. The *S. lacrymans* and *P. strigosozonata* putative HGT was further analyzed at GGPPS level, and the GGPPS phylogenetic tree (Fig. 4) shows that *S. lacrymans* EGO03064 and *P. strigosozonata* EIN09901 are not located within the main clade of Basidiomycota GGPPSs, but within a group of 18 Ascomycota GGPPs. The fungal and plant GGPPS used are listed in Additional file 8 and the entire phylogenetic tree is described in Additional file 9. These Ascomycota species correspond to 11 Eurotiomycetes, 4 Sordariomycetes and 3 Saccharomycetes GGPPs. On the other hand, the *S. lacrymans* and *P. strigosozonata* GGPPSs located outside the di-TPS/GGPPS/P450 cluster (Basidio_63, Basidio_64, Basidio_65) do branch with the main clade of Basidiomycota GGPPSs.

## Discussion

One first question that can be raised concerning an HGT is the time of the transfer and its potential relationship with the radiation of the main fungal lineages. Land plants and fungi diverged long ago and display major metabolic and developmental differences. It is also quite speculative to draw an analogy between plant and fungal evolution, because the origin of fungi and the dating of fungal divergences are highly inconsistent. Some of the difficulties arise from the use of a unique value for the molecular clock, and tying the dating of all tree branches to just one dating event in the fossil record. Studies using several fossils to multi-calibrate the clock tend to favour the origin of fungi around 760 to 1060 million years ago (mya) [35, 36]. Other critical dating suggested by these studies include the diversification of Pezizomycotina at around 520–320 mya, thus coincident with the diversification of land plants (480–435 mya) and the assumed mutualistic associations between Ascomycota and plants [37]. Indeed, it is often suggested that the diversification of both the plant and fungal kingdoms in terrestrial environments was partly dependent on their ancient ecological associations [38–40]. However, the exchange of genetic material between both kingdoms has not been widely studied. The data from Figs. 1 and 2 suggest that the di-TPS HGT may have taken place in an ongoing context of complex mutualistic interactions, involving numerous associations such as mycorrhiza, plant diseases, saprotrophic degradation of plant by fungi and commensal interactions.

The study of Richards et al. [10] highlights five HGT events from fungi to plants and four from plants to fungi.
<table><caption>**Table 1** Partial alignment of fungal di-TPS highlighting catalytic domains</caption><thead><tr><th>Fungus</th><th>Sequence name</th><th>Mg<sup>2+</sup>-PPi binding<br>(D/E enriched region after an alpha helix)</th><th>H<sup>+</sup>-initiated cyclization (DXDD)</th><th>Mg<sup>2+</sup>-ionization/cyclization DEXXE</th></tr></thead><tbody><tr><td>A. take_a</td><td>RDQLGQWHDVEKSNHIGVEL</td><td></td><td>APHTADVDDTAKALLAL</td><td>MLSMFGYQCDEFIETVAGPA</td></tr><tr><td>C. paspali_a</td><td>QQQLSEWNDVEITNHIGVEL</td><td></td><td>APHTADVDDTAKALLAL</td><td>MLSLLGYQTDEFIETVAGSA</td></tr><tr><td>F. fujikuroi_a</td><td>KRQLAVWNDVEDTNHIGVEF</td><td></td><td>APRTADVDDTAKALLAL</td><td>YLSLLGYQTDEYMEAVAGPV</td></tr><tr><td>F. proliferatum_a</td><td>KRQLAVWKDVEETNHIGVEL</td><td></td><td>APRTADVDDTAKALLAL</td><td>YLSLLGYQTDEYMEAVAGPV</td></tr><tr><td>C. purpurea_b</td><td>QLQLSQWDDVEKSNHIGVEL</td><td></td><td>APHTADVDDTAKALLAL</td><td>MLSLLGYQTDEFIEAVAGPA</td></tr><tr><td>M. oryzae_b</td><td>RTQLQNWDDVLSTNHIGIEI</td><td></td><td>APRAPDADDTAKGLLAL</td><td>LVSMLSYQADEFIHKSAAPA</td></tr><tr><td>P. nodorum_a</td><td>EEQLIAWDDVLDTNHIGVEM</td><td></td><td>APRAVDVDDTAKGLLTL</td><td>IISMLGYQIDEFFEAEAAPA</td></tr><tr><td>S. manihoticola_a</td><td>VSQLQQWDDLVESNHIGVEL</td><td></td><td>APRTADVDDTAKGLLAL</td><td>LLSMYGYQNDEFFETHAMAG</td></tr><tr><td>M. graminicola_a</td><td>QSQLATWTDMEATNHIGVEL</td><td></td><td>APGTVDVDDTAKGLLAM</td><td>FMTLINIQVDEFIEGIATQA</td></tr><tr><td>B. sorokiniana_a</td><td>DRQLNAWNDLVSTNHIGVEL</td><td></td><td>VQRAVDVDDTAKGLLAL</td><td>VMSLLGFQTDEFFEAFATPA</td></tr><tr><td>B. maydis_a</td><td>DRQLNAWNDLVSTNHIGVEL</td><td></td><td>IQRAVDVDDTAKGLLAL</td><td>VMSLIGFQTDEFFEAFATPA</td></tr><tr><td>A. benhamiae_a</td><td>RSQLEEWRVSETMHHGFEII</td><td></td><td>PDFQSDMDDTARGMIAL</td><td>IISFLNYQADEFLEAVAGPH</td></tr><tr><td>T. verrucosum_a</td><td>RSQLEEWRVSETMHHGFEII</td><td></td><td>PDFQSDMDDTARGMIAL</td><td>IISFLNYQADEFLEAVAGPH</td></tr><tr><td>T. equinum_a</td><td>RSQLEEWRVSETMHHGFEII</td><td></td><td>PDFQPDMDDTARGMIAL</td><td>VISFLNYQADEFLEAVAGPQ</td></tr><tr><td>T. rubrum_a</td><td>RSQLEEWRVSETMHHGFEII</td><td></td><td>PDFQSDMDDTARGMIAL</td><td>IISFLNYQADEFLEAVAGPH</td></tr><tr><td>T. tonsurans_a</td><td>RSQLEEWRVSETMHHGFEII</td><td></td><td>PDFQPDMDDTARGMIAL</td><td>VISFLNYQADEFLEAVAGPQ</td></tr><tr><td>P. amygdali_b</td><td>RAQLAALDVSTTEHVGFEII</td><td></td><td>PGFQADVDDTAKTISTL</td><td>FIAMLNFQLDEFMEATAGIL</td></tr><tr><td>E. lata_a</td><td>QRQLQNWDVESSTHVGFEVI</td><td></td><td>PFLQADVDDTAKGLISL</td><td>VISFLNFQADEFMEATAGPA</td></tr><tr><td>E. amarillans_a</td><td>RHQLLEWDVAATTHVGFEII</td><td></td><td>PHLPTLDDDTAKGIVSL</td><td>IISFLNYQTDEIMEVVAGAA</td></tr><tr><td>E. typhina_a</td><td>HHQLLEWDVAAATHVGFEII</td><td></td><td>PHLPTLDDDTAKGIVCL</td><td>IISFLNYQTDEFMEWAGAA</td></tr><tr><td>P. amygdali_a</td><td>QRQLSDWDIASTAHVGFEII</td><td></td><td>RSISSDADDTAKAAFAL</td><td>MVSWWDYQADEFMEAVAGLT</td></tr><tr><td>M. oryzae_a</td><td>ERKLSCWNVSETTHVGFEII</td><td></td><td>PSIMADADDTAKTLSAL</td><td>HISMLSYQGDEYMESVAAPA</td></tr><tr><td>P. beta_a</td><td>HSLLEEWDVSTTNNVVGFEII</td><td></td><td>SGIEPDVDDSAKIVTTL</td><td>ALSFFTYQVDEFMEAVAGPA</td></tr><tr><td>A. otae_a</td><td>KDQLEAWDVASADRVGFEIL</td><td></td><td>CNLEADADDTAKALLAL</td><td>ILSFLNYQADEYMEAVGRY</td></tr><tr><td>A. nidulans_c</td><td>SRQLQAWSVKDTHVHGFEII</td><td></td><td>PSIQADADDTAKSLIAV</td><td>VISFLNYQADEFMEAWGRL</td></tr><tr><td>N. fischeri_a</td><td>SRELQAWDVKDTVHVGFEMI</td><td></td><td>PSIQADADDTAKALMTV</td><td>VISFLNYQADEFMEAWGRL</td></tr><tr><td>T. stipitatus_a</td><td>QAKLETWDVEATDHVGSEIL</td><td></td><td>SGVIADADDTARIILTL</td><td>MVSMLNFQVDRFIEAVIGQE</td></tr><tr><td>P. strigosozonata_a</td><td>NTALQTWDVKSTERIALEMI</td><td></td><td>ANVGPDSDDDTAKALTAL</td><td>VISMVNFQVDEFFDIWVQKH</td></tr><tr><td>S. lacrymans_a</td><td>VDALQGWDLISTERIAFEMI</td><td></td><td>ANVGPDSDDDTAKALTAL</td><td>VISMVNFQVDEFFDLWVQKH</td></tr><tr><td>P. anserina_a</td><td>RELKAFWLDATTLPVGFEML</td><td></td><td>DDTAKGITTL</td><td>VISALAYQVDEFIETSVSQL</td></tr><tr><td>G. luxurians_a</td><td>ASQLNSWDSVSKAERVGFEVTI</td><td></td><td>MDTGEDADDSSKVLTAL</td><td>NMCLVTFQIDEFFDSVVQVQ</td></tr><tr><td>G. luxurians_b</td><td>ASQLNSWDSVSKAERVGFEVTI</td><td></td><td>LDVGEDADDSAKAITAL</td><td>AITLINFQVDEFFDSMVQDQ</td></tr><tr><td>G. luxurians_c</td><td>ASQLNSWDSVSKAERVGFEVTI</td><td></td><td>MDTGEDADDSSKVLTAL</td><td>NMCLVTFQIDEFFDSVVQVQ</td></tr><tr><td>L. palustris_a</td><td>ERKLEHWDPACSDFSGVEFL</td><td></td><td>AGMMHDADDTAKALQSL</td><td>VISMLNFQADAYIEKWGKK</td></tr><tr><td>P. teres_a</td><td>DRQLNAWDVTLTLDHVGFEII</td><td></td><td>PYFAPDVDDTSKTITSL</td><td>TIAMLNFQLDEFMEATAGPA</td></tr><tr><td>A. nidulans_a</td><td>QRLLLEAWDLDGTQDVGFEVI</td><td></td><td>PGILADADDTARVLLTL</td><td>HFSLLVYQVDALMESTAIRM</td></tr><tr><td>A. niger_a</td><td>QQALQQWDVASTLQVGSEVL</td><td></td><td>PGFVPDADDTARALLTQ</td><td>VISLLNYQVDEYMESWALL</td></tr><tr><td>A. sojae_a</td><td>SDALNSWNVEATLHVGFEIL</td><td></td><td>PGFVPDADDTARSLLAL</td><td>VISMLNYQADEYMESVWGSL</td></tr><tr><td>T. reesei_a</td><td>QSALNNWDVNQTLHVGFEML</td><td></td><td>PTLLPDADDTAVSLLAM</td><td>VISMLNYQVDEYMESVAYL</td></tr><tr><td>R. rufulum_b</td><td>NRMLKEWNIGSTDRVGFEVL</td><td></td><td>PSCLPDADDTAKTLSAL</td><td>IISMLGFQVDAYMETWAEH</td></tr><tr><td>G. lozoyensis_a</td><td>NHMLQNWDIKGSDRVGFEVL</td><td></td><td>SNFLPDADDTAKANTVL</td><td>WSMLDYQVDGYMEFVLQEK</td></tr><tr><td>A. fumigatus_a</td><td>QTALKGLDRLLATCTLSVL</td><td></td><td>PNACADADDTAKALVAL</td><td>ELSIGLFAQEDDLMERSLASL</td></tr><tr><td>A. oryzae_a</td><td>DAALKGLNGLLRTCTLPVSL</td><td></td><td>PKACPDADDTAKALAAL</td><td>ELSIGIFQEDELMEKSLVNL</td></tr></tbody></table>
<table><caption>**Table 1** Partial alignment of fungal di-TPS highlighting catalytic domains (Continued)</caption><tr><td>A. sojae_b</td><td>DTALKGLNDLLATCTLPVSV</td><td>PKACPDAADDTAKALAAL</td><td>ELSIGIFQDELMEKSLVNL</td></tr><tr><td>A. nidulans_b</td><td>DAALVRLDLLATSTLTVGL</td><td>PKACPDAADDTAKALIAF</td><td>ELSVGIFQEDEEMEKSLENL</td></tr><tr><td>P. marneffei_a</td><td>QKLLDEWNVDSTLHVGFEIL</td><td>PFVGADADDTATTILVL</td><td>SISIHTDHSGSYYHRSDWTT</td></tr><tr><td>T. stipitatus_b</td><td>QKLLDDWDVKSTLHVGFEIL</td><td>PFVGADADDTATTILVL</td><td>CVSIHTDDSDSYYQRSQWTI</td></tr><tr><td>C. purpurea_a</td><td></td><td>APHTADVDDTAKAILAL</td><td>MLSLLGYQTDEFIEAVAGPA</td></tr><tr><td>M. fijiensis_a</td><td>TTQLRHCDJESTNHIGVELV</td><td>GQGATDADDTAKGLLAL</td><td>FMTLINIQVDEFIEAVITPT</td></tr><tr><td>C. globosum_a</td><td>KEMLHGWDVESTDQVGFEIL</td><td>PGSLPDVDDTAKGIEAL</td><td>VLSVLNFHADEYMEGVWEKH</td></tr><tr><td>N. crassa_a</td><td>QSMLLHWDVATTDQVGFELL</td><td>PGGLVDADDTARAIMAL</td><td>VLSMLNYQVDEYMETAVERD</td></tr><tr><td>N. tetrasperma_a</td><td>QSMLLTWDVATTDQVGFELL</td><td>PGGLVDADDTARAIMAL</td><td>VLSMLNYQVDEYMETAVERD</td></tr><tr><td>S. macrospora_a</td><td>KEMLNKWDVASTDQVGFELL</td><td>PGGLVDADDTARAIMAL</td><td>ILSLLNYQVDEYMETAVERD</td></tr><tr><td>R. rufulum_a</td><td>RYLLGCWDASKATPNNDLEIL</td><td>PGSLQDARYIPHAILAL</td><td>YNDYGSATRDVEEGNLNSLD</td></tr><tr><td>B. fuckeliana_a</td><td></td><td>PGVLPDVDDTSKGLEAL</td><td>VLSILNFHADEYMEGIIERH</td></tr><tr><td>X. parietina_a</td><td>RRMLSSWDATSSDDVGEFVL</td><td>PSVLDESDDTGSSIYIL</td><td>VFSALLYDFDHYMEDIAGF</td></tr><tr><td>Z. passerinii_a</td><td>ASQLASWNDVEITNHIGVEL</td><td>PACSRTMMMDQHANGSWE</td><td>YLIRDHSSTMNSRDSA</td></tr></table>

This is of particular interest, as it shows that, on several occasions, barriers of the cell walls as well as incorrect intron splicing between unrelated eukaryotes have been overcome. These cross-species events between plants and fungi occurred at different times during evolution, some of them between basal lineages (prior to fungal *Dikarya* and at the beginning of land plants), and others more recently (during early diversification of Pezizomycotina and following bryophyte differentiation). The number of proposed HGT is clearly incorrectly estimated, since their identification depends not only on the stringency of the method used but also on subsequent events. For example, HGT can occasionally be overestimated due to lineage specific gene loss, as highlighted by the case of the pleiotropic drug resistance ABC transporter family that is found in plants and fungi, but is absent in animals. These are thought to have been present in the ancestors of all three kingdoms, but then to have been lost in the animal lineage. It is thus likely that the ABC transporter genes were not transferred horizontally between plants and fungi [41]. In our study, similarity searches using di-TPSs show that, until now, di-TPS genes have not been detected in most eukaryotic genomes, with the notable exception of plants and fungi. In addition, all higher plants have gibberellins, whereas only a few fungi appear to contain di-TPS genes. This strongly suggests that the fungal di-TPSs we studied were not present in the ancestor of eukaryotes and later lost in most fungal lineages but that they were actively acquired by HGT from other organisms such as plants. The fact that: 1) only a few fungal species contain di-TPSs (Fig. 1); 2) some species contain several di-TPSs (Table 2); 3) some of these di-TPSs are intronless (Additional file 2); and 4) the phylogenetic relationships of fungal di-TPSs are not in agreement with

inferred phylogenetic relationships between species and higher order lineages of fungi (Figs. 1 and 2), indicates that since the original HGT, these genes have been actively transferred and duplicated and have evolved in fungi.

One could therefore ask whether the plant di-TPSs originally transferred to fungi as a functional di-TPS/ GGPPS/P450 gene cluster encoding an end product or whether the selection for production of the diterpene secondary metabolites resulted in the formation of clusters comprising di-TPS, GGPPS and P450 genes. Until recently, it was thought that genes for plant metabolic pathways were not clustered, and it would have been logical to exclude the possibility of a di-TPS cluster HGT from plants to fungi. However, five clusters of genes for plant metabolic pathways have recently been discovered, and all of these are implicated in the synthesis of defence compounds [42]. Two of these clusters, the rice momilactone and phytocassane clusters, even correspond to diterpene biosynthesis genes [43]. The possibility of a past di-TPS cluster HGT should thus be discussed in this context. Part of this scenario can be answered by considering the GGPPS genes found in fungi containing di-TPS. Figure 4 shows that these genes have more similarities to other fungal GGPPSs and that they are not closely related to plant GGPPSs. Another component of this scenario can be found by considering the well-studied gibberellic acid (GA) biosynthetic pathway as a model. In land plants, this pathway is based on further modifications by P450s of ent-kaurene produced by the di-TPS [20, 44]. Plant P450s catalyse a wide variety of monooxygenation/hydroxylation reactions in primary and secondary metabolisms, including GA biosynthesis. Plant genomes contain a large number of P450, and it is
<table><caption>**Table 2** Fungal specie lifestyles and number of GGPPS and P450 genes clustered next to each di-TPS</caption><thead><tr><th>Fungus</th><th></th><th>di-TPS</th><th>GGPPS</th><th>P450</th></tr></thead><tbody><tr><td rowspan="2">Agaricomycetes</td><td>Serpula lacrymans var. lacrymans S7.3<sup>b</sup></td><td>S. lacrymans_a</td><td>1</td><td>4</td></tr><tr><td>Punctularia strigosozonata<sup>b</sup></td><td>P. strigosozonata_a</td><td>1</td><td>3</td></tr><tr><td>Lecaronomycètes</td><td>Xanthoria parietina<sup>c</sup></td><td>X. parietina_a</td><td>1</td><td>3</td></tr><tr><td rowspan="10">Dothideomycetes</td><td>Rhytidhysteron rufulum<sup>a,b</sup></td><td>R. rufulum_a</td><td>0</td><td>0</td></tr><tr><td>Rhytidhysteron rufulum<sup>a,b</sup></td><td>R. rufulum_b</td><td>0</td><td>2</td></tr><tr><td>Pyrenophora teres f. teres 0-1<sup>a</sup></td><td>P. teres_a</td><td>1</td><td>4</td></tr><tr><td>Bipolaris sorokiniana ND90Pr<sup>a</sup></td><td>B. sorokiniana_a</td><td>0</td><td>2</td></tr><tr><td>Bipolaris maydis C5<sup>a</sup></td><td>B. maydis_a</td><td>0</td><td>3</td></tr><tr><td>Phoma betae<sup>a</sup></td><td>P. betae_a</td><td>1</td><td>2</td></tr><tr><td>Phaeosphaeria nodorum SN15<sup>a</sup></td><td>P. nodorum_a</td><td>0</td><td>1</td></tr><tr><td>Mycosphaerella fijiensis<sup>a</sup></td><td>M. fijiensis_a</td><td>0</td><td>0</td></tr><tr><td>Mycosphaerella graminicola IPO323<sup>a</sup></td><td>M. graminicola_a</td><td>1</td><td>3</td></tr><tr><td>Zymoseptoria passerinii SP63<sup>a</sup></td><td>Z. passerinii_a</td><td>0</td><td>0</td></tr><tr><td>Sphaceloma manihoticola<sup>a</sup></td><td>S. manihoticola_a</td><td>1</td><td>3</td></tr><tr><td rowspan="2">Leotiomycetes</td><td>Botryotinia fuckeliana B05.10<sup>a</sup></td><td>B. fuckeliana_a</td><td>0</td><td>2</td></tr><tr><td>Glarea lozoyensis ATCC 20868<sup>b</sup></td><td>G. lozoyensis_a</td><td>0</td><td>2</td></tr><tr><td rowspan="11">Eurotiomycetes</td><td>Aspergillus fumigatus A1163<sup>b,d</sup></td><td>A. fumigatus_a</td><td>0</td><td>1</td></tr><tr><td>Aspergillus nidulans FGSC A4<sup>b</sup></td><td>A. nidulans_a</td><td>0</td><td>2</td></tr><tr><td>Aspergillus nidulans FGSC A4<sup>b</sup></td><td>A. nidulans_b</td><td>0</td><td>1</td></tr><tr><td>Aspergillus nidulans FGSC A4<sup>b</sup></td><td>A. nidulans_c</td><td>1</td><td>2</td></tr><tr><td>Aspergillus niger CBS 513.88<sup>b,d</sup></td><td>A. niger_a</td><td>0</td><td>2</td></tr><tr><td>Aspergillus oryzae RIB40<sup>b</sup></td><td>A. oryzae_a</td><td>1</td><td>1</td></tr><tr><td>Aspergillus sojae NBRC 4239<sup>b</sup></td><td>A. sojae_a</td><td>1</td><td>2</td></tr><tr><td>Aspergillus sojae NBRC 4239<sup>b</sup></td><td>A. sojae_b</td><td>1</td><td>1</td></tr><tr><td>Neosartorya fischeri NRRL 181<sup>b,d</sup></td><td>N. fischeri_a</td><td>0</td><td>0</td></tr><tr><td>Talaromyces stipitatus ATCC 10500<sup>b</sup></td><td>T. stipitatus_a</td><td>0</td><td>2</td></tr><tr><td>Talaromyces stipitatus ATCC 10500<sup>b</sup></td><td>T. stipitatus_b</td><td>0</td><td>3</td></tr><tr><td>Penicillium marneffei ATCC 1822<sup>b,d</sup></td><td>P. marneffei_a</td><td>0</td><td>3</td></tr><tr><td>Arthroderma benhamiae CBS 112371<sup>d</sup></td><td>A. benhamiae_a</td><td>0</td><td>4</td></tr><tr><td>Arthroderma otae CBS 113480<sup>d</sup></td><td>A. otae_a</td><td>0</td><td>4</td></tr><tr><td>Trichophyton equinum CBS 127.97<sup>b,d</sup></td><td>T. equinum_a</td><td>0</td><td>2</td></tr><tr><td>Trichophyton rubrum CBS 118892<sup>b,d</sup></td><td>T. rubrum_a</td><td>0</td><td>3</td></tr><tr><td>Trichophyton tonsurans CBS 112818<sup>b,d</sup></td><td>T. tonsurans_a</td><td>0</td><td>2</td></tr><tr><td>Trichophyton verrucosum HKI 0517<sup>b,d</sup></td><td>T. verrucosum_a</td><td>0</td><td>3</td></tr><tr><td rowspan="11">Sordariomycetes</td><td>Phomopsis amygdali<sup>a</sup></td><td>P. amygdali_a</td><td>1</td><td>5</td></tr><tr><td>Phomopsis amygdali<sup>a</sup></td><td>P. amygdali_b</td><td>1</td><td>1</td></tr><tr><td>Magnaporthe oryzae 70-15<sup>a</sup></td><td>M. oryzae_a</td><td>1</td><td>1</td></tr><tr><td>Magnaporthe oryzae 70-15<sup>a</sup></td><td>M. oryzae_b</td><td>0</td><td>2</td></tr><tr><td>Trichoderma reesei QM6a<sup>b</sup></td><td>T. reesei_a</td><td>0</td><td>2</td></tr><tr><td>Fusarium fujikuroi<sup>a</sup></td><td>F. fujikuroi_a</td><td>1</td><td>4</td></tr><tr><td>Fusarium proliferatum<sup>a</sup></td><td>F. proliferatum_a</td><td>1</td><td>4</td></tr><tr><td>Aciculosporium take<sup>a</sup></td><td>A. take_a</td><td>0</td><td>0</td></tr><tr><td>Claviceps purpurea 20.1<sup>a</sup></td><td>C. purpurea_a</td><td>0</td><td>3</td></tr><tr><td>Claviceps purpurea 20.1<sup>a</sup></td><td>C. purpurea_b</td><td>0</td><td>3</td></tr></tbody></table>
<table><caption>**Table 2** Fungal specie lifestyles and number of GGPPS and P450 genes clustered next to each di-TPS (Continued)</caption><tr><td>Epichloe amarillans E57<sup>a</sup></td><td>E. amarillans_a</td><td>1</td><td>2</td></tr><tr><td>Epichloe typhina E5819<sup>a</sup></td><td>E. typhina_a</td><td>1</td><td>2</td></tr><tr><td>Podospora anserina S mat+<sup>b</sup></td><td>P. anserina_a</td><td>0</td><td>2</td></tr><tr><td>Neurospora crassa OR74A<sup>b</sup></td><td>N. crassa_a</td><td>0</td><td>1</td></tr><tr><td>Neurospora tetrasperma<sup>b</sup></td><td>N. tetrasperma_a</td><td>0</td><td>1</td></tr><tr><td>Sordaria macrospora k-hell<sup>b</sup></td><td>S. macrospora_a</td><td>0</td><td>1</td></tr><tr><td>Chaetomium globosum CBS 148.51<sup>b,d</sup></td><td>C. globosum_a</td><td>0</td><td>3</td></tr><tr><td>Eutypa lata<sup>a</sup></td><td>E. lata_a</td><td>0</td><td>1</td></tr></table>

In superscript: plant parasites<sup>a</sup>, saprophytes<sup>b</sup>, lichens<sup>c</sup>, animal parasites<sup>d</sup>, plant symbionts<sup>e</sup>

estimated that they account for approximately 1 % of genes [44, 45]. It is noteworthy that P450 fungal genes are not closely related to plant P450 genes, since they share only 10-15 % of the amino-acid identity [29]. It is also noteworthy that the GA biosynthetic genes are clustered on a single chromosome in fungi, whereas they are randomly located on multiple chromosomes in plants [20, 44]. Furthermore, the intermediates in GA biosynthesis are different in plants and fungi. This strongly suggests that plants and fungi have evolved their complex GA biosynthesis pathways independently [29] and that di-TPSs were not acquired from the plants as a whole cluster set.

The discovery of most of the di-TPS fungal genes in clusters with P450 and, to a lesser extent, GGPPSs, is not surprising, since a notable characteristic of fungal genomes is that genes involved in successive steps of a metabolic pathway are often physically linked or clustered [46]. However, in the present study Table 2 shows that in Ascomycota these genes are not clustered with a unique organization, suggesting several mobilization modalities of P450 and GGPPS genes following insertion of the di-TPS gene in the fungal genome. A similar event has already been described in fungi, such as in the case of the galactose utilization pathway. It has been demonstrated that in two cases, fungal galactose gene clusters originated independently (once within Ascomycota and once within Basidiomycota) [47].

The presence of di-TPSs in Basidiomycota was not anticipated, since the literature indicates that they appear to be present only in some branches of the Ascomycota. The similarity searches and the phylogeny based on Basidiomycota GGPPS sequences clearly show that these genes share a closer similarity with Ascomycota than with plant GGPPS sequences, suggesting a fungal origin for these di-TPS. The analysis of the genomic context of the di-TPSs in *S. lacrymans* and *P. strigosozonata* highlighted a putative cluster containing one GGPPS and several P450 genes. A comparison between *S. lacrymans* and *P. strigosozonata* genes shows a rather high degree of similarity between both clusters, with the respective

ortholog sequences from *S. lacrymans* and *P. strigosozonata* clusters usually being the first to occur in similarity searches. This cluster was found in *P. strigosozonata* in a more compacted form than in *S. lacrymans*. The GGPPS tree shows that the GGPPS genes from *S. lacrymans* and *P. strigosozonata* clusters branch with none of the Basidiomycota GGPPSs, contrary to the other three GGPPS genes found in *S. lacrymans* and *P. strigosozonata* genomes. This strongly suggests a common, or at least similar, origin for both di-TPS/GGPPS/P450 clusters found in two distantly related Basidiomycota species. It is therefore possible to speculate that the di-TPS genes may have been transferred from Ascomycota to Basidiomycota. An Ascomycota pathogen infecting a wide array of other fungal species appears to be the easiest way to imagine such a gene transfer. There are examples of such mycoparasites in each of the main Ascomycota groups found to contain di-TPS in our study, for example, *Gliocladium roseum*, *Trichoderma* sp. (Sordariomycetes), *Penicillium restrictum* (Eurothiomycetes) [48–50].

It was however surprising to find that the different di-TPS sequences did not group in the tree according to the fungal phylogeny (Lecaronomycetes, Dothideomycetes, Leotiomycetes, Eurotiomycetes, Sordariomycetes) but rather exhibited some similarities depending on lifestyle, with a group comprising mainly plant-associated fungi (pathogens or symbionts) (Fig. 2). On the other hand, sequences from saprophytic and animal pathogenic fungi are scattered into three different paraphyletic groups, and no clear trend can be identified at this level.

Gibberellins have been extensively studied in *F. fujikuroi* since 1926, and this fungus can be considered as a model of a fungal pathogen using a diterpene hormone. A functional gibberellin synthase cluster contains seven genes encoding: a bifunctional *ent*-copalyl and *ent*-kaurene diphosphate synthase, a geranylgeranyl diphosphate synthase, four P450 monooxygenases and a desaturase [51, 52]. Although gibberellin synthase cluster regulation is linked to the nitrogen metabolism and the enhancement of fungal growth, this can also have an impact on infection efficiency. This production of gibberellin induces abnormal
<figure><img src="image_5.png" /><figcaption>**Fig. 3** (See legend on next page.)</figcaption></figure>
(See figure on previous page.)

**Fig. 3** Comparison of gene content and organization in putative diterpene biosynthetic gene clusters, assumed to have undergone HGT between fungi for A. niger_a and T. reesei_a, P. strigosozonata_a and S. lacrymans_a, A. nidulans_a and C. globosum_a, R. rufulum_b and G. lozoyensis_a. The direction and relative size of genes are indicated by arrows; the upper scale indicates the cluster size in bp. The function of putative syntenic genes is: ABH (Alpha/beta hydrolase), di-TPS (diterpene synthase), GGPPS (geranylgeranyl diphosphate synthase), GST (glutathione S-transferase), MSF (major facilitator superfamily transporter), PHT11 (integral membrane protein PTH11-like protein), P450 (Cytochrome P450), UMTAM (S-adenosylmethionine-dependent methyltransferase)

elongation of plants, early dry up of the leaves, lesions on the roots and adventitious roots, and only partially filled, sterile, or empty seeds. Several other fungi belonging to Ascomycota and Basidiomycota have since been identified as able to produce gibberellins involved in pathogenic or, more interestingly, plant growth promoting interactions [29]. Also, our study often detected di-TPS genes in plant fungal pathogens (Fig. 1) that belong not only to the Ascomycota (Leotiomycetes, Sordariomycetes and Dothideomycetes) but also to the Basidiomycota (*M. perniciosa*).

Surprisingly, the potentially hormonal role of fungal diterpenes in a symbiotic or endophytic fungal/plant relationship [53–55], as well as their antimicrobial properties against bacteria and other fungi [56, 57], have rarely been described in the literature. This trend was confirmed in the present study, since a search for di-TPS-containing fungi found only *Xanthoria parietina*, *Fusarium proliferatum*, *Epichloe typhina* and *E. amarillans* as possible examples of symbionts or endophytes. *X. parietina* does not interact with a vascular plant, but lives as a lichen with green algae belonging to the genus *Trebouxia*. *F. proliferatum* is a plant pathogen, although some strains are mutualistic with orchids [58], *E. typhina* and *E. amarillans* produce the choke disease in grasses, but are also considered as endophytes.

It is noteworthy that diterpenes also appear to have found a hormonal role in fungal development, outside the context of a plant/fungal interaction, since several fungi belonging to the Eurotiales, exhibit conidia and cleistothecia formation linked to diterpene synthesis [41, 59, 60]. The saprophyte fungus *Aspergillus nidulans* contains two functional gene clusters involved in the synthesis of the diterpenoid metabolite *ent*-pimara-8(14),15-diene. Overproduction of this metabolite affects the formation of conidia (decrease in number) and cleistothecia (increase in number) [41, 59]. This cluster of eight co-regulated genes identified by Bromann et al. [59] contains a GGPPS, an HMG-CoA reductase, a di-TPS, a translation elongation factor ©, a short-chain dehydrogenase, a hypothetical protein somewhat similar to methyltransferase, a cytochrome P450, and a Zn(II)₂Cys₆-type transcription factor. We found a very similar cluster for the di-TPS containing Eurotiales *Neosartorya fischeri* but not for other Eurotiales (*A. fumigatus*, *A. niger*, *A. oryzae*, *A.*

*sojae*, *Talaromyces stipitatus*, *Penicillium marneffei*). Furthermore, di-TPSs in the latter Eurotiales species did not group in the phylogenetic tree with *A. nidulans* and *N. fischeri* di-TPSs (Fig. 2), thus indicating that they are not as closely related as the enzyme identified by Bromann et al. [59].

Not all di-TPSs synthesize diterpenes with hormonal activity. Toxic diterpenes have also been described in fungi with various lifestyles (plant or animal parasites, endophytes, or saprophytes) and in plants. Fungal toxic diterpenes include various diterpenes such as indole diterpenes (paspaline, paxillin, shearinines, paspalitrems, terpendoles, penitrems, lolitrems, janthitrems, and sulpinines [61–63]), aphidicolins, and fusicoccins. Each species can produce a family of diterpenes with various toxic effects, such as neurotoxins acting: through modulation of the central nervous system and inhibition of the potassium channel in the peripheral nervous system (paxilline from *Penicillium paxilli*, lolitrems from *Neotyphodium lolii*, aflatrems from *Aspergillus flavus*), as an inhibitor of DNA polymerase (aphidicolin from *Phoma betae* [60, 64]), and as an activator of Membrane H⁺-ATPase from higher plants (methyl phomopsenonate and phytotoxin fusicoccins from *Phomopsis amygdali* [65, 66]).

So far, diterpenes have been described as acting either as hormones, or as toxic molecules. In the present study, fungal di-TPSs have been found in species with lifestyles (plant or animal parasites, saprophytes, and symbionts) suiting both of these biological functions. The good conservation of Mg²⁺-PPi binding, H⁺-initiated cyclization and Mg²⁺-ionization/cyclization also suggests that these enzymes are functional in the fungal cell. It is therefore expected that each of the di-TPS genes we studied has been under pressure to evolve in such a way as to better accommodate one of these functions, which could explain why our di-TPS tree (Fig. 2) does not reflect the fungal phylogenetic position but rather the fungal lifestyle. This function-related tree also suggests that in general, when fungi have no di-TPS in their genomes, this is not the result of the loss of an original di-TPS gene. It is however noteworthy that there is evidence that gene loss can also occasionally contribute to this distribution [67]. Fungal species with di-TPS genes are more likely to have obtained their di-TPS gene(s) by an ongoing HGT process,
<figure><img src="image_6.png" /></figure>

**Fig. 4** Phylogenetic tree performed on GGPPS. The fungal and plant GGPPS used are listed in Additional file 8. GGPPS not in di-TPS clusters are black in colour. GGPPS in di-TPS clusters are in red for Basidiomycota, blue for Eurotiomycetes, green for Dothideomycetes and purple for Sordariomycetes. The tree is rooted at the level of the plant clade and sequences belonging to the same monophyletic group are clustered together. The branch lengths that had to be reduced for the readability of the figure are striped. The entire phylogenetic tree is described in Additional file 9

occurring in fungi. One approach to the investigation of this possibility involved searching for gene synteny in unrelated species, resolved within the same clade in the di-TPS tree. The di-TPS tree was indeed found to contain examples of such a potential HGT.

## Conclusions

The analysis of fungal di-TPS traced the evolutionary history of an ancient HGT event of a single di-TPS gene from a plant to an Ascomycota. Genome analysis demonstrated that in fungi, these di-TPSs allowed for the formation of clusters of di-TPS, GGPPS and P450 genes that presumably originated independently and were characterized by different mobilization modalities. We observed that members of these clusters can also be transferred between *Dikarya* fungal species. Our analysis leads to a better understanding of how trabfers of di-TPS have enabled fungi to produce diterpenes acting as hormones or toxins. The variable presence of these compounds has implications for the fungal development and pathogenicity.

## Methods

### Sequence database searches

We used BLAST (BLASTP, TBLASTN and PSI-BLAST [68]) and PipeAlign [69] searches to recover information related to the presence or absence of putative di-TPS genes in the available fungal genomes. We started with the sequence from *Fusarium fujikoroi* as a seed and extended our search with di-TPS sequences from each fungal class using an expectation value threshold of $E \le 10^{-6}$. Searches were performed in the NCBI (http://www.ncbi.nlm.nih.gov/), UniProtKB (http://www.-uniprot.org/), Genoscope (http://www.genoscope.cns.fr/spip/) and JGI (http://www.jgi.doe.gov/) collections [date last accessed 25 January 2014], until no new sequences were found (Additional file 1). Sequences encoding incomplete di-TPSs were discarded. Sesquiterpenoid trichothecene synthase-like genes found in several fungi and having weak similarities to di-TPSs were also discarded. Synteny conservation was determined by BLASTP and TBLASTN searches ($E \le 10^{-6}$), using approximately 10 genes upstream and downstream from the di-TPS genes to find their orthologs. Completed fungal genomes were numbered using NCBI, JGI and 1,000 fungal genome project
databases. For every di-TPS gene accession (Additional file 1), the systematic classification found on the associated web page was used for the fungal phylogenetic tree (Fig. 1).

We used the *Saccharomyces cerevisiae* Geranyl Geranyl Diphosphate Synthase (GGPPS) sequence as a seed to search against the NR database (ncbi http://www.ncbi.nlm.nih.gov), in order to recover additional putative fungal GGPPSs. As a consequence of the large number of Ascomycota sequences available, we selected 100 sequences representative of the diversity of Ascomycota, and used all available GGPPS sequences from the other fungal groups, as well as GGPPS from four plant species (Additional file 3). Putative plant GGPPSs were selected from BLAST searches performed against *Physcomitrella patens*, *Selaginella moellendorffii* and *Zea mays* genome sequences using the *Arabidopsis thaliana* sequence [70] as a seed.

## Gene predictions and intron/exon architecture

Absent or improper gene predictions for the following di-TPS sequences were corrected using the fgenesh SoftBerry software [71]: *Punctularia strigosozonata* (5-scaffold join 1540005-1539226; 1539170-1539050; 1538993-1538707; 1538636-15370098; 1537040-1536985; 1536849-1536666), *Rhystidhysteron rufulum* (NODE_2525 join 7075-7866; 7917-8037; 8091-8371; 8426-10078), *R. rufulum* (NODE_28 join 20355-21152; 21200-21589; 22055-22572; 22828-23059; 23574-23893; 23957-24170; 24225-24307; 24382-24475). The exon and intron positions were formatted into a gff file. The domain architectures (Mg²⁺-PPi binding, H⁺-initiated cyclization and Mg²⁺-ionization/cyclization) of all analysed sequences were predicted from a multiple sequence alignment of the selected sequence using multalin [72] and *Abies grandis* abietadiene synthase as references. Geneious 6.1.6 (Biomatters, www.geneious.com) was used to view the gene structure and generate Table 1.

## Phylogenetic analyses

The same general procedure was used for both datasets (di-TPS and GGPPS). Multiple sequence alignments were computed with Muscle [73], with all default parameters. The alignments were then trimmed using Gblocks [74], with the most relaxed set of parameters. The number of sites retained after trimming was equal to 521 for di-TPS and 186 for GGPPS. Before computing the phylogenies a substitution model selection was realized using ProtTest [75] and the Akaike Information Criterion (AIC). The models selected were LG + $\phi_4$ + I + F in the case of di-TPS and LG + $\phi_4$ in the case of GGPPS. The phylogenies were computed using the SeaView implementation of PhyML [76]. Alpha parameters of the two Gamma distributions, as well as amino acid

invariant frequencies and equilibrium frequencies were estimated using PhyML. Branch supports were determined under PhyML with the Shimodaira-Hasegawa (SH) likelihood test. Branch supports were determined under PhyML with the approximate Likelihood Ratio Test (aLRT) non-parametric branch support based on a Shimodaira-Hasegawa-like procedure.

## Additional files

**Additional file 1: List of di-TPSs for fungi and 11 plant species.**
(DOCX 18 kb)

**Additional file 2: Fungal exons (blue) and introns (black) in putative di-TPS genes.** Fungal species are arranged in the same order as in Table 1 and the subgroups are specified in superscript: Agaricomycetes (Ag), Lecaronomycetes (Lecar), Eurotiomycetes (Eu), Leotiomycetes (Leotio), Sordariomycetes (So), Dothideomycetes (Do). (PDF 2136 kb)

**Additional file 3: Biosynthesis of diterpenoids with symbols indicating the chemical products (rectangles), enzymes (shaded rectangles) and types of catalysed reaction (brackets).** (PDF 45 kb)

**Additional file 4: Best BLAST hits for the genes involved in the putative *A. niger* and *T. reesei* gene cluster** [73]. (DOCX 13 kb)

**Additional file 5: Best BLAST hits for the genes involved in the putative *S. lacrymans* and *P. Strigosozonata* gene cluster.** (DOCX 14 kb)

**Additional file 6: Best BLAST hits for the genes involved in the putative *B. fuckeliana* and *C. globosum* gene cluster.** (DOCX 15 kb)

**Additional file 7: Best BLAST search performed on the genes involved in the putative *R. rufulum* and *G. lozoyensis* gene cluster.**
(DOCX 11 kb)

**Additional file 8: List of GGPPSs containing di-TPS Ascomycota and Basidiomycota species, early diverging lineages of fungi, and three plant species.** (DOCX 23 kb)

**Additional file 9: Entire Phylogenetic tree generated with GGPPS.**
(PDF 11 kb)

## Competing interests

The authors have no competing interests.

## Authors' contributions

MF designed the study and searched databases; MF, VL, CR, GP performed the phylogenetic analyses and interpreted the data. All of the authors contributed to the writing of the manuscript, and reviewed and approved its final version.

## Acknowledgements

We thank Nicolas Lartillot for careful re-reading of the manuscript. All computations have been performed using the LBBE/PRABI cluster.

## Author details

¹Université de Strasbourg, INRA, Inst Natl Recherche Agron, Métab Second Vigne, Unit Mixte Recherche Santé Vigne & Qual Vins, 28 rue de Herrlisheim, F-68021 Colmar, France. ²Université de Strasbourg, CNRS, FRE 2326, Institut de Biologie Moléculaire et Cellulaire du CNRS, UPR 9002, 15 rue René Descartes, F-67084 Strasbourg, France. ³Universite Claude Bernard - Lyon 1, 43 bd. du 11 Novembre 1918, Laboratoire de Biometrie et Biologie Evolutive, UMR CNRS 5558, F-69622 Villeurbanne, France.

Received: 11 September 2014 Accepted: 12 October 2015

Published online: 19 October 2015

## References

1. De Carvalho MO, Loreto ELS. Methods for detection of horizontal transfer of transposable elements in complete genomes. *Gen Mol Biol. 2012;35:1078-84.

2. Fitzpatrick DA. Horizontal gene transfer in fungi. *FEMS Microbiol Lett. 2012;329:1-8.
3. Richards TA. Genome evolution: horizontal movements in the fungi. Curr Biol. 2011;21:R166-8.

4. Hall C, Brachat S, Dietrich FS. Contribution of horizontal gene transfer to the evolution of *Saccharomyces cerevisiae*. Eukaryot Cell. 2005;4:1102-15.

5. Garcia-Vallve S, Romeu A, Palau J. Horizontal gene transfer of glycosyl hydrolases of the rumen fungi. Mol Biol Evol. 2000;17:352-61.

6. 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:e28231.

7. Novo M, Bigey F, Beyne E, Galeote V, Gavory F, Mallet S, et al. Eukaryote-to-eukaryote gene transfer events revealed by the genome sequence of the wine yeast *Saccharomyces cerevisiae* EC1118. Proc Natl Acad Sci USA. 2009;106:16333-8.

8. Friesen TL, Stukenbrock EH, Liu ZH, Meinhardt S, Ling H, Faris JD, et al. Emergence of a new disease as a result of interspecific virulence gene transfer. Nat Genet. 2006;38:953-6.

9. Ma LJ, Van Der Does HC, Borkovich KA, Coleman JJ, Daboussi MJ, Di Pietro A, et al. Comparative genomics reveals mobile pathogenicity chromosomes in *Fusarium*. Nature. 2010;464:367-73.

10. Richards TA, Soanes DM, Foster PG, Leonard G, Thornton CR, Talbot NJ. Phylogenomic analysis demonstrates a pattern of rare and ancient horizontal gene transfer between plants and fungi. Plant Cell. 2009;21:1897-911.

11. Jaramillo VDA, Vargas WA, Sukno SA, Thon MR. Horizontal transfer of a subtilisin gene from plants into an ancestor of the plant pathogenic fungal genus *Colletotrichum*. Plos ONE. 2013;8(3):e59078.

12. Nikolaidis N, Doran N, Cosgrove DJ. Plant expansins in bacteria and fungi: evolution by horizontal gene transfer and independent domain fusion. Mol Biol Evol. 2013;31:376-86.

13. Slot JC, Hibbett DS. Horizontal transfer of a nitrate assimilation gene cluster and ecological transitions in fungi: a phylogeny study. PloS ONE. 2007;2:e1097.

14. MacMilan J. Occurrence of gibberellins in vascular plants, fungi and bacteria. J Plant Growth Regul. 2002;20:387-442.

15. Morrone D, Chambers J, Lowry L, Kim G, Anterola A, Bender K, et al. Gibberellin biosynthesis in bacteria: separate *ent*-copalyl diphosphate and *ent*-kaurene synthases in *Bradyrhizobium japonicum*. FEBS Lett. 2009;583:475-80.

16. Cao R, Zhang Y, Mann FM, Huang C, Mukkamala D, Hudock MP, et al. Diterpene cyclases and the nature of the isoprene fold. Proteins. 2010;15:2417-32.

17. Smanski MJ, Peterson RM, Huang SX, Shen B. Bacterial diterpene synthases: new opportunities for mechanistic enzymology and engineered biosynthesis. Cur Opin Chem Biol. 2012;16:132-41.

18. Aya K, Hiwatashi Y, Kojima M, Sakakibara H, Ueguchi-Tanaka M, Hasebe M, et al. The gibberellin perception system evolved to regulate a pre-existing GAMYB-mediated system during land plant evolution. Nat Commun. 2011;2:544.

19. Hayashi KI, Horie K, Hiwatashi Y, Kawaide H, Yamaguchi S, Hanada A, et al. Endogenous diterpenes derived from *ent*-kaurene, a common gibberellins precursor, regulate protonema differentiation of the moss *Phycomitrella patens*. Plant Physiol. 2010;153:1085-97.

20. Yamagushi S. Gibberellin metabolism and its regulation. Annu Rev Plant Biol. 2008;59:225-51.

21. Keeling CI, Bohlman J. Diterpene resin acids in conifers. Phytochemistry. 2006;67:2415-23.

22. Bohlmann J, Meyer-Gauen G, Croteau R. Plant terpenoid synthases: molecular biology and phylogenetic analysis. Proc Natl Acad Sci USA. 1998;95:4126-33.

23. Withers ST, Keasling JD. Biosynthesis and engineering of isoprenoid small molecules. Appl Microbiol Biotechnol. 2007;73:980-90.

24. Rensing SA, Lang D, Zimmer AD, Terry A, Salamov A, Shapiro H, et al. The *Phycomitrella* genome reveals evolutionary insights into the conquest of land by plants. Science. 2008;319:64-9.

25. Hayashi K, Kawaide H, Notomi M, Sakigi Y, Matsuo A, Nozaki H. Identification and functional analysis of bifunctional *ent*-kaurene synthase from the moss *Phycomitrella patens*. FEBS Lett. 2006;580:6175-81.

26. Chen F, Tholl D, Bohlmann J, Pichersky E. The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J. 2011;66:212-29.

27. Li G, Köllner TG, Yin Y, Jiang Y, Chen H, Xu Y, et al. Nonseed plant *Selaginella moellendorffii* has both seed plant and microbial types of terpene synthases. Proc Natl Acad Sci USA. 2012;109:14711-5.

28. Hedden P, Phillips AL, Rojas MC, Carrera E, Tudzynski B. Gibberellin biosynthesis in plants and fungi: a case of convergent evolution? J Plant Growth Regul. 2002;20:319-31.

29. Bömke C, Tudzynski B. Diversity, regulation, and evolution of the gibberellin biosynthetic pathway in fungi compared to plants and bacteria. Phytochemistry. 2009;70:1876-93.

30. Proctor RH, Mc Cormick SP, Alexander NJ, Desjardins AE. Evidence that a secondary metabolic biosynthetic gene cluster has grown by gene relocation during evolution of the filamentous fungus *Fusarium*. Mol Microbiol. 2009;74:1128-42.

31. Mondego JMC, Carazzolle MF, Costa GGL, Formighieri EF, Parizzi LP, Rincones J, et al. A genome survey of *Moniliophthora perniciosa* gives new insights into Witches' Broom Disease of cacao. BMC Genomics. 2008;9:548.

32. Zerbe P, Hamberger B, Yuen MMS, Chiang A, Sandhu HK, Madilao LL, et al. Gene discovery of modular diterpene metabolism in nonmodel systems. Plant Physiol. 2013;162:1073-91.

33. Peters RJ, Flory JE, Jetter R, Ravn MM, Lee HJ, Coates RM, et al. Abietadiene synthase from grand fir (*Abies grandis*): characterization and mechanism of action of the 'pseudomature' recombinant enzyme. Biochemistry. 2000;39:15592-602.

34. Toyomasu T, Kawaide H, Ishizaki A, Shinoda S, Otsuka M, Mitsuhashi W, et al. Cloning of a full-length cDNA encoding *ent*-kaurene synthase from *Gibberella fujikoroi*: functional analysis of a bifunctional diterpene cyclase. Biosci Biotechnol Biochem. 2000;64:660-4.

35. Taylor JW, Berbee ML. Dating divergences in the fungal tree of life: review and new analyses. Mycologia. 2006;98:838-49.

36. Lücking R, Pfister DH, Plave ER, Lumbsch HT. Fungi evolved right on track. Mycologia. 2009;101:810-22.

37. Sanderson MJ. Molecular data from 27 proteins do not support a precambrian origin of land plants. AM J Bot. 2003;90:954-6.

38. Willis A, Rodrigues BF, Harris PJC. The ecology of arbuscular mycorrhizal fungi. Crit Rev Plant Sci. 2013;32:1-20.

39. Kemen E, Jones JDG. Obligate biotroph parasitism: can we link genomes to lifestyles? Trends Plant Sci. 2012;17:448-547.

40. Bidartondo MI, Read DJ, Trappe JM, Merckx V, Ligrone R, Duckett JG. The dawn of symbiosis between plants and fungi. Biol Lett. 2013;7:574-7.

41. Crouzet J, Trombik T, Fraysse AS, Boutry M. Organization and function of the plant pleiotropic drug resistance ABC transporter family. FEBS Lett. 2005;580:1123-30.

42. Chu HY, Wegel E, Osbourn A. From hormones to secondary metabolism: the emergence of metabolic gene clusters in plants. Plant J. 2011;66:66-79.

43. Swaminathan S, Morrone D, Wang Q, Fulton B, Peters RJ. CYP76M7 is an *ent*-cassadiene C11(-hydrolase defining a second multifunctional diterpenoid biosynthetic gene cluster in rice. Plant Cell. 2009;21:3315-25.

44. Mizutani M, Ohta D. Diversification of P450 genes during land plant evolution. Annu Rev Plant Biol. 2010;61:291-315.

45. Nelson DR, Schuler MA, Paquette SM, Werck-Reichhart D, Bak S. Comparative genomics of rice and *Arabidopsis*. Analysis of 727 cytochrome P450 genes and pseudogenes from a monocot and a dicot. Plant Physiol. 2004;135:756-72.

46. Gacek A, Strauss J. The chromatin code of fungal secondary metabolite gene clusters. Appl Microbiol Biotechnol. 2012;6:1389-404.

47. Slot JC, Rokas A. Multiple GAL pathway gene clusters evolved independently and by different mechanisms in fungi. Proc Natl Acad Sci USA. 2010;107:10136-41.

48. Cota LV, Maffia LA, Mizubuti ESG, Macedo PEF. Biological control by *Clonostachys rosea* as a key component in the integrated management of strawberry gray mold. Biol Control. 2009;50:222-30.

49. Kubicek CP, Herrera-Estrella A, Seidl-Seiboth V, Martinez DA, Druzhinina IS, Thon M, et al. Comparative genome sequence analysis underscores mycoparasitism as the ancestral life style of *Trichoderma*. Genome Biol. 2011;12:R40.

50. Nicoletti R, De Stephano M. *Penicillium restrictum* as an antagonist of plant pathogenic fungi. Dyn Biochem Process Biotech Mol Biol. 2012;6:61-9.

51. Tudzynski B. Biosynthesis of gibberellins in *Gibberella fujikoroi*: biomolecular aspects. Appl Microbiol Biotechnol. 1999;52:198-310.

52. Tudzynski B. Gibberellin biosynthesis in fungi: genes, enzymes, evolution, and impact on biotechnology. Appl Microbiol Biotechnol. 2005;66:597-611.

53. Schardl CL, Young CA, Hesse U, Amyotte SG, Andreeva K, Calie PJ, et al. Plant-Symbiotic Fungi as Chemical Engineers: Multi-Genome Analysis of the Clavicipitaceae Reveals Dynamics of Alkaloid Loci. PloS Genetics. 2013;9, e1003323.
54. Zhang P, Zhou PP, Yu JL. An endophytic taxol-producing fungus from *Taxus x media*, *Aspergillus candidus* MD3. FEMS Microbiol Lett. 2009;293:155-9.

55. Lopes AA, Pupo MT. Biosynthesis of Aphidicolin Proceeds via the Mevalonate Pathway in the Endophytic Fungus *Nigrospora sphaerica*. J Braz Chem Soc. 2011;22:80-5.

56. De Souza JJ, Vieira IJC, Rodrigues-Filho E, Braz-Filho R. Terpenoids from endophytic fungi. Molecules. 2011;16:10604-18.

57. Pongcharoen W, Rukachaisirikul V, Pongpaichit S, Rungjindamai N, Sakayajoj J. Pimarane diterpene and cytochalasin derivatives from the endophytic fungus *Eutypella scoparia* PSU-D44. J Nat Prod. 2006;69:856-8.

58. Tsavkelova EA, Bömke C, Netrusov AI, Weiner J, Tudzynski B. Production of gibberellic acids by an orchid-associated *Fusarium proliferatum* strain. Fungal Genet Biol. 2008;10:1393-403.

59. Bromann K, Toivari M, Viljanen K, Vuoristo A, Ruohonen L, Nakari-Setälä T. Identification and characterization of a novel diterpene gene cluster in *Aspergillus nidulans*. PLoS ONE. 2012;7:e35450.

60. Roncal T, Cordobes S, Sterner O, Ugalde U. Conidiation in *Penicillium cyclopium* is induced by conidiogenone, an endogenous diterpene. Eukaryot Cell. 2002;1:823-9.

61. Parker EJ, Scott DB. Indole-diterpene biosynthesis in ascomycetous fungi. In: An Z, editor. Handbook of industrial mycology. New York, N.Y.: Marcel Dekker; 2004. p. 405-26.

62. Saikia S, Nicholson MJ, Young C, Parker EJ, Scott B. The genetic basis for indole-diterpene chemical diversity in filamentous fungi. Mycol Res. 2008;112:184-99.

63. Zhang S, Monahan BJ, Tkacz JS, Scott B. Indole-diterpene gene cluster from *Aspergillus flavus*. Appl Environ Microbiol. 2004;70:6875-83.

64. Toyomasu T, Nakaminami K, Mie T, Watanabe K, Ito H, Matsui H, et al. Cloning of a gene cluster responsible for the biosynthesis of diterpene aphidicolin, a specific inhibitor of DNA polymerase (λ). Biosci Biotechnol Biochem. 2004;68:146-52.

65. Toyomasu T, Tsukahara M, Kaneko A, Niida R, Litsuhashi W, Dairi T, et al. Fusicoccins are biosynthesized by an unusual chimera diterpene synthase in fungi. Proc Natl Acad Sci USA. 2007;104:3084-8.

66. Toyomasu T, Kaneko A, Tokiwano T, Kanno Y, Kanno Y, Niida R, et al. Biosynthetic gene-based secondary metabolite screening: a new diterpene, methyl phomopsenonate, from the fungus *Phomopsis amygdali*. J Org Chem. 2009;74:1541-8.

67. Boemke C, Rojas MC, Hedden P, Tudzynski B. Loss of Gibberellin Production in *Fusarium verticillioides* (*Gibberella fujikuroi* MP-A) Is Due to a Deletion in the Gibberellinic Acid Gene Cluster. Appl Environ Microbiol. 2008;74:7790-801.

68. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment research tool. J Mol Biol. 1990;215:403-10.

69. Plewniak F, Blanchetti L, Brelivet Y, Carles A, Chalmel F, Lecompte O, et al. PiPeAlign: A new toolkit for protein family analysis. Nucleic Acids Res. 2003;31:3829-32.

70. Beck G, Coman D, Herren E, Ruiz-Sola MA, Rodriguez-Conception M, Gruissem W, et al. Characterization of the GGPP synthase gene family in *Arabidopsis thaliana*. Plant Mol Biol. 2013;82:393-416.

71. Solovyev V, Kosarev P, Seledsov I, Vorobeyev D. Automatic annotation of eukaryotic genes, pseudogenes and promoters. Genome Biol. 2006;7:S10.

72. Corpet F. Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res. 1998;16:10881-90.

73. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792-7.

74. Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540-52.

75. Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27:1164-5.

76. Gouy M, Guindon S, Gascuel O. SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27:221-4.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central