BACKGROUND
secreted proteins are essential for interactions with the environment, particularly for sessile organisms such as fungi.
fungi exhibit a wide diversity of nutritional lifestyles, ranging from strict saprobes, living only of dead organic material, through necrotrophic pathogens to obligate biotrophic mutualists and pathogens, existing within the living tissue of a host.
accordingly, fungal secreted proteins play important roles in the degradation of complex organic compounds, and in modulating directly or indirectly symbiotic and pathogenic interactions with their hosts.
effectors are secreted proteins of plant and animal pathogens, which are involved in manipulating, directly or indirectly, defense reactions and metabolism of their hosts .
for example, the maize pathogen ustilago maydis secretes a chorismate mutase that is translocated into infected maize cells to manipulate the shikimate pathway, which, amongst others, is involved in the production of salicylate, an important signal eliciting host defenses such as the systemic acquired resistance   <cit> .
other effectors prevent the detection of pathogen-associated molecular patterns , which are structurally conserved microbial molecules, by specific receptors of the host  <cit> .
for instance, chitin is an important structural component of fungal cell walls but also acts as a pamp.
however, in compatible plant-pathogen interactions, chitin does not elicit effective host defenses as it is either modified by extracellular deacetylation to chitosan or shielded by chitin-binding effectors .
thus, secreted proteins play essential roles in infection and determine whether a fungus may be an effective pathogen for a given host.
most of the fungal secreted proteins are believed to be transported through the type ii secretion pathway  <cit> .
proteins that are destined for this pathway contain an n-terminal signal peptide that mediates recognition and transport across the endoplasmic reticulum  membrane by a membrane-embedded translocon.
the signal peptide is cleaved off and is thus not part of the mature protein.
during transport through the er and the golgi apparatus the protein is folded, in many cases glycosylated and then travels from the golgi to one of several possible final destinations, one of which is the extracellular space .
secreted proteins have been identified by molecular genetic and biochemical methods, for example the yeast signal sequence trap  <cit> , but these typically identify only a fraction of the presumed total secretome.
the signal peptide that is characteristic of the type ii secretion pathway is highly conserved among eukaryotes, which allowed the development of computer programs such as signalp  <cit> , targetp  <cit> , and phobius  <cit>  for its identification in protein sequences.
this, coupled with the increasing availability of whole genome sequences for fungi, provides us with the resources to identify computationally the predicted set of secreted proteins  in a large number of fungi representing various nutritional lifestyles and phylogenetic lineages.
previous studies of fungal secretomes suggested that the composition of gene families varies greatly among species, and that some gene families appear to be associated with certain phylogenetic lineages and/or lifestyles  <cit> .
using gene model annotations from genome sequences, we performed a computational identification of soluble, secreted proteins, i.e.
proteins containing a signal peptide, but no signal anchor or transmembrane domains, in species representing both a broad phylogenetic diversity as well as diversity in lifestyles.
this allows assessing critically the hypothesis that species sharing a nutritional lifestyle have secretomes that exhibit a similar composition.
in this context, we specifically address the questions  whether secretome size and composition correlates with lifestyle and/or with phylogenetic lineages,  whether such observed differences might be the result of expansion or contraction of certain gene families, and  whether altered occurrences of functional units/domains, independent of gene families, might explain the differences.
RESULTS
prediction and characterization of fungal secretomes
the annotated gene models for  <dig> fungal and three phytopathogenic oomycota genomes were obtained from various databases .
the fungal collection comprised  <dig> ascomycota, eight basidiomycota, two zygomycota, one chytridiomycota, and one microsporidium.
the three oomycota genomes were included as outgroup for the plant pathogens within the mycota.
the selected fungal species did not only cover all phyla of the mycota but also represented a wide range of nutritional lifestyles.
thus,  <dig> species are saprophytes, eight are plant pathogens including one obligate biotroph, eight are vertebrate pathogens , one is an ectomycorrhizal, and one is a polyphagous fungus .
this latter species, i.e.
rhizopus oryzae
<cit> , has been found in association with plants but is also described as a mammalian pathogen and a saprophyte.table  <dig>
species analyzed
ashbya gossypii
aspergillus fumigatus
aspergillus nidulans
batrachochytrium dendrobatidis
botrytis cinerea
candida albicans
candida glabrata
candida lusitaniae
chaetomium globosum
coccidioides immitis
coprinus cinereus
cryptococcus neoformans
debaryomyces hansenii
encephalitozoon cuniculi
fusarium graminearum
kluyveromyces lactis
laccaria bicolor
magnaporthe oryzae
neurospora crassa
phanerochaete chrysosporium
phycomyces blakesleeanus
phytophthora infestans
phytophthora ramorum
phytophthora sojae
pichia stipitis
postia placenta
puccinia graminis
rhizopus oryzae
saccharomyces cerevisiae
schizosaccharomyces pombe
sclerotinia sclerotiorum
sporobolomyces roseus
stagonospora nodorum
trichoderma reesei
ustilago maydis
yarrowia lipolytica
names, codes used in figures, lifestyles, phyla, and subphyla of the species analyzed in this study, ordered alphabetically.
proteins containing a type ii secretion signal sequence were identified using the computer programs signalp  <dig>  and  <dig> , and targetp  <cit> .
a protein was only considered secreted when all three algorithms predicted this feature, according to  <cit> , which yielded an average of  <dig> % signal peptide-containing proteins per species .
as we focused our study on proteins that are soluble, we subtracted all proteins carrying transmembrane domains from this initial set.
these domains were identified by the program phobius that was specifically designed to distinguish between n-terminal transmembrane domains  and “classical” type ii secretion signal sequences  <cit> .
this showed that the soluble secretome constituted on average  <dig> % of the proteome with extreme values of  <dig> % for schizosaccharomyces pombe and  <dig> % for magnaporthe oryzae .
when assessing these soluble secretomes in a phylogenetic context, pezizomycotina and basidiomycota had similar sizes, whereas the average size of saccharomycotina was noticeably smaller .
comparing the soluble secretomes of fungi differing in nutritional lifestyle, saprophytes and animal pathogens had similar sizes, whereas the secretomes of plant pathogens were larger.
these differences held true for both absolute  and relative size .
at the level of absolute size of the soluble secretomes, the  <dig> species were grouped into three classes.
class  <dig> contained  <dig> species with secretomes comprising less than  <dig> proteins .
class  <dig> comprised  <dig> species with secretomes ranging from ca.
<dig> to  <dig> proteins.
class  <dig> included the remaining six species with secretomes between ca.
<dig> and  <dig> proteins.
in parallel, however, the species represented in these absolute secretome size classes also have increasingly large proteomes , which shows relative secretome size to be the more valid comparator for our data set.figure  <dig>
secretome sizes in relation to phylogeny and lifestyle.
absolute  and relative  sizes of secretomes.
a, b, d, e: species from phyla with at least three members.
rare species: the only representative of a phylogenetic lineage in a nutritional lifestyle group.
c, f: all species sharing a nutritional lifestyle except l. bicolor, as it is a symbiont, and r. oryzae, as it is polyphagous.
a, d: species grouped by phylogeny.
b, c, e, f: species grouped by nutritional lifestyle.
diamonds: data points, colored according to classification opposite to boxes.
boxes: lower and upper quartiles.
internal line: median.
lower and upper whiskers: minimum and maximum.
asterisks: means.
greyscales: phylogenetic groups.
colors: lifestyle groups.
comparing relative sizes of the soluble secretomes with phylogeny, within the ascomycota the saccharomycotina and taphrinomycotina have secretomes in the range of  <dig> % to  <dig> % of the proteomes with an average of  <dig> %, whereas the pezizomycotina secretomes range from  <dig> % to  <dig> % of the proteomes with an average of  <dig> % .
the basidiomycota secretomes cover a broad size range from  <dig> % to  <dig> % of the proteomes, with an average of  <dig> %.
overall, ascomycota have relative secretome sizes from  <dig> % to  <dig> % with an average of  <dig> % and thus cover a similar size range as the basidiomycota, with unicellular species at the lower and filamentous species at the higher ends .
in relation to nutritional lifestyle, animal pathogen secretomes range from  <dig> % to  <dig> % of the proteomes with an average of  <dig> %, whereas plant pathogen secretomes range from  <dig> % to  <dig> % of the proteomes with an average of  <dig> % and thus are considerably larger.
the sizes of the secretomes of saprophytes range from  <dig> % to  <dig> % with an average of  <dig> % .
when comparing, in dikarya, the  <dig> unicellular or dimorphic yeasts with the  <dig> filamentous fungi , it is obvious that the former group of taxa typically do not only have small absolute but also small relative secretomes.
the large-proteome filamentous fungi, on the other hand, rather have large secretomes both absolutely and relatively .
the plot of secretome size versus proteome size shows the general trend of increase of the former with increase of the latter .
the similar, small secretome-to-proteome size ratios of the saccharomycotina are clearly seen in the form of a “dip” of the data points in the lower left corner of the graph, while for the other species, these show more variation.
plots of secretome size vs. genome size  and proteome size vs. genome size  show similar, near-linear relationships with increasing variability for the larger genomes, secretomes, and proteomes.
in the latter two comparisons, phytophthora infestans appears as an extreme case, which is based in the massive expansion of its genome, due to increased numbers of repetitive elements  <cit> .figure  <dig>
correlations between the absolute sizes of genomes, proteomes, and soluble secretomes.
absolute sizes of predicted soluble secretomes vs. predicted proteomes.
absolute sizes of predicted soluble secretomes vs. genomes.
absolute sizes of predicted proteomes vs. genomes.
data points are color-coded by phylogeny.
abbreviations for species are given in table  <dig>
for each secretome, protein length distribution was analyzed by depicting the numbers of proteins in defined length intervals  as fractions of the total secretome .
the majority of proteins in each secretome had less than  <dig> aa .
however, proteins of up to  <dig> aa residues were remarkably enriched in four of the five plant pathogenic pezizomycotina analyzed, i.e.
botrytis cinerea, m. oryzae, stagonospora nodorum, and sclerotinia sclerotiorum, making up 10% to 15% of the total secretomes.
a similar enrichment of such small proteins was further observed only in the ectomycorrhiza  fungus laccaria bicolor, in the obligate biotrophic plant pathogen puccinia graminis, both belonging to the basidiomycota, and in the polyphagous zygomycete r. oryzae.figure  <dig>
protein size distributions within soluble secretomes.
number of proteins in five length intervals  as fractions of the total soluble secretomes.
circle areas correspond to the absolute sizes of the soluble secretomes.
species denominators are color-coded by lifestyle: dark blue = saprophytes, dark green = plant pathogens, red = animal pathogens, grey = polyphage , light blue = ectomycorrhiza , light green = obligate biotroph .
abbreviations for species are given in table  <dig>
we further divided the secretomes based on cysteine content and size since previous work indicated that putative effectors are often small cysteine-rich proteins .
small secreted proteins  are defined here to have a mature length of ≤  <dig> aa residues.
if such proteins had a relative cysteine content of ≥5%, as well as ≥ <dig> cysteine residues, we considered them as small secreted cysteine-rich proteins .
using these criteria, we found that cysteine-rich proteins were overall clearly enriched in the ssps , which was not obvious when applying a more relaxed definition for cys-enrichment using only a minimum of four residues but not considering the percentage of cysteine residues .
class  <dig> and  <dig> size secretomes had more than two-fold and three-fold higher sscp-to-secretome ratios than the small class  <dig> size secretomes, respectively, suggesting that the increase in secretome size may connect to a relative increase of sscp content.
pezizomycotina and basidiomycota had close to four-fold higher sscp-to-secretome ratios than saccharomycotina, and notably higher ratios than oomycota.
plant pathogens had a two-and-a-half-fold higher sscp-to-secretome ratio than animal pathogens, while saprophytes were in between and closer to plant pathogens .table  <dig>
relative cysteine contents by secretome size, phylogeny and lifestyle
b) relaxed definition: cysteine-rich ≥ <dig> cysteines
species
c-rich total
ssps
sscps
numbers of species and average contents of total cysteine-rich proteins, ssps and sscps in total secretomes in secretome size-, phylogenetic and lifestyle-groupings containing at least three species; ssp & sscp ≤  <dig> residues without signal peptides.
identification of protein clusters by markov chain clustering
all soluble secreted proteins identified were subjected to an “all-versus-all” blastp and this result was then used by the software tribe-mcl to arrange all sequences into protein clusters  <cit> .
thus,  <dig>  out of  <dig>  proteins were sorted into  <dig>  clusters , leaving  <dig>  singletons .
as similarly seen for whole proteomes  <cit> , certain clusters were either abundant in pezizomycotina or almost specific to them, as visualized by the intensity of red color, e.g.
below the universal clusters in the middle of figure 4a.
other protein clusters were either enriched in saccharomycotina , or practically absent.
furthermore, we identified a large group of oomycota-specific clusters , and a number of clusters specific to individual species.figure  <dig>
tribe-mcl-clustering of proteins and distributions of interpro domains and molecular function and biological process go terms.
heatmap of protein clusters with at least two proteins  and interpro domains , molecular function go terms  and biological process go terms  with at least two occurrences by protein clusters/interpro domains/go terms  and species .
increasing number of proteins/occurrences per cluster/interpro domain/go term and species is indicated by higher intensity of red.
species  are color-coded by lifestyle: dark blue = saprophytes, dark green = plant pathogens, red = animal pathogens, grey = “polyphage , light blue = ectomycorrhiza , light green = obligate biotroph .
species abbreviations are given in table  <dig>  protein clusters  with assigned interpro domains, molecular function or biological process go terms are provided on the right where the intensity of the blue color indicates their abundances.
the heights of the heatmaps correspond to the numbers of individual protein clusters/interpro domains/go terms represented in the respective maps.
to assign putative functions, all sequences were analyzed with interproscan  <cit> .
for the annotation of the secreted proteins, only interpro domains were considered.
on average, 32% of the sequences had at least one interpro  domain assigned.
the highest percentages of annotated proteins were obtained for aspergillus nidulans, aspergillus fumigatus, cryptococcus neoformans and trichoderma reesei , whereas p. graminis, b. dendrobatidis and encephalitozoon cuniculi had the lowest percentages , possibly relating to the intensity of research devoted to different taxa.
to assign interpro domains to the protein clusters determined by tribe-mcl, interpro domains found in members of a given cluster were assigned to the entire cluster if at least 75% of the sequences contained it.
using this approach,  <dig> out of  <dig>  protein clusters  had at least one interpro domain assigned .
moreover, we obtained gene ontology annotations of the types molecular function  and biological process  for the secreted proteins from the uniprotkb database.
on average, 35% of all sequences had at least one molecular function go annotation and 21% had at least one biological process go annotation.
for the molecular function go terms, the most extreme values for ratios of proteins in putative secretomes annotated with at least one term were 68% and 65% in candida albicans and s. pombe, respectively, and 8% and 14% in postia placenta and p. graminis, respectively, whereas for biological process go terms, the most extreme values were 79% and 76% of putative secretome proteins in saccharomyces cerevisiae and s. pombe, respectively, and 9% in both p. graminis and b. dendrobatidis.
as for the interpro domains, we assigned go molecular function and biological process annotations to protein clusters when at least 75% of the proteins in a cluster shared the go term.
thus, of the total of  <dig>  protein clusters,  <dig> were assigned at least one molecular function go term and  <dig> were assigned at least one biological process go term .
we then analyzed whether certain protein clusters determined above by tribe-mcl occur in all species of a given phylogenetic lineage but do not occur in any species outside of this lineage .
no cluster existed, that occurred in all  <dig> fungi.
two clusters contained proteins from all species analyzed except e. cuniculi, which is an intracellular pathogen of mammals exhibiting a strongly reduced genome and hence secretome size .
the proteins grouped in those two clusters were disulphide isomerases , involved in protein folding and maturation in the er  <cit>  and undefined proteins .
interestingly, no cluster existed that was absolutely ascomycota-specific .
the most widely represented cluster within ascomycota contained proteins from  <dig> of the  <dig> species analyzed within this taxon but was not specific to it since it also occurred in other phyla.
the corresponding protein was a mannosyl oligosaccharide glucosidase , the first enzyme in the n-linked oligosaccharide processing pathway, which was found exactly once per species except in chaetomium globosum.
for saccharomycotina, pezizomycotina, and basidiomycota, one specific cluster each was found .
the saccharomycotina-specific cluster  comprised two proteins each in pichia stipitis and c. albicans and one protein in each of the other species.
the cluster had no interpro domain or mf go or bp go term assigned, but  <dig> of  <dig> proteins  contained a cfem  domain according to interproscan that was described as an extracellular n-terminal domain comprising eight cysteines and that is found in certain fungal cell-surface integral membrane proteins  <cit> .
manual analysis showed that ten proteins belonging to cluster no.
<dig> contained an authentic cfem domain, while one sequence had only six cysteines.
only one of these proteins belonged to the sscps conforming to our new stringent definition given above.
the cfem domain is also common in pezizomycotina and basidiomycota, but the corresponding proteins distribute to other clusters, none of which shows full coverage of a specific phylogenetic lineage as cluster  <dig> does.
overall,  <dig> cfem domains were identified in proteins from  <dig> protein clusters and one orphan protein.
the pezizomycotina-specific cluster no.
<dig>  was represented by two proteins in t. reesei, fusarium graminearum, and neurospora crassa, and by one protein in each of the other species.
according to our above definition, all these proteins were sscps.
the single basidiomycota-specific cluster no.
<dig>  comprised one to three members in each of the species analyzed.
we analyzed if certain clusters were exclusively associated with all species exhibiting a given nutritional lifestyle .
we detected one particular cluster  that was found in all plant-pathogenic pezizomycotina but neither in any other species of that subphylum having a different lifestyle nor in plant pathogens outside of the pezizomycotina.
the saprophytic agaricomycetes coprinus cinereus, phanerochaete chrysosporium, and p. placenta had two clusters that encompassed only these species, i.e.
no.
<dig>  and no.
<dig> .
on average 60% of the sequences in each secretome clustered with sequences from other species and another 11% with other sequences exclusively from the same species.
despite this general trend, extreme examples existed where only very small fractions of the secretomes found similar members in other species.
in particular, only 20%, 9% and 14% of the sequences in the secretomes of b. dendrobatidis, e. cuniculi, and p. graminis did so, while 59%, 21%, and 37% clustered exclusively within these species, respectively.
this suggested that large parts of their secretomes are specific for these organisms or the lineages they belong to.
importantly, the heatmap representation of the tribe-mcl results indicated that proteins rather cluster by phylogeny than by nutritional lifestyle.
analysis of occurrence and abundance of interpro domains and molecular function and biological process go terms
above, we assigned interpro domains and molecular function and biological process go terms identified in the proteins of a cluster to that cluster only if these were identified in at least 75% of the proteins forming that cluster.
to analyze the secretomes based on functional traits, irrespective of the clustering of the proteins carrying them, we evaluated the presence and abundance of the interpro domains and go terms in relation to phylogeny and lifestyle.
from totals of  <dig> ,  <dig>  and  <dig>  occurrences of interpro domains, molecular function and biological process go terms, respectively, around 97% of each were derived from multiple occurrences of  <dig>   <dig> and  <dig> individual terms.
for comparison, only  <dig>  of  <dig>  proteins  were members of any of  <dig>  clusters.
moreover, the  <dig>  proteins in clusters accounted for 90% each of all occurrences of the three annotation types .
in the heatmap representations of abundances of interpro domains, molecular function and biological process go terms , distribution patterns were similar to protein distribution patterns  but they were less prominent.
for example, the band that spreads across the entire species spectrum in the middle of the protein cluster heatmap corresponds to the region roughly one third from the top in the interpro domain heatmap.
this band appears broader for the interpro domains  than for the protein clusters , because many of the highly conserved proteins in these clusters contain several different, equally conserved interpro domains.
the saccharomycotina-enriched protein clusters in the upper region of the protein cluster heatmap  are mirrored by saccharomycotina-enriched interpro domains in the lower region of the respective heatmap .
similarly, species-specific protein clusters for b. dendrobatidis, r. oryzae, phycomyces blakesleeanus, p. graminis, and p. placenta in the upper region of figure 4a correspond to specific interpro domains for these species in the lower and upper regions of figure 4b.
similar comparisons can be made for the molecular function and biological process go terms as well .
interestingly, a few species show many more go annotations than all others do, e.g.
s. pombe, s. cerevisiae and c. albicans, caused both by high ratios of the predicted secretomes having at least one annotation and many individual proteins having many different annotations.table  <dig>
distribution of proteins, interpro domains and go terms
total numbers of proteins and of interpro domain, molecular function go and biological process go term occurrences identified in these proteins.
fractions of proteins contained in protein clusters with at least two proteins and in clusters with proteins from at least three species.
fractions of interpro domain, molecular function go and biological process go term occurrences identified in these protein clusters and fractions of total occurrences of these terms derived from individual terms occurring at least twice and in at least three species.
next, we analyzed the interpro domains, i.e.
functional features, for correlation with phylogeny or nutritional lifestyle.
we identified no interpro domain present in all  <dig> species, but five existed in all species except e. cuniculi .
four of these five relate to disulfide-dithiol-mediated redox regulation.
most individual occurrences of the four were identified in the proteins making up the disulphide isomerase protein cluster  <dig> .
the fifth interpro domain present in all species except e. cuniculi represented a domain of glycoside hydrolases .
no other interpro domains were fully conserved within wide phylogenetic lineages .
the two zygomycota exclusively shared two interpro domains, whereas the three oomycota shared  <dig>  comparing interpro domains to lifestyle, no correlation was observed at all .
three molecular function go terms were conserved among all  <dig> species , while by phylogeny, five others were conserved in all species except e. cuniculi .
there existed only a single biological process go term conserved in all species , while two others  were conserved in all species except e. cuniculi .
by lifestyle, no molecular function or biological process go term was conserved .
protein cluster contractions and expansions
so far, we analyzed the absence or presence of protein clusters in all fungi sharing a given phylogenetic lineage or lifestyle.
however, protein families may have contracted or expanded in individual species or parts of evolutionary lineages and such changes may relate to lifestyle.
to reduce noise in quantitative assessment, expansion was defined as the presence of more than twice the median number of proteins for a cluster in a given species, contraction as less than half the median number of proteins.
applying these criteria to the  <dig> protein clusters present in at least three species, that thus contain at least three proteins to be able to calculate medians, we found  <dig> clusters with expansions and/or contractions in at least one species .
within the limits of the applied criteria, no protein cluster was contracted or expanded exclusively in all species belonging to a certain phylogenetic lineage, or exhibiting a certain nutritional lifestyle .
even when scoring contractions or expansions less stringently, i.e.
when allowing that this might have happened in additional species representing other lineages or lifestyles, we only found one protein cluster.
this cluster  was expanded in all five agaricomycotina, being present in five copies in the mammalian pathogen c. neoformans and up to seventeen in the saprophyte p. placenta.
in other basidiomycota, it was expanded in the cereal pathogen p. graminis , in the maize pathogen u. maydis , but not in the saprophyte s. roseus .
this is visualized in figure 5a by a red line just below the empty upper region in the basidiomycota that is gapped only under s. roseus.
beyond, cluster  <dig> was also expanded in the zygomycota p. blakesleeanus  and r. oryzae , which had six and sixteen proteins, respectively.
it had just one copy in four species of saccharomycotina whereas the remaining species entirely lacked it.
except in c. globosum where it was also missing, species in the pezizomycotina had up to three proteins belonging to cluster  <dig> figure  <dig>
expansions and contractions in protein clusters and in interpro domain and molecular function and biological process go term abundances.
heatmap of protein clusters , interpro domains , molecular function go terms  and biological process go terms  represented in at least three species that show an expansion of more than twofold  or a contraction of less than half  of the median of the abundances per species for the respective entry type.
increasing level of expansion or contraction per cluster/interpro domain/go term and species is indicated by higher color intensity.
less than twofold expansion or contraction is shown in monochrome grey.
white indicates no protein in this cluster or occurrence of this interpro domain or go term in this species.
species  are color-coded by lifestyle: dark blue = saprophytes, dark green = plant pathogens, red = animal pathogens, grey = polyphage , light blue = ectomycorrhiza , light green = obligate biotroph .
species abbreviations are given in table  <dig>  protein clusters  with assigned interpro domains, molecular function or biological process go terms are provided on the right where the intensity of the blue color indicates their abundances.
the heights of the heatmaps correspond to the numbers of individual protein clusters/interpro domains/go terms represented in the respective maps.
the lack of general contractions even when using the less stringent condition results probably from the fact that in most clusters the median of the number of proteins was very low  and that the absence of a protein in a species was not scored as a contraction.
however, we did observe in five instances interesting trends of expansions or contractions in certain phylogenetic lineages .
the first of these five protein clusters  was described above as the only one showing a general expansion in a phylogenetic lineage  and some basidiomycota and zygomycota.
the second , apparently arose in the ascomycota since it was not present in any of the other lineages.
the saccharomycotina had between two and five members in each species, whereas in the pezizomycotina the protein number was contracted to one each in f. graminearum, m. oryzae, coccidioides immitis, a. nidulans, and a. fumigatus, visible as a green line in the uppermost part of figure 5a.
t. reesei, c. globosum, s. sclerotiorum, and s. nodorum lacked these proteins completely.
the third protein cluster  showed a trend of contraction in agaricomycetes.
within basidiomycota, members of this cluster were restricted to species of this class with only one copy each, except p. placenta that had five.
in the pezizomycotina, it was found with two to six proteins per species, except in c. immitis that had only one.
all other lineages entirely lacked cluster  <dig>
the fourth cluster  showed trends of expansion in three pezizomycotina.
it contained proteins in each basidiomycota except u. maydis and in each pezizomycotina except f. graminearum, m. oryzae, and c. immitis.
the basidiomycota and t. reesei, n. crassa, c. globosum, s. sclerotiorum and s. nodorum each had one to four of such proteins.
the helotiales b. cinerea, the eurotiales a. nidulans and a. fumigatus showed expansions with five, five, and seven copies, respectively.
all other lineages entirely lacked cluster  <dig> .
finally, cluster  <dig> , was also only found in pezizomycotina and basidiomycota.
in basidiomycota, only u. maydis and c. cinereus had this cluster, containing one and three members, and in ascomycota only seven species within the pezizomycotina with up to six copies in s. nodorum.
the eurotiales a. nidulans and a. fumigatus, but also u. maydis showed a contraction in cluster  <dig>  with one copy each.
when further focusing our analysis on individual species, we observed a number of expanded protein clusters on condition that a given cluster should comprise at least ten members in that species and at least 50% more members than in all other species having this cluster as well .
disregarding the phytophthora species,  <dig> of such cases were found in species of basidiomycota and saccharomycotina, but not of pezizomycotina, and further in b. dendrobatidis and r. oryzae.
in general, members of basidiomycota showed the most extreme cases of such single-species expansions .
contractions and expansions in abundance of interpro domains and of molecular function and biological process go terms
as we analyzed the presence or absence of all interpro domains and molecular function and biological process go terms identified in the  <dig> secretomes independent of the proteins in specific lineages or lifestyles, we also analyzed whether specific contractions or expansions occurred in these data sets.
we applied the same limits for contraction or expansion of abundances of these terms as for protein clusters.
thus, out of  <dig> interpro domains, representing a total of  <dig>  individual interpro domain occurrences,  <dig> showed a contraction or expansion in at least one species .
furthermore, out of  <dig> and  <dig> molecular function and biological process go terms, representing  <dig>  and  <dig>  individual occurrences,  <dig> and  <dig> showed at least one expansion or contraction, respectively .
as for protein clusters, abundance of no interpro domain was contracted or expanded exclusively in all species from any phylogenetic lineage or specific nutritional lifestyle  or subgroup thereof .
nevertheless, we observed trends of contraction or expansion.
these, however, were mainly identified in the context of phylogeny.
for those interpro domains that were identified exclusively  or to at least 75%  in the members of individual protein clusters, expansions and/or contractions of interpro domain abundances in multiple species  typically paralleled the expansions and/or contractions in the corresponding protein clusters .
the remaining, however, did not show such one-to-one associations  and in the following, we highlight some noticeable cases.
ipr <dig>  was the most common domain in the data set and is linked to a large number of gh families.
this domain was identified in all species except e. cuniculi and within fungi was mostly far more abundant in pezizomycotina, with  <dig> to  <dig> occurrences per species, and basidiomycota, with  <dig> to  <dig> occurrences per species, than in saccharomycotina, where it was contracted with three to nine occurrences per species.
interestingly, b. dendrobatidis contained ipr <dig> only once.
the animal pathogens c. immitis  and c. neoformans  and the saprotrophic basidiomycota s. roseus, with five or eight occurrences, also showed contractions for this interpro term within their phylogenetic lineages.
of all the other interpro domains linked to ipr <dig>  only ipr <dig>  and ipr <dig>  showed similar distributions .
the latter is contracted in saccharomycotina.
not surprisingly, molecular function go terms correlating with ipr <dig>  particularly mf go: <dig>  hydrolase activity, hydrolyzing o-glycosyl compounds, the most common molecular function go term in the data set, and mf go: <dig>  hydrolase activity, showed similar patterns of contractions.
these terms, however, also showed expansions in many pezizomycotina and basidiomycota, but no correlation with lifestyle existed for these expansions.
as glycoside hydrolases are enzymes important for such diverse aspects as cell wall remodeling, degradation of host cell walls by plant pathogens to allow spreading through tissues and degradation of polysaccharides to provide carbon sources, all secreted proteins were analyzed against the cazy database by dbcan.
in total, members of  <dig> gh families were identified .
whereas  <dig> and  <dig> gh families were identified in pezizomycotina and basidiomycota, respectively, saccharomycotina contained members in only  <dig> gh families.
in pezizomycotina, ghs ranged from  <dig> proteins in  <dig> families in the mammalian pathogen c. immitis to  <dig> proteins in  <dig> families in the plant pathogen s. nodorum.
in saccharomycotina, gh content ranged from  <dig> proteins in  <dig> families in the mammalian pathogen candida lusitaniae to  <dig> proteins in  <dig> families in the mammalian pathogen c. albicans.
in basidiomycota, ghs ranged from  <dig> proteins in nine families in s. roseus to  <dig> proteins in  <dig> families in c. cinereus.
with respect to nutritional lifestyle, gh content was relatively low in the specialized animal pathogens c. neoformans  and c. immitis , but also in the biotrophic plant pathogen u. maydis  and in the saprophytic yeast s. roseus .
ipr <dig> , the second most abundant interpro domain in the entire data set , is a domain that is found in pectate lyase, iota-carrageenase and ghs from family  <dig>  including various galacturonases.
within dikarya, it was absent from saccharomycotina except for contracted single occurrences in s. cerevisiae and ashbya gossypii, from s. pombe, and from s. roseus.
in pezizomycotina, it was typically more abundant than in basidiomycota.
expansion was seen in eurotiales, in helotiales and in f. graminearum.
ipr <dig>  was present in all ascomycota except debaryomyces hansenii, but only in four of eight basidiomycota .
within pezizomycotina, it was expanded in all plant pathogens.
this domain is contained in zinc-dependent metalloproteases from family m <dig> that are typically soluble secreted proteins, synthesized as zymogenic proproteins, that degrade proteins by hydrolyzing single residues from the c-termini .
ipr <dig>  describes an n-acetylglucosamine-binding domain that occurs in chitinases of fungi and plants.
it is ubiquitous in pezizomycotina except c. immitis that lacks it, and is expanded in all plant pathogens.
in both saccharomycotina and basidiomycota, however, only half of the species share it and in each, it is contracted with one occurrence only.
as expected, the corresponding molecular function go term, go: <dig>  chitin binding, shows practically the same pattern of abundances, contractions and expansions.
interestingly, though, the expansions of chitin-binding domains in a number of mostly plant pathogenic pezizomycotina was not mirrored by the terms mf go: <dig>  chitinase activity and bp go: <dig>  chitin catabolic process, relating to enzymatic activity.
either these species may have some chitin-active enzymes containing multiple chitin-binding domains, or these domains may occur independently of enzymes, for example in effectors shielding the fungi from recognition by their host plants.
the last example of expansions in abundance of interpro domains are ipr <dig>  gmc oxidoreductase, n-terminal and ipr <dig>  gmc oxidoreductase, c-terminal.
glucose methanol choline oxidoreductases are flavoenzymes utilizing various substrates.
within fungi, these interpro domains were only found in pezizomycotina and basidiomycota.
interestingly, while these domains tend towards expansion in the plant pathogenic pezizomycotina they are lacking in the specialized animal pathogens c. immitis and c. neoformans, respectively.
c. cinereus has at least thrice as many occurrences of both gmc oxidoreductase interpro domains as any other basidiomycota species.
as for protein clusters, we also observed a number of interpro domains, molecular function and biological process go terms that showed major expansions specifically in a single species, with more than ten occurrences and at least 50% more occurrences than in any other species having that interpro domain.
many of these expansions reflected expansions in the protein clusters, annotated with these interpro domains or go terms .
however, some prominent examples existed of expanded terms not linked to expanded protein clusters, such as ipr <dig>  in the pezizomycotina m. oryzae with  <dig> occurrences, and ipr <dig>  in the basidiomycota c. cinereus with  <dig> occurrences.
discussion
the completion of annotated proteomes allows for in depth analyses of secretomes by bioinformatics.
our comparative analyses of fungal and oomycota secretomes, representing a broad range of taxa and nutritional lifestyles, revealed that primarily phylogenetic relationships determine the size and the composition of the secretomes, and to a lesser extent the nutritional lifestyle.
we identified protein clusters that are specific for or expanded in certain lineages.
others are expanded with respect to lifestyle, but only within individual lineages.
secretome size and composition
secretome size correlates with phylogenetic lineages.
saccharomycotina have smaller secretomes, both absolute and relative, than pezizomycotina or basidiomycota, which have similar sized secretomes in both respects.
this may be because yeasts are less capable of utilizing complex organic polymers as substrates, since they often live in niches providing access to sugars and other compounds that are easier to metabolize.
secretome size is also associated with lifestyle.
plant pathogens have the largest, saprophytes have intermediate, and animal pathogens have the smallest secretomes.
to some extent, this finding may be explained by the assumption that animal hosts represent a nutritionally simpler environment than plant hosts or litter.
furthermore, plant pathogens, particularly obligate biotrophs such as p. graminis, may secrete many effectors that would not provide benefits to saprophytes.
it was suggested that portions of their secretomes participate in the coevolutionary arms race with host factors thus leading to sequence diversification  <cit> .
besides, p. graminis initially infects its alternate host berberis to switch thereafter to a cereal host.
it must be presumed that the sets of secreted proteins involved in infection of both hosts do only partially overlap which will further increase the total size of the rusts’ secretome.
moreover, the typically small size of animal pathogen secretomes may reflect an adaptation to the flexibility of the animal immune system.
against a system effectively recognizing non-self molecules and mounting system-wide defense reactions in a circulatory system, a “stealth mode” is likely an efficient way of evasion.
in such a scenario, it appears as a selective advantage to minimize the numbers of putative pamps and thus to reduce the size of the secretome.
comparison of the mammalian pathogens within ascomycota points to another link between environmental influences and secretome size.
the three candida species, c. immitis and a. fumigatus have secretomes of 203- <dig>   <dig> and  <dig> proteins, respectively.
besides causing invasive mycoses in immunocompromised hosts, candida species are commensals on mammalian skin and mucosa.
c. immitis is able to infect immunocompetent hosts and beyond lives as a saprophyte on host carcasses, with a concomitant lack of many enzymes required for decomposition of plant materials  <cit> .
aspergillus species typically are environmental saprophytes and only a few of them can in addition act as opportunistic pathogens of humans and animals  <cit> .
comparison of the two aspergillus species in our data sets shows that a. nidulans, having only the saprophytic lifestyle, has absolutely and relatively a larger secretome, and secretes  <dig> more enzymes than a. fumigatus.
these examples provide support for the notion that dominance of the saprophytic phase is associated with a larger secretome.
furthermore, the higher number of secreted enzymes allows a. nidulans to utilize a wider variety of substrates present in soil or plant debris, which are likely more difficult to degrade than those available in animal hosts.
in addition, competition with other microorganisms is conceivably less severe in animal hosts than in soil or litter, which may affect the number of secreted proteins involved in combatting competitors.
the secretomes in the oomycota are considerably larger than in most pezizomycotina and basidiomycota and larger than most fungal plant pathogen secretomes in absolute, but not in relative numbers.
the high absolute numbers of secreted proteins are derived, at least in part, from expansions in the rxlr- and crinkler-effector families which comprise a considerable fraction of these secretomes  <cit> .
we analyzed the secretomes for the distribution of size classes and cysteine richness to determine whether sscps, previously discussed to harbor many effectors, are enriched in certain phylogenetic lineages or nutritional lifestyles.
examples for this type of effectors are some well characterized proteins from the biotrophic plant pathogen cladosporium fulvum
<cit> .
other sscps are hydrophobins covering aerial mycelia and spores  <cit> , the membrane-bound cfem-proteins  <cit> , antimicrobial proteins important to fight off competitors  <cit>  and peptides such as ptrtoxb of pyrenophora tritici-repentis targeting the host  <cit> .
using our definition of sscps  our analysis shows a correlation of both phylogeny and lifestyle with the percentage of sscps in secretomes.
by phylogeny, the relative sscp content is similar in pezizomycotina and in basidiomycota, while it is almost four-fold lower in saccharomycotina.
as for lifestyle, the sscp fraction of the total secretomes is 30% higher in plant pathogens than in saprophytes, presumably because of the considerable numbers of potential effectors.
in animal pathogens, this fraction is only half of that of saprophytes.
the fact that the highest individual ssp and sscp fractions of the total secretomes occur in p. graminis and the hemibiotrophic pezizomycotina m. oryzae fits with the assumption that effectors would play eminent roles particularly in biotrophy.
our results correspond well with the large secretomes and high numbers of ssps recently reported for obligate biotrophic phytopathogens melampsora larici-populina and p. graminis
<cit> .
presence/absence of and contractions/expansions in protein clusters, interpro domains and go terms
we clustered all proteins into families, assigned interpro domain and molecular function and biological process go annotations, and identified clusters with expansions or contractions in specific taxonomic groups.
it is important to keep in mind that what we define here as an expansion may be a recent lineage specific series of gene duplications or could be due to ancient gene duplications followed by gene loss in the other lineages.
to distinguish between these two alternative hypotheses will require in depth phylogenetic analyses of each family.
however, in either case, we consider that an expansion represents the presence of lineage specific selective pressure for those species to maintain the gene, regardless of whether the gene duplication events were recent or ancient.
the secretomes of the saccharomycotina lack many protein clusters completely, which may relate to their unicellular lifestyle and overall reduced genomes.
some of the clusters missing in saccharomycotina are expanded or contracted in basidiomycota compared to pezizomycotina or vice versa.
apart from reflecting the specific phylogenetic histories of the corresponding species, expansion of protein families may correlate with their abilities to cope with more complex and/or variable environments.
assessment of distinct phylogenetic lineages revealed that saccharomycotina, pezizomycotina, and basidiomycota each contain a single lineage-specific protein cluster.
the saccharomycotina-specific cluster contains the cfem motif, but narrowly misses the threshold for interpro domain assignment and has no go terms assigned.
the cfem domain is not saccharomycotina-specific, as we identified it in  <dig> protein clusters across the dikarya.
this example illustrates the intrinsic difficulties in analysis of protein functions by assigning interpro domains or go terms to protein clusters.
some clusters have not been assigned terms because the assignment threshold was missed or proteins containing the same terms were clustered in different protein clusters.
for this reason, we also performed analyses for the interpro domains and molecular function and biological process go terms themselves.
four of five interpro domains absent only in e. cuniculi relate to er- or golgi-localized proteins, which may reflect the major reduction in secretion in this intracellular parasite  <cit> .
the fifth represents the catalytic core of glycoside hydrolases.
no other interpro domain showed absolute specificities either by phylogeny or by lifestyle.
by our definition, we found no protein cluster or interpro domain either contracted or expanded exclusively within all species representing a particular lifestyle or phylogenetic lineage.
nevertheless, five protein clusters were expanded or contracted in some species in certain phylogenetic lineages.
for interpro domain abundances, we found two terms to be contracted in saccharomycotina and further observed trends of expansions, as well as contractions, in subsets of species within phylogenetic lineages.
in saccharomycotina, the glycoside hydrolase catalytic domain  and glycoside hydrolase, family  <dig>  were contracted according to the definition applied, consistent with the go- and particularly the cazy-data.
this may be indicative of a narrower spectrum of exploitable polysaccharides in this taxon.
two domains are expanded in all plant pathogenic pezizomycotina.
some of the zinc-metalloproteases from family m <dig> , that typically degrade extracellular matrix proteins, may be involved in degradation of plant tissue, others however could mediate the suppression of plant defense by degradation of host proteins involved in pathogen recognition in the apoplast or playing other important roles in defense.
the other expanded ipr domain is the chitin-binding domain, type  <dig>  that occurs in single or multiple copies in fungal chitinases and other glycoside hydrolases.
in plant pathogenic fungi, it may mediate binding of the containing proteins to the fungal cell wall to avoid host recognition by shielding or modifying chitin, which is a pamp as discussed above.
in addition, chitinases or other glycoside hydrolases containing this domain could serve host cell wall degradation or fungal morphogenesis.
several interpro domains for glycoside hydrolases and other enzymes involved in degradation, as well as go terms relating to such activities, show considerable tendencies towards expansion in saprophytes and plant pathogens in the pezizomycotina and basidiomycota.
for glycoside hydrolases, the cazy-data confirmed this tendency.
the likely explanation is that these species can utilize a wider variety of complex carbohydrates than yeasts or animal pathogens.
examples are iprs  <dig> and  <dig>  that are particularly expanded in plant pathogenic pezizomycotina plus the saprotrophic basidiomycota c. cinereus.
proteins containing these domains, including various alcohol and sugar peroxidases, produce hydrogen peroxide, which may amongst others modify target molecules directly or serve as substrate for oxygen radical generation by fenton chemistry.
radicals may serve the breakdown of lignin in white rot fungi or of polysaccharides as in cellulose degradation by brown rot fungi  <cit> .
the role of peroxide generated by plant pathogens is diverse  <cit> .
for example, it may interfere with plant defense compounds by direct reduction, as was recently proposed for a fad-dependent glucose dehydrogenase from the pezizomycotina glomerella cingulata that causes anthracnose on various fruits  <cit> .
apparently, only few protein clusters specific for lineages exist, which may be explained by the broad phylogenetic distance covered by the fungi analyzed thus obstructing the clustering of related proteins.
in addition, particularly in the case of plant pathogenic fungi, proteins required for host interaction may be under positive selection due to the “arms-race” between host and pathogen  <cit> .
thus they may not be recognized as originating from the same ancestral gene anymore, even between relatively closely related species as shown for three species of powdery mildew fungi that only shared less than  <dig> effectors  <cit> .
finally, losses or gains of protein families need not be equivalent with losses or gains of biological functions, as these may be provided by other protein families.
it has been shown that changes in lifestyle  have occurred often in fungal evolution  and such events need not be accompanied by the same losses or gains of proteins.
with respect to expansion of protein clusters in individual species only, we observed for example an expansion of peptidases  and serine endopeptidases of family s <dig> that cleave target proteins close to alanine-rich c-terminal tripeptides ) in b. dendrobatidis, confirming previous observations  <cit> .
this may reflect an adaptation to the animal pathogenic lifestyle and similar expansions were observed for peptidases in other fungi exhibiting this lifestyle.
a study comparing closely related animal pathogens and saprophytes showed that in c. immitis and c. posadasii families of subtilisins and other family s <dig> serine peptidases, which typically are non-specific endopeptidases , and mentioned m <dig> metalloproteases are expanded  <cit> .
these families are also expanded in the non-pathogen uncinocarpus reesii, but not in the more distantly related histoplasma capsulatum.
in contrast to these onygenales, within the eurotiales these peptidase families are expanded neither in the saprophyte a. nidulans nor in the opportunistic pathogens a. fumigatus and a. terreus, illustrating different extents of adaptation to the pathogenic lifestyle.
complementary, whereas the onygenales show strong contractions in, or even lack of a number of plant cell wall degrading enzyme families , the helotiales all contain similar numbers of these enzymes  <cit> .
this suggests that there is a considerable phylogenetic signal in the data of species comparisons irrespective of the degree of kinships and furthermore, that lifestyles are not linked to the presence, absence or abundance of certain protein families of functional groups in a straight-forward causal relationship.
CONCLUSIONS
we identified and analyzed the soluble secretomes of fungi to test the hypothesis that the secretomes of species belonging to different phylogenetic lineages but sharing the same nutritional lifestyle also share a similar composition.
secretome size correlates with lifestyle, with plant pathogens having the largest and animal pathogens having the smallest secretomes.
in contrast, the signal of lifestyle on the specific composition of secretomes is restricted, for which we propose two explanations.
first, many of the species analyzed, except the obligate biotrophs, do not follow only one lifestyle throughout their life cycle.
for example, even a sophisticated plant pathogen such as f. graminearum has an extended saprophytic growth phase during which it has to compete with soil microorganisms that will not play a role during the pathogenic phase.
thus, the secretome of such a fungus will show the imprint of adaptation to several habitats.
second, expansions and contractions of gene families and functional domains reflect a strong phylogenetic signal that may overrule that of lifestyle.
however, within each lifestyle, we find expansions in gene families with important roles.
in summary, we found only limited evidence for converging evolution to shape the overall composition of fungal secretomes.

 secreted proteins are essential for interactions with the environment, particularly for sessile organisms such as fungi. fungi exhibit a wide diversity of nutritional lifestyles, ranging from strict saprobes, living only of dead organic material, through necrotrophic pathogens to obligate biotrophic mutualists and pathogens, existing within the living tissue of a host. accordingly, fungal secreted proteins play important roles in the degradation of complex organic compounds, and in modulating directly or indirectly symbiotic and pathogenic interactions with their hosts. effectors are secreted proteins of plant and animal pathogens, which are involved in manipulating, directly or indirectly, defense reactions and metabolism of their hosts. for example, the maize pathogen ustilago maydis secretes a chorismate mutase that is translocated into infected maize cells to manipulate the shikimate pathway, which, amongst others, is involved in the production of salicylate, an important signal eliciting host defenses such as the systemic acquired resistance. other effectors prevent the detection of pathogen-associated molecular patterns, which are structurally conserved microbial molecules, by specific receptors of the host. for instance, chitin is an important structural component of fungal cell walls but also acts as a pamp. however, in compatible plant-pathogen interactions, chitin does not elicit effective host defenses as it is either modified by extracellular deacetylation to chitosan or shielded by chitin-binding effectors. thus, secreted proteins play essential roles in infection and determine whether a fungus may be an effective pathogen for a given host. most of the fungal secreted proteins are believed to be transported through the type ii secretion pathway. proteins that are destined for this pathway contain an n-terminal signal peptide that mediates recognition and transport across the endoplasmic reticulum membrane by a membrane-embedded translocon. the signal peptide is cleaved off and is thus not part of the mature protein. during transport through the er and the golgi apparatus the protein is folded, in many cases glycosylated and then travels from the golgi to one of several possible final destinations, one of which is the extracellular space. secreted proteins have been identified by molecular genetic and biochemical methods, for example the yeast signal sequence trap, but these typically identify only a fraction of the presumed total secretome. the signal peptide that is characteristic of the type ii secretion pathway is highly conserved among eukaryotes, which allowed the development of computer programs such as signalp, targetp, and phobius for its identification in protein sequences. this, coupled with the increasing availability of whole genome sequences for fungi, provides us with the resources to identify computationally the predicted set of secreted proteins in a large number of fungi representing various nutritional lifestyles and phylogenetic lineages. previous studies of fungal secretomes suggested that the composition of gene families varies greatly among species, and that some gene families appear to be associated with certain phylogenetic lineages and/or lifestyles. using gene model annotations from genome sequences, we performed a computational identification of soluble, secreted proteins, i.e. proteins containing a signal peptide, but no signal anchor or transmembrane domains, in species representing both a broad phylogenetic diversity as well as diversity in lifestyles. this allows assessing critically the hypothesis that species sharing a nutritional lifestyle have secretomes that exhibit a similar composition. in this context, we specifically address the questions whether secretome size and composition correlates with lifestyle and/or with phylogenetic lineages, whether such observed differences might be the result of expansion or contraction of certain gene families, and whether altered occurrences of functional units/domains, independent of gene families, might explain the differences.  prediction and characterization of fungal secretomes the annotated gene models for fungal and three phytopathogenic oomycota genomes were obtained from various databases. the fungal collection comprised ascomycota, eight basidiomycota, two zygomycota, one chytridiomycota, and one microsporidium. the three oomycota genomes were included as outgroup for the plant pathogens within the mycota. the selected fungal species did not only cover all phyla of the mycota but also represented a wide range of nutritional lifestyles. thus, species are saprophytes, eight are plant pathogens including one obligate biotroph, eight are vertebrate pathogens, one is an ectomycorrhizal, and one is a polyphagous fungus. this latter species, i.e. rhizopus oryzae, has been found in association with plants but is also described as a mammalian pathogen and a saprophyte.table species analyzed ashbya gossypii aspergillus fumigatus aspergillus nidulans batrachochytrium dendrobatidis botrytis cinerea candida albicans candida glabrata candida lusitaniae chaetomium globosum coccidioides immitis coprinus cinereus cryptococcus neoformans debaryomyces hansenii encephalitozoon cuniculi fusarium graminearum kluyveromyces lactis laccaria bicolor magnaporthe oryzae neurospora crassa phanerochaete chrysosporium phycomyces blakesleeanus phytophthora infestans phytophthora ramorum phytophthora sojae pichia stipitis postia placenta puccinia graminis rhizopus oryzae saccharomyces cerevisiae schizosaccharomyces pombe sclerotinia sclerotiorum sporobolomyces roseus stagonospora nodorum trichoderma reesei ustilago maydis yarrowia lipolytica names, codes used in figures, lifestyles, phyla, and subphyla of the species analyzed in this study, ordered alphabetically. proteins containing a type ii secretion signal sequence were identified using the computer programs signalp and, and targetp. a protein was only considered secreted when all three algorithms predicted this feature, according to, which yielded an average of % signal peptide-containing proteins per species. as we focused our study on proteins that are soluble, we subtracted all proteins carrying transmembrane domains from this initial set. these domains were identified by the program phobius that was specifically designed to distinguish between n-terminal transmembrane domains and “classical” type ii secretion signal sequences. this showed that the soluble secretome constituted on average % of the proteome with extreme values of % for schizosaccharomyces pombe and % for magnaporthe oryzae. when assessing these soluble secretomes in a phylogenetic context, pezizomycotina and basidiomycota had similar sizes, whereas the average size of saccharomycotina was noticeably smaller. comparing the soluble secretomes of fungi differing in nutritional lifestyle, saprophytes and animal pathogens had similar sizes, whereas the secretomes of plant pathogens were larger. these differences held true for both absolute and relative size. at the level of absolute size of the soluble secretomes, the species were grouped into three classes. class contained species with secretomes comprising less than proteins. class comprised species with secretomes ranging from ca. to proteins. class included the remaining six species with secretomes between ca. and proteins. in parallel, however, the species represented in these absolute secretome size classes also have increasingly large proteomes, which shows relative secretome size to be the more valid comparator for our data set.figure secretome sizes in relation to phylogeny and lifestyle. absolute and relative sizes of secretomes. a, b, d, e: species from phyla with at least three members. rare species: the only representative of a phylogenetic lineage in a nutritional lifestyle group. c, f: all species sharing a nutritional lifestyle except l. bicolor, as it is a symbiont, and r. oryzae, as it is polyphagous. a, d: species grouped by phylogeny. b, c, e, f: species grouped by nutritional lifestyle. diamonds: data points, colored according to classification opposite to boxes. boxes: lower and upper quartiles. internal line: median. lower and upper whiskers: minimum and maximum. asterisks: means. greyscales: phylogenetic groups. colors: lifestyle groups. comparing relative sizes of the soluble secretomes with phylogeny, within the ascomycota the saccharomycotina and taphrinomycotina have secretomes in the range of % to % of the proteomes with an average of %, whereas the pezizomycotina secretomes range from % to % of the proteomes with an average of %. the basidiomycota secretomes cover a broad size range from % to % of the proteomes, with an average of %. overall, ascomycota have relative secretome sizes from % to % with an average of % and thus cover a similar size range as the basidiomycota, with unicellular species at the lower and filamentous species at the higher ends. in relation to nutritional lifestyle, animal pathogen secretomes range from % to % of the proteomes with an average of %, whereas plant pathogen secretomes range from % to % of the proteomes with an average of % and thus are considerably larger. the sizes of the secretomes of saprophytes range from % to % with an average of %. when comparing, in dikarya, the unicellular or dimorphic yeasts with the filamentous fungi, it is obvious that the former group of taxa typically do not only have small absolute but also small relative secretomes. the large-proteome filamentous fungi, on the other hand, rather have large secretomes both absolutely and relatively. the plot of secretome size versus proteome size shows the general trend of increase of the former with increase of the latter. the similar, small secretome-to-proteome size ratios of the saccharomycotina are clearly seen in the form of a “dip” of the data points in the lower left corner of the graph, while for the other species, these show more variation. plots of secretome size vs. genome size and proteome size vs. genome size show similar, near-linear relationships with increasing variability for the larger genomes, secretomes, and proteomes. in the latter two comparisons, phytophthora infestans appears as an extreme case, which is based in the massive expansion of its genome, due to increased numbers of repetitive elements.figure correlations between the absolute sizes of genomes, proteomes, and soluble secretomes. absolute sizes of predicted soluble secretomes vs. predicted proteomes. absolute sizes of predicted soluble secretomes vs. genomes. absolute sizes of predicted proteomes vs. genomes. data points are color-coded by phylogeny. abbreviations for species are given in table  for each secretome, protein length distribution was analyzed by depicting the numbers of proteins in defined length intervals as fractions of the total secretome. the majority of proteins in each secretome had less than aa. however, proteins of up to aa residues were remarkably enriched in four of the five plant pathogenic pezizomycotina analyzed, i.e. botrytis cinerea, m. oryzae, stagonospora nodorum, and sclerotinia sclerotiorum, making up 10% to 15% of the total secretomes. a similar enrichment of such small proteins was further observed only in the ectomycorrhiza fungus laccaria bicolor, in the obligate biotrophic plant pathogen puccinia graminis, both belonging to the basidiomycota, and in the polyphagous zygomycete r. oryzae.figure protein size distributions within soluble secretomes. number of proteins in five length intervals as fractions of the total soluble secretomes. circle areas correspond to the absolute sizes of the soluble secretomes. species denominators are color-coded by lifestyle: dark blue = saprophytes, dark green = plant pathogens, red = animal pathogens, grey = polyphage, light blue = ectomycorrhiza, light green = obligate biotroph. abbreviations for species are given in table  we further divided the secretomes based on cysteine content and size since previous work indicated that putative effectors are often small cysteine-rich proteins. small secreted proteins are defined here to have a mature length of ≤  aa residues. if such proteins had a relative cysteine content of ≥5%, as well as ≥ cysteine residues, we considered them as small secreted cysteine-rich proteins. using these criteria, we found that cysteine-rich proteins were overall clearly enriched in the ssps, which was not obvious when applying a more relaxed definition for cys-enrichment using only a minimum of four residues but not considering the percentage of cysteine residues. class and size secretomes had more than two-fold and three-fold higher sscp-to-secretome ratios than the small class size secretomes, respectively, suggesting that the increase in secretome size may connect to a relative increase of sscp content. pezizomycotina and basidiomycota had close to four-fold higher sscp-to-secretome ratios than saccharomycotina, and notably higher ratios than oomycota. plant pathogens had a two-and-a-half-fold higher sscp-to-secretome ratio than animal pathogens, while saprophytes were in between and closer to plant pathogens.table relative cysteine contents by secretome size, phylogeny and lifestyle b) relaxed definition: cysteine-rich ≥ cysteines species c-rich total ssps sscps numbers of species and average contents of total cysteine-rich proteins, ssps and sscps in total secretomes in secretome size-, phylogenetic and lifestyle-groupings containing at least three species; ssp & sscp ≤  residues without signal peptides. identification of protein clusters by markov chain clustering all soluble secreted proteins identified were subjected to an “all-versus-all” blastp and this result was then used by the software tribe-mcl to arrange all sequences into protein clusters. thus, out of proteins were sorted into clusters, leaving singletons. as similarly seen for whole proteomes, certain clusters were either abundant in pezizomycotina or almost specific to them, as visualized by the intensity of red color, e.g. below the universal clusters in the middle of figure 4a. other protein clusters were either enriched in saccharomycotina, or practically absent. furthermore, we identified a large group of oomycota-specific clusters, and a number of clusters specific to individual species.figure tribe-mcl-clustering of proteins and distributions of interpro domains and molecular function and biological process go terms. heatmap of protein clusters with at least two proteins and interpro domains, molecular function go terms and biological process go terms with at least two occurrences by protein clusters/interpro domains/go terms and species. increasing number of proteins/occurrences per cluster/interpro domain/go term and species is indicated by higher intensity of red. species are color-coded by lifestyle: dark blue = saprophytes, dark green = plant pathogens, red = animal pathogens, grey = “polyphage, light blue = ectomycorrhiza, light green = obligate biotroph. species abbreviations are given in table  protein clusters with assigned interpro domains, molecular function or biological process go terms are provided on the right where the intensity of the blue color indicates their abundances. the heights of the heatmaps correspond to the numbers of individual protein clusters/interpro domains/go terms represented in the respective maps. to assign putative functions, all sequences were analyzed with interproscan. for the annotation of the secreted proteins, only interpro domains were considered. on average, 32% of the sequences had at least one interpro domain assigned. the highest percentages of annotated proteins were obtained for aspergillus nidulans, aspergillus fumigatus, cryptococcus neoformans and trichoderma reesei, whereas p. graminis, b. dendrobatidis and encephalitozoon cuniculi had the lowest percentages, possibly relating to the intensity of research devoted to different taxa. to assign interpro domains to the protein clusters determined by tribe-mcl, interpro domains found in members of a given cluster were assigned to the entire cluster if at least 75% of the sequences contained it. using this approach, out of protein clusters had at least one interpro domain assigned. moreover, we obtained gene ontology annotations of the types molecular function and biological process for the secreted proteins from the uniprotkb database. on average, 35% of all sequences had at least one molecular function go annotation and 21% had at least one biological process go annotation. for the molecular function go terms, the most extreme values for ratios of proteins in putative secretomes annotated with at least one term were 68% and 65% in candida albicans and s. pombe, respectively, and 8% and 14% in postia placenta and p. graminis, respectively, whereas for biological process go terms, the most extreme values were 79% and 76% of putative secretome proteins in saccharomyces cerevisiae and s. pombe, respectively, and 9% in both p. graminis and b. dendrobatidis. as for the interpro domains, we assigned go molecular function and biological process annotations to protein clusters when at least 75% of the proteins in a cluster shared the go term. thus, of the total of protein clusters, were assigned at least one molecular function go term and were assigned at least one biological process go term. we then analyzed whether certain protein clusters determined above by tribe-mcl occur in all species of a given phylogenetic lineage but do not occur in any species outside of this lineage. no cluster existed, that occurred in all fungi. two clusters contained proteins from all species analyzed except e. cuniculi, which is an intracellular pathogen of mammals exhibiting a strongly reduced genome and hence secretome size. the proteins grouped in those two clusters were disulphide isomerases, involved in protein folding and maturation in the er and undefined proteins. interestingly, no cluster existed that was absolutely ascomycota-specific. the most widely represented cluster within ascomycota contained proteins from of the species analyzed within this taxon but was not specific to it since it also occurred in other phyla. the corresponding protein was a mannosyl oligosaccharide glucosidase, the first enzyme in the n-linked oligosaccharide processing pathway, which was found exactly once per species except in chaetomium globosum. for saccharomycotina, pezizomycotina, and basidiomycota, one specific cluster each was found. the saccharomycotina-specific cluster comprised two proteins each in pichia stipitis and c. albicans and one protein in each of the other species. the cluster had no interpro domain or mf go or bp go term assigned, but of proteins contained a cfem domain according to interproscan that was described as an extracellular n-terminal domain comprising eight cysteines and that is found in certain fungal cell-surface integral membrane proteins. manual analysis showed that ten proteins belonging to cluster no. contained an authentic cfem domain, while one sequence had only six cysteines. only one of these proteins belonged to the sscps conforming to our new stringent definition given above. the cfem domain is also common in pezizomycotina and basidiomycota, but the corresponding proteins distribute to other clusters, none of which shows full coverage of a specific phylogenetic lineage as cluster does. overall, cfem domains were identified in proteins from protein clusters and one orphan protein. the pezizomycotina-specific cluster no. was represented by two proteins in t. reesei, fusarium graminearum, and neurospora crassa, and by one protein in each of the other species. according to our above definition, all these proteins were sscps. the single basidiomycota-specific cluster no. comprised one to three members in each of the species analyzed. we analyzed if certain clusters were exclusively associated with all species exhibiting a given nutritional lifestyle. we detected one particular cluster that was found in all plant-pathogenic pezizomycotina but neither in any other species of that subphylum having a different lifestyle nor in plant pathogens outside of the pezizomycotina. the saprophytic agaricomycetes coprinus cinereus, phanerochaete chrysosporium, and p. placenta had two clusters that encompassed only these species, i.e. no. and no.. on average 60% of the sequences in each secretome clustered with sequences from other species and another 11% with other sequences exclusively from the same species. despite this general trend, extreme examples existed where only very small fractions of the secretomes found similar members in other species. in particular, only 20%, 9% and 14% of the sequences in the secretomes of b. dendrobatidis, e. cuniculi, and p. graminis did so, while 59%, 21%, and 37% clustered exclusively within these species, respectively. this suggested that large parts of their secretomes are specific for these organisms or the lineages they belong to. importantly, the heatmap representation of the tribe-mcl results indicated that proteins rather cluster by phylogeny than by nutritional lifestyle. analysis of occurrence and abundance of interpro domains and molecular function and biological process go terms above, we assigned interpro domains and molecular function and biological process go terms identified in the proteins of a cluster to that cluster only if these were identified in at least 75% of the proteins forming that cluster. to analyze the secretomes based on functional traits, irrespective of the clustering of the proteins carrying them, we evaluated the presence and abundance of the interpro domains and go terms in relation to phylogeny and lifestyle. from totals of, and occurrences of interpro domains, molecular function and biological process go terms, respectively, around 97% of each were derived from multiple occurrences of and individual terms. for comparison, only of proteins were members of any of clusters. moreover, the proteins in clusters accounted for 90% each of all occurrences of the three annotation types. in the heatmap representations of abundances of interpro domains, molecular function and biological process go terms, distribution patterns were similar to protein distribution patterns but they were less prominent. for example, the band that spreads across the entire species spectrum in the middle of the protein cluster heatmap corresponds to the region roughly one third from the top in the interpro domain heatmap. this band appears broader for the interpro domains than for the protein clusters, because many of the highly conserved proteins in these clusters contain several different, equally conserved interpro domains. the saccharomycotina-enriched protein clusters in the upper region of the protein cluster heatmap are mirrored by saccharomycotina-enriched interpro domains in the lower region of the respective heatmap. similarly, species-specific protein clusters for b. dendrobatidis, r. oryzae, phycomyces blakesleeanus, p. graminis, and p. placenta in the upper region of figure 4a correspond to specific interpro domains for these species in the lower and upper regions of figure 4b. similar comparisons can be made for the molecular function and biological process go terms as well. interestingly, a few species show many more go annotations than all others do, e.g. s. pombe, s. cerevisiae and c. albicans, caused both by high ratios of the predicted secretomes having at least one annotation and many individual proteins having many different annotations.table distribution of proteins, interpro domains and go terms total numbers of proteins and of interpro domain, molecular function go and biological process go term occurrences identified in these proteins. fractions of proteins contained in protein clusters with at least two proteins and in clusters with proteins from at least three species. fractions of interpro domain, molecular function go and biological process go term occurrences identified in these protein clusters and fractions of total occurrences of these terms derived from individual terms occurring at least twice and in at least three species. next, we analyzed the interpro domains, i.e. functional features, for correlation with phylogeny or nutritional lifestyle. we identified no interpro domain present in all species, but five existed in all species except e. cuniculi. four of these five relate to disulfide-dithiol-mediated redox regulation. most individual occurrences of the four were identified in the proteins making up the disulphide isomerase protein cluster. the fifth interpro domain present in all species except e. cuniculi represented a domain of glycoside hydrolases. no other interpro domains were fully conserved within wide phylogenetic lineages. the two zygomycota exclusively shared two interpro domains, whereas the three oomycota shared comparing interpro domains to lifestyle, no correlation was observed at all. three molecular function go terms were conserved among all species, while by phylogeny, five others were conserved in all species except e. cuniculi. there existed only a single biological process go term conserved in all species, while two others were conserved in all species except e. cuniculi. by lifestyle, no molecular function or biological process go term was conserved. protein cluster contractions and expansions so far, we analyzed the absence or presence of protein clusters in all fungi sharing a given phylogenetic lineage or lifestyle. however, protein families may have contracted or expanded in individual species or parts of evolutionary lineages and such changes may relate to lifestyle. to reduce noise in quantitative assessment, expansion was defined as the presence of more than twice the median number of proteins for a cluster in a given species, contraction as less than half the median number of proteins. applying these criteria to the protein clusters present in at least three species, that thus contain at least three proteins to be able to calculate medians, we found clusters with expansions and/or contractions in at least one species. within the limits of the applied criteria, no protein cluster was contracted or expanded exclusively in all species belonging to a certain phylogenetic lineage, or exhibiting a certain nutritional lifestyle. even when scoring contractions or expansions less stringently, i.e. when allowing that this might have happened in additional species representing other lineages or lifestyles, we only found one protein cluster. this cluster was expanded in all five agaricomycotina, being present in five copies in the mammalian pathogen c. neoformans and up to seventeen in the saprophyte p. placenta. in other basidiomycota, it was expanded in the cereal pathogen p. graminis, in the maize pathogen u. maydis, but not in the saprophyte s. roseus. this is visualized in figure 5a by a red line just below the empty upper region in the basidiomycota that is gapped only under s. roseus. beyond, cluster was also expanded in the zygomycota p. blakesleeanus and r. oryzae, which had six and sixteen proteins, respectively. it had just one copy in four species of saccharomycotina whereas the remaining species entirely lacked it. except in c. globosum where it was also missing, species in the pezizomycotina had up to three proteins belonging to cluster figure expansions and contractions in protein clusters and in interpro domain and molecular function and biological process go term abundances. heatmap of protein clusters, interpro domains, molecular function go terms and biological process go terms represented in at least three species that show an expansion of more than twofold or a contraction of less than half of the median of the abundances per species for the respective entry type. increasing level of expansion or contraction per cluster/interpro domain/go term and species is indicated by higher color intensity. less than twofold expansion or contraction is shown in monochrome grey. white indicates no protein in this cluster or occurrence of this interpro domain or go term in this species. species are color-coded by lifestyle: dark blue = saprophytes, dark green = plant pathogens, red = animal pathogens, grey = polyphage, light blue = ectomycorrhiza, light green = obligate biotroph. species abbreviations are given in table  protein clusters with assigned interpro domains, molecular function or biological process go terms are provided on the right where the intensity of the blue color indicates their abundances. the heights of the heatmaps correspond to the numbers of individual protein clusters/interpro domains/go terms represented in the respective maps. the lack of general contractions even when using the less stringent condition results probably from the fact that in most clusters the median of the number of proteins was very low and that the absence of a protein in a species was not scored as a contraction. however, we did observe in five instances interesting trends of expansions or contractions in certain phylogenetic lineages. the first of these five protein clusters was described above as the only one showing a general expansion in a phylogenetic lineage and some basidiomycota and zygomycota. the second, apparently arose in the ascomycota since it was not present in any of the other lineages. the saccharomycotina had between two and five members in each species, whereas in the pezizomycotina the protein number was contracted to one each in f. graminearum, m. oryzae, coccidioides immitis, a. nidulans, and a. fumigatus, visible as a green line in the uppermost part of figure 5a. t. reesei, c. globosum, s. sclerotiorum, and s. nodorum lacked these proteins completely. the third protein cluster showed a trend of contraction in agaricomycetes. within basidiomycota, members of this cluster were restricted to species of this class with only one copy each, except p. placenta that had five. in the pezizomycotina, it was found with two to six proteins per species, except in c. immitis that had only one. all other lineages entirely lacked cluster the fourth cluster showed trends of expansion in three pezizomycotina. it contained proteins in each basidiomycota except u. maydis and in each pezizomycotina except f. graminearum, m. oryzae, and c. immitis. the basidiomycota and t. reesei, n. crassa, c. globosum, s. sclerotiorum and s. nodorum each had one to four of such proteins. the helotiales b. cinerea, the eurotiales a. nidulans and a. fumigatus showed expansions with five, five, and seven copies, respectively. all other lineages entirely lacked cluster. finally, cluster, was also only found in pezizomycotina and basidiomycota. in basidiomycota, only u. maydis and c. cinereus had this cluster, containing one and three members, and in ascomycota only seven species within the pezizomycotina with up to six copies in s. nodorum. the eurotiales a. nidulans and a. fumigatus, but also u. maydis showed a contraction in cluster with one copy each. when further focusing our analysis on individual species, we observed a number of expanded protein clusters on condition that a given cluster should comprise at least ten members in that species and at least 50% more members than in all other species having this cluster as well. disregarding the phytophthora species, of such cases were found in species of basidiomycota and saccharomycotina, but not of pezizomycotina, and further in b. dendrobatidis and r. oryzae. in general, members of basidiomycota showed the most extreme cases of such single-species expansions. contractions and expansions in abundance of interpro domains and of molecular function and biological process go terms as we analyzed the presence or absence of all interpro domains and molecular function and biological process go terms identified in the secretomes independent of the proteins in specific lineages or lifestyles, we also analyzed whether specific contractions or expansions occurred in these data sets. we applied the same limits for contraction or expansion of abundances of these terms as for protein clusters. thus, out of interpro domains, representing a total of individual interpro domain occurrences, showed a contraction or expansion in at least one species. furthermore, out of and molecular function and biological process go terms, representing and individual occurrences, and showed at least one expansion or contraction, respectively. as for protein clusters, abundance of no interpro domain was contracted or expanded exclusively in all species from any phylogenetic lineage or specific nutritional lifestyle or subgroup thereof. nevertheless, we observed trends of contraction or expansion. these, however, were mainly identified in the context of phylogeny. for those interpro domains that were identified exclusively or to at least 75% in the members of individual protein clusters, expansions and/or contractions of interpro domain abundances in multiple species typically paralleled the expansions and/or contractions in the corresponding protein clusters. the remaining, however, did not show such one-to-one associations and in the following, we highlight some noticeable cases. ipr was the most common domain in the data set and is linked to a large number of gh families. this domain was identified in all species except e. cuniculi and within fungi was mostly far more abundant in pezizomycotina, with to occurrences per species, and basidiomycota, with to occurrences per species, than in saccharomycotina, where it was contracted with three to nine occurrences per species. interestingly, b. dendrobatidis contained ipr only once. the animal pathogens c. immitis and c. neoformans and the saprotrophic basidiomycota s. roseus, with five or eight occurrences, also showed contractions for this interpro term within their phylogenetic lineages. of all the other interpro domains linked to ipr only ipr and ipr showed similar distributions. the latter is contracted in saccharomycotina. not surprisingly, molecular function go terms correlating with ipr particularly mf go: hydrolase activity, hydrolyzing o-glycosyl compounds, the most common molecular function go term in the data set, and mf go: hydrolase activity, showed similar patterns of contractions. these terms, however, also showed expansions in many pezizomycotina and basidiomycota, but no correlation with lifestyle existed for these expansions. as glycoside hydrolases are enzymes important for such diverse aspects as cell wall remodeling, degradation of host cell walls by plant pathogens to allow spreading through tissues and degradation of polysaccharides to provide carbon sources, all secreted proteins were analyzed against the cazy database by dbcan. in total, members of gh families were identified. whereas and gh families were identified in pezizomycotina and basidiomycota, respectively, saccharomycotina contained members in only gh families. in pezizomycotina, ghs ranged from proteins in families in the mammalian pathogen c. immitis to proteins in families in the plant pathogen s. nodorum. in saccharomycotina, gh content ranged from proteins in families in the mammalian pathogen candida lusitaniae to proteins in families in the mammalian pathogen c. albicans. in basidiomycota, ghs ranged from proteins in nine families in s. roseus to proteins in families in c. cinereus. with respect to nutritional lifestyle, gh content was relatively low in the specialized animal pathogens c. neoformans and c. immitis, but also in the biotrophic plant pathogen u. maydis and in the saprophytic yeast s. roseus. ipr, the second most abundant interpro domain in the entire data set, is a domain that is found in pectate lyase, iota-carrageenase and ghs from family including various galacturonases. within dikarya, it was absent from saccharomycotina except for contracted single occurrences in s. cerevisiae and ashbya gossypii, from s. pombe, and from s. roseus. in pezizomycotina, it was typically more abundant than in basidiomycota. expansion was seen in eurotiales, in helotiales and in f. graminearum. ipr was present in all ascomycota except debaryomyces hansenii, but only in four of eight basidiomycota. within pezizomycotina, it was expanded in all plant pathogens. this domain is contained in zinc-dependent metalloproteases from family m that are typically soluble secreted proteins, synthesized as zymogenic proproteins, that degrade proteins by hydrolyzing single residues from the c-termini. ipr describes an n-acetylglucosamine-binding domain that occurs in chitinases of fungi and plants. it is ubiquitous in pezizomycotina except c. immitis that lacks it, and is expanded in all plant pathogens. in both saccharomycotina and basidiomycota, however, only half of the species share it and in each, it is contracted with one occurrence only. as expected, the corresponding molecular function go term, go: chitin binding, shows practically the same pattern of abundances, contractions and expansions. interestingly, though, the expansions of chitin-binding domains in a number of mostly plant pathogenic pezizomycotina was not mirrored by the terms mf go: chitinase activity and bp go: chitin catabolic process, relating to enzymatic activity. either these species may have some chitin-active enzymes containing multiple chitin-binding domains, or these domains may occur independently of enzymes, for example in effectors shielding the fungi from recognition by their host plants. the last example of expansions in abundance of interpro domains are ipr gmc oxidoreductase, n-terminal and ipr gmc oxidoreductase, c-terminal. glucose methanol choline oxidoreductases are flavoenzymes utilizing various substrates. within fungi, these interpro domains were only found in pezizomycotina and basidiomycota. interestingly, while these domains tend towards expansion in the plant pathogenic pezizomycotina they are lacking in the specialized animal pathogens c. immitis and c. neoformans, respectively. c. cinereus has at least thrice as many occurrences of both gmc oxidoreductase interpro domains as any other basidiomycota species. as for protein clusters, we also observed a number of interpro domains, molecular function and biological process go terms that showed major expansions specifically in a single species, with more than ten occurrences and at least 50% more occurrences than in any other species having that interpro domain. many of these expansions reflected expansions in the protein clusters, annotated with these interpro domains or go terms. however, some prominent examples existed of expanded terms not linked to expanded protein clusters, such as ipr in the pezizomycotina m. oryzae with occurrences, and ipr in the basidiomycota c. cinereus with occurrences. discussion the completion of annotated proteomes allows for in depth analyses of secretomes by bioinformatics. our comparative analyses of fungal and oomycota secretomes, representing a broad range of taxa and nutritional lifestyles, revealed that primarily phylogenetic relationships determine the size and the composition of the secretomes, and to a lesser extent the nutritional lifestyle. we identified protein clusters that are specific for or expanded in certain lineages. others are expanded with respect to lifestyle, but only within individual lineages. secretome size and composition secretome size correlates with phylogenetic lineages. saccharomycotina have smaller secretomes, both absolute and relative, than pezizomycotina or basidiomycota, which have similar sized secretomes in both respects. this may be because yeasts are less capable of utilizing complex organic polymers as substrates, since they often live in niches providing access to sugars and other compounds that are easier to metabolize. secretome size is also associated with lifestyle. plant pathogens have the largest, saprophytes have intermediate, and animal pathogens have the smallest secretomes. to some extent, this finding may be explained by the assumption that animal hosts represent a nutritionally simpler environment than plant hosts or litter. furthermore, plant pathogens, particularly obligate biotrophs such as p. graminis, may secrete many effectors that would not provide benefits to saprophytes. it was suggested that portions of their secretomes participate in the coevolutionary arms race with host factors thus leading to sequence diversification. besides, p. graminis initially infects its alternate host berberis to switch thereafter to a cereal host. it must be presumed that the sets of secreted proteins involved in infection of both hosts do only partially overlap which will further increase the total size of the rusts’ secretome. moreover, the typically small size of animal pathogen secretomes may reflect an adaptation to the flexibility of the animal immune system. against a system effectively recognizing non-self molecules and mounting system-wide defense reactions in a circulatory system, a “stealth mode” is likely an efficient way of evasion. in such a scenario, it appears as a selective advantage to minimize the numbers of putative pamps and thus to reduce the size of the secretome. comparison of the mammalian pathogens within ascomycota points to another link between environmental influences and secretome size. the three candida species, c. immitis and a. fumigatus have secretomes of 203- and proteins, respectively. besides causing invasive mycoses in immunocompromised hosts, candida species are commensals on mammalian skin and mucosa. c. immitis is able to infect immunocompetent hosts and beyond lives as a saprophyte on host carcasses, with a concomitant lack of many enzymes required for decomposition of plant materials. aspergillus species typically are environmental saprophytes and only a few of them can in addition act as opportunistic pathogens of humans and animals. comparison of the two aspergillus species in our data sets shows that a. nidulans, having only the saprophytic lifestyle, has absolutely and relatively a larger secretome, and secretes more enzymes than a. fumigatus. these examples provide support for the notion that dominance of the saprophytic phase is associated with a larger secretome. furthermore, the higher number of secreted enzymes allows a. nidulans to utilize a wider variety of substrates present in soil or plant debris, which are likely more difficult to degrade than those available in animal hosts. in addition, competition with other microorganisms is conceivably less severe in animal hosts than in soil or litter, which may affect the number of secreted proteins involved in combatting competitors. the secretomes in the oomycota are considerably larger than in most pezizomycotina and basidiomycota and larger than most fungal plant pathogen secretomes in absolute, but not in relative numbers. the high absolute numbers of secreted proteins are derived, at least in part, from expansions in the rxlr- and crinkler-effector families which comprise a considerable fraction of these secretomes. we analyzed the secretomes for the distribution of size classes and cysteine richness to determine whether sscps, previously discussed to harbor many effectors, are enriched in certain phylogenetic lineages or nutritional lifestyles. examples for this type of effectors are some well characterized proteins from the biotrophic plant pathogen cladosporium fulvum. other sscps are hydrophobins covering aerial mycelia and spores, the membrane-bound cfem-proteins, antimicrobial proteins important to fight off competitors and peptides such as ptrtoxb of pyrenophora tritici-repentis targeting the host. using our definition of sscps our analysis shows a correlation of both phylogeny and lifestyle with the percentage of sscps in secretomes. by phylogeny, the relative sscp content is similar in pezizomycotina and in basidiomycota, while it is almost four-fold lower in saccharomycotina. as for lifestyle, the sscp fraction of the total secretomes is 30% higher in plant pathogens than in saprophytes, presumably because of the considerable numbers of potential effectors. in animal pathogens, this fraction is only half of that of saprophytes. the fact that the highest individual ssp and sscp fractions of the total secretomes occur in p. graminis and the hemibiotrophic pezizomycotina m. oryzae fits with the assumption that effectors would play eminent roles particularly in biotrophy. our results correspond well with the large secretomes and high numbers of ssps recently reported for obligate biotrophic phytopathogens melampsora larici-populina and p. graminis. presence/absence of and contractions/expansions in protein clusters, interpro domains and go terms we clustered all proteins into families, assigned interpro domain and molecular function and biological process go annotations, and identified clusters with expansions or contractions in specific taxonomic groups. it is important to keep in mind that what we define here as an expansion may be a recent lineage specific series of gene duplications or could be due to ancient gene duplications followed by gene loss in the other lineages. to distinguish between these two alternative hypotheses will require in depth phylogenetic analyses of each family. however, in either case, we consider that an expansion represents the presence of lineage specific selective pressure for those species to maintain the gene, regardless of whether the gene duplication events were recent or ancient. the secretomes of the saccharomycotina lack many protein clusters completely, which may relate to their unicellular lifestyle and overall reduced genomes. some of the clusters missing in saccharomycotina are expanded or contracted in basidiomycota compared to pezizomycotina or vice versa. apart from reflecting the specific phylogenetic histories of the corresponding species, expansion of protein families may correlate with their abilities to cope with more complex and/or variable environments. assessment of distinct phylogenetic lineages revealed that saccharomycotina, pezizomycotina, and basidiomycota each contain a single lineage-specific protein cluster. the saccharomycotina-specific cluster contains the cfem motif, but narrowly misses the threshold for interpro domain assignment and has no go terms assigned. the cfem domain is not saccharomycotina-specific, as we identified it in protein clusters across the dikarya. this example illustrates the intrinsic difficulties in analysis of protein functions by assigning interpro domains or go terms to protein clusters. some clusters have not been assigned terms because the assignment threshold was missed or proteins containing the same terms were clustered in different protein clusters. for this reason, we also performed analyses for the interpro domains and molecular function and biological process go terms themselves. four of five interpro domains absent only in e. cuniculi relate to er- or golgi-localized proteins, which may reflect the major reduction in secretion in this intracellular parasite. the fifth represents the catalytic core of glycoside hydrolases. no other interpro domain showed absolute specificities either by phylogeny or by lifestyle. by our definition, we found no protein cluster or interpro domain either contracted or expanded exclusively within all species representing a particular lifestyle or phylogenetic lineage. nevertheless, five protein clusters were expanded or contracted in some species in certain phylogenetic lineages. for interpro domain abundances, we found two terms to be contracted in saccharomycotina and further observed trends of expansions, as well as contractions, in subsets of species within phylogenetic lineages. in saccharomycotina, the glycoside hydrolase catalytic domain and glycoside hydrolase, family were contracted according to the definition applied, consistent with the go- and particularly the cazy-data. this may be indicative of a narrower spectrum of exploitable polysaccharides in this taxon. two domains are expanded in all plant pathogenic pezizomycotina. some of the zinc-metalloproteases from family m, that typically degrade extracellular matrix proteins, may be involved in degradation of plant tissue, others however could mediate the suppression of plant defense by degradation of host proteins involved in pathogen recognition in the apoplast or playing other important roles in defense. the other expanded ipr domain is the chitin-binding domain, type that occurs in single or multiple copies in fungal chitinases and other glycoside hydrolases. in plant pathogenic fungi, it may mediate binding of the containing proteins to the fungal cell wall to avoid host recognition by shielding or modifying chitin, which is a pamp as discussed above. in addition, chitinases or other glycoside hydrolases containing this domain could serve host cell wall degradation or fungal morphogenesis. several interpro domains for glycoside hydrolases and other enzymes involved in degradation, as well as go terms relating to such activities, show considerable tendencies towards expansion in saprophytes and plant pathogens in the pezizomycotina and basidiomycota. for glycoside hydrolases, the cazy-data confirmed this tendency. the likely explanation is that these species can utilize a wider variety of complex carbohydrates than yeasts or animal pathogens. examples are iprs and that are particularly expanded in plant pathogenic pezizomycotina plus the saprotrophic basidiomycota c. cinereus. proteins containing these domains, including various alcohol and sugar peroxidases, produce hydrogen peroxide, which may amongst others modify target molecules directly or serve as substrate for oxygen radical generation by fenton chemistry. radicals may serve the breakdown of lignin in white rot fungi or of polysaccharides as in cellulose degradation by brown rot fungi. the role of peroxide generated by plant pathogens is diverse. for example, it may interfere with plant defense compounds by direct reduction, as was recently proposed for a fad-dependent glucose dehydrogenase from the pezizomycotina glomerella cingulata that causes anthracnose on various fruits. apparently, only few protein clusters specific for lineages exist, which may be explained by the broad phylogenetic distance covered by the fungi analyzed thus obstructing the clustering of related proteins. in addition, particularly in the case of plant pathogenic fungi, proteins required for host interaction may be under positive selection due to the “arms-race” between host and pathogen. thus they may not be recognized as originating from the same ancestral gene anymore, even between relatively closely related species as shown for three species of powdery mildew fungi that only shared less than effectors. finally, losses or gains of protein families need not be equivalent with losses or gains of biological functions, as these may be provided by other protein families. it has been shown that changes in lifestyle have occurred often in fungal evolution and such events need not be accompanied by the same losses or gains of proteins. with respect to expansion of protein clusters in individual species only, we observed for example an expansion of peptidases and serine endopeptidases of family s that cleave target proteins close to alanine-rich c-terminal tripeptides ) in b. dendrobatidis, confirming previous observations. this may reflect an adaptation to the animal pathogenic lifestyle and similar expansions were observed for peptidases in other fungi exhibiting this lifestyle. a study comparing closely related animal pathogens and saprophytes showed that in c. immitis and c. posadasii families of subtilisins and other family s serine peptidases, which typically are non-specific endopeptidases, and mentioned m metalloproteases are expanded. these families are also expanded in the non-pathogen uncinocarpus reesii, but not in the more distantly related histoplasma capsulatum. in contrast to these onygenales, within the eurotiales these peptidase families are expanded neither in the saprophyte a. nidulans nor in the opportunistic pathogens a. fumigatus and a. terreus, illustrating different extents of adaptation to the pathogenic lifestyle. complementary, whereas the onygenales show strong contractions in, or even lack of a number of plant cell wall degrading enzyme families, the helotiales all contain similar numbers of these enzymes. this suggests that there is a considerable phylogenetic signal in the data of species comparisons irrespective of the degree of kinships and furthermore, that lifestyles are not linked to the presence, absence or abundance of certain protein families of functional groups in a straight-forward causal relationship.  we identified and analyzed the soluble secretomes of fungi to test the hypothesis that the secretomes of species belonging to different phylogenetic lineages but sharing the same nutritional lifestyle also share a similar composition. secretome size correlates with lifestyle, with plant pathogens having the largest and animal pathogens having the smallest secretomes. in contrast, the signal of lifestyle on the specific composition of secretomes is restricted, for which we propose two explanations. first, many of the species analyzed, except the obligate biotrophs, do not follow only one lifestyle throughout their life cycle. for example, even a sophisticated plant pathogen such as f. graminearum has an extended saprophytic growth phase during which it has to compete with soil microorganisms that will not play a role during the pathogenic phase. thus, the secretome of such a fungus will show the imprint of adaptation to several habitats. second, expansions and contractions of gene families and functional domains reflect a strong phylogenetic signal that may overrule that of lifestyle. however, within each lifestyle, we find expansions in gene families with important roles. in summary, we found only limited evidence for converging evolution to shape the overall composition of fungal secretomes. 
