BACKGROUND
transmembrane transporters of unicellular organisms, like yeast, are one of the primary links between the outer world and the metabolic pathways inside the cell.
the importance of these genes is seen in the substantial proportion of transporter genes within the yeast genome   <cit> .
in saccharomyces cerevisiae, for example, over  <dig> genes encoding transporter proteins have been found  <cit> .
different species of yeast require different substrates to be transported, and the number and kind of transporters present in the genome therefore vary between species.
the presence or absence of a transporter may provide insight into the niche preference and metabolic ability of a yeast  <cit> .
an example is seen in the loss or inactivation of the galactose transporter gal2p  in several hemiascomycete species  <cit> .
similarly, expansion or contraction of gene numbers in various transporter subfamilies can indicate metabolic abilities: for example, peroxisomal and long chain fatty acid transporters have undergone amplification in y. lipolytica, which is known to grow on fatty acids  <cit> .
more unexpectedly, genes belonging to two heavy metal transporter subfamilies, sit and ct <dig>  have been amplified  <dig> and  <dig> times respectively in y. lipolytica  <cit> , possibly as the result of a natural symbiosis with bacteria providing iron-siderophores and/or other natural chelators.
yeast transporters have been classified into a number of families and subfamilies, based on functional and phylogenetic criteria  <cit> .
within the hemiascomycete phylum,  <dig> small phylogenetic transporter subfamilies have been identified comprising a total of  <dig> transporters named according to their evolutionary patterns  <cit> .
the subfamily which is most variable in gene number across yeast species is the anion:cation symporter subfamily  of the major facilitator superfamily  <cit> .
these transporters are all putative weak acid permeases, which take up anionic vitamins  in symport with protons or cations.
in saccharomyces cerevisiae the allantoate permease family consists of eight proteins where the allantoate permease yjr152w  has given its name to the subfamily.
the other members of this family are ylr004c , yll055w , ygr260w , yil166c, yal067c , ygr065c , ycr028c   <cit>  .
to investigate the evolution of this gene family, and track its expansions and contractions across hemiascomycetes, we have performed a phylogenetic analysis of these genes across a number of hemiascomycete species for which whole genomes are available.
the dal <dig> subfamily is a ubiquitous transporter found in all hemiascomycete species, but with a very complex evolution involving repeated gene losses and duplications.
RESULTS
the relationship between the members of the dal <dig> subfamily
the allantoate permease transporters are known to be an ancient subfamily  <cit> , and as expected we found almost all members of the dal <dig> group represented in the hemiascomycete, euascomycete and basidiomycete taxa included in our analysis.
we identified three major clusters in the phylogenetic tree of all known members of the dal <dig> subfamily: one including the three genes seo <dig>  vht <dig> and fen <dig> and their orthologs, another including yct <dig>  dal <dig> and th <dig>  and the third including yil166c, yol163-2w and tna <dig>  for more details .
the relationships within and between these clades are consistent with those in earlier papers  with the addition of yol163-2w, which is a pseudogene in s. cerevisiae, in the last group.
the two genes seo <dig> and vht <dig> appear to be the result of an ancient duplication in an ancestor of the saccharomycotina group.
members from the pezizomycotina group cluster outside of the saccharomycotina seo <dig> and vht <dig> clade.
for most genes within this subfamily, the phylogenetic relationships we observed were similar to known taxonomic relationships .
with euascomycete, basidiomycete and hemiascomycete species each clustering together, and basidiomycetes at the base of each tree.
within the hemiascomycetes, we expect to find y. lipolytica at the base of the clade, with kluyveromyces and candida species diverging subsequently, followed by the saccharomyces sensu stricto species  <cit> .
we observed this phylogenetic pattern overall, but with some exceptions.
the relationship of the dal <dig> genes
the dal <dig> gene tree provides one exception to the expected phylogenetic relationships.
s. pombe has two copies of this gene: one is basal to the hemiascomycetes, as expected, but the second clusters with k. lactis .
this clustering of a s. pombe gene with k. lactis has very low bootstrap support, however, and is probably artifactual.
c. guilliermondi and d. hansenii both possess multiple paralogs of this gene, and each has one representative that clustered outside the hemiascomycetes.
these genes may be ancient paralogs, which may have evolved functions distinct from those described for dal <dig>
the relationship of the tna <dig> genes
another unexpected relationship was found in the gene tree of tna <dig>  in which four species k. waltii, c. tropicalis, y. lipolytica and s. pombe had several gene duplicates, both recent and ancient, and spread unevenly in the gene tree.
in our analyses we found five copies of this gene from c. tropicalis and s. pombe clustering at three different places in the tree and three copies from k. waltii in two clusters.
for y. lipolytica we found  <dig> gene copies clustering at five different places .
a distant group of species, group c, include only species from the pre-whole genome duplication event in the yeast clade and could constitute another gene member with a different function or specialization from tna <dig>  the paralog was obviously lost in an ancestor of the whole genome duplication lineage.
this might also be the case in the groups of a and b where y. lipolytica has several recent duplications.
the gene tree indicates ancient duplication of this gene, since many species have gene copies in several clusters  with reciprocal losses in the other clusters which result in a non exact match of the taxonomic phylogeny.
the relationship of the yil166c and yol163-2w genes
the evolution of the group of genes which we originally believed to be orthologs of yil166c is more complex than we anticipated.
these genes appear to be members of two paralogous clades, one consisting of yil166c from s. cerevisiae and closely related genes in the sensu stricto species, c. guilliermondi and d. bruxellensis.
gene members from u. maydis and a. oryzae were also found in this, the smallest of the clades of the genes in the dal <dig> subfamily .
the second paralogous clade consists of genes from only hemiascomycete species, including the contiguous s. cerevisiae genes yol163w and yol162w, constituting a putative single frame-shifted pseudogene.
s. cerevisiae and s. bayanus were the only representatives from the saccharomyces sensu stricto group that we found.
more than half the members of this clade have gene duplications.
besides y. lipolytica, c. guillerimondi, c. albicans, c. tropicalis, a. gossipii and d. hansenii each have two or more paralogous copies, indicating multiple recent gene duplications after species divergence.
these two clades, yil166c and yol163-2w, have a common base represented by two euascomycete species, a. nidulans and a. fumigatus, and three gene copies from c. guillermondi .
there are two further points of interest in this region of the tree .
the first is the presence of a u. maydis  gene at the base of the hemiascomycete yil166c clade.
although the node has only 60% bootstrap support, the node below this is 100% supported, which results in a confident clade in which a basidiomycete is nested within ascomycetes.
this relationship suggests that there has been a horizontal gene transfer event from either euascomycetes or hemiascomycetes to this basidiomycete species.
it will be interesting to see if a similar sequence can be detected in any other basidiomycete species.
secondly, while most species are represented in most gene clades, the species represented in each of these two clades  are substantially different.
the yil166c group contains the sensu stricto species, d. bruxellensis and c. guilliermondi, while the other clade contains the other species from the hemiascomycetes.
the saccharomyces sensu stricto species, c. guilliermondi and d. bruxellensis have all been found in alcoholic beverages like wine  <cit> .
the presence of these orthologous transporters in these taxa could be due to selection for retention of a transporter specially adapted to an environment of low oxygen and high ethanol concentration.
it seems that the genes from these species are true orthologs with a common origin, and have either been lost from the other hemiascomycetes or horizontally transferred into the "wine species".
the gene transfer is either from a species close to u. maydis, which branched off before the euascomycetes, or the gene in u. maydis is also a horizontally transferred gene, and the horizontal gene transfer occurred from a species in the euascomycetes clade to u. maydis, and the "wine species".
these hypotheses are currently speculative, and much further work would be required to investigate them further.
another unusual behaviour of the hemiascomycete orthologs of yil166c is that this transporter appears to have undergone species-specific duplication in every species in our analysis.
some have been duplicated prior to speciation like c. tropicalis and c. albicans.
the intraspecific duplication has resulted in two or three copies for most of the species in the trees except for y. lipolytica.
the pattern of gene duplication in this species is unusual: for each member of the subfamily, y. lipolytica has either three or more paralogs or no copy of the gene at all.
we found no orthologs of yil166c, yct <dig>  th <dig> and vht <dig> in y. lipolytica, but for fen <dig> we found three copies, yol163w four copies, dal <dig> five copies, seo <dig> six copies and for tna <dig>  <dig> copies, that met our criteria.
the selective advantage of these amplifications for y. lipolytica is not understood, but it could be a compensation for the loss of the other dal <dig> transporters or a selective advantage coupled to a sudden change in the environment and/or the hydrocarbon diet.
the amplification we observe of anion:cation symporters in y. lipolytica has also been seen in drug:h+antiporter transporters and quinate:h+symporters  <cit> .
how have the gene expansions taken place?
duplications of genes could be on the gene level, segmental, chromosomal or whole genome duplications.
the duplications offer the opportunity for copies to evolve different functions, either broader or more specialized.
gene copies that provide an advantage for the organism may be preserved and go to fixation under selection.
others might lose any function, become pseudogenes and ultimately be lost from the genome.
in the dal <dig> family the number of copies various tremendously within and between species for the different gene family members.
this expansion and shrinkage is most evident in y. lipolytica where we found three to twenty paralogs for five of the nine subfamilies and in the others, no copy of the gene at all.
the other species has also expanded into two to four copies for most genes in this subfamily.
in some cases it is evident that the duplications occurred before speciation, for example see c. albicans and c. tropicalis for yol163w, while in other cases there have been recent independent duplications in several species, see y. lipolytica and a. gossypii in yol163w .
a similar pattern can be observed in the lower part of the tree dal <dig>  where a few lineages  seem to have undergone a recent duplication and at the same time an ancestor of s. cerevisiae and the sensu stricto group deleted their copy.
most species-specific genes are located near telomeres in saccharomyces sensu stricto species  <cit> .
four  of the eight genes  in this subfamily are subtelomeric in s. cerevisiae.
it is possible that the complexity that we see in many trees, especially for y. lipolytica and d. hansenii, is a result of horizontal gene transfer coupled with intraspecific duplications.
maybe the duplications that we see are a result of specialization for different substrates or maybe even for a different function.
most of the permeases are complex in nature, with many different substrates, and for many of the species investigated here the transporter proteins are poorly studied.
CONCLUSIONS
yeast is one of the most important model organisms in biology research including evolution.
the yeast lineage, as part of the kingdom of fungi, has revealed both horizontal transfers and whole genome duplication and is a very important group regarding new insights into evolutionary pathways  <cit> .
crucial genes often show an interesting evolution and transmembrane proteins are often essential for the survival of the yeast cells.
the transmembrane dal <dig> subfamily includes genes that seem to be very important due to its many duplications in almost all yeast species investigated, but it also includes genes that obviously can be deleted without any harm done to the organism.
when we analysed species from the whole kingdom of fungi we found that the dal <dig> subfamily has a very complex evolution in yeast, not seen before, with intra- and interspecific duplications and unusual relationships indicating specialization, deletions and maybe even horizontal gene transfer.
we believe that the dal <dig> subfamily will be important in future investigations of evolution in fungi, especially the evolution of transmembrane proteins and their specialization.

 transmembrane transporters of unicellular organisms, like yeast, are one of the primary links between the outer world and the metabolic pathways inside the cell. the importance of these genes is seen in the substantial proportion of transporter genes within the yeast genome. in saccharomyces cerevisiae, for example, over genes encoding transporter proteins have been found. different species of yeast require different substrates to be transported, and the number and kind of transporters present in the genome therefore vary between species. the presence or absence of a transporter may provide insight into the niche preference and metabolic ability of a yeast. an example is seen in the loss or inactivation of the galactose transporter gal2p in several hemiascomycete species. similarly, expansion or contraction of gene numbers in various transporter subfamilies can indicate metabolic abilities: for example, peroxisomal and long chain fatty acid transporters have undergone amplification in y. lipolytica, which is known to grow on fatty acids. more unexpectedly, genes belonging to two heavy metal transporter subfamilies, sit and ct have been amplified and times respectively in y. lipolytica, possibly as the result of a natural symbiosis with bacteria providing iron-siderophores and/or other natural chelators. yeast transporters have been classified into a number of families and subfamilies, based on functional and phylogenetic criteria. within the hemiascomycete phylum, small phylogenetic transporter subfamilies have been identified comprising a total of transporters named according to their evolutionary patterns. the subfamily which is most variable in gene number across yeast species is the anion:cation symporter subfamily of the major facilitator superfamily. these transporters are all putative weak acid permeases, which take up anionic vitamins in symport with protons or cations. in saccharomyces cerevisiae the allantoate permease family consists of eight proteins where the allantoate permease yjr152w has given its name to the subfamily. the other members of this family are ylr004c, yll055w, ygr260w, yil166c, yal067c, ygr065c, ycr028c. to investigate the evolution of this gene family, and track its expansions and contractions across hemiascomycetes, we have performed a phylogenetic analysis of these genes across a number of hemiascomycete species for which whole genomes are available. the dal subfamily is a ubiquitous transporter found in all hemiascomycete species, but with a very complex evolution involving repeated gene losses and duplications.  the relationship between the members of the dal subfamily the allantoate permease transporters are known to be an ancient subfamily, and as expected we found almost all members of the dal group represented in the hemiascomycete, euascomycete and basidiomycete taxa included in our analysis. we identified three major clusters in the phylogenetic tree of all known members of the dal subfamily: one including the three genes seo vht and fen and their orthologs, another including yct dal and th and the third including yil166c, yol163-2w and tna for more details. the relationships within and between these clades are consistent with those in earlier papers with the addition of yol163-2w, which is a pseudogene in s. cerevisiae, in the last group. the two genes seo and vht appear to be the result of an ancient duplication in an ancestor of the saccharomycotina group. members from the pezizomycotina group cluster outside of the saccharomycotina seo and vht clade. for most genes within this subfamily, the phylogenetic relationships we observed were similar to known taxonomic relationships. with euascomycete, basidiomycete and hemiascomycete species each clustering together, and basidiomycetes at the base of each tree. within the hemiascomycetes, we expect to find y. lipolytica at the base of the clade, with kluyveromyces and candida species diverging subsequently, followed by the saccharomyces sensu stricto species. we observed this phylogenetic pattern overall, but with some exceptions. the relationship of the dal genes the dal gene tree provides one exception to the expected phylogenetic relationships. s. pombe has two copies of this gene: one is basal to the hemiascomycetes, as expected, but the second clusters with k. lactis. this clustering of a s. pombe gene with k. lactis has very low bootstrap support, however, and is probably artifactual. c. guilliermondi and d. hansenii both possess multiple paralogs of this gene, and each has one representative that clustered outside the hemiascomycetes. these genes may be ancient paralogs, which may have evolved functions distinct from those described for dal the relationship of the tna genes another unexpected relationship was found in the gene tree of tna in which four species k. waltii, c. tropicalis, y. lipolytica and s. pombe had several gene duplicates, both recent and ancient, and spread unevenly in the gene tree. in our analyses we found five copies of this gene from c. tropicalis and s. pombe clustering at three different places in the tree and three copies from k. waltii in two clusters. for y. lipolytica we found gene copies clustering at five different places. a distant group of species, group c, include only species from the pre-whole genome duplication event in the yeast clade and could constitute another gene member with a different function or specialization from tna the paralog was obviously lost in an ancestor of the whole genome duplication lineage. this might also be the case in the groups of a and b where y. lipolytica has several recent duplications. the gene tree indicates ancient duplication of this gene, since many species have gene copies in several clusters with reciprocal losses in the other clusters which result in a non exact match of the taxonomic phylogeny. the relationship of the yil166c and yol163-2w genes the evolution of the group of genes which we originally believed to be orthologs of yil166c is more complex than we anticipated. these genes appear to be members of two paralogous clades, one consisting of yil166c from s. cerevisiae and closely related genes in the sensu stricto species, c. guilliermondi and d. bruxellensis. gene members from u. maydis and a. oryzae were also found in this, the smallest of the clades of the genes in the dal subfamily. the second paralogous clade consists of genes from only hemiascomycete species, including the contiguous s. cerevisiae genes yol163w and yol162w, constituting a putative single frame-shifted pseudogene. s. cerevisiae and s. bayanus were the only representatives from the saccharomyces sensu stricto group that we found. more than half the members of this clade have gene duplications. besides y. lipolytica, c. guillerimondi, c. albicans, c. tropicalis, a. gossipii and d. hansenii each have two or more paralogous copies, indicating multiple recent gene duplications after species divergence. these two clades, yil166c and yol163-2w, have a common base represented by two euascomycete species, a. nidulans and a. fumigatus, and three gene copies from c. guillermondi. there are two further points of interest in this region of the tree. the first is the presence of a u. maydis gene at the base of the hemiascomycete yil166c clade. although the node has only 60% bootstrap support, the node below this is 100% supported, which results in a confident clade in which a basidiomycete is nested within ascomycetes. this relationship suggests that there has been a horizontal gene transfer event from either euascomycetes or hemiascomycetes to this basidiomycete species. it will be interesting to see if a similar sequence can be detected in any other basidiomycete species. secondly, while most species are represented in most gene clades, the species represented in each of these two clades are substantially different. the yil166c group contains the sensu stricto species, d. bruxellensis and c. guilliermondi, while the other clade contains the other species from the hemiascomycetes. the saccharomyces sensu stricto species, c. guilliermondi and d. bruxellensis have all been found in alcoholic beverages like wine. the presence of these orthologous transporters in these taxa could be due to selection for retention of a transporter specially adapted to an environment of low oxygen and high ethanol concentration. it seems that the genes from these species are true orthologs with a common origin, and have either been lost from the other hemiascomycetes or horizontally transferred into the "wine species". the gene transfer is either from a species close to u. maydis, which branched off before the euascomycetes, or the gene in u. maydis is also a horizontally transferred gene, and the horizontal gene transfer occurred from a species in the euascomycetes clade to u. maydis, and the "wine species". these hypotheses are currently speculative, and much further work would be required to investigate them further. another unusual behaviour of the hemiascomycete orthologs of yil166c is that this transporter appears to have undergone species-specific duplication in every species in our analysis. some have been duplicated prior to speciation like c. tropicalis and c. albicans. the intraspecific duplication has resulted in two or three copies for most of the species in the trees except for y. lipolytica. the pattern of gene duplication in this species is unusual: for each member of the subfamily, y. lipolytica has either three or more paralogs or no copy of the gene at all. we found no orthologs of yil166c, yct th and vht in y. lipolytica, but for fen we found three copies, yol163w four copies, dal five copies, seo six copies and for tna copies, that met our criteria. the selective advantage of these amplifications for y. lipolytica is not understood, but it could be a compensation for the loss of the other dal transporters or a selective advantage coupled to a sudden change in the environment and/or the hydrocarbon diet. the amplification we observe of anion:cation symporters in y. lipolytica has also been seen in drug:h+antiporter transporters and quinate:h+symporters. how have the gene expansions taken place? duplications of genes could be on the gene level, segmental, chromosomal or whole genome duplications. the duplications offer the opportunity for copies to evolve different functions, either broader or more specialized. gene copies that provide an advantage for the organism may be preserved and go to fixation under selection. others might lose any function, become pseudogenes and ultimately be lost from the genome. in the dal family the number of copies various tremendously within and between species for the different gene family members. this expansion and shrinkage is most evident in y. lipolytica where we found three to twenty paralogs for five of the nine subfamilies and in the others, no copy of the gene at all. the other species has also expanded into two to four copies for most genes in this subfamily. in some cases it is evident that the duplications occurred before speciation, for example see c. albicans and c. tropicalis for yol163w, while in other cases there have been recent independent duplications in several species, see y. lipolytica and a. gossypii in yol163w. a similar pattern can be observed in the lower part of the tree dal where a few lineages seem to have undergone a recent duplication and at the same time an ancestor of s. cerevisiae and the sensu stricto group deleted their copy. most species-specific genes are located near telomeres in saccharomyces sensu stricto species. four of the eight genes in this subfamily are subtelomeric in s. cerevisiae. it is possible that the complexity that we see in many trees, especially for y. lipolytica and d. hansenii, is a result of horizontal gene transfer coupled with intraspecific duplications. maybe the duplications that we see are a result of specialization for different substrates or maybe even for a different function. most of the permeases are complex in nature, with many different substrates, and for many of the species investigated here the transporter proteins are poorly studied.  yeast is one of the most important model organisms in biology research including evolution. the yeast lineage, as part of the kingdom of fungi, has revealed both horizontal transfers and whole genome duplication and is a very important group regarding new insights into evolutionary pathways. crucial genes often show an interesting evolution and transmembrane proteins are often essential for the survival of the yeast cells. the transmembrane dal subfamily includes genes that seem to be very important due to its many duplications in almost all yeast species investigated, but it also includes genes that obviously can be deleted without any harm done to the organism. when we analysed species from the whole kingdom of fungi we found that the dal subfamily has a very complex evolution in yeast, not seen before, with intra- and interspecific duplications and unusual relationships indicating specialization, deletions and maybe even horizontal gene transfer. we believe that the dal subfamily will be important in future investigations of evolution in fungi, especially the evolution of transmembrane proteins and their specialization. 
