BACKGROUND
over the last decades, the organization of the secretory and endocytic pathways has been studied in a variety of different eukaryotic organisms.
these studies have revealed that the principal organization of the endomembrane system and the molecular machineries involved in vesicular trafficking are conserved among all eukaryotes.
in general, the transport between different intracellular organelles is mediated by cargo-laden vesicles that bud from a donor and then specifically fuse with an acceptor compartment.
key players during the final fusion step are the so-called snare proteins.
these proteins are a large family of small cytoplasmically orientated membrane proteins that are typically tail-anchored.
their key characteristic is the so-called snare domain, an extended stretch of heptad repeats that is usually connected to a single transmembrane domain by a short linker.
snare domains of heterologous sets of snare proteins have the ability to assemble into tight, parallel four-helix bundle complexes.
it is thought that complex formation between snare proteins from opposing membranes provide the energy that drives membrane fusion .
although the concept of snare-mediated membrane fusion emerged in the 1990s  <cit> , several snare proteins had been discovered earlier.
the field of membrane trafficking was greatly inspired by ingenious genetic screens in the single-cell fungi saccharomyces cerevisiae for mutants with defects in the secretory pathway in the late 1970s  <cit> .
in these mutants, secretion was blocked at higher temperatures.
these so-called sec mutants accumulate secretory proteins at the point in the secretory pathway that is blocked and often, they could be distinguished by their phenotype.
in the following years, it became clear that sec genes generally encode for key components of the machineries that mediate trafficking of a transport vesicle, including a few snare proteins.
all at once, straightforward genetic tools had catapulted baker's yeast to the front row for studying vesicle trafficking pathways.
in fact, other genetic screens for mutants showing defects in other transport processes  have turned up novel yeast snares.
the last yeast snare was tracked down in  <dig>  yielding an overall repertoire of  <dig> different snare genes in baker's yeast  <cit>  .
for most trafficking steps in baker's yeast, distinct units of four interacting snare domains have now been assigned .
however, not all details of the trafficking routes in yeast are resolved.
generally, different snare family members are localized on distinct organelles and vesicles that demarcate particular membranes.
however, it is becoming clear that some trafficking steps cannot be identified by a unique set of snares, since some snares are involved in more than one trafficking step.
in addition, they may participate in different snare complexes .
to appraise s. cerevisiae as model organism for studying vesicle trafficking pathways better, it is also necessary to understand the limits of the organism, since baker's yeast appears to be a secondarily reduced organism.
for example, s. cerevisiae has only a limited ability to produce hyphae .
the capability to produce polarized hyphae, however, is a hallmark of the fungal kingdom, allowing them to colonize and exploit new substrates efficiently.
in addition, whereas most fungi are multicellular organisms, s. cerevisiae most of the times grows as single cells that reproduce asexually.
it thus seems possible that adaptation to a relatively simple lifestyle was accompanied by a degeneration of the intracellular trafficking itinerary together with the involved machineries.
a comparison of the snare gene repertoires of s. cerevisiae with that of five other fungal species, including three species of filamentous fungi, revealed that, in general, the members of the snare family are largely conserved in fungi  <cit> .
this notion has been corroborated by subsequent inspections of the snare repertoire of few other fungal species  <cit> .
together, these studies put forward that the yeast lifestyle did not entail a radical change in the intracellular trafficking pathways.
the bioinformatic strategies used to identify the homologs of the yeast snares, however, did not provide a universal classification scheme.
consequently, it was impossible to compare the snare sets of different fungi with each other in an unambiguous manner, let alone comparing them with snare sets from other eukaryotes like animals.
hence it remains unclear whether snare types, which are possibly linked to a novel transport step, evolved or degenerated in particular fungal lineages.
fungi are thought to be more closely related to animals than to plants and, currently, are placed with animals and several protistan taxa into the monophyletic group of opisthokonta.
notably, morphological observations show that the appearance of analogous intracellular organelles such as the vacuoles/lysosomes and the golgi apparatus are markedly different between yeast and animal cells.
it is unclear, however, whether these differences are reflected in their respective snare inventory.
recently, using hidden markov models  we have established a precise classification for all eukaryotic snare proteins  <cit> .
based on this classification, we were able to deduce that the emergence of multicellularity in animals went along with an expansion of the set of snare proteins, in particular the expansion of the snares involved in endosomal trafficking  <cit> .
in order to analyze the differences between the snare repertoires of fungi and animals and the evolutionary changes in the snare inventory of different fungal lineages, we have collected and classified the snares from a broad range of different fungi according to our hmms, making use of the substantial number of available genomic sequences.
to identify specification and duplication events, and to pinpoint the differences between the set of snare proteins in fungi and animals, we used phylogenetic information to reveal the history of snare genes.
afterwards, we mapped features of the fungal snare evolution onto the species tree.
methods
sequences
recently, we published a classification for snare sequences, based on more than  <dig> protein sequences from  <dig> eukaryotic species, including  <dig> sequences from  <dig> species of fungi.
we were able to identify twenty basic subtypes within the snare protein family and developed specific as well as sensitive hmms for each subtype  <cit> .
we used these models to identify snare proteins in fungal genomes that became available recently.
the models yielded  <dig> additional sequences from the nr-database at ncbi and various genome projects , baylor college of medicine, j. craig venter institute and broad institute) and several est databases   <cit> .
the newly found protein sequences underwent a visual verification process to remove redundancy and misassembled sequences.
overall, our data contained  <dig> unique snare sequences from  <dig> fungal species.
the dataset was incorporated into the snare database  <cit>  and can be accessed on the project homepage .
the snare sets of all fungal species with complete genomes sequenced can be found in .
for a more detailed inspection of the evolutionary changes within fungi, we made use of the observation that within a subtype, sequence regions besides the snare motif were also well conserved.
consequently, we constructed alignments of the entire sequences for each subgroup using muscle  <cit> .
sites with more than  <dig> percent gaps were removed from the alignments.
in addition, we excluded all sequences which contained more than  <dig> percent gaps.
these alignments are also available via the supplementary section of our projects homepage.
based on these alignments, we constructed phylogenetic trees for each of the fungal snares .
we also generated sequence alignments in which we graphically incorporated information about the domain structure based on the available high resolution structure information .
phylogeny
the phylogenetic reconstruction was composed of two different analytical approaches and was basically done as described in  <cit> .
the first approach used iqpnni   <cit>  to construct phylogenetic trees from the curated alignments.
we used a gamma distribution as a model for rate heterogeneity with four rate categories for the estimation of the gamma distribution parameter.
the proportion of invariable sites was estimated from the data and the jones, taylor and thornton  distance matrix  <cit>  served as a substitution matrix.
we used the stopping rule of iqpnni, but the calculation had to run for at least the suggested number of iterations.
the default values were used for the remaining parameter.
in addition, likelihood-mapping was applied to determine the confidence of the edges in the calculated trees.
the second approach used the phylip package  <cit>  to apply a distance-based bootstrap analysis with  <dig> replicates to each of the curated alignments.
standard settings were used for seqboot, the jtt distance matrix and also a gamma distribution  for protdist, as were standard options for neighbour.
if required, the random seed was set to nine.
we used the almost unbiased  test  <cit>  to address the systematically biased bootstrap values.
we obtained the sitewise log-likelihoods needed for the au test using a modified version of phyml  <cit>  and the test was performed using consel  <cit> .
the reconstructed iqpnni trees served as starting points to join the results of both calculations.
the inner edges of the trees were labeled with their likelihood-mapping and corrected bootstrap support values.
we provide nexus files for each reconstructed tree on our project homepage .
the trees can be interactively explored using, for example, splitstree  <dig>  <cit> .
to gain initial insights on the placement of fungi within the opisthokonta, especially the basal fungal species batrachochytrium dendrobatidis and blastocladiella emersonii, we first performed an extensive phylogenetic analysis based only on the sequence of the snare domain.
this analysis was carried out to ensure that the evolutionary history of the individual snare groups indeed reflects the monophyletic origin of fungi.
for this analysis, we used all fungal species together with the choanoflagellate monosiga brevicollis, the slime mold dictyostelium discoideum, and ten selected animal species .
we also included the snares found in the animal pathogen encephalitozoon cuniculi, which belongs to the microsporidia.
in earlier times, this clade was considered to be of protist origin, but it is now believed to have diverged from early fungi.
however, we purposefully did not include their sequences into the general phylogenetic analysis, as e. cuniculi appears to have a markedly reduced set of snare proteins  <cit>  with clearly deviating sequences, making them sensitive to long branching artefacts.
the reconstructed phylogenetic trees demonstrate that fungi  are well separated from other opisthokonts .
with this corroborated, we then performed a more detailed phylogenetic analysis of the snares in fungi using the entire sequence.
RESULTS
in order to catalog the snare inventory of fungi, we increased our previous collection  <cit>  with newly published genomes, and complemented this set with ests from different fungal species, obtaining an overall set of more than  <dig> snare sequences from  <dig> fungal species.
the dataset comprises about  <dig> species with nearly complete sets of snares .
all sequences and species used in this study can also be found in our snare database on our projects homepage .
the collected sequences were then classified according to our previously established hmm models  <cit>  that are based only on the highly conserved snare motif.
the four main classes were then further subdivided into twenty distinct groups according to our hmm models.
each of these subgroups represents a distinct conserved trafficking step.
for phylogenetic analysis of the fungal lineage, we aligned entire sequences for each subgroup, as evolutionary changes have usually occurred by paralogous expansion of particular snare types.
based on these alignments, we constructed phylogenetic trees for each of the snare types in fungi .
fungi contain a compact set of snare proteins
overall, the fungal genomes examined in this work comprise a remarkably consistent and homologous set of snare proteins.
in each species with a complete genome, we only found a little more than twenty different snares, whereas in metazoans, for example species with more than  <dig> snares are common.
a marked expansion of the metazoan set of snares occurred during the rise of multicellularity  <cit> .
in earlier studies it was put forward that the development of multicellularity is generally linked to an expansion of the snares  <cit> .
in contrast to this idea, we usually only found members of the basic snare proteins to be present in each fungi species.
this implies that fungi, although they also have evolved multicellularity, have a relatively unchanged set of snares and have not significantly expanded the repertoire of the proto-eukaryotic cell.
it was not the objective of the present study to discuss the function of each snare protein coded for in fungal genomes in detail, as the sequence information alone is insufficient for recognizing the exact trafficking step the protein is involved in .
instead, we will mainly focus on highlighting the idiosyncrasies that we came across.
notably, the snare proteins involved in secretory trafficking routes between the endoplasmatic reticulum , the golgi apparatus  and the trans-golgi network  were generally present as singletons and were highly conserved in all fungi.
noticeable changes affected the set of snare proteins involved in secretion and in the endosomal/vacuolar trafficking pathway.
an interesting observation is that the set of snare proteins in fungi is very similar to the set of the single-cell choanoflagellate monosiga brevicollis, which is thought to be closely related to animals.
in contrast to m. brevicollis, lower animals possess an enlarged set of snare proteins  <cit> .
a surprising aspect of the chytrids batrachochytrium dendrobatidis and blastocladiella emersonii is that they possess a snare that was classified as qb.iii.d-type  <cit> .
snares belonging to this subclass were first found in green plants and termed novel plant snares   <cit> .
our previous analysis had shown that npsn can be found in plants and several protist lineages, but not in choanoflagellates and animals  <cit> .
the fact that basal fungi retained a bona fide npsn suggests that late-diverging fungi have lost npsn independently.
interestingly, the chytrids are the most primitive of the fungi; many are aquatic and possess a flagellum.
consistent with their more basal position, sequences from the chytrids batrachochytrium dendrobatidis and blastocladiella emersonii are usually placed close to the root of the fungal lineage in our analysis.
these results suggest that npsn was a member of the original snare repertoire of the eukaryotic cenancestor but was then lost in the choanoflagellate/animal and the fungal lineage.
another protein family has been widely used as a marker for the different trafficking steps of the endomembrane system.
these are rab proteins, a family of small gtp-binding proteins that regulate different vesicular trafficking steps by cycling between active gtp and inactive gdp forms  <cit> .
it has been observed that fungi contain a relatively small number of rab proteins  <cit> , whereas their number has expanded in several other lineages such as animals  <cit> , land plants <cit> , and some protists  <cit> .
it would therefore be interesting to investigate whether the pattern of duplications and diversifications found for the rab family in fungi  <cit>  correlates with that of the snare family  <cit> .
the snare protein vam <dig> is a molecular "invention" of fungi
our analysis shows that fungi contain three different snares of the qc.iii-type: tlg <dig>  syntaxin  <dig>  and vam <dig>  <cit> .
the choanoflagellate m. brevicollis and basal animals contain only two qc.iii-snares, syx <dig> and syx <dig>  which are homologous to tlg <dig> and syx <dig>  respectively.
remarkably, we found vam <dig> homologs to be present in all complete fungal genomes, even in the basal fungus b. dendrobatidis.
however, no vam <dig> homolog was discovered in other eukaryotes, suggesting that the qc.iii-snare vam <dig> is a new innovation in fungi.
of note, this status of vam <dig> had been recognized earlier by yoshizawa and colleagues  <cit> .
however, one should keep in mind that their analysis was based on an inferior cluster approach that led, from a biological point of view, to a number of wrongly assigned clusters.
in addition, they have missed several well-known snare proteins.
although the snare motif of vam <dig> is not greatly diverged, it remains unclear whether an ancestral tlg <dig> or syx <dig> gave rise to vam <dig>  nonetheless, considering its unique domain structure, vam <dig> is the most divergent snare protein in the repertoire of fungi.
in contrast to most other snare proteins, vam <dig> does not carry a transmembrane region at the c-terminal end of its snare motif.
instead, the protein seems to be anchored to a membrane by means of an n-terminal extension that harbors a phox homology  domain, which acts as phosphoinositide binding module.
in fact, vam <dig> has been shown to interact specifically with phosphatidylinositol 3-phosphate p)  <cit> , a lipid that is specific for endosomal and vacuolar membranes.
a large body of work carried out in yeast implies that vam <dig> catalyzes homotypic vacuole fusion with its partner snares, vam <dig> , vti <dig>  and nyv <dig>   <cit> .
interestingly, however, the homologous factor in the basidiomycete ustilago maydis, termed yup <dig>  was found to be present on vacuoles and rapidly moving endosomal organelles.
moreover, a temperature sensitive mutation in yup <dig> blocks sorting to vacuoles and a defect in endosomal recycling  <cit>  suggests that the function of vam <dig> is not restricted to homotypic vacuole fusion, a process that, as we will discuss below, may be restricted to saccharomycotina.
changes in the endosomal and vacuolar set of snare proteins in saccharomycotina
the yeast s. cerevisiae belongs to the saccharomycotina clade, the members of which are sometimes referred to as hemiascomycetes.
s. cerevisiae possesses two endosomal syntaxins , pep <dig> and vam <dig>  which are believed to be involved in trafficking to late endosomes and to vacuoles, respectively.
in agreement with the earlier assessment by gupta and heath  <cit> , our phylogenetic tree shows that these two snares have been duplicated in the lineage of saccharomycotina only .
in addition, our analysis also demonstrates that the duplication of ancestral qa.iii.b-snare into pep <dig> and vam <dig> is independent of the duplications of the endosomal syntaxins in animals, i.e.
syntaxin  <dig>  syntaxin  <dig>  syntaxin  <dig> and syntaxin  <dig>
the bifurcation of an ancestral endosomal syntaxin within the saccharomycotina clade into proteins with distinct roles, namely pep <dig> and vam <dig>  may suggest that the specialized machinery involved in homotypic vacuole fusion arose only within the hemiascomycetes.
homotypic vacuole fusion is a process that occurs during cell division in baker's yeast.
during budding, the few large vacuoles maintained in s. cerevisiae are fragmented.
multiple small vacuoles then segregate into the budding daughter cell and subsequently fuse again into larger vacuoles  <cit> .
this phenomenon has been exploited by william wickner and colleagues  <cit>  to study the role of snares and other factors in the fusion process.
this process, however, does not take place in, for example, s. pombe, a yeast that recently was proposed to belong to a third ascomycetes class, the archiascomycotina  <cit> .
this suggests that homotypic vacuole fusion, described first in baker's yeast, is a process specific to hemiascomycetes.
notably, the two qa.iii.b-snares of yarrowia lipolytica are much less diverged.
consequently, it seems interesting to find out at which point in the evolution of the saccharomycotina lineage the process of homotypic vacuole fusion evolved.
in addition, the sets of proteins catalyzing trafficking to the vacuole  and maintaining a low vacuole copy number  genes have only partial functional overlap  <cit> .
it would thus be interesting to investigate whether the factors specifically involved in the process of homotypic vacuole fusion are restricted to the saccharomycotina.
previously, it has been debated whether the r-snare involved in trafficking towards vacuoles in yeast, nyv <dig>  is indeed homologous to the animal r-snare vamp <dig> .
in fact, the snare motif of nyv <dig> from baker's yeast is only distantly related to vamp <dig> in animals  <cit> ; however, in our analysis, both clearly fall into the same group .
it has even been thought that nyv <dig> does not contain a canonical profilin-like n-terminal extension, often also referred to as longin domain  <cit> .
a recent structural analysis has demonstrated that the n-terminal extension of nyv <dig> indeed contains a profilin-like fold of this type  <cit> , similar to the ones found in sec <dig>   <cit>  and ykt <dig>   <cit> .
in general, all fungi possess a nyv <dig> homolog.
the only fungal lineages that do not appear to contain a nyv <dig> homolog are members of schizosaccharomyces .
our phylogenetic analysis shows that nyv <dig> sequences are remarkably derived in saccharomycetales, bringing to light why it has been difficult before to acknowledge nyv <dig> as a bona fide endosomal r-snare.
in s. cerevisiae, nyv <dig> seems to interact only with the q.iii.b-type syntaxin vam <dig> but not with pep <dig>  which is thought to interact preferentially with the r-snare ykt <dig>  hence, it seems possible that nyv <dig> and vam <dig> coevolved in saccharomycotina.
the larger set of secretory snares in s. cerevisiae can be attributed to a whole genome duplication
in contrast to the notion that baker's yeast is a reduced organism, no fungal genome inspected contained more snare genes than s. cerevisiae with its  <dig> different snare factors .
these  <dig> genes are clearly more than the 20/ <dig> genuine snares of the proto-eukaryotic cell.
two of the added snares vam <dig>  a novel, fungi-specific snare, and vam <dig>  a saccharomycotina-specific snare were discussed above.
in addition, baker's yeast contains four closely related pairs of secretory snare proteins, sso <dig> and sso <dig> , snc <dig> and snc <dig> , sec <dig> and spo <dig> , and sro <dig> and sro <dig> .
our analysis demonstrates that these doubled set of secretory factors very probably arose from a well established whole genome duplication  in s. cerevisiae  <cit> .
indeed, we found the same duplicated set of snare proteins in the closely related species s. paradoxus, s. mikatae, s. bayanus and s. kudriavzevii, all of which are collectively referred to as saccharomyces senso-stricto species .
remarkably, some but not all of the same secretory snare genes were also found in duplicates in s. castellii, candida glabrata, and vanderwaltozyma polyspora .
the ancestor of this group of species  and of the saccharomyces sensu stricto species are believed to have undergone the same wgd followed by major secondary reductions, while several closely related species like eremothecium gossypii , kluyveromyces waltii, kluyveromyces lactis, and saccharomyces kluyveri have kept a single copy of the genome.
these species are grouped with the ones that underwent a wgd as, "saccharomyces complex" <cit> .
indeed, we found only one copy of each of the secretory snares in the latter species.
interestingly, the duplicates found in s. castellii, c. glabrata and v. polyspora usually split at clearly different positions in our phylogenetic trees as compared to the duplicates of the saccharomyces sensu stricto species.
this is particularly noticeable for the two sec <dig> proteins from s. castellii and also from v. polyspora that are highly similar, whereas the two qbc-snares, sec <dig> and spo <dig>  exhibit marked asymmetrical rates of evolution in the saccharomyces sensu stricto species, pointing to a different fate of the duplicated secretory snares in the two lineages.
similarly, the two sso proteins of c. glabrata are clearly more diverged than the duplicated set of the saccharomyces sensu stricto species.
in fact, it has been proposed that these lineages must have diverged and adapted to different living conditions very soon after the wgd, so that the genes retained in duplicates diverged largely independently  <cit> .
although only about 14% of all duplicated genes were retained after the wgd, all duplicated snare genes are involved in the process of secretion in baker's yeast.
this suggests that these factors presented a selective advantage for yeast.
the selective advantage of the duplicated secretory set is not entirely clear, since so far, only subtle functional differences of the duplicated snares have been found.
as mentioned above, only one homologous pair, the two qbc-snares  is highly diverged in the saccharomyces sensu-stricto species.
sec <dig> and spo <dig> possess only approximately 37% identity, whereas the homologs of the secretory syntaxins  and r-snares  are 79% and 74% identical, respectively.
whereas sec <dig>  particularly in saccharomycotina, possesses a rather long n-terminal extension, this region was reduced in spo <dig>  <cit> .
interestingly, sec <dig> and spo <dig> are specialized to different developmental stages in baker's yeast.
whereas sec <dig> interacts with both secretory syntaxins and synaptobrevins in secretion during vegetative growth, the more deviated spo <dig> is required only for sporulation.
during sporulation in s. cerevisiae, the prospore membrane, which envelops each daughter nucleus during meiosis, is generated de novo adjacent to the spindle pole bodies by the secretory snare machinery  <cit> .
nevertheless, sec <dig> and spo <dig> show some partial functional redundancy during sporulation, because only the spo20/sec <dig> double mutants exhibit a complete loss of prospore membranes, while spo <dig> mutants show a milder sporulation phenotype and sec <dig> mutants have no sporulation phenotype  <cit> .
furthermore, the two secretory syntaxins sso <dig> and sso <dig> seem to have some functional differences, because only sso <dig> seems to be necessary during sporulation  <cit> .
degeneration of the tomosyn snare motif in saccharomycotina
in our previous study, we established that tomosyn , a regulatory r-snare without a membrane anchor, is a member of the basic snare repertoire of the proto-eukaryotic cell.
structural studies on the yeast tomosyn homolog, sro <dig>  have demonstrated that this factor carries a large n-terminal domain composed of two consecutive seven-bladed β-propeller domains  <cit> .
however, tomosyn cannot serve as a fusogenic r-snare because it does not have a transmembrane anchor at its c-terminal end  <cit> .
in a molecular process that is not completely understood to date, tomosyn is believed to function in establishing polarized secretion by controlling the accessibility of the t-snares acceptor complexes  <cit> .
as mentioned above, in baker's yeast, two close homologs of tomosyn exist, sro <dig> and sro <dig>   <cit> .
they are often confused with the factor lethal giant larvae , which has a function in establishing epithelial cell polarity in animals  <cit> .
the main basis for this confusion is that sro <dig> and sro <dig>  like lgl, do not possess the canonical r-snare motif found in tomosyns.
in fact, lgl derived independently from tomosyn in animals  <cit> , where it lost its snare domain  <cit> .
interestingly, the canonical r-snare motif of tomosyn is present in other fungal tomosyns  <cit> .
our current analysis demonstrates that the r-snare motif degenerated in the saccharomycetales clade.
two markedly diverged secretory syntaxins within pezizomycotina
our analysis also demonstrates that secretory snares, in particular the secretory syntaxin , have been duplicated recurrently in several fungal lineages, but no clear pattern between duplications and lifestyle or taxonomical grouping has been detected.
usually, duplicated snares are not very greatly diverged in fungi.
however, as gupta and heath have already had pointed out, the sodariomycete neurospora crassa contains two clearly diverged secretory syntaxins, nsyn <dig> and nsyn <dig>  <cit> .
a similar set of diverged secretory syntaxins has been described recently in trichoderma reesei  <cit> , which also belongs to the pezizomycotina.
the pezizomycotina are filamentous fungi, which are sometimes also called euascomycetes.
besides the saccharomycotina, they are the second major subphylum of the ascomycotina, which are defined by the production of a specialized structure, the ascus, that surrounds the meiotic spores.
most pezizomycotina species are hyphal fungi.
hyphae grow at their tips, a process that involves highly polarized secretion of cell wall material.
the two syntaxins of t. reesei, sso <dig> and sso <dig>  were found to be involved in different secretory processes.
whereas sso <dig> was found in the apical compartments of actively growing hyphae, sso <dig> was found in older, non-growing hyphae.
in vivo fret studies also suggest that the two proteins exert their secretory function at spatially segregated areas of the plasma membrane  <cit> .
our analysis now shows that two diverged secretory syntaxins in n. crassa and t. reesei represent a split that arose within the lineage of pezizomycotina .
notably, some species of the eurotiomycetes, a lineage of the pezizomycotina including the genus aspergillus, only possess one sso, suggesting that their lineage split before the duplication of secretory syntaxins  <cit> .
some aspergillus species, a. oryzae and a. flavus, do possess two different secretory syntaxins.
however, this duplication appears to have occurred independently within the eurotiomycetes.
in order to avoid confusion between the well-known duplicated set of secretory syntaxins in s. cerevisiae  we suggest to name these two clearly diverged secretory syntaxins in pezizomycotina ssoa and ssob.
we noted that the duplication of the secretory syntaxin was not accompanied by a duplication of the interacting snare partners of sso, snc and sec <dig>  these usually exist only as one copy in the respective species of pezizomycotina, suggesting that, as shown for t. reesei  <cit> , snc and sec <dig> can interact with two different syntaxins.
interestingly, sec <dig> in a few sodariomycetes, namely magnaporthe grisea, neurospora crassa, podospora anserina and chaetomium globosum, is clearly diverged.
snares recapitulate the phylogeny of fungi
multi-gene phylogenetic studies have started to change our view of the phylogenetic relationships within the fungi kingdom.
in particular, the clades generated by those studies  <cit>  differ markedly from previous classifications of yeasts that were usually defined mostly by morphological analysis and growth responses.
we have already noted that the individual gene trees built from orthologous snares generally recapitulate the newer classification.
in fact, the sequences of species that belong to the same clade usually congregate in the individual trees.
this indicates that snare proteins have diversified only slowly during evolution, but fast enough to reflect species evolution within the fungi kingdom.
for a better resolution of species evolution, we concatenated the snare sequences that are present as singletons.
the resulting phylogenetic tree agrees very well with the recently refined classification of the major fungal lineages .
the tree clearly resolves the major fungi phyla basidiomycotina and ascomycotina.
the latter clearly split into two major subphyla, the pezizomycotina and the saccharomycotina, whereas the two schizosaccharomyces species sit outside the subphyla clades.
our data corroborate the notion that within the pezizomycotina, the eurotiomycetidae split first.
as discussed above, the eurotiomycetidae possess only one secretory syntaxin, sso, whereas we found two clearly diverged sso sequences in all other pezizomycotina, suggesting that this trait emerged after the split of the eurotiomycetidae.
within the saccharomycotina, two major clusters were observed, one being the candida clade, which comprises the organisms that translate the codon ctg as serine instead of leucine  <cit>  and the other being the saccharomyces complex mentioned above.
generally, these two clusters are clearly separated in all individual trees built from orthologous snares, whereas yarrowia lipolytica, as the sole sequenced member of a third cluster of the saccharomycotina, splits earlier.
furthermore, our analysis sheds some initial light on the evolutionary position of the few basal fungi available.
the concatenated tree includes the two mucoromycotina, rhizopus oryzae and phycomyces blakesleeanus, and the chytrid fungus batrachochytrium dendrobatidis, all of which are clearly separated from basidiomycotina.
CONCLUSIONS
our study confirms that fungi, in contrast to animals and green plants  <cit> , generally possess a relatively simple set of snare proteins with only minimal variation between different lineages.
in addition, our analysis reveals that single-cell yeasts, by and large, did not reduce their primordial set of snare proteins, suggesting that they have maintained the basic transport pathways of the proto-eukaryotic cell.
in fact, one can observe quite the contrary: for example, the hemiascomycete s. cerevisiae has gained a few additional snares during evolution.
taken together, these points imply that a basic set of membrane fusion proteins is able to sustain and is necessary for the lifestyle of multicellular fungi as well as of single-cell yeasts.
hence, multicellularity per se does not need an expansion of the repertoire of snare proteins.
moreover, the snares of hemiascomycetes are often more deviated than the ones of basal fungi.
these facts highlight that the trafficking pathways of baker's yeast are not only different to those in animal cells but also are different to those of many other fungi.
since the extent of the differences is unknown to date, these ought to be investigated in more detail by comparing the trafficking pathways of yeast with those of more basal fungi.
it is well established that clear differences can be seen in the appearance of organelles of the endocytic system of yeast and animal cells, most strikingly between the yeast vacuole and the animal lysosomes.
yet, is challenging to compare their endocytic systems as both consist of an interconnected and complex network of organelles.
the current view is that the endosomal system in yeast is simpler than that in an animal cell, but that they also have common principles and pathways  <cit> .
this view is now corroborated by our inspection of the evolution of the snare repertoire in fungi.
originally, the endocytic systems of the progenitors of fungal and animal cells probably operated with the same basic snare collection inherited from the proto-eukaryotic cell.
however, all fungi possess a novel, fungi-specific snare factor, vam <dig>  which can bind to pip-containing membranes by means of its n-terminal px-domain.
although the exact function of this snare protein is not clear at the moment, it seems to play a role in endosomal trafficking, revealing a clear difference in the endosomal snare collection of fungi and animals.
during evolution, the lineage that gave rise to baker's yeast has gained few additional snare proteins.
the largest transformation of the snare collection of baker's yeast was caused by a relatively recent whole genome duplication within the saccharomyces lineage, resulting in an expansion of the secretory set of snares.
although the consequences of this duplication are not well understood, it seems that at least one of the new secretory snares, spo <dig>  has specialized for the process of sporulation.
another gene duplication took place in the endosomal syntaxins.
earlier in evolution, an ancestral endosomal syntaxin  was duplicated, giving rise to pep <dig> and vam <dig> in the lineage of saccharomycotina, while other fungi only possess one endosomal syntaxin.
it is unclear whether one of the duplicated endosomal syntaxins in hemiascomycetous yeasts retained the ancestral function while the second gained a new function  or whether each one carries out part of the ancestral function .
it seems possible that one of the two syntaxins specialized to the process of homotypic vacuole fusion.
this hypothesis could be tested by investigating the function of the ancestral endosomal syntaxin in endosomal trafficking of more basal fungal species like ustilago maydis  <cit> .
the presence of two diverged endosomal syntaxins in hemiascomycetes yeast has also to be taken into consideration when comparing endosomal trafficking in yeast with that in animal cells.
in addition, one has to bear in mind that independent and more substantial expansions of the snare repertoire occurred in animals, possibly establishing more complex endosomal trafficking pathways  <cit> .
for example, during the evolution of animals, an ancestral endosomal syntaxin  gave rise to three different syntaxins, syx <dig>  syx <dig> and syx <dig>  in addition, several other novel endosomal snares were added.
when comparing yeast with animal cells, it is also important not to combine different animal species, as the original animal snare collection was modified in many different ways in different animal lineages.
for instance, vertebrates possess few additional endosomal snares, e.g.
syx <dig>  although the rise of vertebrates mainly resulted in a robust expansion of the secretory set of snares.
abbreviations
au: almost unbiased; er: endoplasmatic reticulum; est: expressed sequence tags; fret: fluorescence resonance energy transfer; ga: golgi apparatus; hmm: hidden markov model; iqpnni: important quartet puzzling and nearest neighbor interchange; jtt: jones: taylor and thornton; lgl: lethal giant larvae; nr: non redundant; npsn: new plant snare; nsf: n-ethylmaleimide-sensitive factor; px: phox homology; snare: soluble nsf attachment receptor; tgn: trans golgi network; wgd: whole genome duplication.
authors' contributions
together the three authors curated the snare data base.
nk performed the sequence and phylogenetic analysis.
thk participated in the phylogenetic analyis and helped to draft the manuscript.
thk and df conceived the study and participated in its design and coordination.
df drafted the manuscript.
all authors read and approved the final manuscript.
supplementary material
additional file 1
classification of the complete set of snare proteins of saccharomyces cerevisiae according to our hmm models.
the different types are listed according to their intracellular distribution and are allocated into interacting snare units based on biological knowledge.
in addition, the uniprotkb/swiss-prot accession number and the initial descriptions are given for each snare protein.
click here for file
additional file 2
list of fungal snares.
list of snare repertoires from different fungi with complete genomes sequenced.
entire sequences were used to calculate phylogenetic trees.
click here for file
additional file 3
alignment and tree files.
detailed information to all sequences and trees are contained in a zipped file, which contains all phylogenetic trees in nexus-format  and all alignments used to calculate the trees in fasta-format .
in addition, a text file is provided  that allows to identify the sequences used for our analysis in the ncbi data base: listed is the sequence number & name used in our alignments & trees and the respective gi string.
click here for file
additional file 4
structure-based sequence alignments.
structure-based sequence alignments of several fungi snares.
click here for file
additional file 5
list of snares used as outgroups.
list of snares from species used as outgroups for calculating phylogenetic trees based on sequence alignments of the snare motif region only.
click here for file
additional file 6
additional sequence peculiarities of fungal snares.
a text describing few additional sequence peculiarities, in particular of layer residues, of fungal snares.
click here for file
acknowledgements
we are greatly indebted to d. h. huson for generously supporting this project.
this research was supported by the max-planck-society .

 over the last decades, the organization of the secretory and endocytic pathways has been studied in a variety of different eukaryotic organisms. these studies have revealed that the principal organization of the endomembrane system and the molecular machineries involved in vesicular trafficking are conserved among all eukaryotes. in general, the transport between different intracellular organelles is mediated by cargo-laden vesicles that bud from a donor and then specifically fuse with an acceptor compartment. key players during the final fusion step are the so-called snare proteins. these proteins are a large family of small cytoplasmically orientated membrane proteins that are typically tail-anchored. their key characteristic is the so-called snare domain, an extended stretch of heptad repeats that is usually connected to a single transmembrane domain by a short linker. snare domains of heterologous sets of snare proteins have the ability to assemble into tight, parallel four-helix bundle complexes. it is thought that complex formation between snare proteins from opposing membranes provide the energy that drives membrane fusion. although the concept of snare-mediated membrane fusion emerged in the 1990s, several snare proteins had been discovered earlier. the field of membrane trafficking was greatly inspired by ingenious genetic screens in the single-cell fungi saccharomyces cerevisiae for mutants with defects in the secretory pathway in the late 1970s. in these mutants, secretion was blocked at higher temperatures. these so-called sec mutants accumulate secretory proteins at the point in the secretory pathway that is blocked and often, they could be distinguished by their phenotype. in the following years, it became clear that sec genes generally encode for key components of the machineries that mediate trafficking of a transport vesicle, including a few snare proteins. all at once, straightforward genetic tools had catapulted baker's yeast to the front row for studying vesicle trafficking pathways. in fact, other genetic screens for mutants showing defects in other transport processes have turned up novel yeast snares. the last yeast snare was tracked down in yielding an overall repertoire of different snare genes in baker's yeast. for most trafficking steps in baker's yeast, distinct units of four interacting snare domains have now been assigned. however, not all details of the trafficking routes in yeast are resolved. generally, different snare family members are localized on distinct organelles and vesicles that demarcate particular membranes. however, it is becoming clear that some trafficking steps cannot be identified by a unique set of snares, since some snares are involved in more than one trafficking step. in addition, they may participate in different snare complexes. to appraise s. cerevisiae as model organism for studying vesicle trafficking pathways better, it is also necessary to understand the limits of the organism, since baker's yeast appears to be a secondarily reduced organism. for example, s. cerevisiae has only a limited ability to produce hyphae. the capability to produce polarized hyphae, however, is a hallmark of the fungal kingdom, allowing them to colonize and exploit new substrates efficiently. in addition, whereas most fungi are multicellular organisms, s. cerevisiae most of the times grows as single cells that reproduce asexually. it thus seems possible that adaptation to a relatively simple lifestyle was accompanied by a degeneration of the intracellular trafficking itinerary together with the involved machineries. a comparison of the snare gene repertoires of s. cerevisiae with that of five other fungal species, including three species of filamentous fungi, revealed that, in general, the members of the snare family are largely conserved in fungi. this notion has been corroborated by subsequent inspections of the snare repertoire of few other fungal species. together, these studies put forward that the yeast lifestyle did not entail a radical change in the intracellular trafficking pathways. the bioinformatic strategies used to identify the homologs of the yeast snares, however, did not provide a universal classification scheme. consequently, it was impossible to compare the snare sets of different fungi with each other in an unambiguous manner, let alone comparing them with snare sets from other eukaryotes like animals. hence it remains unclear whether snare types, which are possibly linked to a novel transport step, evolved or degenerated in particular fungal lineages. fungi are thought to be more closely related to animals than to plants and, currently, are placed with animals and several protistan taxa into the monophyletic group of opisthokonta. notably, morphological observations show that the appearance of analogous intracellular organelles such as the vacuoles/lysosomes and the golgi apparatus are markedly different between yeast and animal cells. it is unclear, however, whether these differences are reflected in their respective snare inventory. recently, using hidden markov models we have established a precise classification for all eukaryotic snare proteins. based on this classification, we were able to deduce that the emergence of multicellularity in animals went along with an expansion of the set of snare proteins, in particular the expansion of the snares involved in endosomal trafficking. in order to analyze the differences between the snare repertoires of fungi and animals and the evolutionary changes in the snare inventory of different fungal lineages, we have collected and classified the snares from a broad range of different fungi according to our hmms, making use of the substantial number of available genomic sequences. to identify specification and duplication events, and to pinpoint the differences between the set of snare proteins in fungi and animals, we used phylogenetic information to reveal the history of snare genes. afterwards, we mapped features of the fungal snare evolution onto the species tree. methods sequences recently, we published a classification for snare sequences, based on more than protein sequences from eukaryotic species, including sequences from species of fungi. we were able to identify twenty basic subtypes within the snare protein family and developed specific as well as sensitive hmms for each subtype. we used these models to identify snare proteins in fungal genomes that became available recently. the models yielded additional sequences from the nr-database at ncbi and various genome projects, baylor college of medicine, j. craig venter institute and broad institute) and several est databases. the newly found protein sequences underwent a visual verification process to remove redundancy and misassembled sequences. overall, our data contained unique snare sequences from fungal species. the dataset was incorporated into the snare database and can be accessed on the project homepage. the snare sets of all fungal species with complete genomes sequenced can be found in. for a more detailed inspection of the evolutionary changes within fungi, we made use of the observation that within a subtype, sequence regions besides the snare motif were also well conserved. consequently, we constructed alignments of the entire sequences for each subgroup using muscle. sites with more than percent gaps were removed from the alignments. in addition, we excluded all sequences which contained more than percent gaps. these alignments are also available via the supplementary section of our projects homepage. based on these alignments, we constructed phylogenetic trees for each of the fungal snares. we also generated sequence alignments in which we graphically incorporated information about the domain structure based on the available high resolution structure information. phylogeny the phylogenetic reconstruction was composed of two different analytical approaches and was basically done as described in. the first approach used iqpnni to construct phylogenetic trees from the curated alignments. we used a gamma distribution as a model for rate heterogeneity with four rate categories for the estimation of the gamma distribution parameter. the proportion of invariable sites was estimated from the data and the jones, taylor and thornton distance matrix served as a substitution matrix. we used the stopping rule of iqpnni, but the calculation had to run for at least the suggested number of iterations. the default values were used for the remaining parameter. in addition, likelihood-mapping was applied to determine the confidence of the edges in the calculated trees. the second approach used the phylip package to apply a distance-based bootstrap analysis with replicates to each of the curated alignments. standard settings were used for seqboot, the jtt distance matrix and also a gamma distribution for protdist, as were standard options for neighbour. if required, the random seed was set to nine. we used the almost unbiased test to address the systematically biased bootstrap values. we obtained the sitewise log-likelihoods needed for the au test using a modified version of phyml and the test was performed using consel. the reconstructed iqpnni trees served as starting points to join the results of both calculations. the inner edges of the trees were labeled with their likelihood-mapping and corrected bootstrap support values. we provide nexus files for each reconstructed tree on our project homepage. the trees can be interactively explored using, for example, splitstree. to gain initial insights on the placement of fungi within the opisthokonta, especially the basal fungal species batrachochytrium dendrobatidis and blastocladiella emersonii, we first performed an extensive phylogenetic analysis based only on the sequence of the snare domain. this analysis was carried out to ensure that the evolutionary history of the individual snare groups indeed reflects the monophyletic origin of fungi. for this analysis, we used all fungal species together with the choanoflagellate monosiga brevicollis, the slime mold dictyostelium discoideum, and ten selected animal species. we also included the snares found in the animal pathogen encephalitozoon cuniculi, which belongs to the microsporidia. in earlier times, this clade was considered to be of protist origin, but it is now believed to have diverged from early fungi. however, we purposefully did not include their sequences into the general phylogenetic analysis, as e. cuniculi appears to have a markedly reduced set of snare proteins with clearly deviating sequences, making them sensitive to long branching artefacts. the reconstructed phylogenetic trees demonstrate that fungi are well separated from other opisthokonts. with this corroborated, we then performed a more detailed phylogenetic analysis of the snares in fungi using the entire sequence.  in order to catalog the snare inventory of fungi, we increased our previous collection with newly published genomes, and complemented this set with ests from different fungal species, obtaining an overall set of more than snare sequences from fungal species. the dataset comprises about species with nearly complete sets of snares. all sequences and species used in this study can also be found in our snare database on our projects homepage. the collected sequences were then classified according to our previously established hmm models that are based only on the highly conserved snare motif. the four main classes were then further subdivided into twenty distinct groups according to our hmm models. each of these subgroups represents a distinct conserved trafficking step. for phylogenetic analysis of the fungal lineage, we aligned entire sequences for each subgroup, as evolutionary changes have usually occurred by paralogous expansion of particular snare types. based on these alignments, we constructed phylogenetic trees for each of the snare types in fungi. fungi contain a compact set of snare proteins overall, the fungal genomes examined in this work comprise a remarkably consistent and homologous set of snare proteins. in each species with a complete genome, we only found a little more than twenty different snares, whereas in metazoans, for example species with more than snares are common. a marked expansion of the metazoan set of snares occurred during the rise of multicellularity. in earlier studies it was put forward that the development of multicellularity is generally linked to an expansion of the snares. in contrast to this idea, we usually only found members of the basic snare proteins to be present in each fungi species. this implies that fungi, although they also have evolved multicellularity, have a relatively unchanged set of snares and have not significantly expanded the repertoire of the proto-eukaryotic cell. it was not the objective of the present study to discuss the function of each snare protein coded for in fungal genomes in detail, as the sequence information alone is insufficient for recognizing the exact trafficking step the protein is involved in. instead, we will mainly focus on highlighting the idiosyncrasies that we came across. notably, the snare proteins involved in secretory trafficking routes between the endoplasmatic reticulum, the golgi apparatus and the trans-golgi network were generally present as singletons and were highly conserved in all fungi. noticeable changes affected the set of snare proteins involved in secretion and in the endosomal/vacuolar trafficking pathway. an interesting observation is that the set of snare proteins in fungi is very similar to the set of the single-cell choanoflagellate monosiga brevicollis, which is thought to be closely related to animals. in contrast to m. brevicollis, lower animals possess an enlarged set of snare proteins. a surprising aspect of the chytrids batrachochytrium dendrobatidis and blastocladiella emersonii is that they possess a snare that was classified as qb.iii.d-type. snares belonging to this subclass were first found in green plants and termed novel plant snares. our previous analysis had shown that npsn can be found in plants and several protist lineages, but not in choanoflagellates and animals. the fact that basal fungi retained a bona fide npsn suggests that late-diverging fungi have lost npsn independently. interestingly, the chytrids are the most primitive of the fungi; many are aquatic and possess a flagellum. consistent with their more basal position, sequences from the chytrids batrachochytrium dendrobatidis and blastocladiella emersonii are usually placed close to the root of the fungal lineage in our analysis. these results suggest that npsn was a member of the original snare repertoire of the eukaryotic cenancestor but was then lost in the choanoflagellate/animal and the fungal lineage. another protein family has been widely used as a marker for the different trafficking steps of the endomembrane system. these are rab proteins, a family of small gtp-binding proteins that regulate different vesicular trafficking steps by cycling between active gtp and inactive gdp forms. it has been observed that fungi contain a relatively small number of rab proteins, whereas their number has expanded in several other lineages such as animals, land plants, and some protists. it would therefore be interesting to investigate whether the pattern of duplications and diversifications found for the rab family in fungi correlates with that of the snare family. the snare protein vam is a molecular "invention" of fungi our analysis shows that fungi contain three different snares of the qc.iii-type: tlg syntaxin and vam. the choanoflagellate m. brevicollis and basal animals contain only two qc.iii-snares, syx and syx which are homologous to tlg and syx respectively. remarkably, we found vam homologs to be present in all complete fungal genomes, even in the basal fungus b. dendrobatidis. however, no vam homolog was discovered in other eukaryotes, suggesting that the qc.iii-snare vam is a new innovation in fungi. of note, this status of vam had been recognized earlier by yoshizawa and colleagues. however, one should keep in mind that their analysis was based on an inferior cluster approach that led, from a biological point of view, to a number of wrongly assigned clusters. in addition, they have missed several well-known snare proteins. although the snare motif of vam is not greatly diverged, it remains unclear whether an ancestral tlg or syx gave rise to vam nonetheless, considering its unique domain structure, vam is the most divergent snare protein in the repertoire of fungi. in contrast to most other snare proteins, vam does not carry a transmembrane region at the c-terminal end of its snare motif. instead, the protein seems to be anchored to a membrane by means of an n-terminal extension that harbors a phox homology domain, which acts as phosphoinositide binding module. in fact, vam has been shown to interact specifically with phosphatidylinositol 3-phosphate p), a lipid that is specific for endosomal and vacuolar membranes. a large body of work carried out in yeast implies that vam catalyzes homotypic vacuole fusion with its partner snares, vam, vti and nyv. interestingly, however, the homologous factor in the basidiomycete ustilago maydis, termed yup was found to be present on vacuoles and rapidly moving endosomal organelles. moreover, a temperature sensitive mutation in yup blocks sorting to vacuoles and a defect in endosomal recycling suggests that the function of vam is not restricted to homotypic vacuole fusion, a process that, as we will discuss below, may be restricted to saccharomycotina. changes in the endosomal and vacuolar set of snare proteins in saccharomycotina the yeast s. cerevisiae belongs to the saccharomycotina clade, the members of which are sometimes referred to as hemiascomycetes. s. cerevisiae possesses two endosomal syntaxins, pep and vam which are believed to be involved in trafficking to late endosomes and to vacuoles, respectively. in agreement with the earlier assessment by gupta and heath, our phylogenetic tree shows that these two snares have been duplicated in the lineage of saccharomycotina only. in addition, our analysis also demonstrates that the duplication of ancestral qa.iii.b-snare into pep and vam is independent of the duplications of the endosomal syntaxins in animals, i.e. syntaxin syntaxin syntaxin and syntaxin the bifurcation of an ancestral endosomal syntaxin within the saccharomycotina clade into proteins with distinct roles, namely pep and vam may suggest that the specialized machinery involved in homotypic vacuole fusion arose only within the hemiascomycetes. homotypic vacuole fusion is a process that occurs during cell division in baker's yeast. during budding, the few large vacuoles maintained in s. cerevisiae are fragmented. multiple small vacuoles then segregate into the budding daughter cell and subsequently fuse again into larger vacuoles. this phenomenon has been exploited by william wickner and colleagues to study the role of snares and other factors in the fusion process. this process, however, does not take place in, for example, s. pombe, a yeast that recently was proposed to belong to a third ascomycetes class, the archiascomycotina. this suggests that homotypic vacuole fusion, described first in baker's yeast, is a process specific to hemiascomycetes. notably, the two qa.iii.b-snares of yarrowia lipolytica are much less diverged. consequently, it seems interesting to find out at which point in the evolution of the saccharomycotina lineage the process of homotypic vacuole fusion evolved. in addition, the sets of proteins catalyzing trafficking to the vacuole and maintaining a low vacuole copy number genes have only partial functional overlap. it would thus be interesting to investigate whether the factors specifically involved in the process of homotypic vacuole fusion are restricted to the saccharomycotina. previously, it has been debated whether the r-snare involved in trafficking towards vacuoles in yeast, nyv is indeed homologous to the animal r-snare vamp. in fact, the snare motif of nyv from baker's yeast is only distantly related to vamp in animals; however, in our analysis, both clearly fall into the same group. it has even been thought that nyv does not contain a canonical profilin-like n-terminal extension, often also referred to as longin domain. a recent structural analysis has demonstrated that the n-terminal extension of nyv indeed contains a profilin-like fold of this type, similar to the ones found in sec and ykt. in general, all fungi possess a nyv homolog. the only fungal lineages that do not appear to contain a nyv homolog are members of schizosaccharomyces. our phylogenetic analysis shows that nyv sequences are remarkably derived in saccharomycetales, bringing to light why it has been difficult before to acknowledge nyv as a bona fide endosomal r-snare. in s. cerevisiae, nyv seems to interact only with the q.iii.b-type syntaxin vam but not with pep which is thought to interact preferentially with the r-snare ykt hence, it seems possible that nyv and vam coevolved in saccharomycotina. the larger set of secretory snares in s. cerevisiae can be attributed to a whole genome duplication in contrast to the notion that baker's yeast is a reduced organism, no fungal genome inspected contained more snare genes than s. cerevisiae with its different snare factors. these genes are clearly more than the 20/ genuine snares of the proto-eukaryotic cell. two of the added snares vam a novel, fungi-specific snare, and vam a saccharomycotina-specific snare were discussed above. in addition, baker's yeast contains four closely related pairs of secretory snare proteins, sso and sso, snc and snc, sec and spo, and sro and sro. our analysis demonstrates that these doubled set of secretory factors very probably arose from a well established whole genome duplication in s. cerevisiae. indeed, we found the same duplicated set of snare proteins in the closely related species s. paradoxus, s. mikatae, s. bayanus and s. kudriavzevii, all of which are collectively referred to as saccharomyces senso-stricto species. remarkably, some but not all of the same secretory snare genes were also found in duplicates in s. castellii, candida glabrata, and vanderwaltozyma polyspora. the ancestor of this group of species and of the saccharomyces sensu stricto species are believed to have undergone the same wgd followed by major secondary reductions, while several closely related species like eremothecium gossypii, kluyveromyces waltii, kluyveromyces lactis, and saccharomyces kluyveri have kept a single copy of the genome. these species are grouped with the ones that underwent a wgd as, "saccharomyces complex". indeed, we found only one copy of each of the secretory snares in the latter species. interestingly, the duplicates found in s. castellii, c. glabrata and v. polyspora usually split at clearly different positions in our phylogenetic trees as compared to the duplicates of the saccharomyces sensu stricto species. this is particularly noticeable for the two sec proteins from s. castellii and also from v. polyspora that are highly similar, whereas the two qbc-snares, sec and spo exhibit marked asymmetrical rates of evolution in the saccharomyces sensu stricto species, pointing to a different fate of the duplicated secretory snares in the two lineages. similarly, the two sso proteins of c. glabrata are clearly more diverged than the duplicated set of the saccharomyces sensu stricto species. in fact, it has been proposed that these lineages must have diverged and adapted to different living conditions very soon after the wgd, so that the genes retained in duplicates diverged largely independently. although only about 14% of all duplicated genes were retained after the wgd, all duplicated snare genes are involved in the process of secretion in baker's yeast. this suggests that these factors presented a selective advantage for yeast. the selective advantage of the duplicated secretory set is not entirely clear, since so far, only subtle functional differences of the duplicated snares have been found. as mentioned above, only one homologous pair, the two qbc-snares is highly diverged in the saccharomyces sensu-stricto species. sec and spo possess only approximately 37% identity, whereas the homologs of the secretory syntaxins and r-snares are 79% and 74% identical, respectively. whereas sec particularly in saccharomycotina, possesses a rather long n-terminal extension, this region was reduced in spo. interestingly, sec and spo are specialized to different developmental stages in baker's yeast. whereas sec interacts with both secretory syntaxins and synaptobrevins in secretion during vegetative growth, the more deviated spo is required only for sporulation. during sporulation in s. cerevisiae, the prospore membrane, which envelops each daughter nucleus during meiosis, is generated de novo adjacent to the spindle pole bodies by the secretory snare machinery. nevertheless, sec and spo show some partial functional redundancy during sporulation, because only the spo20/sec double mutants exhibit a complete loss of prospore membranes, while spo mutants show a milder sporulation phenotype and sec mutants have no sporulation phenotype. furthermore, the two secretory syntaxins sso and sso seem to have some functional differences, because only sso seems to be necessary during sporulation. degeneration of the tomosyn snare motif in saccharomycotina in our previous study, we established that tomosyn, a regulatory r-snare without a membrane anchor, is a member of the basic snare repertoire of the proto-eukaryotic cell. structural studies on the yeast tomosyn homolog, sro have demonstrated that this factor carries a large n-terminal domain composed of two consecutive seven-bladed β-propeller domains. however, tomosyn cannot serve as a fusogenic r-snare because it does not have a transmembrane anchor at its c-terminal end. in a molecular process that is not completely understood to date, tomosyn is believed to function in establishing polarized secretion by controlling the accessibility of the t-snares acceptor complexes. as mentioned above, in baker's yeast, two close homologs of tomosyn exist, sro and sro. they are often confused with the factor lethal giant larvae, which has a function in establishing epithelial cell polarity in animals. the main basis for this confusion is that sro and sro like lgl, do not possess the canonical r-snare motif found in tomosyns. in fact, lgl derived independently from tomosyn in animals, where it lost its snare domain. interestingly, the canonical r-snare motif of tomosyn is present in other fungal tomosyns. our current analysis demonstrates that the r-snare motif degenerated in the saccharomycetales clade. two markedly diverged secretory syntaxins within pezizomycotina our analysis also demonstrates that secretory snares, in particular the secretory syntaxin, have been duplicated recurrently in several fungal lineages, but no clear pattern between duplications and lifestyle or taxonomical grouping has been detected. usually, duplicated snares are not very greatly diverged in fungi. however, as gupta and heath have already had pointed out, the sodariomycete neurospora crassa contains two clearly diverged secretory syntaxins, nsyn and nsyn. a similar set of diverged secretory syntaxins has been described recently in trichoderma reesei, which also belongs to the pezizomycotina. the pezizomycotina are filamentous fungi, which are sometimes also called euascomycetes. besides the saccharomycotina, they are the second major subphylum of the ascomycotina, which are defined by the production of a specialized structure, the ascus, that surrounds the meiotic spores. most pezizomycotina species are hyphal fungi. hyphae grow at their tips, a process that involves highly polarized secretion of cell wall material. the two syntaxins of t. reesei, sso and sso were found to be involved in different secretory processes. whereas sso was found in the apical compartments of actively growing hyphae, sso was found in older, non-growing hyphae. in vivo fret studies also suggest that the two proteins exert their secretory function at spatially segregated areas of the plasma membrane. our analysis now shows that two diverged secretory syntaxins in n. crassa and t. reesei represent a split that arose within the lineage of pezizomycotina. notably, some species of the eurotiomycetes, a lineage of the pezizomycotina including the genus aspergillus, only possess one sso, suggesting that their lineage split before the duplication of secretory syntaxins. some aspergillus species, a. oryzae and a. flavus, do possess two different secretory syntaxins. however, this duplication appears to have occurred independently within the eurotiomycetes. in order to avoid confusion between the well-known duplicated set of secretory syntaxins in s. cerevisiae we suggest to name these two clearly diverged secretory syntaxins in pezizomycotina ssoa and ssob. we noted that the duplication of the secretory syntaxin was not accompanied by a duplication of the interacting snare partners of sso, snc and sec these usually exist only as one copy in the respective species of pezizomycotina, suggesting that, as shown for t. reesei, snc and sec can interact with two different syntaxins. interestingly, sec in a few sodariomycetes, namely magnaporthe grisea, neurospora crassa, podospora anserina and chaetomium globosum, is clearly diverged. snares recapitulate the phylogeny of fungi multi-gene phylogenetic studies have started to change our view of the phylogenetic relationships within the fungi kingdom. in particular, the clades generated by those studies differ markedly from previous classifications of yeasts that were usually defined mostly by morphological analysis and growth responses. we have already noted that the individual gene trees built from orthologous snares generally recapitulate the newer classification. in fact, the sequences of species that belong to the same clade usually congregate in the individual trees. this indicates that snare proteins have diversified only slowly during evolution, but fast enough to reflect species evolution within the fungi kingdom. for a better resolution of species evolution, we concatenated the snare sequences that are present as singletons. the resulting phylogenetic tree agrees very well with the recently refined classification of the major fungal lineages. the tree clearly resolves the major fungi phyla basidiomycotina and ascomycotina. the latter clearly split into two major subphyla, the pezizomycotina and the saccharomycotina, whereas the two schizosaccharomyces species sit outside the subphyla clades. our data corroborate the notion that within the pezizomycotina, the eurotiomycetidae split first. as discussed above, the eurotiomycetidae possess only one secretory syntaxin, sso, whereas we found two clearly diverged sso sequences in all other pezizomycotina, suggesting that this trait emerged after the split of the eurotiomycetidae. within the saccharomycotina, two major clusters were observed, one being the candida clade, which comprises the organisms that translate the codon ctg as serine instead of leucine and the other being the saccharomyces complex mentioned above. generally, these two clusters are clearly separated in all individual trees built from orthologous snares, whereas yarrowia lipolytica, as the sole sequenced member of a third cluster of the saccharomycotina, splits earlier. furthermore, our analysis sheds some initial light on the evolutionary position of the few basal fungi available. the concatenated tree includes the two mucoromycotina, rhizopus oryzae and phycomyces blakesleeanus, and the chytrid fungus batrachochytrium dendrobatidis, all of which are clearly separated from basidiomycotina.  our study confirms that fungi, in contrast to animals and green plants, generally possess a relatively simple set of snare proteins with only minimal variation between different lineages. in addition, our analysis reveals that single-cell yeasts, by and large, did not reduce their primordial set of snare proteins, suggesting that they have maintained the basic transport pathways of the proto-eukaryotic cell. in fact, one can observe quite the contrary: for example, the hemiascomycete s. cerevisiae has gained a few additional snares during evolution. taken together, these points imply that a basic set of membrane fusion proteins is able to sustain and is necessary for the lifestyle of multicellular fungi as well as of single-cell yeasts. hence, multicellularity per se does not need an expansion of the repertoire of snare proteins. moreover, the snares of hemiascomycetes are often more deviated than the ones of basal fungi. these facts highlight that the trafficking pathways of baker's yeast are not only different to those in animal cells but also are different to those of many other fungi. since the extent of the differences is unknown to date, these ought to be investigated in more detail by comparing the trafficking pathways of yeast with those of more basal fungi. it is well established that clear differences can be seen in the appearance of organelles of the endocytic system of yeast and animal cells, most strikingly between the yeast vacuole and the animal lysosomes. yet, is challenging to compare their endocytic systems as both consist of an interconnected and complex network of organelles. the current view is that the endosomal system in yeast is simpler than that in an animal cell, but that they also have common principles and pathways. this view is now corroborated by our inspection of the evolution of the snare repertoire in fungi. originally, the endocytic systems of the progenitors of fungal and animal cells probably operated with the same basic snare collection inherited from the proto-eukaryotic cell. however, all fungi possess a novel, fungi-specific snare factor, vam which can bind to pip-containing membranes by means of its n-terminal px-domain. although the exact function of this snare protein is not clear at the moment, it seems to play a role in endosomal trafficking, revealing a clear difference in the endosomal snare collection of fungi and animals. during evolution, the lineage that gave rise to baker's yeast has gained few additional snare proteins. the largest transformation of the snare collection of baker's yeast was caused by a relatively recent whole genome duplication within the saccharomyces lineage, resulting in an expansion of the secretory set of snares. although the consequences of this duplication are not well understood, it seems that at least one of the new secretory snares, spo has specialized for the process of sporulation. another gene duplication took place in the endosomal syntaxins. earlier in evolution, an ancestral endosomal syntaxin was duplicated, giving rise to pep and vam in the lineage of saccharomycotina, while other fungi only possess one endosomal syntaxin. it is unclear whether one of the duplicated endosomal syntaxins in hemiascomycetous yeasts retained the ancestral function while the second gained a new function or whether each one carries out part of the ancestral function. it seems possible that one of the two syntaxins specialized to the process of homotypic vacuole fusion. this hypothesis could be tested by investigating the function of the ancestral endosomal syntaxin in endosomal trafficking of more basal fungal species like ustilago maydis. the presence of two diverged endosomal syntaxins in hemiascomycetes yeast has also to be taken into consideration when comparing endosomal trafficking in yeast with that in animal cells. in addition, one has to bear in mind that independent and more substantial expansions of the snare repertoire occurred in animals, possibly establishing more complex endosomal trafficking pathways. for example, during the evolution of animals, an ancestral endosomal syntaxin gave rise to three different syntaxins, syx syx and syx in addition, several other novel endosomal snares were added. when comparing yeast with animal cells, it is also important not to combine different animal species, as the original animal snare collection was modified in many different ways in different animal lineages. for instance, vertebrates possess few additional endosomal snares, e.g. syx although the rise of vertebrates mainly resulted in a robust expansion of the secretory set of snares. abbreviations au: almost unbiased; er: endoplasmatic reticulum; est: expressed sequence tags; fret: fluorescence resonance energy transfer; ga: golgi apparatus; hmm: hidden markov model; iqpnni: important quartet puzzling and nearest neighbor interchange; jtt: jones: taylor and thornton; lgl: lethal giant larvae; nr: non redundant; npsn: new plant snare; nsf: n-ethylmaleimide-sensitive factor; px: phox homology; snare: soluble nsf attachment receptor; tgn: trans golgi network; wgd: whole genome duplication. authors' contributions together the three authors curated the snare data base. nk performed the sequence and phylogenetic analysis. thk participated in the phylogenetic analyis and helped to draft the manuscript. thk and df conceived the study and participated in its design and coordination. df drafted the manuscript. all authors read and approved the final manuscript. supplementary material additional file 1 classification of the complete set of snare proteins of saccharomyces cerevisiae according to our hmm models. the different types are listed according to their intracellular distribution and are allocated into interacting snare units based on biological knowledge. in addition, the uniprotkb/swiss-prot accession number and the initial descriptions are given for each snare protein. click here for file additional file 2 list of fungal snares. list of snare repertoires from different fungi with complete genomes sequenced. entire sequences were used to calculate phylogenetic trees. click here for file additional file 3 alignment and tree files. detailed information to all sequences and trees are contained in a zipped file, which contains all phylogenetic trees in nexus-format and all alignments used to calculate the trees in fasta-format. in addition, a text file is provided that allows to identify the sequences used for our analysis in the ncbi data base: listed is the sequence number & name used in our alignments & trees and the respective gi string. click here for file additional file 4 structure-based sequence alignments. structure-based sequence alignments of several fungi snares. click here for file additional file 5 list of snares used as outgroups. list of snares from species used as outgroups for calculating phylogenetic trees based on sequence alignments of the snare motif region only. click here for file additional file 6 additional sequence peculiarities of fungal snares. a text describing few additional sequence peculiarities, in particular of layer residues, of fungal snares. click here for file acknowledgements we are greatly indebted to d. h. huson for generously supporting this project. this research was supported by the max-planck-society.