expansions and contractions of gene families and in the number of occurrences of functional domains are largely lineage specific, e.g. while the secretomes of fungi with similar lifestyles share certain characteristics, the expansion and contraction of gene families is largely lineage specific, and not shared among all fungi of a given lifestyle. secretome size correlates with phylogeny and to a lesser extent with lifestyle. it is generally assumed that fungi with similar nutritional lifestyles have similar secretome compositions. in this study, we test this hypothesis by annotating and comparing the soluble secretomes, defined as the sets of proteins containing classical signal peptides but lacking transmembrane domains of fungi representing a broad diversity of nutritional lifestyles. however, within a given lifestyle a few general trends exist, such as the expansion of secreted family m metallopeptidases and chitin-binding proteins in plant pathogenic pezizomycotina. small secreted cysteine-rich proteins, which may comprise many effectors important for the interaction of plant pathogens with their hosts, are defined here to have a mature length of ≤  aa residues, at least four cysteines, and a total cysteine content of ≥5%. plant pathogens and saprophytes have larger secretomes than animal pathogens. relative sscp content is noticeably higher in plant pathogens than in animal pathogens, while saprophytes were in between and closer to plant pathogens. contraction of glycoside hydrolases in saccharomycotina, and are only weakly correlated with lifestyle.
 since the first fungal genome sequences became available, investigators have been employing comparative genomics to understand how fungi have evolved to occupy diverse ecological niches. the secretome, i.e. the entirety of all proteins secreted by an organism, is of particular importance, as by these proteins fungi acquire nutrients and communicate with their surroundings.  it is generally assumed that fungi with similar nutritional lifestyles have similar secretome compositions. in this study, we test this hypothesis by annotating and comparing the soluble secretomes, defined as the sets of proteins containing classical signal peptides but lacking transmembrane domains of fungi representing a broad diversity of nutritional lifestyles. secretome size correlates with phylogeny and to a lesser extent with lifestyle. plant pathogens and saprophytes have larger secretomes than animal pathogens. small secreted cysteine-rich proteins, which may comprise many effectors important for the interaction of plant pathogens with their hosts, are defined here to have a mature length of ≤  aa residues, at least four cysteines, and a total cysteine content of ≥5%. sscps are found enriched in the secretomes of the pezizomycotina and basidiomycota in comparison to saccharomycotina. relative sscp content is noticeably higher in plant pathogens than in animal pathogens, while saprophytes were in between and closer to plant pathogens. expansions and contractions of gene families and in the number of occurrences of functional domains are largely lineage specific, e.g. contraction of glycoside hydrolases in saccharomycotina, and are only weakly correlated with lifestyle. however, within a given lifestyle a few general trends exist, such as the expansion of secreted family m metallopeptidases and chitin-binding proteins in plant pathogenic pezizomycotina.  while the secretomes of fungi with similar lifestyles share certain characteristics, the expansion and contraction of gene families is largely lineage specific, and not shared among all fungi of a given lifestyle. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords fungisecretomecysteine-richsspsscpclustercontractionexpansionphylogenynutritional lifestyleissue-copyright-statement© the author 2014 