we also observed a close relationship between the gene yil166c and its homologs in the saccharomyces sensu stricto species and two "wine spoiler" yeasts, dekkera bruxellensis and candida guilliermondi, which could possibly be the result of horizontal gene transfer between these distantly related species. the transmembrane dal subfamily was found to have a very complex evolution in yeast with intra- and interspecific duplications and unusual relationships indicating specialization, specific deletions and maybe even horizontal gene transfer. all are putative anion:cation symporters of vitamins. we believe that this group will be important in future investigations of evolution in fungi and especially the evolution of transmembrane proteins and their specialization. we studied the complex evolution of this group in species from the kingdom of fungi particularly the ascomycetes. the dal subfamily is an old and important group since it is represented in both basidiomycetes and ascomycetes. the dal subfamily of transmembrane transporter genes has eight known members in s. cerevisiae. in the analyses we detect several well defined groups without s. cerevisiae representation suggesting new gene members in this subfamily with perhaps altered specialization or function.
 genes continuously duplicate and the duplicated copies remain in the genome or get deleted. the dal subfamily of transmembrane transporter genes has eight known members in s. cerevisiae. all are putative anion:cation symporters of vitamins. the dal subfamily is an old and important group since it is represented in both basidiomycetes and ascomycetes. we studied the complex evolution of this group in species from the kingdom of fungi particularly the ascomycetes.  we identified numerous gene duplications creating sets of orthologous and paralogous genes. in different lineages the dal subfamily members expanded or contracted and in some lineages a specific member could not be found at all. we also observed a close relationship between the gene yil166c and its homologs in the saccharomyces sensu stricto species and two "wine spoiler" yeasts, dekkera bruxellensis and candida guilliermondi, which could possibly be the result of horizontal gene transfer between these distantly related species. in the analyses we detect several well defined groups without s. cerevisiae representation suggesting new gene members in this subfamily with perhaps altered specialization or function.  the transmembrane dal subfamily was found to have a very complex evolution in yeast with intra- and interspecific duplications and unusual relationships indicating specialization, specific deletions and maybe even horizontal gene transfer. we believe that this group will be important in future investigations of evolution in fungi and especially the evolution of transmembrane proteins and their specialization. 