BACKGROUND
since their original discovery  <cit> , pseudogenes have been found in all studied genomes so far, within the three kingdoms of life  <cit> .
their proportions vary greatly, however, from one organism to another, depending on lifestyle  and on genome properties .
pseudogenes correspond to ca.
3% of the gene repertoire of drosophila melanogaster, while there are approximately as many pseudogenes as functional genes in the human genome .
pseudogenes are often species-specific and, within small genomes, tend to accumulate in chromosomal regions such as subtelomeres or heterochromatin, minimizing possible deleterious effects  <cit> .
the historical definition of a pseudogene is a dna sequence that looks like an active gene but has lost its ability to code for a functional product, due to more or less extensive mutational disablements.
pseudogene formation is frequently observed in pathogenic organisms undergoing reductive evolution while benefiting from host functions  <cit> .
pseudogenes can also correspond to the non-functionalization of a duplicated gene copy  <cit>  that originated either from dna duplication or from retro-transposition  <cit> .
sequence degradation may range from a single disabling mutation, such as a frameshift or an in-frame stop codon in a protein-coding gene, to extensive changes including numerous insertions and deletions  <cit> .
in all genomes, most pseudogenes are likely to disappear with time by the accumulation of successive mutations  <cit> .
those with limited alterations, however, may be repaired by reverse mutations, gene conversion or may be reactivated by translational recoding events  <cit> .
some pseudogenes can also acquire a new functional role, such as the control of another gene expression, or the generation of genetic diversity , and give rise to new genes, such as the human xist non-coding rna gene, which evolved from a pseudogene of a protein-coding sequence  <cit> .
in fungi, only pseudogenes of particular genes or pathways of interest have been described  <cit> .
in saccharomyces cerevisiae, two systematic analyses performed by different approaches  <cit>  both concluded to the paucity of pseudogenes of anciently protein-coding sequences in the yeast genome.
in order to determine the mechanisms of pseudogene formation across several related species, their age and the functions affected by pseudogenization, we performed the first systematic comparative analysis of the pseudogene repertoire in a set of eight hemiascomycetous yeasts spanning a large evolutionary range, similar or larger than the phylum of chordates  <cit> .
we show that these genomes also contain a limited number of pseudogenes, independently of their global level of gene redundancy.
most pseudogenes originate from duplicated gene copies resulting from previous dna duplication events, but a few could correspond to retro-processed sequences.
some pseudogenes are formed from single-copy gene and hence correspond to a functional loss in the yeast species.
the relative paucity of pseudogenes suggests that functionally inactive genes are rapidly eliminated in hemiascomycetous genomes because their sequence diverges rapidly or because they are successively truncated, or entirely deleted.
methods
genomes and protein sequences
sequences and annotations were taken from the génolevures database  for candida glabrata, zygosaccharomyces rouxii, kluyveromyces lactis, kluyveromyces thermotolerans, saccharomyces kluyveri, debaryomyces hansenii and yarrowia lipolytica, and from the saccharomyces genome database  for saccharomyces cerevisiae.
all these sequenced strains are haploids.
note that the annotation files already contained  <dig> protein-coding pseudogenes.
the protein databank used for comparison contains the translational products of the  <dig>  annotated cds from all eight studied genomes .
we used the protein family classification constructed by the génolevures consortium  <cit> .
whole genome shotgun sequences of s. cerevisiae haploid strains awri <dig>  jay <dig>  m <dig>  rm11-1a, yjm <dig>  yps <dig> and ec <dig> were retrieved at embl , as well as the genomic scaffolds of the strain ec <dig> .
identification of the pseudogenes
we developed a set of automatic procedures to obtain an exhaustive list of the potential pseudogenes of anciently protein-coding genes in a given genome.
each genome sequence is given to fasty  <cit>  and translated into the  <dig> frames for detection of amino acid similarity against the translational products of the  <dig>  annotated cds from the studied genomes .
the matches selected have at least 25% identity and are considered statistically significant if their estimated z-score  <cit>  is greater or equal to  <dig> .
this cut-off was chosen after examining the distribution of the z-scores for all obtained results.
when several matches overlap, the similarity region is delimited from the left-most to the right-most aligned sequences.
protein domain conservation
all matches similar to a protein of unknown function were queried against the pfam database  <cit>  of conserved domains with the hmmer algorithm <cit> .
identification of the pseudogenes in s. cerevisiae strains
the regions corresponding to the pseudogenes in the s. cerevisiae reference strain s288c were retrieved in the whole genome shotgun sequences of the strains awri <dig>  ec <dig>  jay <dig>  m <dig>  rm11-1a, yjm <dig> and yps <dig> by a fasta search at the nucleotide level.
these identified regions were then queried against the protein bestmatches of the corresponding pseudogenes in s288c with fasty.
we concluded to the presence of a pseudogene in the considered strain if evidence for coding sequence degradation is proposed in the fasty alignment .
the pseudogenes present in strains m <dig> and yps <dig> but not in s288c were identified based on the work by  <cit> .
selection criteria
pseudogenes are defined from the similarity regions using the following criteria: 1) no overlap or partial overlap with an already annotated functional genetic element on the same dna strand; 2) possible overlap of less than  <dig> nucleotides with an already annotated functional genetic element on the opposite dna strand; 3) evidence for coding sequence degradation: at least one in-frame stop codon, one frameshift mutation, or a truncation of at least 30% relatively to the cognate cds .
each selected potential pseudogene is aligned to its bestmatch with genewise  <cit>  to predict the corresponding coding sequence.
fifteen candidates were discarded after this step because genewise predicted an intact cds while fasty introduced a frameshift mutation in the alignment.
the pseudogenes identified in this study are available at the génolevures website .
the algorithm for pseudogene detection was written for our own analysis, considering that annotations of the genomes studied are of high quality.
it is available upon request to the authors.
pseudogene nomenclature
each detected pseudogene is numbered serially from left to right of the chromosome, based on the nomenclature proposed in  <cit> .
the name indicates the species , the project or strain number , the chromosome , and the pseudogene nature of the sequence  followed by the serial number .
sequence divergence analysis
we assumed that substitutions at the third codon position of a protein coding sequence  evolve according to a neutral molecular clock.
estimations of the evolutionary distances by several substitution models  are saturated for half of our data set and, therefore, were not retained for analysis.
p-distances are then computed for amino acid and for nucleotide sequences on the third codon position.
to conserve the frame of the coding sequences, the nucleotide alignment is derived from the amino acid alignment  of the translated products of the pseudogene and its cognate gene obtained by genewise.
the nucleotide alignments between functional paralogs are derived from the amino acid alignments obtained by muscle  <cit> .
synteny conservation
pairwise synteny was examined in windows of five cds upstream and five cds downstream of the pseudogene and any of its homolog in another genome.
we confirmed conservation of synteny if at least two pairs of homologous neighbors are found within the window.
detection of processed pseudogenes
to determine if a pseudogene could arose from a retro-transposition event, we systematically checked for classical retro-transposition hallmarks: lack of intron with respect to their paralog cds, retrotransposon-related sequences within the first flanking annotated elements of the pseudogene , and polya-tail at the 3'-end of the pseudogene.
for this last check, we analyzed the a content of the 3'-end flanking region  of each pseudogene.
we considered as a potential poly-tail, a window of  <dig> nucleotides containing at least  <dig> adenines, with at least one stretch of  <dig> adenines, and less than  <dig> thymines.
to avoid fortuitous signals due to at-rich sequences, the same procedure was performed on  <dig> random sequences .
only sequences for which no signal was detected in the random sequences were considered as potential poly-tail.
statistical analysis
correlation tests were done with the method of spearman implemented in r  <cit> .
we accepted a correlation between data sets when the p-value was lower than  <dig> .
RESULTS
strategy
we systematically searched for all potential pseudogenes of ancient protein-coding genes within  <dig> completely sequenced and annotated genomes of hemiascomycetes: saccharomyces cerevisiae and candida glabrata , zygosaccharomyces rouxii, kluyveromyces lactis, kluyveromyces thermotolerans, saccharomyces kluyveri , and debaryomyces hansenii and yarrowia lipolytica, members of the «ctg» group and of the dipodascaceae, respectively .
to do this, we compared each genome sequence to the set of  <dig>  protein-coding sequences  annotated in these genomes .
we considered here as pseudogene any sequence that simultaneously -i- does not overlap an already annotated genetic element, -ii- shares sequence similarity with an annotated cds  and -iii- shows disabling mutations in the reading frame: in-frame stop codon, frameshift mutation or truncation of more than 30% relatively to the cds.
we chose this limit of 30% because, among the functional members of a given protein coding gene family, the length variation does not exceed 30% in the majority of cases .
inactive pseudogenes resulting solely from mutations in promoters, as well as pseudogenes without detectable similarity among the eight studied genomes are, therefore, excluded from our analysis.
similarly, we do not question the genome annotations: an annotated gene that appears truncated compared to other homologous sequences will not be considered as a pseudogene because truncation does not necessarily imply inactivation .
complete set of detected pseudogenes and sequence degradation patterns
we found a total of  <dig> pseudogenes among the  <dig> yeast genomes .
they correspond to  <dig> distinct protein families or singletons.
the total number of pseudogenes varies significantly between species: from  <dig> in c. glabrata to  <dig> in y. lipolytica, corresponding to only  <dig> % and  <dig> % of the annotated cds, respectively .
a. in megabases, except rdna.
b. pseudogenes of protein coding sequences only.
the percentage of pseudogenes relative to cds is indicated in parenthesis.
c. proportion of « full-size » pseudogenes, i.e.
extending over more than 70% of their bestmatch length.
d. proportion of pseudogenes extensively truncated at their 3'-end.
e. proportion of pseudogenes extensively truncated at their 5'-end.
the pseudogenes identified exhibit a large panel of sequence degradation, ranging from only a few mutational disablements to extensive truncations.
we distinguished here the «full-size» pseudogenes, i.e.
those extending over more than 70% of their closest functional homolog, from the «truncated» pseudogenes .
in s. cerevisiae, s. kluyveri, k. thermotolerans and k. lactis, the number of «full-size» pseudogenes exceeds truncated ones.
the opposite is true for c. glabrata, z. rouxii, d. hansenii and y. lipolytica.
there are as many pseudogenes truncated from their 3'-end, as those truncated from their 5'-end , consistent with the idea that truncation results from deletion at the dna level rather than incomplete cdna formation, characteristic of retro-processed pseudogenes.
«full-size» pseudogenes contain more disabling mutations  than truncated pseudogenes  .
we found four pseudogenes  inactivated by insertion of a dna fragment of mitochondrial origin , as previously described  <cit> .
sequence divergence between pseudogenes and their best functional homologs
in the absence of the actual ancestral coding sequences of the pseudogenes, we estimated the distances between pseudogenes and their closest functional homologs  to measure their degree of mutational decay .
distributions of the p-distances were calculated for each yeast species .
they range from  <dig> to  <dig>  globally.
the distributions are different if one considers pseudogenes whose bestmatch is not in the same genome and pseudogenes whose bestmatch is a paralog in the same genome.
in the first case, distances are high, consistent with the generally long evolutionary distances between studied species.
in the second case, distances are much lower than that observed between functional paralogs in the same genome, suggesting that many pseudogenes arose after recent gene duplication events.
the genomes of s. kluyveri and k. thermotolerans contain the highest proportions of highly diverged pseudogenes , while the genomes of s. cerevisiae and y. lipolytica contain the highest proportions of less diverged pseudogenes .
note that functional paralogs  in all these genomes are issued from both ancient  and recent  duplications  <cit> .
pseudogenes across species, conservation or independent formation
a pseudogene should be free from any functional pressure and, therefore, should not be conserved over long evolutionary periods, unless it acquires a functional role.
we, therefore, searched for pseudogenes conserved in several yeast species as a possible functional signature.
a total of  <dig> pseudogenes are shared between at least two species but  <dig> of them are located outside of any conserved synteny block and probably correspond to independent pseudogenization events of homologous genes in different lineages.
the remaining  <dig> pseudogenes are found among five conserved synteny blocks, defining five groups of orthologous pseudogenes .
note that we did not consider the syntenic orthologs of yil009c-a , ypl052w  and yor239w  as pseudogenes, because they all bear the programmed translational frameshift present in the functional genes of s. cerevisiae  <cit> .
according to the species phylogeny and following parsimony criterion, three of these sets of orthologous pseudogenes also correspond to independent pseudogenization events of the same gene in distinct lineages .
pseudogenes sakl0gp <dig> and klth0gp <dig>  corresponding to ydr285w  of s. cerevisiae are conserved between s. kluyveri and k. thermotolerans.
zip <dig> is a transverse filament protein of the synaptonemal complex that is required for normal levels of meiotic recombination and pairing between homologous chromosomes during meiosis.
the  <dig> orthologous pseudogenes are «full-size» and share a single common frameshift mutation .
the most parsimonious hypothesis is that the frameshift mutation appeared in the common ancestor of these  <dig> species.
rec <dig> is a meiotic gene in s. cerevisiae, which has been lost several times in hemiascomycetes  <cit> .
this gene is also pseudogenized in three of the newly sequenced genomes studied.
sakl0gp <dig> and klth0bp <dig> are conserved between s. kluyveri and k. thermotolerans, suggesting a pseudogenization event in their common ancestor.
the orthologous gene is also pseudogenized in z. rouxii , and probably corresponds to another independent event.
unlike in the previous case, the pseudogenes of this block of conserved synteny are truncated and each copy accumulated several species-specific disabling mutations.
conservation of pseudogenes across s. cerevisiae strains
the genome sequences for the s. cerevisiae non-reference strains are draft assemblies, not suitable for the pseudogene detection procedure developed in this work.
however, to get a first estimation of pseudogene conservation among strains, we examined the conservation of the pseudogenes identified in the s. cerevisiae reference strain s288c  among  <dig> other sequenced strains .
we found that  <dig> of the pseudogenes identified in s288c are also pseudogenized in all other strains, suggesting that the pseudogenization event occurred in their common ancestor.
the  <dig> other pseudogenes of s288c sometimes correspond to intact coding sequences in one or a few other strains, indicating more recent pseudogenization events .
<dig> of these  <dig> pseudogenes contain at least two degrading mutations  and  <dig> of them are truncated at their 3' or 5' end.
the only pseudogene specific to s288c  contains only one internal stop codon.
p: homologous pseudogene, intact: homolog with non-degraded coding sequence.
distribution of the pseudogenes among protein families and corresponding functions
based on sequence similarity, pseudogenes can be attributed to the functional gene families classified according to their predicted translational products  <cit> .
such families contain groups of orthologs between species, as well as groups of paralogs resulting from gene duplications and losses .
in each yeast species, the majority of pseudogenes  belong to gene families with a functionally characterized s. cerevisiae member, which can be used to infer the probable function of the gene that was pseudogenized.
interestingly, we observed a bias towards transporters, proteins acting at the periphery of the cell, and enzymes .
when normalized against the ratios of these functional categories among active genes, the bias is conserved for transporters and proteins acting at the periphery of the cell, but not for enzymes.
among pseudogenes with no ascribed function, there is no over-representation of domains related to these three categories.
although the number of pseudogenes is not directly correlated to the number of its active paralogs in a given family, we observed frequent formation of pseudogenes in these functional categories.
for example, there are <dig> pseudogenes in d. hansenii for the sugar transporter family,  <dig> pseudogenes in c. glabrata for the lectin-like protein family, and  <dig> pseudogenes in z. rouxii for the dup <dig> gene family coding for membrane proteins .
s.c. is for s. cerevisiae.
anumber of pseudogenes with a homolog in s. cerevisiae coding for proteins involved in transport and/or acting at the periphery of the cell.
bnumber of pseudogenes without homolog in s. cerevisiae coding for enzymes.
a total of  <dig> pseudogenes are similar to singleton genes, and  <dig>  belong to gene families specific to a given species for which, unless in s. cerevisiae, there is usually no functional indication .
clustering of the pseudogenes in subtelomeres
gene densities are nearly constant along yeast chromosomes, except for subtelomeres where the number of active genes is reduced  <cit> , and where genes can be transcriptionally silenced  <cit> .
we, therefore, examined the distribution of the pseudogenes along the chromosomes, separating the subtelomeres  from the central regions .
in all yeast species, pseudogenes exceed active genes in number  in the subtelomeres.
in s. cerevisiae, c. glabrata and z. rouxii, more than half of the pseudogenes are found in subtelomeres .
this proportion is lower in k. thermotolerans, d. hansenii, s. kluyveri and k. lactis , and drops to only 5% in y. lipolytica.
the presence of rdna loci in the subtelomeres of y. lipolytica  might prevent some pseudogenes accumulating in these regions, as well as in d. hansenii  and in c. glabrata   <cit> .
a. percentage of pseudogenes in subtelomeric regions .
b. percentage of active genes in subtelomeric regions.
c. number of pseudogenes without annotated functional paralog in the genome.
d. number of pseudogenes with annotated functional paralog whose closest homolog is in the same genome  or in another genome .
*.
species for which all chromosomes are fully sequenced, including their telomeric repeats.
possible origin of the pseudogenes
we tried to define the origin of yeast pseudogenes based on the presence or absence of a paralog in the same genome and on the conservation of synteny between species .
the presence of an active paralog in the genome reveals the occurrence of a previous duplication event, followed by mutational inactivation of one of the  <dig> copies.
the absence of any paralogs reveals the degradation of a single copy gene, hence the loss of the corresponding function in the species.
the conservation of synteny is a signature for mutational sequence degradation at the origin of the pseudogene.
the non-conservation of synteny is compatible with a species-specific duplication, either a retro-transposition event or a segmental duplication.
pseudogenes originating after species-specific gene duplications
more than half of the pseudogenes arose after a gene duplication event specific to their own genome .
this corresponds to the vast majority of pseudogenes in s. cerevisiae, d. hansenii and y. lipolytica , ca.
70% of the pseudogenes in c. glabrata and z. rouxii, and less than 60% in k. lactis, k. thermotolerans and s. kluyveri.
the majority of these pseudogenes probably arose after segmental duplication, but some may be processed pseudogenes, as suggested by a poly tract at their 3'-end or their location next to a retrotransposon-related sequence  <cit> .
this would correspond to only  <dig> % of the pseudogenes in y. lipolytica, and 20% of the pseudogenes in k. thermotolerans.
only one pseudogene  was identified by the lack of intron compared to its closest functional homolog.
a. number of pseudogenes originating from mutational inactivation of a duplicated gene copy formed after speciation.
b. number of pseudogenes originating from mutational inactivation of a duplicated gene copy formed before speciation.
c. number of pseudogenes originating from mutational inactivation of a single copy gene.
d. number of duplicated pseudogenes among the first category .
e. number of pseudogenes being part of a duplicated segment involving other adjacent genes among the first category .
f. number of retro-processed pseudogenes, among the first category , identified by: either the presence of a 3' poly-tail  or the proximity of retrotransposon-related sequence .
in each species, the candidates identified by these  <dig> criteria are different.
interestingly, we found duplications of pseudogenes, as evidenced by their common pattern of degrading mutations.
for example, in y. lipolytica, we found five pseudogenes similar to yali0a14927g , which, based on phylogenetic analysis, suggest two original pseudogenization events in the history of this family.
note that these pseudogenes are not located in subtelomeric regions and are not part of larger duplicated regions including other genes.
these pseudogenes have thus not been maintained in the genome by some selective pressure on the duplication of adjacent genes.
similarly, two such cases of pseudogene multiplication were encountered in k. lactis, as well as  <dig> cases of pseudogene duplication .
most pseudogenes located at the chromosome ends arose after a species-specific gene duplication : from 60% in k. thermotolerans to 100% in y. lipolytica.
see table  <dig> for legend.
pseudogenes originating from ancestral gene duplications
pseudogenes formed after ancestral gene duplications are found in all yeast genomes .
four pseudogenes, sace0bp <dig>  cagl0dp <dig>  caglohp <dig> and cagl0mp <dig> were formed after the whole-genome duplication event that occurred in the common ancestor of c. glabrata and s. cerevisiae  <cit>  .
this is a very small number compared to the extensive loss of duplicated genes that occurred by deletion  <cit> .
strikingly, there are more pseudogenes originated from other ancestral duplications in protoploid saccharomycetaceae and in d. hansenii and y. lipolytica, suggesting that selective pressure on duplicated genes are different after whole-genome duplication and other duplication events, such as segmental duplications.
pseudogenes originating without previous gene duplication and putative loss of function
pseudogenes with no functional paralog  are found in all the studied genomes and could, therefore, correspond to a functional loss.
the genomes of k. lactis, k. thermotolerans, s. kluyveri and c. glabrata contain the highest proportions of such pseudogenes .
the lowest proportions are found in d. hansenii and in y. lipolytica , consistent with the higher gene redundancy in these genomes  <cit> .
there is no bias towards any functional category .
pseudogenes homologous to essential genes in s. cerevisiae are found in all species, except in k. lactis.
repeated pseudogenization of the same gene is also encountered in this category.
for example, pseudogenes cagl0mp <dig> and zyro0gp <dig>  similar to the essential chromosomal passenger gene yjr089w , have no homolog in their respective genome.
finally, we noticed three interesting pseudogenes in s. cerevisiae that correspond to horizontally acquired genes.
sace0fp <dig> and sace0fp <dig> have no homolog in s. cerevisiae but were identified by similarity to deha2d01122g  and deha2e07282g , respectively.
absence of these genes among hemiascomycetes suggests that four independent horizontal gene transfer events have occurred, the two genes in s. cerevisiae being secondarily pseudogenized.
the third one, sace0np <dig> is a duplicated and identical copy of sace0fp <dig>  revealing the expansion of horizontally transferred inactivated genes.
the transfers must have occurred in an ancestor of all the studied strains of s. cerevisiae because the homologs of sace0fp <dig>  sace0fp <dig> and sace0np <dig> are pseudogenized in the  <dig> other strains examined.
discussion
we present here the first systematic search for pseudogenes in the genomes of eight distinct yeast species spanning the whole evolutionary spectrum of hemiascomycetes.
because our method of detection relies on sequence similarity with annotated protein coding genes, the total number of pseudogenes identified  represents a minimal estimate .
among the genomes analyzed here, only s. cerevisiae has been the subject of previous systematic analysis of the pseudogene content  <cit> .
our analysis revealed a smaller number of pseudogenes in this species  because we deliberately ignored partial overlaps with annotated genetic elements.
on the contrary, additional pseudogenes were identified based on our multi-species sequence comparisons.
these differences illustrate the difficulties in listing all pseudogenes in a given genome.
despite this fact and numerical variations between yeast species, our results demonstrate that the proportion of pseudogenes compared to active genes remains low in all hemiascomycetous genomes studied .
our comparative analysis of pseudogenes across an entire yeast phylum provides a unique data set to examine their origin and evolutionary conservation.
we found that pseudogenes in yeasts are formed either by disabling mutations  or by extensive truncations.
the general absence of conservation of pseudogenes between yeast species is consistent with their large evolutionary distances  <cit>  and indicates that new pseudogenes were formed within each lineage.
however, intra-species conservation of pseudogenes is high: about 80% of the pseudogenes in the s. cerevisiae reference strain are old enough to be conserved among  <dig> other strains of this species.
poor conservation of pseudogenes is also observed between mammalian species  <cit> .
most of the pseudogenes correspond to duplicated gene copies, illustrating the extensive dynamics of gene duplications in the yeast genomes, most probably segmental duplications, as observed in s. cerevisiae  <cit> .
only few pseudogenes correspond to duplicated gene copies formed by the ancestral whole-genome duplication common to c. glabrata and s. cerevisiae  <cit> .
this is consistent with the idea that most duplicated copies were lost by complete deletion, as previously proposed  <cit> , or with the possibility that pseudogenes have been degraded beyond recognition, given the time elapsed since the whole-genome duplication event.
among the protoploid saccharomycetaceae, the presence of species-specific pseudogenes within conserved synteny blocks indicates their relatively recent formation, to the notable exception of two ancestral pseudogenes conserved in s. kluyveri and k. thermotolerans, the more closely related pairs of species studied.
the general absence of conserved pseudogenes confirms that most of the observed pseudogenes are on their way to complete degradation.
according to doniger et al.,  <cit>   <dig> pseudogenes containing an internal in-frame stop codon  and  <dig> pseudogenes containing a frameshift mutation are present in either m <dig> or yps <dig>  while they correspond to intact coding sequences in s288c.
from the alignment of the three genomes within coding regions, it appears that the loss of a gene by deletion is much more frequent ) than the creation of a pseudogene by a substitution event  or by an indel event ).
according to lynch et al.
<cit> , the rate of gene loss,  <dig>  × 10- <dig> per gene per cell division, is much higher than the base-substitutional rate,  <dig>  × 10- <dig> per site per cell division.
it should then be much more frequent to lose a gene by deletion than to create a pseudogene.
unlike their abundance in mammalian genomes  <cit> , only few pseudogenes may have originated from retro-transposition events in yeasts.
but their presence suggests that, as experimentally demonstrated in s. cerevisiae  <cit> , retro-transposition occurred in these genomes, although a very small number of active retrotransposons are usually present  <cit> .
the number of retro-processed pseudogenes in yeasts is, however, probably underestimated because their identification is difficult given the small number of intron-containing genes.
the few detected pseudogenes homologous to intron-containing genes are all truncated and do not span the intron insertion-site, except in one case .
pseudogenes without functional paralog suggest functional loss in the corresponding species, unless a non-homologous gene encodes a similar function  <cit> .
in s. cerevisiae, the fact that about half of the pseudogenes correspond to unknown functions is striking given that there remain only 17% of the genes not yet functionally characterized in this species  <cit> .
this suggests that the divergence of functional repertoire between yeast species primarily concerns functions not yet identified.
the frequent occurrence of pseudogenes corresponding to transporters or protein acting at the periphery of the cell may be correlated with the tendency of such genes to cluster in the subtelomeric regions, which are highly dynamic in terms of gene duplications and losses  <cit>  and often concern functions involved in the adaptation of the species to its environment  <cit> .
the highly dynamical behavior of such families is supported by the fact that  <dig> to  <dig> times more pseudogenes than genes are found in the subtelomeric regions.
the vast majority of pseudogenes in the subtelomeres originated from a species-specific duplication.
this bias could be correlated to the accelerated base-pair substitution observed in the subtelomeres, which probably also contributes to adaptive evolution  <cit> .
it also suggests that pseudogenes issued from ancestral duplication and function loss are not preferentially maintained in subtelomeres with respect to central regions of chromosomes.
in s. cerevisiae, pseudogenes of highly connected genes are significantly under-represented, possibly reflecting the lower propensity of gene loss among these genes: 6% of the pseudogenes of genes coding for protein in complexes against 28% of the active genes  <cit>  and 59% of pseudogenes of genes with a genetic interaction profile against 75% of the active genes  <cit> .
whether pseudogenes correspond to intermediate gene states before complete erasure or steady states conferring selective advantages remains an open question and may be case-specific.
our data provide no indication as to whether pseudogenes may be transcribed or not.
however, according to available experimental data in s. cerevisiae ,  <dig> of the  <dig> detected pseudogenes appear to be transcribed .
evidence of transcription also exists in c. glabrata and y. lipolytica .
the transcription products of pseudogenes may be directly targeted to the degradation machinery such as nmd  <cit> , but they could also play some role in the cell by interfering with the expression of functional genes.
moreover, as previously suggested  <cit> , pseudogenes with few disabling mutations could constitute a reservoir of functional protein products if recoding events occur, such as programmed frameshift , or bypass of the stop-codon  <cit> .
the number of disabling mutations within each pseudogene is given in columns  <dig> and  <dig>  the relative length  of the pseudogene with respect to its closest functional homolog is indicated in column  <dig>  last column indicates the reference of the data set where evidence of transcription is found .
for a given chromosome, all identified pseudogenes are separated by several genes.
there is no bias of any kind among these pseudogenes.
CONCLUSIONS
pseudogenes are found in all yeast genomes, albeit in limited number compared to genomes of multi-cellular eukaryotes.
they mostly result from lineage-specific mutational degradations that may correlate with species adaptation to their environment.
yeast pseudogenes show a wide range of mutational alterations, consistent with their rapid evolution, hence their absence of conservation between species.
along with complete gene deletion, pseudogene formation contributes to the rapid genome evolution by gene duplication and loss in yeasts.
the paucity of observed pseudogenes across the entire phylum of hemiascomycetes suggests that pseudogene formation is not the main mechanism of gene loss within these genomes.
this could be explained by the low estimated rate for pseudogene formation across s. cerevisiae strains, compared to the estimated rate of gene deletion.
however, despite their unlikely occurrence, pseudogenes do exist in the yeast genomes.
they appear mainly by species-specific duplications and testify for the adaptation of the cell to its environment.
their poor conservation across species suggests that most of them are likely to disappear.
authors' contributions
il designed and conducted the analysis, implemented the algorithms and wrote the manuscript.
bd contributed to the interpretation of the results and to the writing of the manuscript.
both authors read and approved the final manuscript.
supplementary material
additional file 1
table s <dig> and table s <dig>  - table s <dig>  list and characteristics of the detected pseudogenes in the  <dig> studied yeast genomes.
- table s <dig>  list of gene families  with detected pseudogenes.
click here for file
additional file 2
table s <dig> and descriptions of tables s <dig> and s <dig>  descriptions of tables s <dig> and s <dig>  - table s <dig>  lost functions among the eight studied species.
a list of the pseudogenes with no paralog in the different genomes
click here for file
acknowledgements
we thank gilles fischer, bertrand llorente and eduardo rocha for helpful comments on this work and our colleagues from the "unité de génétique moléculaire des levures" and from the génolevures consortium for fruitful discussions.
this work was supported by grant anr-05-blan- <dig> from the agence nationale de la recherche .
bd is a member of the institut universitaire de france.

 since their original discovery, pseudogenes have been found in all studied genomes so far, within the three kingdoms of life. their proportions vary greatly, however, from one organism to another, depending on lifestyle and on genome properties. pseudogenes correspond to ca. 3% of the gene repertoire of drosophila melanogaster, while there are approximately as many pseudogenes as functional genes in the human genome. pseudogenes are often species-specific and, within small genomes, tend to accumulate in chromosomal regions such as subtelomeres or heterochromatin, minimizing possible deleterious effects. the historical definition of a pseudogene is a dna sequence that looks like an active gene but has lost its ability to code for a functional product, due to more or less extensive mutational disablements. pseudogene formation is frequently observed in pathogenic organisms undergoing reductive evolution while benefiting from host functions. pseudogenes can also correspond to the non-functionalization of a duplicated gene copy that originated either from dna duplication or from retro-transposition. sequence degradation may range from a single disabling mutation, such as a frameshift or an in-frame stop codon in a protein-coding gene, to extensive changes including numerous insertions and deletions. in all genomes, most pseudogenes are likely to disappear with time by the accumulation of successive mutations. those with limited alterations, however, may be repaired by reverse mutations, gene conversion or may be reactivated by translational recoding events. some pseudogenes can also acquire a new functional role, such as the control of another gene expression, or the generation of genetic diversity, and give rise to new genes, such as the human xist non-coding rna gene, which evolved from a pseudogene of a protein-coding sequence. in fungi, only pseudogenes of particular genes or pathways of interest have been described. in saccharomyces cerevisiae, two systematic analyses performed by different approaches both concluded to the paucity of pseudogenes of anciently protein-coding sequences in the yeast genome. in order to determine the mechanisms of pseudogene formation across several related species, their age and the functions affected by pseudogenization, we performed the first systematic comparative analysis of the pseudogene repertoire in a set of eight hemiascomycetous yeasts spanning a large evolutionary range, similar or larger than the phylum of chordates. we show that these genomes also contain a limited number of pseudogenes, independently of their global level of gene redundancy. most pseudogenes originate from duplicated gene copies resulting from previous dna duplication events, but a few could correspond to retro-processed sequences. some pseudogenes are formed from single-copy gene and hence correspond to a functional loss in the yeast species. the relative paucity of pseudogenes suggests that functionally inactive genes are rapidly eliminated in hemiascomycetous genomes because their sequence diverges rapidly or because they are successively truncated, or entirely deleted. methods genomes and protein sequences sequences and annotations were taken from the génolevures database for candida glabrata, zygosaccharomyces rouxii, kluyveromyces lactis, kluyveromyces thermotolerans, saccharomyces kluyveri, debaryomyces hansenii and yarrowia lipolytica, and from the saccharomyces genome database for saccharomyces cerevisiae. all these sequenced strains are haploids. note that the annotation files already contained protein-coding pseudogenes. the protein databank used for comparison contains the translational products of the annotated cds from all eight studied genomes. we used the protein family classification constructed by the génolevures consortium. whole genome shotgun sequences of s. cerevisiae haploid strains awri jay m rm11-1a, yjm yps and ec were retrieved at embl, as well as the genomic scaffolds of the strain ec. identification of the pseudogenes we developed a set of automatic procedures to obtain an exhaustive list of the potential pseudogenes of anciently protein-coding genes in a given genome. each genome sequence is given to fasty and translated into the frames for detection of amino acid similarity against the translational products of the annotated cds from the studied genomes. the matches selected have at least 25% identity and are considered statistically significant if their estimated z-score is greater or equal to. this cut-off was chosen after examining the distribution of the z-scores for all obtained results. when several matches overlap, the similarity region is delimited from the left-most to the right-most aligned sequences. protein domain conservation all matches similar to a protein of unknown function were queried against the pfam database of conserved domains with the hmmer algorithm. identification of the pseudogenes in s. cerevisiae strains the regions corresponding to the pseudogenes in the s. cerevisiae reference strain s288c were retrieved in the whole genome shotgun sequences of the strains awri ec jay m rm11-1a, yjm and yps by a fasta search at the nucleotide level. these identified regions were then queried against the protein bestmatches of the corresponding pseudogenes in s288c with fasty. we concluded to the presence of a pseudogene in the considered strain if evidence for coding sequence degradation is proposed in the fasty alignment. the pseudogenes present in strains m and yps but not in s288c were identified based on the work by. selection criteria pseudogenes are defined from the similarity regions using the following criteria: 1) no overlap or partial overlap with an already annotated functional genetic element on the same dna strand; 2) possible overlap of less than nucleotides with an already annotated functional genetic element on the opposite dna strand; 3) evidence for coding sequence degradation: at least one in-frame stop codon, one frameshift mutation, or a truncation of at least 30% relatively to the cognate cds. each selected potential pseudogene is aligned to its bestmatch with genewise to predict the corresponding coding sequence. fifteen candidates were discarded after this step because genewise predicted an intact cds while fasty introduced a frameshift mutation in the alignment. the pseudogenes identified in this study are available at the génolevures website. the algorithm for pseudogene detection was written for our own analysis, considering that annotations of the genomes studied are of high quality. it is available upon request to the authors. pseudogene nomenclature each detected pseudogene is numbered serially from left to right of the chromosome, based on the nomenclature proposed in. the name indicates the species, the project or strain number, the chromosome, and the pseudogene nature of the sequence followed by the serial number. sequence divergence analysis we assumed that substitutions at the third codon position of a protein coding sequence evolve according to a neutral molecular clock. estimations of the evolutionary distances by several substitution models are saturated for half of our data set and, therefore, were not retained for analysis. p-distances are then computed for amino acid and for nucleotide sequences on the third codon position. to conserve the frame of the coding sequences, the nucleotide alignment is derived from the amino acid alignment of the translated products of the pseudogene and its cognate gene obtained by genewise. the nucleotide alignments between functional paralogs are derived from the amino acid alignments obtained by muscle. synteny conservation pairwise synteny was examined in windows of five cds upstream and five cds downstream of the pseudogene and any of its homolog in another genome. we confirmed conservation of synteny if at least two pairs of homologous neighbors are found within the window. detection of processed pseudogenes to determine if a pseudogene could arose from a retro-transposition event, we systematically checked for classical retro-transposition hallmarks: lack of intron with respect to their paralog cds, retrotransposon-related sequences within the first flanking annotated elements of the pseudogene, and polya-tail at the 3'-end of the pseudogene. for this last check, we analyzed the a content of the 3'-end flanking region of each pseudogene. we considered as a potential poly-tail, a window of nucleotides containing at least adenines, with at least one stretch of adenines, and less than thymines. to avoid fortuitous signals due to at-rich sequences, the same procedure was performed on random sequences. only sequences for which no signal was detected in the random sequences were considered as potential poly-tail. statistical analysis correlation tests were done with the method of spearman implemented in r. we accepted a correlation between data sets when the p-value was lower than.  strategy we systematically searched for all potential pseudogenes of ancient protein-coding genes within completely sequenced and annotated genomes of hemiascomycetes: saccharomyces cerevisiae and candida glabrata, zygosaccharomyces rouxii, kluyveromyces lactis, kluyveromyces thermotolerans, saccharomyces kluyveri, and debaryomyces hansenii and yarrowia lipolytica, members of the «ctg» group and of the dipodascaceae, respectively. to do this, we compared each genome sequence to the set of protein-coding sequences annotated in these genomes. we considered here as pseudogene any sequence that simultaneously -i- does not overlap an already annotated genetic element, -ii- shares sequence similarity with an annotated cds and -iii- shows disabling mutations in the reading frame: in-frame stop codon, frameshift mutation or truncation of more than 30% relatively to the cds. we chose this limit of 30% because, among the functional members of a given protein coding gene family, the length variation does not exceed 30% in the majority of cases. inactive pseudogenes resulting solely from mutations in promoters, as well as pseudogenes without detectable similarity among the eight studied genomes are, therefore, excluded from our analysis. similarly, we do not question the genome annotations: an annotated gene that appears truncated compared to other homologous sequences will not be considered as a pseudogene because truncation does not necessarily imply inactivation. complete set of detected pseudogenes and sequence degradation patterns we found a total of pseudogenes among the yeast genomes. they correspond to distinct protein families or singletons. the total number of pseudogenes varies significantly between species: from in c. glabrata to in y. lipolytica, corresponding to only % and % of the annotated cds, respectively. a. in megabases, except rdna. b. pseudogenes of protein coding sequences only. the percentage of pseudogenes relative to cds is indicated in parenthesis. c. proportion of « full-size » pseudogenes, i.e. extending over more than 70% of their bestmatch length. d. proportion of pseudogenes extensively truncated at their 3'-end. e. proportion of pseudogenes extensively truncated at their 5'-end. the pseudogenes identified exhibit a large panel of sequence degradation, ranging from only a few mutational disablements to extensive truncations. we distinguished here the «full-size» pseudogenes, i.e. those extending over more than 70% of their closest functional homolog, from the «truncated» pseudogenes. in s. cerevisiae, s. kluyveri, k. thermotolerans and k. lactis, the number of «full-size» pseudogenes exceeds truncated ones. the opposite is true for c. glabrata, z. rouxii, d. hansenii and y. lipolytica. there are as many pseudogenes truncated from their 3'-end, as those truncated from their 5'-end, consistent with the idea that truncation results from deletion at the dna level rather than incomplete cdna formation, characteristic of retro-processed pseudogenes. «full-size» pseudogenes contain more disabling mutations than truncated pseudogenes. we found four pseudogenes inactivated by insertion of a dna fragment of mitochondrial origin, as previously described. sequence divergence between pseudogenes and their best functional homologs in the absence of the actual ancestral coding sequences of the pseudogenes, we estimated the distances between pseudogenes and their closest functional homologs to measure their degree of mutational decay. distributions of the p-distances were calculated for each yeast species. they range from to globally. the distributions are different if one considers pseudogenes whose bestmatch is not in the same genome and pseudogenes whose bestmatch is a paralog in the same genome. in the first case, distances are high, consistent with the generally long evolutionary distances between studied species. in the second case, distances are much lower than that observed between functional paralogs in the same genome, suggesting that many pseudogenes arose after recent gene duplication events. the genomes of s. kluyveri and k. thermotolerans contain the highest proportions of highly diverged pseudogenes, while the genomes of s. cerevisiae and y. lipolytica contain the highest proportions of less diverged pseudogenes. note that functional paralogs in all these genomes are issued from both ancient and recent duplications. pseudogenes across species, conservation or independent formation a pseudogene should be free from any functional pressure and, therefore, should not be conserved over long evolutionary periods, unless it acquires a functional role. we, therefore, searched for pseudogenes conserved in several yeast species as a possible functional signature. a total of pseudogenes are shared between at least two species but of them are located outside of any conserved synteny block and probably correspond to independent pseudogenization events of homologous genes in different lineages. the remaining pseudogenes are found among five conserved synteny blocks, defining five groups of orthologous pseudogenes. note that we did not consider the syntenic orthologs of yil009c-a, ypl052w and yor239w as pseudogenes, because they all bear the programmed translational frameshift present in the functional genes of s. cerevisiae. according to the species phylogeny and following parsimony criterion, three of these sets of orthologous pseudogenes also correspond to independent pseudogenization events of the same gene in distinct lineages. pseudogenes sakl0gp and klth0gp corresponding to ydr285w of s. cerevisiae are conserved between s. kluyveri and k. thermotolerans. zip is a transverse filament protein of the synaptonemal complex that is required for normal levels of meiotic recombination and pairing between homologous chromosomes during meiosis. the orthologous pseudogenes are «full-size» and share a single common frameshift mutation. the most parsimonious hypothesis is that the frameshift mutation appeared in the common ancestor of these species. rec is a meiotic gene in s. cerevisiae, which has been lost several times in hemiascomycetes. this gene is also pseudogenized in three of the newly sequenced genomes studied. sakl0gp and klth0bp are conserved between s. kluyveri and k. thermotolerans, suggesting a pseudogenization event in their common ancestor. the orthologous gene is also pseudogenized in z. rouxii, and probably corresponds to another independent event. unlike in the previous case, the pseudogenes of this block of conserved synteny are truncated and each copy accumulated several species-specific disabling mutations. conservation of pseudogenes across s. cerevisiae strains the genome sequences for the s. cerevisiae non-reference strains are draft assemblies, not suitable for the pseudogene detection procedure developed in this work. however, to get a first estimation of pseudogene conservation among strains, we examined the conservation of the pseudogenes identified in the s. cerevisiae reference strain s288c among other sequenced strains. we found that of the pseudogenes identified in s288c are also pseudogenized in all other strains, suggesting that the pseudogenization event occurred in their common ancestor. the other pseudogenes of s288c sometimes correspond to intact coding sequences in one or a few other strains, indicating more recent pseudogenization events. of these pseudogenes contain at least two degrading mutations and of them are truncated at their 3' or 5' end. the only pseudogene specific to s288c contains only one internal stop codon. p: homologous pseudogene, intact: homolog with non-degraded coding sequence. distribution of the pseudogenes among protein families and corresponding functions based on sequence similarity, pseudogenes can be attributed to the functional gene families classified according to their predicted translational products. such families contain groups of orthologs between species, as well as groups of paralogs resulting from gene duplications and losses. in each yeast species, the majority of pseudogenes belong to gene families with a functionally characterized s. cerevisiae member, which can be used to infer the probable function of the gene that was pseudogenized. interestingly, we observed a bias towards transporters, proteins acting at the periphery of the cell, and enzymes. when normalized against the ratios of these functional categories among active genes, the bias is conserved for transporters and proteins acting at the periphery of the cell, but not for enzymes. among pseudogenes with no ascribed function, there is no over-representation of domains related to these three categories. although the number of pseudogenes is not directly correlated to the number of its active paralogs in a given family, we observed frequent formation of pseudogenes in these functional categories. for example, there are pseudogenes in d. hansenii for the sugar transporter family, pseudogenes in c. glabrata for the lectin-like protein family, and pseudogenes in z. rouxii for the dup gene family coding for membrane proteins. s.c. is for s. cerevisiae. anumber of pseudogenes with a homolog in s. cerevisiae coding for proteins involved in transport and/or acting at the periphery of the cell. bnumber of pseudogenes without homolog in s. cerevisiae coding for enzymes. a total of pseudogenes are similar to singleton genes, and belong to gene families specific to a given species for which, unless in s. cerevisiae, there is usually no functional indication. clustering of the pseudogenes in subtelomeres gene densities are nearly constant along yeast chromosomes, except for subtelomeres where the number of active genes is reduced, and where genes can be transcriptionally silenced. we, therefore, examined the distribution of the pseudogenes along the chromosomes, separating the subtelomeres from the central regions. in all yeast species, pseudogenes exceed active genes in number in the subtelomeres. in s. cerevisiae, c. glabrata and z. rouxii, more than half of the pseudogenes are found in subtelomeres. this proportion is lower in k. thermotolerans, d. hansenii, s. kluyveri and k. lactis, and drops to only 5% in y. lipolytica. the presence of rdna loci in the subtelomeres of y. lipolytica might prevent some pseudogenes accumulating in these regions, as well as in d. hansenii and in c. glabrata. a. percentage of pseudogenes in subtelomeric regions. b. percentage of active genes in subtelomeric regions. c. number of pseudogenes without annotated functional paralog in the genome. d. number of pseudogenes with annotated functional paralog whose closest homolog is in the same genome or in another genome. *. species for which all chromosomes are fully sequenced, including their telomeric repeats. possible origin of the pseudogenes we tried to define the origin of yeast pseudogenes based on the presence or absence of a paralog in the same genome and on the conservation of synteny between species. the presence of an active paralog in the genome reveals the occurrence of a previous duplication event, followed by mutational inactivation of one of the copies. the absence of any paralogs reveals the degradation of a single copy gene, hence the loss of the corresponding function in the species. the conservation of synteny is a signature for mutational sequence degradation at the origin of the pseudogene. the non-conservation of synteny is compatible with a species-specific duplication, either a retro-transposition event or a segmental duplication. pseudogenes originating after species-specific gene duplications more than half of the pseudogenes arose after a gene duplication event specific to their own genome. this corresponds to the vast majority of pseudogenes in s. cerevisiae, d. hansenii and y. lipolytica, ca. 70% of the pseudogenes in c. glabrata and z. rouxii, and less than 60% in k. lactis, k. thermotolerans and s. kluyveri. the majority of these pseudogenes probably arose after segmental duplication, but some may be processed pseudogenes, as suggested by a poly tract at their 3'-end or their location next to a retrotransposon-related sequence. this would correspond to only % of the pseudogenes in y. lipolytica, and 20% of the pseudogenes in k. thermotolerans. only one pseudogene was identified by the lack of intron compared to its closest functional homolog. a. number of pseudogenes originating from mutational inactivation of a duplicated gene copy formed after speciation. b. number of pseudogenes originating from mutational inactivation of a duplicated gene copy formed before speciation. c. number of pseudogenes originating from mutational inactivation of a single copy gene. d. number of duplicated pseudogenes among the first category. e. number of pseudogenes being part of a duplicated segment involving other adjacent genes among the first category. f. number of retro-processed pseudogenes, among the first category, identified by: either the presence of a 3' poly-tail or the proximity of retrotransposon-related sequence. in each species, the candidates identified by these criteria are different. interestingly, we found duplications of pseudogenes, as evidenced by their common pattern of degrading mutations. for example, in y. lipolytica, we found five pseudogenes similar to yali0a14927g, which, based on phylogenetic analysis, suggest two original pseudogenization events in the history of this family. note that these pseudogenes are not located in subtelomeric regions and are not part of larger duplicated regions including other genes. these pseudogenes have thus not been maintained in the genome by some selective pressure on the duplication of adjacent genes. similarly, two such cases of pseudogene multiplication were encountered in k. lactis, as well as cases of pseudogene duplication. most pseudogenes located at the chromosome ends arose after a species-specific gene duplication: from 60% in k. thermotolerans to 100% in y. lipolytica. see table for legend. pseudogenes originating from ancestral gene duplications pseudogenes formed after ancestral gene duplications are found in all yeast genomes. four pseudogenes, sace0bp cagl0dp caglohp and cagl0mp were formed after the whole-genome duplication event that occurred in the common ancestor of c. glabrata and s. cerevisiae. this is a very small number compared to the extensive loss of duplicated genes that occurred by deletion. strikingly, there are more pseudogenes originated from other ancestral duplications in protoploid saccharomycetaceae and in d. hansenii and y. lipolytica, suggesting that selective pressure on duplicated genes are different after whole-genome duplication and other duplication events, such as segmental duplications. pseudogenes originating without previous gene duplication and putative loss of function pseudogenes with no functional paralog are found in all the studied genomes and could, therefore, correspond to a functional loss. the genomes of k. lactis, k. thermotolerans, s. kluyveri and c. glabrata contain the highest proportions of such pseudogenes. the lowest proportions are found in d. hansenii and in y. lipolytica, consistent with the higher gene redundancy in these genomes. there is no bias towards any functional category. pseudogenes homologous to essential genes in s. cerevisiae are found in all species, except in k. lactis. repeated pseudogenization of the same gene is also encountered in this category. for example, pseudogenes cagl0mp and zyro0gp similar to the essential chromosomal passenger gene yjr089w, have no homolog in their respective genome. finally, we noticed three interesting pseudogenes in s. cerevisiae that correspond to horizontally acquired genes. sace0fp and sace0fp have no homolog in s. cerevisiae but were identified by similarity to deha2d01122g and deha2e07282g, respectively. absence of these genes among hemiascomycetes suggests that four independent horizontal gene transfer events have occurred, the two genes in s. cerevisiae being secondarily pseudogenized. the third one, sace0np is a duplicated and identical copy of sace0fp revealing the expansion of horizontally transferred inactivated genes. the transfers must have occurred in an ancestor of all the studied strains of s. cerevisiae because the homologs of sace0fp sace0fp and sace0np are pseudogenized in the other strains examined. discussion we present here the first systematic search for pseudogenes in the genomes of eight distinct yeast species spanning the whole evolutionary spectrum of hemiascomycetes. because our method of detection relies on sequence similarity with annotated protein coding genes, the total number of pseudogenes identified represents a minimal estimate. among the genomes analyzed here, only s. cerevisiae has been the subject of previous systematic analysis of the pseudogene content. our analysis revealed a smaller number of pseudogenes in this species because we deliberately ignored partial overlaps with annotated genetic elements. on the contrary, additional pseudogenes were identified based on our multi-species sequence comparisons. these differences illustrate the difficulties in listing all pseudogenes in a given genome. despite this fact and numerical variations between yeast species, our results demonstrate that the proportion of pseudogenes compared to active genes remains low in all hemiascomycetous genomes studied. our comparative analysis of pseudogenes across an entire yeast phylum provides a unique data set to examine their origin and evolutionary conservation. we found that pseudogenes in yeasts are formed either by disabling mutations or by extensive truncations. the general absence of conservation of pseudogenes between yeast species is consistent with their large evolutionary distances and indicates that new pseudogenes were formed within each lineage. however, intra-species conservation of pseudogenes is high: about 80% of the pseudogenes in the s. cerevisiae reference strain are old enough to be conserved among other strains of this species. poor conservation of pseudogenes is also observed between mammalian species. most of the pseudogenes correspond to duplicated gene copies, illustrating the extensive dynamics of gene duplications in the yeast genomes, most probably segmental duplications, as observed in s. cerevisiae. only few pseudogenes correspond to duplicated gene copies formed by the ancestral whole-genome duplication common to c. glabrata and s. cerevisiae. this is consistent with the idea that most duplicated copies were lost by complete deletion, as previously proposed, or with the possibility that pseudogenes have been degraded beyond recognition, given the time elapsed since the whole-genome duplication event. among the protoploid saccharomycetaceae, the presence of species-specific pseudogenes within conserved synteny blocks indicates their relatively recent formation, to the notable exception of two ancestral pseudogenes conserved in s. kluyveri and k. thermotolerans, the more closely related pairs of species studied. the general absence of conserved pseudogenes confirms that most of the observed pseudogenes are on their way to complete degradation. according to doniger et al., pseudogenes containing an internal in-frame stop codon and pseudogenes containing a frameshift mutation are present in either m or yps while they correspond to intact coding sequences in s288c. from the alignment of the three genomes within coding regions, it appears that the loss of a gene by deletion is much more frequent ) than the creation of a pseudogene by a substitution event or by an indel event ). according to lynch et al., the rate of gene loss, × 10- per gene per cell division, is much higher than the base-substitutional rate, × 10- per site per cell division. it should then be much more frequent to lose a gene by deletion than to create a pseudogene. unlike their abundance in mammalian genomes, only few pseudogenes may have originated from retro-transposition events in yeasts. but their presence suggests that, as experimentally demonstrated in s. cerevisiae, retro-transposition occurred in these genomes, although a very small number of active retrotransposons are usually present. the number of retro-processed pseudogenes in yeasts is, however, probably underestimated because their identification is difficult given the small number of intron-containing genes. the few detected pseudogenes homologous to intron-containing genes are all truncated and do not span the intron insertion-site, except in one case. pseudogenes without functional paralog suggest functional loss in the corresponding species, unless a non-homologous gene encodes a similar function. in s. cerevisiae, the fact that about half of the pseudogenes correspond to unknown functions is striking given that there remain only 17% of the genes not yet functionally characterized in this species. this suggests that the divergence of functional repertoire between yeast species primarily concerns functions not yet identified. the frequent occurrence of pseudogenes corresponding to transporters or protein acting at the periphery of the cell may be correlated with the tendency of such genes to cluster in the subtelomeric regions, which are highly dynamic in terms of gene duplications and losses and often concern functions involved in the adaptation of the species to its environment. the highly dynamical behavior of such families is supported by the fact that to times more pseudogenes than genes are found in the subtelomeric regions. the vast majority of pseudogenes in the subtelomeres originated from a species-specific duplication. this bias could be correlated to the accelerated base-pair substitution observed in the subtelomeres, which probably also contributes to adaptive evolution. it also suggests that pseudogenes issued from ancestral duplication and function loss are not preferentially maintained in subtelomeres with respect to central regions of chromosomes. in s. cerevisiae, pseudogenes of highly connected genes are significantly under-represented, possibly reflecting the lower propensity of gene loss among these genes: 6% of the pseudogenes of genes coding for protein in complexes against 28% of the active genes and 59% of pseudogenes of genes with a genetic interaction profile against 75% of the active genes. whether pseudogenes correspond to intermediate gene states before complete erasure or steady states conferring selective advantages remains an open question and may be case-specific. our data provide no indication as to whether pseudogenes may be transcribed or not. however, according to available experimental data in s. cerevisiae, of the detected pseudogenes appear to be transcribed. evidence of transcription also exists in c. glabrata and y. lipolytica. the transcription products of pseudogenes may be directly targeted to the degradation machinery such as nmd, but they could also play some role in the cell by interfering with the expression of functional genes. moreover, as previously suggested, pseudogenes with few disabling mutations could constitute a reservoir of functional protein products if recoding events occur, such as programmed frameshift, or bypass of the stop-codon. the number of disabling mutations within each pseudogene is given in columns and the relative length of the pseudogene with respect to its closest functional homolog is indicated in column last column indicates the reference of the data set where evidence of transcription is found. for a given chromosome, all identified pseudogenes are separated by several genes. there is no bias of any kind among these pseudogenes.  pseudogenes are found in all yeast genomes, albeit in limited number compared to genomes of multi-cellular eukaryotes. they mostly result from lineage-specific mutational degradations that may correlate with species adaptation to their environment. yeast pseudogenes show a wide range of mutational alterations, consistent with their rapid evolution, hence their absence of conservation between species. along with complete gene deletion, pseudogene formation contributes to the rapid genome evolution by gene duplication and loss in yeasts. the paucity of observed pseudogenes across the entire phylum of hemiascomycetes suggests that pseudogene formation is not the main mechanism of gene loss within these genomes. this could be explained by the low estimated rate for pseudogene formation across s. cerevisiae strains, compared to the estimated rate of gene deletion. however, despite their unlikely occurrence, pseudogenes do exist in the yeast genomes. they appear mainly by species-specific duplications and testify for the adaptation of the cell to its environment. their poor conservation across species suggests that most of them are likely to disappear. authors' contributions il designed and conducted the analysis, implemented the algorithms and wrote the manuscript. bd contributed to the interpretation of the results and to the writing of the manuscript. both authors read and approved the final manuscript. supplementary material additional file 1 table s and table s - table s list and characteristics of the detected pseudogenes in the studied yeast genomes. - table s list of gene families with detected pseudogenes. click here for file additional file 2 table s and descriptions of tables s and s descriptions of tables s and s - table s lost functions among the eight studied species. a list of the pseudogenes with no paralog in the different genomes click here for file acknowledgements we thank gilles fischer, bertrand llorente and eduardo rocha for helpful comments on this work and our colleagues from the "unité de génétique moléculaire des levures" and from the génolevures consortium for fruitful discussions. this work was supported by grant anr-05-blan- from the agence nationale de la recherche. bd is a member of the institut universitaire de france.
