BACKGROUND
hypocreales is home to a wide array of ecologically diverse fungi.
some are devastating plant pathogens , while others form numerous lineages of both insect pathogens  and mycoparasites   <cit> .
at the divergence of the four most derived families  of hypocreales, there was a major shift away from plant-based nutrition to either insect pathogenesis or fungal parasitism, i.e., mycoparasites  <cit> .
two major lineages of mycoparasites are found within the order, although other mycoparasites exist .
the first and larger of these two lineages is the family hypocreaceae, most notable for mycoparasitic trichoderma spp.
used in biological control of plant pathogenic fungi, and trichoderma reesei e.g.
simmons  <cit> , the industrial workhorse for cellulase production  <cit> .
the second major lineage of mycoparasites, the genus tolypocladium, is nested within the insect pathogenic family, ophiocordycipitaceae  <cit> .
most species of tolypocladium parasitize the truffles of elaphomyces , ectomycorrhizal fungi closely related to aspergillus and penicillium  <cit> .
tolypocladium ophioglossoides  quandt, kepler & spatafora is a commonly collected truffle parasite with a broad geographic distribution throughout many parts of the northern hemisphere  <cit> .
there are, however, a few tolypocladium species that attack insects and rotifers, and based on current multigene phylogenies some of these are inferred to be reversals to insect pathogenesis  <cit> .
one of these is a beetle pathogen, t. inflatum, which was the first source of the immunosuppressant drug, cyclosporin a  <cit> .
evidence from multigene studies has also shown a close phylogenetic relationship between t. ophioglossoides and t. inflatum  <cit> .
secondary metabolism is defined as the synthesis of often bioactive, small molecules that are not essential to the growth of an organism.
genes related to production of secondary metabolites are often clustered together in close proximity within a genome and coregulated  <cit> .
a wide variety of secondary metabolites including the ergot alkaloids, fumonisins, and destruxins, is produced by species of hypocreales .
many of these metabolites are produced by nonribosomal peptide synthetases , which are often large, multi-modular proteins that produce short peptides frequently incorporating non-standard amino acids.
nrps modules are composed of three primary functional domains including adenylation , thiolation , and condensation  domains  <cit> .
due to their high level of amino acid and nucleotide conservation, the a-domains are frequently used to reconstruct the evolutionary histories of these genes  <cit> .
polyketide synthases  are another class of secondary metabolite producing enzymes that are common in fungi and are also modular in nature.
they are related to fatty acid synthases  <cit> , and assemble small bioactive molecules based on acetyl-coa or malonyl-coa subunits  <cit> .
the other major classes of secondary metabolite-producing enzymes are terpene synthases and dimethylallyltryptophan  synthases, both of which have been reported from hypocrealean taxa.
fungal secondary metabolites clusters often include genes required for regulation of expression of the gene cluster and decoration, epimerization, and transport of the mature secondary metabolite  <cit> .
peptaibols, or peptaibiotics, are antibiotic secondary metabolites products produced by very large nrps enzymes .
their name is a derivative of their structure as they are peptides containing the uncommon non-proteinogenic amino acid, α-amino isobutryic acid , and a c-terminal amino ethanol  <cit> .
the presence of aib residues promotes helix formation, and several of these helices form multimeric units that in turn form voltage gated ion channels capable of inserting into cell membranes where they disrupt membrane potential causing leakiness  <cit> .
peptaibols are produced by trichoderma spp.
and other members of hypocreaceae , leading to the proposition that they may play a role in mycoparasitism.
there is at least one empirical study to support this in trichoderma  <cit> .
further studies have found that peptaibols function, along with cell wall degrading enzymes, to synergistically inhibit new cell wall synthesis in fungal prey of tr.
harzianum .
wiest et al.
<cit>  identified and characterized the first peptaibol synthetase nrps modular structure  from tr.
virens along with the  <dig> residue peptaibol product.
since that time, several other peptaibol nrps genes have been identified in tr.
virens and other species  <cit> .
efrapeptins are another class of peptaibiotics described almost exclusively from tolypocladium spp.
that have antifungal and insecticidal properties .
they differ from orthodox peptaibols by the presence of a mitochondrial atpase inhibiting c-terminal “blocking group” n-peptido-1-isobutyl-2-ethylamine  <cit> .
eight of these, named efrapeptins a, and c-i, have been isolated from t. inflatum  <cit> .
numerous genomes from species of the mycoparasitic genus trichoderma  have been sequenced  <cit> , and more recently the genomes of several insect pathogens in hypocreales have been completed  , including t. inflatum, the beetle pathogenic congener of t. ophioglossoides  <cit> .
the genomes of all of these species are rich in secondary metabolite genes and clusters, ranging from  <dig> to  <dig> secondary metabolite gene clusters per genome.
comparisons of gene content and expression are beginning to shed light on mechanisms underlying host specificity and the evolution of primary and secondary metabolism.
in this study, the draft genome of the truffle parasite t. ophioglossoides was generated to compare the gene content and secondary metabolite content of this truffle parasite to those of closely related insect pathogens and more distantly related mycoparasites.
the secondary metabolite potential of t. ophioglossoides is characterized with a focus on understanding the evolution of gene clusters encoding for peptaibiotics.
RESULTS
genome assembly and structure
the draft genome assembly of the t. ophioglossoides cbs  <dig> is approximately  <dig>  megabases  and is assembled on  <dig> scaffolds.
the assembly, with an n <dig> of  <dig>  base pairs , captures the majority of gene space and synteny , in spite of the exclusive use of short read technology.
in addition, the core eukaryotic genes mapping approach   <cit>  identified  <dig> complete and  <dig> complete and partial core eukaryotic genes, estimating the assembly completeness at  <dig>  % or  <dig>  %.
we predict  <dig>  protein-coding genes resulting in  <dig>  protein models, of which  <dig>  have support from rna, making the size and number of protein models similar to that of the beetle pathogen, t. inflatum, which is  <dig>  mb and has  <dig>  protein models.
they both also share high gc contents,  <dig>  % and 58 %, for t. ophioglossoides and t. inflatum, respectively.
these species are very closely related, which is reflected in the genome scale phylogeny  and the large regions of shared synteny .
there are, however, two large scale inversions and some small rearrangements with disagreement, suggesting these species have diverged despite the relatively short branch lengths separating them in the species phylogeny .
the large number of small non-syntenic points in the mummerplot alignment could be due to the reliance on short read technology and the presence of short contigs in both assemblies.table  <dig> genome statistics for tolypocladium ophioglossoides compared to t. inflatum
t. ophioglossoides
t. inflatum
fig.
<dig> maximum likelihood phylogeny of hypocrealean taxa analyzed in this study.
raxml tree created using concatenated alignment of  <dig>  protein clusters identified in the hal pipeline  <cit>  with  <dig>  amino acid positions and  <dig>  distinct alignment patterns.
branches are colored based on host/nutritional association.
both pathogens and saprobes of plant material are colored green.
species that possess peptaibol a-domains are denoted by a blue ‘α’
fig.
<dig> whole genome nucleotide alignment of the tolypocladium ophioglossoides and the t. inflatum genomes.
mummerplot visualization of nucmer alignment of all  <dig> unmasked t. ophioglossoides scaffolds against the reference 101 t. inflatum scaffolds.
red lines represent sequence aligning in the same direction while blue lines represent inversions.
the black arrows point out the major inversions between the two genomes
an overview of secondary metabolites in t. ophioglossoides
the t. ophioglossoides and t. inflatum genomes harbor  <dig> and  <dig> core secondary metabolite genes – nrps, pks, terpene synthase, and dmats – spread across  <dig> and  <dig> secondary metabolite gene clusters, respectively.
similar to the contrast observed between closely related metarhizium species , t. ophioglossoides and t. inflatum differ in the types of core secondary metabolite genes they possess with  <dig> shared between the two species,  <dig> unique to t. inflatum, and  <dig> unique to t. ophioglossoides .
notably, t. ophioglossoides does not contain the nrps gene, or any of the other genes in the sima cluster responsible for the production of cyclosporin a in t. inflatum  <cit> .
however, t. ophioglossoides does share a hypocreales-conserved core set of genes that flank the sima region with t. inflatum.
cyclosporin production has been reported from several species of tolypocladium  <cit> , however, none of the truffle parasites, including t. ophioglossoides, have been demonstrated to produce the compound.
whether possession of the sima gene is a derived character state, or whether there has been a single or multiple losses within the genus, remains unknown and requires further sampling.
secondary metabolites core genes predicted from the t. ophioglossoides genome include  <dig> pkss, two pks-like genes,  <dig> nrpss, six nrps-like genes, three hybrid nrps-pks genes, and four terpenes .
no dmat synthases were identified in the t. ophioglossoides genome.
based on a-domain homology, two putative siderophore synthetases  were among the predicted nrpss .
this is in contrast to t. inflatum, which possesses three putative siderophore synthetases   <cit> .
the entire pseurotin-a precursor synthetase hybrid nrps-pks cluster  was identified in the t. ophioglossoides genome.
pseurotin-a, an antifungal compound described from several aspergillus spp.
, was also recently identified in the genome of m. robertsii  <cit> .
the disjunct distribution of this secondary metabolite cluster raises several questions about the evolutionary mechanisms  that may have led to this distribution.
t. ophioglossoides possesses the destruxins synthetase nrps gene  .
destruxins are known for their insecticidal properties in metarhizium spp.
<cit> , and the entire destruxins synthetase cluster was characterized by wang et al.
in m. robertsii.
homologs of the other essential genes in the destruxins cluster are present in the t. ophioglossoides toph_ <dig> cluster, except for dtxs <dig> , an aspartic decarboxylase responsible for producing β-alanine, one of the amino acids incorporated into destruxins.
there are inversions in this cluster between m. robertsii and t. ophioglossoides as well, including the dtxs <dig> aldo-keto reductase homolog  and an abc transporter , which was not found to be essential in destruxins production in m. robertsii  <cit> .
for these reasons, it is not likely that t. ophioglossoides produces destruxins, but possibly produces another group of related compounds.
the sequenced strain of tr.
virens shares a homolog of the destruxins nrps gene  , but destruxins have not been reported to be produced by that species either  <cit> .
to date, only two secondary metabolites have been reported to be produced by t. ophioglossoides: ophiocordin  and ophiosetin .
ophiosetin is structurally similar to equisetin, an antibiotic with inhibitory activity of hiv- <dig> integrase, and both are produced by nrps-pks hybrid genes.
based on phylogenetic analysis of a-domains, cluster synteny with the equisetin cluster  <cit> , and sequence homology, this study identifies the putative ophiosetin synthetase cluster around the hybrid nrps-pks, toph_ <dig> .
further studies involving transformations and chemical verification and characterization of this cluster will be necessary to confirm this genotype-chemotype linkage.
ophiocordin is a polyketide and no putative gene or gene cluster related to its production was identified here.
except for the two peptaibiotic clusters discussed below, the remaining  <dig> secondary metabolite gene clusters are not yet associated with a specific gene product.
peptaibiotics of tolypocladium
among the wide assortment of secondary metabolite genes and gene clusters, the draft genome of t. ophioglossoides possesses four peptaibiotic nrps genes located within three gene clusters.
phylogenetic analyses of fungal a-domains from a variety of nrps genes known to produce specific products and from whole genome mining of the hypocrealean species in fig.
<dig>  revealed that trichoderma peptaibol a-domains group into three clades, one of which  is well supported in all analyses , the other two clades have strong support with the exception of a few a-domains from one of the t. ophioglossoides peptaibiotic genes  .
all a-domains from the four peptaibiotic genes in t. ophioglossoides fall within these three clades , and representation in the three peptaibol clades is exclusive to a-domains from trichoderma and tolypocladium.
t. inflatum also possesses three peptaibiotic nrps genes.
strikingly, this limits the presence of peptaibiotic a-domains and genes to the two sampled mycoparasitic lineages in hypocreales.fig.
<dig> phylogeny of the three peptaibol a-domain clades.
peptaibol focused, a-domain phylogeny created using raxml, showing the  <dig> major clades.
green branches represent those supported by ≥70 % bootstrap support.
tolypocladium spp.
a-domains are colored in blue and teal, and trichoderma spp.
a-domains are colored in yellow, red, or orange
there are differences between species and gene membership within a-domain clades, however.
for instance, clade  <dig> is enriched in a-domains from the trichoderma peptaibols , while clade  <dig> contains more a-domains from tolypocladium spp.
.
importantly clade  <dig> contains a-domains that encode for incorporation of aib, as well as other a-domains that encode for incorporation of isovaline, leucine, isoleucine, alanine, glycine, valine, and serine .
clade  <dig> is known to include a-domains that encode for valine, glutamine, asparagine, leucine, and isoleucine .
clade  <dig>  which is known only to incorporate a single amino acid, proline , has a relatively equal distribution of both trichoderma  and tolypocladium  a-domains.
prolines are proposed to play an important structural role in peptaibols by creating a kink in the peptaibol chain, and clade  <dig> occupies a long branch within the tree, suggesting it is highly diverged from the other a-domains .
peptaibols were reported from t. geodes  <cit>  based on chemical isolation, but all previous reports of the compounds were identified from fungi in hypocreaceae  <cit>  .
a gene cluster responsible for the production of these peptaibols has not been identified, and to date no genomic sequence data have been produced for t. geodes from which to predict which genes or clusters may be responsible for its production.
efrapeptins, which are peptaibiotics originally described from t. inflatum, have been reported from several species of tolypocladium  <cit> , and it remains unknown whether the products produced by these clusters in t. ophioglossoides are of the efrapeptin class, or more traditional class of peptaibols.
regardless, the phylogenetic diversity of peptaibiotic nrpss, as revealed by phylogenomic analyses of a-domains, supports a greater chemical diversity of peptaibiotics than currently known from chemical analyses.
to interpret the evolutionary history of the peptaibol a-domains within hypocreales, the peptaibol a-domain clade tree was reconciled with the species tree .
because clade  <dig> is inferred to have an independent origin from clades  <dig> and  <dig> within the a-domain phylogeny, it was reconciled with the species tree separately .
we have included closely related a-domains from a complete a-domain phylogeny of hypocreales that are not a part of the three peptaibol a-domain clades in order to root the peptaibol a-domain clades.
the deep coalescence of peptaibol a-domains from all three clades at the common ancestor of trichoderma and tolypocladium  suggests that the presence of peptaibiotics in tolypocladium is an ancient attribute of tolypocladium genomes and not a product of a more recent horizontal transfer.
analyzing clades  <dig> and  <dig>  only  <dig> of the a-domains in t. ophioglossoides and t. inflatum coalesce at their most recent common ancestor, whereas  <dig> coalesce more deeply at the divergence of the four most derived families in hypocreales .fig.
<dig> a-domain tree/species tree reconciliation.
reconciliation of the peptaibiotic a-domain  tree with the species tree.
see additional file  <dig> for clade  <dig> reconciliation.
to root the tree, closely related outgroup a-domains from cochliobolus carbonum , fusarium graminearum , cordyceps militaris , and epichloë festucae  were included.
other abbreviations include: neurospora crassa , verticillium dahliae , f. verticillioides , f. oxysporum , nectria haematococcum , tr.
atroviride , t. reesei , tr.
virens , beauveria bassiana , claviceps purpurea , metarhizium robertsii , m. acridum , t. ophioglossoides , t. inflatum , ophiocordyceps sinensis
there exist at least two possible explanations to these findings in tolypocladium spp.
first, the divergence of the a-domains could be significant enough so as to distort the evolutionary history as represented in the phylogenetic tree, especially since many of the branches of the tolypocladium a-domain tree do not have strong bootstrap support.
this explanation is not supported by the data, because based on the species phylogeny amino acid divergence levels between tolypocladium species is less than that of trichoderma species .
second, incomplete lineage sorting could lead to this pattern of coalescence, suggesting that the ancestor to the genus possessed a multitude of these a-domains within one or several peptaibiotic genes that have undergone a complex history of ancient gains  and losses.
in contrast, five of the a-domains inferred to be present in the common ancestor of trichoderma spp.
deeply coalesce, and  <dig> are shared in the common ancestor.
this pattern indicates a higher degree of domain tree – species tree congruence, which could be explained by two different mechanisms including 1) vertical descent with expansion and maintenance or 2) horizontal gene transfer from tolypocladium to common ancestor of trichoderma.
to test for signatures of horizontal gene transfer, the a-domain tree was reconciled, with a modified species tree in which trichoderma was sister to tolypocladium, and outgroup a-domains were included to root the domain tree.
this produced a smaller deep coalescent cost , but this is due to the fewer number of extinctions  required , because the number of duplications remained the same .
in this simulated reconciliation,  <dig> a-domains were inferred to deeply coalesce, and of these, seven were inherited in each lineage.
taken together, this suggests that the diversity of trichoderma a-domains cannot solely be characterized as the product of horizontal transfer from tolypocladium a-domains.
instead it suggests that the common ancestor of trichoderma possessed a small number of peptaibiotic nrpss and a-domains that largely diversified in a manner consistent with speciation of the genus.
t. ophioglossoides peptaibiotic gene clusters
the four peptaibiotic genes  are located in three gene clusters on three different scaffolds .
genes toph_ <dig>  and toph_ <dig>  are located within the same gene cluster and separated by only nine genes, including one pks.
in total, the cluster contains  <dig> genes, many of which are typically found within secondary metabolite clusters including two multidrug transporters , a decarboxylase , an esterase , an epimerase , and a leucine zipper transcription factor  among others .
at more than  <dig>  amino acids in length, toph_ <dig> encodes for the largest peptaibiotic nrps produced by t. ophioglossoides, and is the largest gene  in the genome.
gene toph_ <dig> , located on a separate scaffold, is only  <dig>  amino acids in length and shares a high degree of amino acid identity with modules  <dig> and  <dig> of toph_ <dig> .
the genes surrounding toph_ <dig> include those typically found in an intact secondary metabolite cluster, including an abc transporter, transcription factors, and an esterase, and there are no indications that toph_ <dig> is a nonfunctional pseudogene.
however, no short peptaibiotics have been described so far, but the two a-domains of this gene both fall within peptaibiotic clade  <dig> , and it is clearly orthologous to part of toph_ <dig>  so while it may not be producing a peptaibiotic, it is analyzed here only within the evolutionary context of these genes.fig.
<dig> large peptaibiotic gene clusters of t. ophioglossoides.
panel a shows t. ophioglossoides scaffold  <dig> that includes two peptaibiotic nrps genes, while panel b shows the other large nrps gene which lies scaffold  <dig>  these clusters were predicted by the antismash pipeline, and putative annotation of genes within clusters are given.
proteins without putative functions are colored gray, and sizes of genes and spacing are approximate in scale
fig.
<dig> a-domain synteny map of for peptaibiotic genes in tolypocladium and trichoderma.
black lines connect orthologous a-domains as determined by mlbp ≥50 % in the reduced peptaibol a-domain tree
the remaining t. ophioglossoides peptaibiotic nrps gene, toph_ <dig> , is located within a cluster containing two pks genes , several abc transporters , an esterase , an epimerase , a hydrolase , two cytochrome p450s , and a radr transcription regulator  .
it remains to be seen if the products of the pkss are incorporated into the peptide created by toph_ <dig>  some of the a-domains within toph_ <dig> are divergent , especially those that group  as the earliest diverging lineages of clade  <dig>  and the presence of these a-domains in the tree, causes the support for this clade to weaken substantially.
peptaibiotic cluster synteny between mycoparasitic and insect pathogenic tolypocladium spp.
despite also possessing three large peptaibiotic genes in two clusters, t. inflatum peptaibiotic genes are highly divergent from t. ophioglossoides and the gene clusters are not located in the same regions of the genome .
alignment of t. ophioglossoides scaffold containing two peptaibiotic nrps genes  with the t. inflatum scaffold containing the single peptaibiotic nrps gene , revealed high synteny in the t. ophioglossoides two-peptaibiotic gene cluster region, except for the absence of the two peptaibiotic genes, themselves .
no other genes are present on the t. inflatum scaffold in the positions where the peptaibiotic nrps genes are located in the t. ophioglossoides genome.
this alignment did reveal, however, the presence of a truncated n terminal portion of an nrps  that aligns with the n terminal region of toph_ <dig>  when included in the entire hypocrealean a-domain phylogeny from additional file  <dig>  this t. inflatum nrps “relic”, tinf <dig>  is most closely related to the first a-domain of toph_ <dig> .
while retrotransposon relics are reported from secondary metabolite gene clusters  <cit>  and secondary metabolite gene modules may be truncated  <cit> , this is the first report of a relic nrps domain that remains within an otherwise intact secondary metabolite cluster including the pks located within the cluster.
due to the truncated nature of this protein model  and lack of other functional domains, it is unlikely that this relic produces a secondary metabolite peptide.
there is no evidence of transposable elements in or near this cluster in either of the tolypocladium genomes.
as mentioned above, this t. inflatum scaffold which is syntenous with the gene cluster in t. ophioglossoides containing two peptaibiotic nrps genes, also contains the nrps gene tinf <dig>  but this gene is located approximately  <dig> kb downstream of the relic cluster .fig.
<dig> analysis of cluster containing two peptaibiotic genes in t. ophioglossoides and corresponding region of the t. inflatum genome.
a. alignment showing high degree of synteny between t. ophioglossoides scaffold containing two peptaibiotic genes,toph_ <dig> and toph_ <dig>  and scaffold  <dig> in t. inflatum, which also contains the peptaibiotic gene tinf <dig>  that is over 380 kb upstream of the syntenic region.
b. subsection of the large a-domain tree showing the relationship between the tinf <dig> relic nrps and the first a domain of toph_ <dig>  c. map of major peptaibiotic genes within the tolypocladium genomes with the number of modules per gene shown below the scaffold.
gray area represents the syntenic regions between these clusters
the region of the scaffold containing the t. ophioglossoides third large peptaibiotic gene, toph_ <dig>  does not align well with any portion of the t. inflatum genome.
similarly, the final peptaibol cluster in t. inflatum, containing two peptaibiotic nrps genes , does not align well to any portion of the t. ophioglossoides genome.
the lack of synteny between these clusters in tolypocladium spp.
highlights the significant amount of genomic rearrangements between these closely related taxa.
campbell et al.
<cit>  observed patterns of differential gene loss in botrytis spp.
within an ancient, horizontally-transferred, secondary metabolite gene cluster, leading to a patchy distribution of the genes within the clusters.
this is not the pattern seen in the peptaibiotic clusters in tolypocladium spp., in which the protein models are not reciprocal best blast hits .
thus, despite the fact that their products may have similar functions, these peptaibiotic nrps genes are divergent and located within nonhomologous gene clusters.
mixed homology of peptaibol a-domains in hypocreales
using the moderate to strongly supported nodes  in the a-domain phylogeny as a guide for module homology, the peptaibol nrps genes are more conserved among the trichoderma species examined as compared to species of tolypocladium , a finding reflected in the domain tree – species tree reconciliation analyses .
using whole genome data of the sampled species of trichoderma, a-domains from the tr.
virens three peptaibol nrps genes   <cit>  were identified for the phylogenetic analyses; tr.
virens peptaibols contain  <dig>   <dig>  and  <dig> modules respectively.
tr.
reesei has two peptaibol nrpss  which possess  <dig> and  <dig> modules.
in the annotation of tr.
atroviride imi  <dig>  the a-domain hmm identified one  <dig> module peptaibol nrps , and several single a-domain protein models that group within the three peptabioitic clades and are all located on scaffold  <dig>  further examination of this tr.
atroviride gene region using the jgi genome browser  revealed that all of the ab initio gene predictions of that region predict a single protein model that is the approximate length of the  <dig> module peptaibol genes in tr.
virens and tr.
reesei.
degenkolb et al.
<cit>  identified a homolog of the  <dig> modular peptaibol gene from a different strain of tr.
atroviride, and thus this is likely a mis-annotation of tr.
atroviride scaffold  <dig>  alignment of this scaffold in tr.
atroviride  to those of tr.
virens and tr.
reesei  revealed high nucleotide homology.
however, the flanking regions did not align well, and a blast search of the tr.
atroviride genome using the nucleotide sequences from tr.
virens and tr.
reesei revealed that the  <dig> modular peptaibol gene in tr.
atroviride is located within a different portion of the genome than in the other two trichoderma spp.
.
comparing the trichoderma a-domains, each species possesses one ortholog of the large  <dig> or  <dig> module nrps gene , and at the nucleotide level these regions of their genomes also align and are syntenous .
the a-domains are syntenous in their arrangement within the large peptaibol nrps genes across the species, except for:  the insertion of a clade  <dig> domain at the third module, and  a duplication of either the seventeenth or eighteenth a-domains, which are most closely related.
within the two  <dig> module nrps genes in tr.
virens and tr.
reesei, there is complete synteny of the a-domains.
in tr.
virens, it has been demonstrated that this  <dig> module nrps, tex <dig>  is responsible for two different sizes of peptaibols   <cit> .
due to differences in annotation , the  <dig> module peptaibol gene in tr.
atroviride is not compared in this analysis.
the short  <dig> module peptaibol synthetase from trichoderma spp.
is found only in tr.
virens.
between these three groups of peptaibol nrpss, the terminating residues are all orthologous, as well as the initiating residues in the larger classes of peptaibol nrpss.
in tolypocladium, there is a very different pattern of homology and synteny between the peptaibiotic nrpss of the two species.
only a few of the a-domain relationships within tolypocladium are statistically supported  .
the first, third, and last a-domains of the largest nrps in both species  are orthologous, but not the other domains within those two genes.
there are several instances of intragenic module duplications which are known to occur within nrps genes and have been proposed to play a role in the evolution of novel metabolites  <cit> .
within toph_ <dig>  for example, there is strong support for a shared ancestry between modules  <dig>   <dig> and  <dig> , indicating that these modules are the product of lineage specific duplications .
this indicates a more complicated evolutionary history of these genes in tolypocladium.
the lack of module synteny and orthology between tolypocladium peptaibiotic gene modules is comparable to the lack of genomic synteny observed between their clusters.
part of this is due to the deep coalescence of the tolypocladium a-domains.
this evidence indicates that tolypocladium peptaibiotic genes are not highly similar but are the products of more ancient divergences.
this is notable, because in contrast to trichoderma spp., all of which exhibit some degree of mycoparasitism , t. ophioglossoides and t. inflatum have different ecologies, which are characterized by mycoparasitism and insect pathogenicity, respectively.
thus, if peptaibols are important in successful mycoparasitism  nrps genes in more ecologically diverse lineages.
future genome-scale studies sampling other hypocrealean lineages containing both mycoparasites and insect pathogens, including the genus polycephalomyces, could further test this hypothesis.
CONCLUSIONS
the genome of t. ophioglossoides is rich in secondary metabolite gene clusters, and  <dig> out of  <dig> of these clusters have no putative product.
given this potential and its life history as a mycoparasite, this species should be targeted for future studies to discover novel natural compounds with potential antibiosis, including antifungal, activity.
the sima nrps gene cluster, responsible for the production of the immunosuppressant cyclosporin, is not present in the t. ophioglossoides genome, but three large peptaibiotic genes are present within two clusters.
these are the first data to suggest the potential for peptaibiotic production from a mycoparasitic species of tolypocladium.
this study confirms the presence of three phylogenetic clades of peptaibiotic nrps a-domains from tolypocladium and trichoderma spp., and that peptaibiotics in general are limited to the mycoparasitic lineages of hypocreales, based on current sampling.
reconciliation of the a-domain tree with the organismal phylogeny reveals that the peptaibiotic nrpss of trichoderma and tolypocladium are likely the product of different mechanisms of diversification.
trichoderma is characterized by a-domain diversification that is largely consistent with speciation whereas tolypocladium is characterized by a-domain diversity that shows patterns of deep coalescence.
deep coalescence is inconsistent with peptaibiotic nrps diversity being the product of hgt to tolypocladium, rather it is the product of complex patterns of lineage sorting and gains and losses of a-domains from hypocrealean ancestors.
while the diversity of peptaibiotic nrpss in trichoderma could possibly be explained by hgt in the common ancestor of the three species, none of the tolypocladium peptaibiotic nrpss analyzed here are candidates for hgt.
further research is required to identify the structures of specific metabolites of the tolypocladium gene clusters and to determine if these peptaibiotics are produced during mycoparasitism by t. ophioglossoides or if these genes are present in other mycoparasitic lineages of hypocreales.

 hypocreales is home to a wide array of ecologically diverse fungi. some are devastating plant pathogens, while others form numerous lineages of both insect pathogens and mycoparasites. at the divergence of the four most derived families of hypocreales, there was a major shift away from plant-based nutrition to either insect pathogenesis or fungal parasitism, i.e., mycoparasites. two major lineages of mycoparasites are found within the order, although other mycoparasites exist. the first and larger of these two lineages is the family hypocreaceae, most notable for mycoparasitic trichoderma spp. used in biological control of plant pathogenic fungi, and trichoderma reesei e.g. simmons, the industrial workhorse for cellulase production. the second major lineage of mycoparasites, the genus tolypocladium, is nested within the insect pathogenic family, ophiocordycipitaceae. most species of tolypocladium parasitize the truffles of elaphomyces, ectomycorrhizal fungi closely related to aspergillus and penicillium. tolypocladium ophioglossoides quandt, kepler & spatafora is a commonly collected truffle parasite with a broad geographic distribution throughout many parts of the northern hemisphere. there are, however, a few tolypocladium species that attack insects and rotifers, and based on current multigene phylogenies some of these are inferred to be reversals to insect pathogenesis. one of these is a beetle pathogen, t. inflatum, which was the first source of the immunosuppressant drug, cyclosporin a. evidence from multigene studies has also shown a close phylogenetic relationship between t. ophioglossoides and t. inflatum. secondary metabolism is defined as the synthesis of often bioactive, small molecules that are not essential to the growth of an organism. genes related to production of secondary metabolites are often clustered together in close proximity within a genome and coregulated. a wide variety of secondary metabolites including the ergot alkaloids, fumonisins, and destruxins, is produced by species of hypocreales. many of these metabolites are produced by nonribosomal peptide synthetases, which are often large, multi-modular proteins that produce short peptides frequently incorporating non-standard amino acids. nrps modules are composed of three primary functional domains including adenylation, thiolation, and condensation domains. due to their high level of amino acid and nucleotide conservation, the a-domains are frequently used to reconstruct the evolutionary histories of these genes. polyketide synthases are another class of secondary metabolite producing enzymes that are common in fungi and are also modular in nature. they are related to fatty acid synthases, and assemble small bioactive molecules based on acetyl-coa or malonyl-coa subunits. the other major classes of secondary metabolite-producing enzymes are terpene synthases and dimethylallyltryptophan synthases, both of which have been reported from hypocrealean taxa. fungal secondary metabolites clusters often include genes required for regulation of expression of the gene cluster and decoration, epimerization, and transport of the mature secondary metabolite. peptaibols, or peptaibiotics, are antibiotic secondary metabolites products produced by very large nrps enzymes. their name is a derivative of their structure as they are peptides containing the uncommon non-proteinogenic amino acid, α-amino isobutryic acid, and a c-terminal amino ethanol. the presence of aib residues promotes helix formation, and several of these helices form multimeric units that in turn form voltage gated ion channels capable of inserting into cell membranes where they disrupt membrane potential causing leakiness. peptaibols are produced by trichoderma spp. and other members of hypocreaceae, leading to the proposition that they may play a role in mycoparasitism. there is at least one empirical study to support this in trichoderma. further studies have found that peptaibols function, along with cell wall degrading enzymes, to synergistically inhibit new cell wall synthesis in fungal prey of tr. harzianum. wiest et al. identified and characterized the first peptaibol synthetase nrps modular structure from tr. virens along with the residue peptaibol product. since that time, several other peptaibol nrps genes have been identified in tr. virens and other species. efrapeptins are another class of peptaibiotics described almost exclusively from tolypocladium spp. that have antifungal and insecticidal properties. they differ from orthodox peptaibols by the presence of a mitochondrial atpase inhibiting c-terminal “blocking group” n-peptido-1-isobutyl-2-ethylamine. eight of these, named efrapeptins a, and c-i, have been isolated from t. inflatum. numerous genomes from species of the mycoparasitic genus trichoderma have been sequenced, and more recently the genomes of several insect pathogens in hypocreales have been completed, including t. inflatum, the beetle pathogenic congener of t. ophioglossoides. the genomes of all of these species are rich in secondary metabolite genes and clusters, ranging from to secondary metabolite gene clusters per genome. comparisons of gene content and expression are beginning to shed light on mechanisms underlying host specificity and the evolution of primary and secondary metabolism. in this study, the draft genome of the truffle parasite t. ophioglossoides was generated to compare the gene content and secondary metabolite content of this truffle parasite to those of closely related insect pathogens and more distantly related mycoparasites. the secondary metabolite potential of t. ophioglossoides is characterized with a focus on understanding the evolution of gene clusters encoding for peptaibiotics.  genome assembly and structure the draft genome assembly of the t. ophioglossoides cbs is approximately megabases and is assembled on scaffolds. the assembly, with an n of base pairs, captures the majority of gene space and synteny, in spite of the exclusive use of short read technology. in addition, the core eukaryotic genes mapping approach identified complete and complete and partial core eukaryotic genes, estimating the assembly completeness at  % or  %. we predict protein-coding genes resulting in protein models, of which have support from rna, making the size and number of protein models similar to that of the beetle pathogen, t. inflatum, which is  mb and has protein models. they both also share high gc contents,  % and 58 %, for t. ophioglossoides and t. inflatum, respectively. these species are very closely related, which is reflected in the genome scale phylogeny and the large regions of shared synteny. there are, however, two large scale inversions and some small rearrangements with disagreement, suggesting these species have diverged despite the relatively short branch lengths separating them in the species phylogeny. the large number of small non-syntenic points in the mummerplot alignment could be due to the reliance on short read technology and the presence of short contigs in both assemblies.table genome statistics for tolypocladium ophioglossoides compared to t. inflatum t. ophioglossoides t. inflatum fig. maximum likelihood phylogeny of hypocrealean taxa analyzed in this study. raxml tree created using concatenated alignment of protein clusters identified in the hal pipeline with amino acid positions and distinct alignment patterns. branches are colored based on host/nutritional association. both pathogens and saprobes of plant material are colored green. species that possess peptaibol a-domains are denoted by a blue ‘α’ fig. whole genome nucleotide alignment of the tolypocladium ophioglossoides and the t. inflatum genomes. mummerplot visualization of nucmer alignment of all unmasked t. ophioglossoides scaffolds against the reference 101 t. inflatum scaffolds. red lines represent sequence aligning in the same direction while blue lines represent inversions. the black arrows point out the major inversions between the two genomes an overview of secondary metabolites in t. ophioglossoides the t. ophioglossoides and t. inflatum genomes harbor and core secondary metabolite genes – nrps, pks, terpene synthase, and dmats – spread across and secondary metabolite gene clusters, respectively. similar to the contrast observed between closely related metarhizium species, t. ophioglossoides and t. inflatum differ in the types of core secondary metabolite genes they possess with shared between the two species, unique to t. inflatum, and unique to t. ophioglossoides. notably, t. ophioglossoides does not contain the nrps gene, or any of the other genes in the sima cluster responsible for the production of cyclosporin a in t. inflatum. however, t. ophioglossoides does share a hypocreales-conserved core set of genes that flank the sima region with t. inflatum. cyclosporin production has been reported from several species of tolypocladium, however, none of the truffle parasites, including t. ophioglossoides, have been demonstrated to produce the compound. whether possession of the sima gene is a derived character state, or whether there has been a single or multiple losses within the genus, remains unknown and requires further sampling. secondary metabolites core genes predicted from the t. ophioglossoides genome include pkss, two pks-like genes, nrpss, six nrps-like genes, three hybrid nrps-pks genes, and four terpenes. no dmat synthases were identified in the t. ophioglossoides genome. based on a-domain homology, two putative siderophore synthetases were among the predicted nrpss. this is in contrast to t. inflatum, which possesses three putative siderophore synthetases. the entire pseurotin-a precursor synthetase hybrid nrps-pks cluster was identified in the t. ophioglossoides genome. pseurotin-a, an antifungal compound described from several aspergillus spp., was also recently identified in the genome of m. robertsii. the disjunct distribution of this secondary metabolite cluster raises several questions about the evolutionary mechanisms that may have led to this distribution. t. ophioglossoides possesses the destruxins synthetase nrps gene. destruxins are known for their insecticidal properties in metarhizium spp., and the entire destruxins synthetase cluster was characterized by wang et al. in m. robertsii. homologs of the other essential genes in the destruxins cluster are present in the t. ophioglossoides toph_ cluster, except for dtxs, an aspartic decarboxylase responsible for producing β-alanine, one of the amino acids incorporated into destruxins. there are inversions in this cluster between m. robertsii and t. ophioglossoides as well, including the dtxs aldo-keto reductase homolog and an abc transporter, which was not found to be essential in destruxins production in m. robertsii. for these reasons, it is not likely that t. ophioglossoides produces destruxins, but possibly produces another group of related compounds. the sequenced strain of tr. virens shares a homolog of the destruxins nrps gene, but destruxins have not been reported to be produced by that species either. to date, only two secondary metabolites have been reported to be produced by t. ophioglossoides: ophiocordin and ophiosetin. ophiosetin is structurally similar to equisetin, an antibiotic with inhibitory activity of hiv- integrase, and both are produced by nrps-pks hybrid genes. based on phylogenetic analysis of a-domains, cluster synteny with the equisetin cluster, and sequence homology, this study identifies the putative ophiosetin synthetase cluster around the hybrid nrps-pks, toph_. further studies involving transformations and chemical verification and characterization of this cluster will be necessary to confirm this genotype-chemotype linkage. ophiocordin is a polyketide and no putative gene or gene cluster related to its production was identified here. except for the two peptaibiotic clusters discussed below, the remaining secondary metabolite gene clusters are not yet associated with a specific gene product. peptaibiotics of tolypocladium among the wide assortment of secondary metabolite genes and gene clusters, the draft genome of t. ophioglossoides possesses four peptaibiotic nrps genes located within three gene clusters. phylogenetic analyses of fungal a-domains from a variety of nrps genes known to produce specific products and from whole genome mining of the hypocrealean species in fig.  revealed that trichoderma peptaibol a-domains group into three clades, one of which is well supported in all analyses, the other two clades have strong support with the exception of a few a-domains from one of the t. ophioglossoides peptaibiotic genes. all a-domains from the four peptaibiotic genes in t. ophioglossoides fall within these three clades, and representation in the three peptaibol clades is exclusive to a-domains from trichoderma and tolypocladium. t. inflatum also possesses three peptaibiotic nrps genes. strikingly, this limits the presence of peptaibiotic a-domains and genes to the two sampled mycoparasitic lineages in hypocreales.fig. phylogeny of the three peptaibol a-domain clades. peptaibol focused, a-domain phylogeny created using raxml, showing the major clades. green branches represent those supported by ≥70 % bootstrap support. tolypocladium spp. a-domains are colored in blue and teal, and trichoderma spp. a-domains are colored in yellow, red, or orange there are differences between species and gene membership within a-domain clades, however. for instance, clade is enriched in a-domains from the trichoderma peptaibols, while clade contains more a-domains from tolypocladium spp.. importantly clade contains a-domains that encode for incorporation of aib, as well as other a-domains that encode for incorporation of isovaline, leucine, isoleucine, alanine, glycine, valine, and serine. clade is known to include a-domains that encode for valine, glutamine, asparagine, leucine, and isoleucine. clade which is known only to incorporate a single amino acid, proline, has a relatively equal distribution of both trichoderma and tolypocladium a-domains. prolines are proposed to play an important structural role in peptaibols by creating a kink in the peptaibol chain, and clade occupies a long branch within the tree, suggesting it is highly diverged from the other a-domains. peptaibols were reported from t. geodes based on chemical isolation, but all previous reports of the compounds were identified from fungi in hypocreaceae. a gene cluster responsible for the production of these peptaibols has not been identified, and to date no genomic sequence data have been produced for t. geodes from which to predict which genes or clusters may be responsible for its production. efrapeptins, which are peptaibiotics originally described from t. inflatum, have been reported from several species of tolypocladium, and it remains unknown whether the products produced by these clusters in t. ophioglossoides are of the efrapeptin class, or more traditional class of peptaibols. regardless, the phylogenetic diversity of peptaibiotic nrpss, as revealed by phylogenomic analyses of a-domains, supports a greater chemical diversity of peptaibiotics than currently known from chemical analyses. to interpret the evolutionary history of the peptaibol a-domains within hypocreales, the peptaibol a-domain clade tree was reconciled with the species tree. because clade is inferred to have an independent origin from clades and within the a-domain phylogeny, it was reconciled with the species tree separately. we have included closely related a-domains from a complete a-domain phylogeny of hypocreales that are not a part of the three peptaibol a-domain clades in order to root the peptaibol a-domain clades. the deep coalescence of peptaibol a-domains from all three clades at the common ancestor of trichoderma and tolypocladium suggests that the presence of peptaibiotics in tolypocladium is an ancient attribute of tolypocladium genomes and not a product of a more recent horizontal transfer. analyzing clades and only of the a-domains in t. ophioglossoides and t. inflatum coalesce at their most recent common ancestor, whereas coalesce more deeply at the divergence of the four most derived families in hypocreales.fig. a-domain tree/species tree reconciliation. reconciliation of the peptaibiotic a-domain tree with the species tree. see additional file for clade reconciliation. to root the tree, closely related outgroup a-domains from cochliobolus carbonum, fusarium graminearum, cordyceps militaris, and epichloë festucae were included. other abbreviations include: neurospora crassa, verticillium dahliae, f. verticillioides, f. oxysporum, nectria haematococcum, tr. atroviride, t. reesei, tr. virens, beauveria bassiana, claviceps purpurea, metarhizium robertsii, m. acridum, t. ophioglossoides, t. inflatum, ophiocordyceps sinensis there exist at least two possible explanations to these findings in tolypocladium spp. first, the divergence of the a-domains could be significant enough so as to distort the evolutionary history as represented in the phylogenetic tree, especially since many of the branches of the tolypocladium a-domain tree do not have strong bootstrap support. this explanation is not supported by the data, because based on the species phylogeny amino acid divergence levels between tolypocladium species is less than that of trichoderma species. second, incomplete lineage sorting could lead to this pattern of coalescence, suggesting that the ancestor to the genus possessed a multitude of these a-domains within one or several peptaibiotic genes that have undergone a complex history of ancient gains and losses. in contrast, five of the a-domains inferred to be present in the common ancestor of trichoderma spp. deeply coalesce, and are shared in the common ancestor. this pattern indicates a higher degree of domain tree – species tree congruence, which could be explained by two different mechanisms including 1) vertical descent with expansion and maintenance or 2) horizontal gene transfer from tolypocladium to common ancestor of trichoderma. to test for signatures of horizontal gene transfer, the a-domain tree was reconciled, with a modified species tree in which trichoderma was sister to tolypocladium, and outgroup a-domains were included to root the domain tree. this produced a smaller deep coalescent cost, but this is due to the fewer number of extinctions required, because the number of duplications remained the same. in this simulated reconciliation, a-domains were inferred to deeply coalesce, and of these, seven were inherited in each lineage. taken together, this suggests that the diversity of trichoderma a-domains cannot solely be characterized as the product of horizontal transfer from tolypocladium a-domains. instead it suggests that the common ancestor of trichoderma possessed a small number of peptaibiotic nrpss and a-domains that largely diversified in a manner consistent with speciation of the genus. t. ophioglossoides peptaibiotic gene clusters the four peptaibiotic genes are located in three gene clusters on three different scaffolds. genes toph_ and toph_ are located within the same gene cluster and separated by only nine genes, including one pks. in total, the cluster contains genes, many of which are typically found within secondary metabolite clusters including two multidrug transporters, a decarboxylase, an esterase, an epimerase, and a leucine zipper transcription factor among others. at more than amino acids in length, toph_ encodes for the largest peptaibiotic nrps produced by t. ophioglossoides, and is the largest gene in the genome. gene toph_, located on a separate scaffold, is only amino acids in length and shares a high degree of amino acid identity with modules and of toph_. the genes surrounding toph_ include those typically found in an intact secondary metabolite cluster, including an abc transporter, transcription factors, and an esterase, and there are no indications that toph_ is a nonfunctional pseudogene. however, no short peptaibiotics have been described so far, but the two a-domains of this gene both fall within peptaibiotic clade, and it is clearly orthologous to part of toph_ so while it may not be producing a peptaibiotic, it is analyzed here only within the evolutionary context of these genes.fig. large peptaibiotic gene clusters of t. ophioglossoides. panel a shows t. ophioglossoides scaffold that includes two peptaibiotic nrps genes, while panel b shows the other large nrps gene which lies scaffold these clusters were predicted by the antismash pipeline, and putative annotation of genes within clusters are given. proteins without putative functions are colored gray, and sizes of genes and spacing are approximate in scale fig. a-domain synteny map of for peptaibiotic genes in tolypocladium and trichoderma. black lines connect orthologous a-domains as determined by mlbp ≥50 % in the reduced peptaibol a-domain tree the remaining t. ophioglossoides peptaibiotic nrps gene, toph_, is located within a cluster containing two pks genes, several abc transporters, an esterase, an epimerase, a hydrolase, two cytochrome p450s, and a radr transcription regulator. it remains to be seen if the products of the pkss are incorporated into the peptide created by toph_ some of the a-domains within toph_ are divergent, especially those that group as the earliest diverging lineages of clade and the presence of these a-domains in the tree, causes the support for this clade to weaken substantially. peptaibiotic cluster synteny between mycoparasitic and insect pathogenic tolypocladium spp. despite also possessing three large peptaibiotic genes in two clusters, t. inflatum peptaibiotic genes are highly divergent from t. ophioglossoides and the gene clusters are not located in the same regions of the genome. alignment of t. ophioglossoides scaffold containing two peptaibiotic nrps genes with the t. inflatum scaffold containing the single peptaibiotic nrps gene, revealed high synteny in the t. ophioglossoides two-peptaibiotic gene cluster region, except for the absence of the two peptaibiotic genes, themselves. no other genes are present on the t. inflatum scaffold in the positions where the peptaibiotic nrps genes are located in the t. ophioglossoides genome. this alignment did reveal, however, the presence of a truncated n terminal portion of an nrps that aligns with the n terminal region of toph_ when included in the entire hypocrealean a-domain phylogeny from additional file this t. inflatum nrps “relic”, tinf is most closely related to the first a-domain of toph_. while retrotransposon relics are reported from secondary metabolite gene clusters and secondary metabolite gene modules may be truncated, this is the first report of a relic nrps domain that remains within an otherwise intact secondary metabolite cluster including the pks located within the cluster. due to the truncated nature of this protein model and lack of other functional domains, it is unlikely that this relic produces a secondary metabolite peptide. there is no evidence of transposable elements in or near this cluster in either of the tolypocladium genomes. as mentioned above, this t. inflatum scaffold which is syntenous with the gene cluster in t. ophioglossoides containing two peptaibiotic nrps genes, also contains the nrps gene tinf but this gene is located approximately kb downstream of the relic cluster.fig. analysis of cluster containing two peptaibiotic genes in t. ophioglossoides and corresponding region of the t. inflatum genome. a. alignment showing high degree of synteny between t. ophioglossoides scaffold containing two peptaibiotic genes,toph_ and toph_ and scaffold in t. inflatum, which also contains the peptaibiotic gene tinf that is over 380 kb upstream of the syntenic region. b. subsection of the large a-domain tree showing the relationship between the tinf relic nrps and the first a domain of toph_ c. map of major peptaibiotic genes within the tolypocladium genomes with the number of modules per gene shown below the scaffold. gray area represents the syntenic regions between these clusters the region of the scaffold containing the t. ophioglossoides third large peptaibiotic gene, toph_ does not align well with any portion of the t. inflatum genome. similarly, the final peptaibol cluster in t. inflatum, containing two peptaibiotic nrps genes, does not align well to any portion of the t. ophioglossoides genome. the lack of synteny between these clusters in tolypocladium spp. highlights the significant amount of genomic rearrangements between these closely related taxa. campbell et al. observed patterns of differential gene loss in botrytis spp. within an ancient, horizontally-transferred, secondary metabolite gene cluster, leading to a patchy distribution of the genes within the clusters. this is not the pattern seen in the peptaibiotic clusters in tolypocladium spp., in which the protein models are not reciprocal best blast hits. thus, despite the fact that their products may have similar functions, these peptaibiotic nrps genes are divergent and located within nonhomologous gene clusters. mixed homology of peptaibol a-domains in hypocreales using the moderate to strongly supported nodes in the a-domain phylogeny as a guide for module homology, the peptaibol nrps genes are more conserved among the trichoderma species examined as compared to species of tolypocladium, a finding reflected in the domain tree – species tree reconciliation analyses. using whole genome data of the sampled species of trichoderma, a-domains from the tr. virens three peptaibol nrps genes were identified for the phylogenetic analyses; tr. virens peptaibols contain and modules respectively. tr. reesei has two peptaibol nrpss which possess and modules. in the annotation of tr. atroviride imi the a-domain hmm identified one module peptaibol nrps, and several single a-domain protein models that group within the three peptabioitic clades and are all located on scaffold further examination of this tr. atroviride gene region using the jgi genome browser revealed that all of the ab initio gene predictions of that region predict a single protein model that is the approximate length of the module peptaibol genes in tr. virens and tr. reesei. degenkolb et al. identified a homolog of the modular peptaibol gene from a different strain of tr. atroviride, and thus this is likely a mis-annotation of tr. atroviride scaffold alignment of this scaffold in tr. atroviride to those of tr. virens and tr. reesei revealed high nucleotide homology. however, the flanking regions did not align well, and a blast search of the tr. atroviride genome using the nucleotide sequences from tr. virens and tr. reesei revealed that the modular peptaibol gene in tr. atroviride is located within a different portion of the genome than in the other two trichoderma spp.. comparing the trichoderma a-domains, each species possesses one ortholog of the large or module nrps gene, and at the nucleotide level these regions of their genomes also align and are syntenous. the a-domains are syntenous in their arrangement within the large peptaibol nrps genes across the species, except for: the insertion of a clade domain at the third module, and a duplication of either the seventeenth or eighteenth a-domains, which are most closely related. within the two module nrps genes in tr. virens and tr. reesei, there is complete synteny of the a-domains. in tr. virens, it has been demonstrated that this module nrps, tex is responsible for two different sizes of peptaibols. due to differences in annotation, the module peptaibol gene in tr. atroviride is not compared in this analysis. the short module peptaibol synthetase from trichoderma spp. is found only in tr. virens. between these three groups of peptaibol nrpss, the terminating residues are all orthologous, as well as the initiating residues in the larger classes of peptaibol nrpss. in tolypocladium, there is a very different pattern of homology and synteny between the peptaibiotic nrpss of the two species. only a few of the a-domain relationships within tolypocladium are statistically supported. the first, third, and last a-domains of the largest nrps in both species are orthologous, but not the other domains within those two genes. there are several instances of intragenic module duplications which are known to occur within nrps genes and have been proposed to play a role in the evolution of novel metabolites. within toph_ for example, there is strong support for a shared ancestry between modules and, indicating that these modules are the product of lineage specific duplications. this indicates a more complicated evolutionary history of these genes in tolypocladium. the lack of module synteny and orthology between tolypocladium peptaibiotic gene modules is comparable to the lack of genomic synteny observed between their clusters. part of this is due to the deep coalescence of the tolypocladium a-domains. this evidence indicates that tolypocladium peptaibiotic genes are not highly similar but are the products of more ancient divergences. this is notable, because in contrast to trichoderma spp., all of which exhibit some degree of mycoparasitism, t. ophioglossoides and t. inflatum have different ecologies, which are characterized by mycoparasitism and insect pathogenicity, respectively. thus, if peptaibols are important in successful mycoparasitism nrps genes in more ecologically diverse lineages. future genome-scale studies sampling other hypocrealean lineages containing both mycoparasites and insect pathogens, including the genus polycephalomyces, could further test this hypothesis.  the genome of t. ophioglossoides is rich in secondary metabolite gene clusters, and out of of these clusters have no putative product. given this potential and its life history as a mycoparasite, this species should be targeted for future studies to discover novel natural compounds with potential antibiosis, including antifungal, activity. the sima nrps gene cluster, responsible for the production of the immunosuppressant cyclosporin, is not present in the t. ophioglossoides genome, but three large peptaibiotic genes are present within two clusters. these are the first data to suggest the potential for peptaibiotic production from a mycoparasitic species of tolypocladium. this study confirms the presence of three phylogenetic clades of peptaibiotic nrps a-domains from tolypocladium and trichoderma spp., and that peptaibiotics in general are limited to the mycoparasitic lineages of hypocreales, based on current sampling. reconciliation of the a-domain tree with the organismal phylogeny reveals that the peptaibiotic nrpss of trichoderma and tolypocladium are likely the product of different mechanisms of diversification. trichoderma is characterized by a-domain diversification that is largely consistent with speciation whereas tolypocladium is characterized by a-domain diversity that shows patterns of deep coalescence. deep coalescence is inconsistent with peptaibiotic nrps diversity being the product of hgt to tolypocladium, rather it is the product of complex patterns of lineage sorting and gains and losses of a-domains from hypocrealean ancestors. while the diversity of peptaibiotic nrpss in trichoderma could possibly be explained by hgt in the common ancestor of the three species, none of the tolypocladium peptaibiotic nrpss analyzed here are candidates for hgt. further research is required to identify the structures of specific metabolites of the tolypocladium gene clusters and to determine if these peptaibiotics are produced during mycoparasitism by t. ophioglossoides or if these genes are present in other mycoparasitic lineages of hypocreales. 