although the snare repertoire of baker's yeast is highly conserved, our analysis reveals that it is more deviated than the ones of basal fungi. moreover, since different snare family members are thought to demarcate different organelles and vesicles, this approach allowed us to compare the organization of the endomembrane systems of yeast and animal cells. over the past years, the yeast saccharomyces cerevisiae has served as a powerful model organism for studying the organization of the secretory and endocytic pathways, and a few years ago, its entire set of snares was compiled. here, we make use of the increasing amount of genomic data to investigate the history of the snare family during fungi evolution. however, all fungi contain a novel soluble snare protein, vam which carries an n-terminal px-domain that acts as a phosphoinositide binding module. in eukaryotic cells, directional transport between different compartments of the endomembrane system is mediated by vesicles that bud from a donor organelle and then fuse with an acceptor organelle. for instance, the endosomal syntaxins pep and vam arose from a gene duplication that occurred within the saccharomycotina clade. this highlights that the trafficking pathways of baker's yeast are not only different to those in animal cells but also are somewhat different to those of many other fungi. our data corroborate the notion that fungi generally encompass a relatively simple set of snare proteins, mostly comprising the snares of the proto-eukaryotic cell.
 in eukaryotic cells, directional transport between different compartments of the endomembrane system is mediated by vesicles that bud from a donor organelle and then fuse with an acceptor organelle. a family of integral membrane proteins, termed soluble n-ethylmaleimide-sensitive factor attachment receptor proteins, constitute the key machineries of these different membrane fusion events. over the past years, the yeast saccharomyces cerevisiae has served as a powerful model organism for studying the organization of the secretory and endocytic pathways, and a few years ago, its entire set of snares was compiled.  here, we make use of the increasing amount of genomic data to investigate the history of the snare family during fungi evolution. moreover, since different snare family members are thought to demarcate different organelles and vesicles, this approach allowed us to compare the organization of the endomembrane systems of yeast and animal cells. our data corroborate the notion that fungi generally encompass a relatively simple set of snare proteins, mostly comprising the snares of the proto-eukaryotic cell. however, all fungi contain a novel soluble snare protein, vam which carries an n-terminal px-domain that acts as a phosphoinositide binding module. in addition, the points in fungal evolution, at which lineage-specific duplications and diversifications occurred, could be determined. for instance, the endosomal syntaxins pep and vam arose from a gene duplication that occurred within the saccharomycotina clade.  although the snare repertoire of baker's yeast is highly conserved, our analysis reveals that it is more deviated than the ones of basal fungi. this highlights that the trafficking pathways of baker's yeast are not only different to those in animal cells but also are somewhat different to those of many other fungi. 