although a variety of functions are affected by pseudogenization, there is a bias towards functions involved in the adaptation of the yeasts to their environment, and towards genes of unknown functions. we found that most pseudogenes in yeasts originate from mutational degradation of gene copies formed after species-specific duplications but duplications of pseudogenes themselves are also encountered. in all yeasts, except in y. lipolytica, pseudogenes tend to cluster in subtelomeric regions where they can outnumber the number of functional genes from to times. reiterated pseudogenization of some genes is often observed in different lineages and may affect functions essential in s. cerevisiae, which are, therefore, lost in other species. our work illustrates for the first time the formation of pseudogenes in different branches of hemiascomycetous yeasts, showing their limited conservation and how they testify for the adaptation of the yeasts functional repertoires. pseudogenes are ubiquitous genetic elements that derive from functional genes after mutational inactivation. among yeasts, only the genome of the s. cerevisiae reference strain has been analyzed so far for pseudogenes. we present here the first comparative analysis of pseudogenes within the fully sequenced and annotated genomes of eight yeast species, spanning the entire phylogenetic range of hemiascomycetes. pseudogenes are generally not conserved between the yeast species studied, consistent with their large evolutionary distances, but tend to be conserved among s. cerevisiae strains.
 pseudogenes are ubiquitous genetic elements that derive from functional genes after mutational inactivation. characterization of pseudogenes is important to understand genome dynamics and evolution, and its significance increases when several genomes of related organisms can be compared. among yeasts, only the genome of the s. cerevisiae reference strain has been analyzed so far for pseudogenes.  we present here the first comparative analysis of pseudogenes within the fully sequenced and annotated genomes of eight yeast species, spanning the entire phylogenetic range of hemiascomycetes. a total of pseudogenes were found, out of which mutational degradation patterns and consequences on the genetic repertoire of each species could be identified. we found that most pseudogenes in yeasts originate from mutational degradation of gene copies formed after species-specific duplications but duplications of pseudogenes themselves are also encountered. in all yeasts, except in y. lipolytica, pseudogenes tend to cluster in subtelomeric regions where they can outnumber the number of functional genes from to times. pseudogenes are generally not conserved between the yeast species studied, consistent with their large evolutionary distances, but tend to be conserved among s. cerevisiae strains. reiterated pseudogenization of some genes is often observed in different lineages and may affect functions essential in s. cerevisiae, which are, therefore, lost in other species. although a variety of functions are affected by pseudogenization, there is a bias towards functions involved in the adaptation of the yeasts to their environment, and towards genes of unknown functions.  our work illustrates for the first time the formation of pseudogenes in different branches of hemiascomycetous yeasts, showing their limited conservation and how they testify for the adaptation of the yeasts functional repertoires. 