using nrps adenylation domain phylogenies and reconciliation of the domain tree with the organismal phylogeny, it is demonstrated that the distribution of these domains is likely not the product of horizontal gene transfer between mycoparasitic lineages, but represents independent losses in insect pathogenic lineages. peptaibiotics are a group of secondary metabolites almost exclusively described from trichoderma species of hypocreaceae. two major mycoparasitic lineages, the family hypocreaceae and the genus tolypocladium, exist within the fungal order, hypocreales. tolypocladium species are mainly truffle parasites, but a few species are insect pathogens. peptaibiotic genes are less conserved between species of tolypocladium and are the product of complex patterns of lineage sorting and module duplication. the draft genome sequence of the truffle parasite tolypocladium ophioglossoides was generated and numerous secondary metabolite clusters were discovered, many of which have no known putative product. phylogenomics and comparative genomics can provide insights into the evolution of secondary metabolite genes, their distribution across a broader range of taxa, and their possible function related to host specificity. peptaibiotics are produced by nonribosomal peptide synthetases and have antibiotic and antifungal activities. however, three large peptaibiotic gene clusters were identified using phylogenetic analyses. peptaibiotic nrps genes are restricted to mycoparasitic lineages of hypocreales, based on current sampling. peptaibiotic genes are absent from the predominantly plant and insect pathogenic lineages of hypocreales, and are therefore exclusive to the largely mycoparasitic lineages.
 two major mycoparasitic lineages, the family hypocreaceae and the genus tolypocladium, exist within the fungal order, hypocreales. peptaibiotics are a group of secondary metabolites almost exclusively described from trichoderma species of hypocreaceae. peptaibiotics are produced by nonribosomal peptide synthetases and have antibiotic and antifungal activities. tolypocladium species are mainly truffle parasites, but a few species are insect pathogens.  the draft genome sequence of the truffle parasite tolypocladium ophioglossoides was generated and numerous secondary metabolite clusters were discovered, many of which have no known putative product. however, three large peptaibiotic gene clusters were identified using phylogenetic analyses. peptaibiotic genes are absent from the predominantly plant and insect pathogenic lineages of hypocreales, and are therefore exclusive to the largely mycoparasitic lineages. using nrps adenylation domain phylogenies and reconciliation of the domain tree with the organismal phylogeny, it is demonstrated that the distribution of these domains is likely not the product of horizontal gene transfer between mycoparasitic lineages, but represents independent losses in insect pathogenic lineages. peptaibiotic genes are less conserved between species of tolypocladium and are the product of complex patterns of lineage sorting and module duplication. in contrast, these genes are more conserved within the genus trichoderma and consistent with diversification through speciation.  peptaibiotic nrps genes are restricted to mycoparasitic lineages of hypocreales, based on current sampling. phylogenomics and comparative genomics can provide insights into the evolution of secondary metabolite genes, their distribution across a broader range of taxa, and their possible function related to host specificity. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords secondary metabolismhypocrealesmycoparasiteslineage sortingissue-copyright-statement© the author 2015 