Combined phylogeny and neighborhood analysis of the evolution of the ABC transporters conferring multiple drug resistance in hemiascomycete yeasts
BMC Genomics volume 10, Article number: 459 (2009)
Pleiotropic Drug Resistant transporters (PDR) are members of the ATP-Binding Cassette (ABC) subfamily which export antifungals and other xenobiotics in fungi and plants. This subfamily of transmembrane transporters has nine known members in Saccharomyces cerevisiae.
We have analyzed the complex evolution of the pleiotropic drug resistance proteins (Pdrp) subfamily where gene duplications and deletions occur independently in individual genomes.
This study was carried out on 62 Pdrp from nine hemiascomycetous species, seven of which span 6 of the 14 clades of the Saccharomyces complex while the two others species, Debaryomyces hansenii and Yarrowia lipolytica, are further apart from an evolutive point of view.
Combined phylogenetic and neighborhood analyses enabled us to identify five Pdrp clusters in the Saccharomyces complex. Three of them comprise orthologs of the Pdrp sensu stricto, Pdr5p, Pdr10p, Pdr12p, Pdr15p, Snq2p and YNR070wp. The evolutive pathway of the orthologs of Snq2 and YNR070w is particularly complex due to a tandem gene array in Eremothecium gossypii, Kluyveromyces lactis and Saccharomyces (Lachancea) kluyveri. This pathway and different cases of duplications and deletions were clarified by using a neighborhood analysis based on synteny.
For the two distant species, Yarrowia lipolytica and Debaryomyces hansenii, no neighborhood evidence is available for these clusters and many homologs of Pdr5 and Pdr15 are phylogenetically assigned to species-based clusters.
Two other clusters comprise the orthologs of the sensu lato Pdrp, Aus1p/Pdr11p and YOL075cp respectively. The evolutionary pathway of these clusters is simpler. Nevertheless, orthologs of these genes are missing in some species.
Numerous duplications were traced among the Hemiascomycetous Pdrp studied. The role of the Whole Genome Duplication (WGD) is sorted out and our analyses confirm the common ancestrality of Pdr5p and Pdr15p. A tandem gene array is observed in Eremothecium gossypii. One of the copies is the ortholog of Snq2 while the other one is lost in the post-WGD species.
The neighborhood analysis provides an efficient method to trace the history of genes and disentangle the orthology and paralogy relationships.
The phenotype for pleiotropic drug resistance (PDR) in Saccharomyces cerevisiae (SACE) was discovered in 1973 by RANK and BECH-HANSEN, who reported that single gene mutations were responsible for resistance to multiple drugs of different chemical structures and different targets. The genetic and molecular mapping revealed that these mutations are located in the transcription factors Pdr1p and Pdr3p [2–4] and increased the expression of a series of target genes  including those encoding the ABC efflux pumps Pdr5p [6–8], Snq2p , Pdr10p and Pdr15p [10, 11] or Pdr11p .
Analysis of the full genome sequence of Saccharomyces cerevisiae revealed that these transporters are members of a large phylogenetic subfamily of ABC transporters named Pdrp (for pleiotropic drug resistance proteins), where nuclear binding folds (NBF) alternate with domains comprising six predicted transmembrane spans (TMS) to form a NBF-TMS-NBF-TMS pattern . NBF includes a distinctive series of amino acid sequence motifs, including Walker A, Walker B and ABC signatures involved in ATP binding and hydrolysis.
The inventory of all full-sized Saccharomyces cerevisiae ABC transporters reveals six Pdrp sensu stricto namely: Pdr5p, Pdr10p, Pdr12p, Pdr15p, Snq2p and YNR070wp , characterized by three criteria: (1) The alternation NBF-TMS-NBF-TMS, which is different from the classical TMS-NBD-TMS-NBD topology reported for the mammalian multidrug resistance MDR and MRP drug efflux pumps [14, 15]; (2) The presence of a cysteine residue (C) instead of the lysine residue (K) in N-terminal Walker A motifs as well as the specific NVEQ motif in the C-terminal ABC signature [6, 13]; (3) The efflux of multiple drugs [16–18].
The structural determinants for the broad substrate specificity of Pdrp pumps have not yet been identified even though dozens of publications have reported hundreds of amphiphilic substrates for Pdr5p and Snq2p. These two pumps share many substrates even though some drugs are preferentially handled by either one or the other [17, 19]. In contrast, Pdr12p exhibits clear substrate specificity for short chain weak organic acids .
Fully sequenced Ascomycota and Basidiomycota were all examined and found to contain several Pdr proteins  but their evolutionary lineage is unclear. A previous phylogenetic classification of Pdrp from five species spanning the full Hemiascomycetes phylum identified seven subclusters . GBELSKA et al.  showed that in contrast to the other ABC subfamilies which remain as single copy in each species during evolution, a series of gene duplications occurs independently in the PDR subfamily in five hemiascomycetous species. GAUR et al.  compared the NBF domains of S. cerevisiae and Candida albicans Pdrp. They concluded that C. albicans contains a cluster of four proteins homologs to Pdr5p and a single protein homolog to Snq2p.
In addition to the Pdrp clusters sensu stricto, the S. cerevisiae genome contains two other ABC phylogenetic clusters exhibiting the alternation of NBF-TMS-NBF-TMS domains . One cluster comprising Aus1p and Pdr11p are sterol-influx pumps  while another cluster comprising YOL075cp is of undetermined function. As the members of these two phylogenetic clusters show neither drug efflux properties nor the K/C substitution in the Walker A1 motif, we classified them as Pdrp sensu lato.
The aim of the present paper is to trace the evolution of the Pdrp subfamily in nine species belonging to the Hemiascomycetes phylum, believed to have separated from the filamentous fungi 300 to 400 million years ago [26, 27]. Seven of these species span the 14 clades constituting the Saccharomyces complex . Among those, three novel genome sequences from the Lachancea or Zygosaccharomyces clades were analyzed: Saccharomyces (Lachancea) kluyveri (SAKL), Zygosaccharomyces rouxii (ZYRO) and Kluyveromyces (Lachancea) thermotolerans (KLTH) . Two hemiascomycetous sequences of older lineage; Debaryomyces hansenii (DEHA) and Yarrowia lipolytica (YALI) were also explored. An exhaustive database of all the chromosomal neighbors of the PDR genes from all nine hemiascomycetous species was developed (SERET et al., in preparation) and allowed us to carry out combined phylogenetic and neighborhood analysis. Gene neighborhood analysis has been pioneered by BYRNE and WOLFE for tracing the ancestral gene blocks produced by the Whole Genome Duplication (WGD). We will illustrate the extent and the limits of the information provided by analysis of gene neighborhood data for deciphering the evolution of the large Pdrp subfamily in addition to classical phylogenetic analysis based only on amino acid sequence similarity.
The Pdrp transporters are a well characterized family of genes implied in drug resistance. All genes share typical features, such as Walker A and B and an ABC signature, and are very similar in terms of sequence. In the Génolevures database, the family GL3C0025 which contains all SACE Pdrp comprises 62 members (Table 1) of which 6 are considered as pseudogenes or fragments (Additional file 1). The phylogenetic tree in Figure 1, based exclusively on full-length proteins to improve multiple alignment, shows that the 56 full-Pdr proteins of the GL3C0025 family (Additional file 2) branch into five clusters labeled A to E. Each cluster is named according to its S. cerevisiae member(s): A, the cluster Pdr12p (8 members); B, the cluster Snq2p/YNR070wp (11 members); C, the cluster Pdr5p/10p/15p (25 members) which ramifies into 4 subclusters labeled C1, C2, C3, C4; D, the cluster YOL075cp (7 members); E, the cluster Aus1p/Pdr11p (5 members).
Neighborhood analysis (Additional file 3) was carried out to delineate Subsets of Orthologs by Neighborhood and Similarity (SONS). Seven SONS labeled "a" to "g" (Figure 2) provide information on the evolution of 40 proteins displaying informative neighborhood among the 62 PDR from the GL3C0025 family.
For most genes, the phylogenetic clusters are confirmed by the neighborhood analysis. But in some cases, divergences are observed. In the next sections, each phylogenetic cluster of genes is discussed in an evolutionary perspective integrating the evidence from the neighborhood analysis.
Cluster A (Pdr12p)
This cluster comprises 8 members in SACE, CAGL (Candida glabrata), ZYRO, SAKL, KLLA (Kluyveromyces lactis) and YALI (Figure 1). No full-Pdrp orthologs are detected in KLTH, even though the neighborhood is well maintained on chromosome E, ERGO (Eremothecium gossypii) and DEHA. The neighborhood analysis (Figure 2 SONS a) shows a clear connection (9 neighbors) between SAKL0H10670g and the PDR belonging to the other species, except YALI. A second PDR copy is present in SAKL on chromosome C. SAKL0CO5654g only shares one neighbor, belonging to the GL3R2315 family. The very high amino acid identity (98%, 17 substitutions) between the two paralogs suggests a very recent duplication. A tandem gene array occurred in ZYRO where the two copies diverged a little more (87% amino acid identity). None of these duplications (on different chromosomes in SAKL and within the same chromosome in ZYRO) were maintained in other species. Pdr12p and Snq2p are in closely related clusters as they share a high similarity rate. Despite this proximity of sequence, neighborhood evidence shows that they evolved separately, at least since the speciation time of KLLA.
Cluster B (Snq2p/YNR070wp)
The phylogenetic Snq2p cluster comprises 11 members that are represented in all the examined species except YALI. A shared neighborhood is present from ERGO to SACE (Figure 2 SONS b). It comprises a single ortholog of Snq2p in each species except in SACE where YNR070wp is a paralog of Snq2p. Moreover, neighborhood analysis shows an apparent tandem gene array in ERGO chromosome B. However, the two copies are different (38% amino acid identity). ERGO0B8140p has 63% identity to Snq2p and only 41% identity to Pdr5p (Figure 3). In contrast, ERGO0B8162p has a higher amino acid identity to Pdr5p (53%) than to Snq2p (36%). A similar heterologous tandem in a similar neighborhood is observed for KLLA and SAKL with, however, the insertion of an increasing number of genes between the two members of the repeat. ZYRO also presents two genes in this neighborhood but they are no longer in tandem gene array. It is noteworthy that KLTH only has one gene, an ortholog of Snq2, in this neighborhood.
The phylogenetic tree in Figure 1 also shows that paralogs of Snq2 are observed in SACE (YNR070w), ZYRO (ZYRO0B14762g) and SAKL (SAKL0H21010g) in non-shared chromosomal environments. YNR070w and ZYRO0B14762g probably derive from a species specific duplication event. The case of SAKL0H21010g suggests a duplication followed by a rapid evolution of the copy. In terms of identity, the DEHA0A03696g gene, which does not share any neighbor with any Snq2 or Pdr12's homologs, is closer to Snq2p (56%) than to Pdr12p (44%).
Cluster C (Pdr5p/10p/15p)
The phylogenetic C cluster contains 25 members including the SACE Pdr5p/10p/15p. Figure 2 SONS c shows that in the pre-WGD species ZYRO, three copies of Pdr5/15 orthologs (ZYRO0D11836g, ZYRO0D11858g, ZYRO0D11880g) are organized in tandem on chromosome D. The triplicated ortholog copies of Pdr5/15 from ZYRO are inserted in a neighborhood that can be traced back in KLLA, SAKL and KLTH but none of these chromosomal environments contain a Pdr5/15 homolog. However, in each of these species a Pdr5/15-like copy is present in a neighborhood different from that of the ZYRO triplicate. The KLTH, SAKL and KLLA copies show a strong phylogenetic link with Pdr5p/15p: nearly 67% amino acid identity for KLTH0G19448p, nearly 70% amino acid identity for SAKL0C08008p and nearly 65% amino acid identity for KLLA0F21692p. KLLA0F21692g and SAKL0C08008g share a common neighborhood (see Figure 2, SONS f) which is not shared with KLTH0G19448g (see Figure 2, SONS g).
The ancestor of the ZYRO copies was duplicated by the WGD and produced the ohnologs Pdr5 and Pdr15. In fact, Pdr15, which is located on chromosome D, possesses two common neighbors with ZYRO while Pdr5, located on chromosome F, shares five common neighbors with ZYRO. Only the Pdr15 copy was retained in CAGL in this neighborhood. A second PDR, CAGL0M01760g, with a high amino acid identity with Pdr5p/15p is present in CAGL in another neighborhood.
Pdr10p shows high amino acid identity to Pdr5p/15p (65% and 62% respectively). Segment alignment reveals that Pdr10p can be differentiated from other SACE Pdrp by a 22 amino acid insert (SDAAIMGNDKTVAKEHYSSPSS), located at the intersection of the two half PDR-regions (NBF-TMS-insert-NBF-TMS). This marker means that Pdr10p can be detected in S. cerevisiae, Saccharomyces bayanus, Saccharomyces paradoxus and Saccharomyces mikatae but no ortholog was found in CAGL and in the pre-WGD species (data not shown, inferred from Broad Institute ). The event that gave rise to Pdr10p is thus post-WGD.
In addition to the major subcluster C1 (Figure 1) comprising Pdr5/10/15p, cluster C contains three additional phylogenetic subclusters. Two of these subclusters show a spectacular species-based clustering. Indeed, very similar proteins tend to cluster with proteins belonging to the same species rather than with orthologs in other species. One of the subcluster (C3) comprises five DEHA members. The neighborhood analysis suggests a putative link between DEHA2F16478g and the orthologs of YOL075c. The other subcluster (C4) comprises four YALI members. Neighborhood analysis of this subcluster is not conclusive but phenotype analysis reveals that one of them (YALI0C20265p) is an azole efflux transporter with properties similar to those of Pdr5p/Snq2p  while another (YALI0E14729p) has the unique property to export C16 hydrocarbon residues .
The last subcluster (C2) contains five members:SAKL0C06996g, SAKL0C11704g, KLLA0D03476g as well as ERGO0B08162g and ERGO0G05126g which are very similar (73% amino acid identity). Three of them, ERGO0B08162g, KLLA0D03476g and SAKL0C11704g are in tandem gene array. Interestingly, they belong to the C cluster (Pdr5p/10p/15p) when using phylogenetic criteria but share some neighbors with the B cluster (Snq2p). This is also the case of ZYRO0D17710g of the subcluster C1.
Cluster D (YOL075cp)
This cluster comprises a total of 7 members; one per species except in DEHA and ERGO. A spectacular shared neighborhood presence is observed for the orthologs of YOL075c (see Figure 2, SONS d). Eight common neighbors are, for instance, identified in the post-WGD species SACE and CAGL, four common neighbors are identified in the pre-WGD species and four common neighbors link ZYRO, one of the pre-WGD species to the post-WGD species. The YALI homologue (YALI0D25828p) has no conclusive neighborhood but it exhibits high amino acid identity (41%) to the SACE YOL075cp. This suggests that all hemiascomycetous orthologs of YOL075c originate from an ancestor existent prior to the speciation event of YALI.
The neighborhood analysis links a tandem gene array made up of DEHA2F16478g and the fragment DEHA2F16500g to cluster D. However, the DEHA2F16478g is classed in subcluster C3 by phylogeny and does not display the specific motifs of YOL075cp orthologs.
Cluster E (Aus1p/Pdr11p)
This cluster comprises members of four species only: SACE, CAGL, SAKL and ERGO. Using neighborhood analysis, an ortholog of Aus1p was only identified in CAGL. There is no evidence that the two SACE paralogs, Aus1p on chromosome O and Pdr11p on chromosome I, issue from the WGD as there is no shared neighborhood with pre-WGD species. In SACE, both Aus1p and Pdr11p are essential for anaerobic growth as they contribute to the uptake of ergosterol under conditions which prevent the aerobic synthesis of sterols . The sterol transporter function was recently demonstrated for the CAGL ortholog of Aus1p . The two characteristic Pdrp motifs, N-terminal Walker A and the C-terminal ABC signature, remain, but only partially (Additional file 1).
Six fragments of 72 to 1099 amino acids belonging to the GL3C0025 protein family have been identified. Despite their subtelomeric location, that we define as the last 30 kb of each chromosome following FAIRHEAD and DUJON, the two fragments KLTH0E17138g and KLTH0H00110g share the neighborhood of KLTH0G19448g, the ortholog of Pdr5/15, which may suggest duplications followed by pseudogenization. Because phylogenetic analysis of fragments may be error-prone, their evolutionary significance will not be analyzed.
Orthologous PDR transporters with typical NBF-TMS-NBF-TMS topology are found in all higher fungi such as S. cerevisiae[6, 13], the hemiascomycetous human pathogen C. albicans, Euascomycetes, Archiascomycetes[37, 38], Basidiomycetes and all plants  including Arabidopsis thaliana and Nicotiana plumbaginifolia. In contrast, the PDR topology has not been detected in protozoa  or chordates  where multiple drug resistance results from the overexpression of MDR/MRP efflux transporters of opposite TMS-NBF-TMS-NBF topology.
Our analysis of 56 full-sized Pdrp from 9 different Hemiascomycetes shows that the Génolevures subfamily GL3C0025  is subdivided into five phylogenetic clusters labeled A to E (Figure 1). Each of these clusters can be characterized by its SACE members originally identified by DECOTTIGNIES and GOFFEAU. Within the phylogenetic Saccharomyces complex defined by KURTZMAN, a total of 40 Pdrp homologs sharing significant chromosomal neighborhood were classified in seven SONS, labeled "a" to "g" (Figure 2), in which phylogenetic and neighborhood data are combined.
It now appears that the Pdrp sensu stricto phenotype can be allocated only to cluster A (Pdr12p), to cluster B (Snq2p/YNR070wp) and to cluster C (Pdr5p/10p/15p), given that all three share not only the NBF-TMS-NBF-TMS topology but also the two additional Pdrp traits: efflux drug pumping and typical Walker A1 and ABC Signature 1 motifs.
On the other hand, no experimental information is available concerning the unknown substrates of the members of cluster D (YOL075cp), while the SACE members of cluster E (Aus1p/Pdr11p) are clearly sterol influx pumps . Moreover, as neither cluster D nor cluster E contains members with the typical Pdrp motifs, they may be considered to be Pdrp sensu lato.
Pdr12p is a Pdr sensu stricto, according to our three criteria. Its origin remains questionable as no Pdr12p neighbors are shared either with Pdr5p/15p or Snq2p. In SACE, the function of Pdr12p, which effluxes food preservatives such as short chain weak acids , is different from that of Snq2p and Pdr5p which share the function to efflux a series of antifungal azoles and other hydrophobic substrates . Nevertheless, the SACE Pdr12p sequence is closer to SACE Snq2p (46%) than to SACE Pdr5p (38%). This may suggest a common origin of both Pdr12p and Snq2p clusters. The sequence analysis of new yeast species, phylogenetically situated between DEHA and ERGO, should make it possible to disentangle the exact evolution lineage.
The unlinked ancestor of Pdr5p/15p ohnologs of KLTH jumped and triplicated on ZYRO chromosome D. An ancestor of the three copies, already in this neighborhood, was duplicated through WGD and produced Pdr5p and Pdr15p in all Saccharomycetes.
The phylogenetic tree of the Pdr5p/15p clusters shows several subclusters (C3 and C4) which aggregates members belonging to specific species such as DEHA or YALI. These are cases where species-based clustering is suspected of hindering a functional inference based on sequence similarity only. For example, the subcluster C4 contains four YALI members belonging to four different chromosomal fragments in unshared environments. Despite being inside the same subcluster, theses genes seem to have very different functions. One of them, YALI0E14729g, exhibits the function of alkane efflux, unique as yet, while another, YALI0C20265g, seems to control the classic azole resistance function of the SACE Pdr5p/15p pumps.
The evolutionary patterns of the Pdrp sensu lato are simpler and totally independent from that of the Pdr sensu stricto. The members of Pdr sensu lato Aus1p/Pdr11p cluster are present in four different species only: ERGO, SAKL, CAGL and SACE. These genes may be lost in other species. Indeed, sterol influx pumps are essential only under anaerobiosis conditions. However, no clear relation between the presence of homologs of Aus1p or Pdr11p sterol-influx pumps and facultative anaerobic growth could be detected among the species composing the Saccharomyces complex .
The members of the Pdr sensu lato YOL075c cluster are of unknown function and belong to species of both ancestral clade (YALI0D25828g) and recent clades. No homolog of YOL075c is found in DEHA and ERGO. Extensive neighborhood preservation is observed for all Kluyveromycetes and Saccharomycetes species. Even though the substrates of this cluster of putative transporters are unknown, its continued presence in many pre- and post-WGD species from the Saccharomyces complex indicates that these gene products exert an important physiological function. This function and the subcellular localization of the transporter have not yet been determined. Curiously all members of the YOL075c SONS have at least 100 less amino acid residues (1247 to 1328) than all other members of the GL3C0025 family. Moreover, the cysteine of SGCT motif of the N-Walker A is replaced by the "classic" lysine residue and the NVEQ motif of the C-Walker Signature is replaced by SGGE (Additional file 1).
In three cases only, we observed a divergence between the neighborhood and phylogenetic analyses. ERGO0B08162g, KLLA0D03476g and SAKL0C11704g are allocated to the C cluster (Pdr5p/15p) by phylogeny and to the SONS b (Snq2p/YNR070wp) by neighborhood. This divergence can be explained by the fact that they belong to a tandem gene array in which a copy has been subjected to neofunctionalization. In another case, ZYRO0D17710g (subcluster C1) is allocated to SONS b (Snq2p/YNR070wp) by neighborhood. Its high amino acid similarity with Pdr5p (70%) and significant difference of sequence with SAKL0C11704g and ZYRO0A04114g suggest an evolution based on gene conversion with one of the member of the ZYRO0D11836g, ZYRO0D11858g and ZYRO0D11880g tandem gene array. Finally, the divergence observed for DEHA2F16478g (subcluster C3 and SONS d) may be linked to a species-based clustering of DEHA genes within the subcluster C3.
In order to assess the evolution of genes using a neighbourhood analysis, we considered four more species in addition to those examined in the study by Gbelska et al. . Our classification is consistent for the 24 genes in the five species we share with Gbelska et al. Moreover, the additional species enabled us to follow more accurately the evolution of the chromosomal environments of the Pdrp and to distinguish between paralogous and orthologous status, and, in so doing, clarify the fate of some ERGO and KLLA genes. We found that some orthology relationships are confused by SYNERGY [23, 46] (see Additional file 4). Our analysis converged with the Yeast Gene Order Browser , except for the orthology relationship of the two members of the KLLA tandem gene array which are inverted in the YGOB.
We demonstrate that, combined with phylogenetic analysis, neighborhood analysis allows us to trace back the detailed evolutionary history of most members of the phylogenetic subfamilies of Pdrp sensu stricto from SACE, CAGL, ZYRO, KLTH, SAKL, KLLA and ERGO. These species span 6 out of the 14 clades from the phylogenic Saccharomyces complex  that may have separated from the other yeasts around 150 million years ago. We may therefore conclude that Pdr5p and Pdr15p issued from a single ancestral gene present before the WGD. In contrast, both Pdr10p (a Pdr5p paralog) and YNR070wp (a Snq2p paralog) originated from independent post-WGD duplication events.
We show that for a large subfamily such as Pdrp, whose members share high amino acid sequence similarity, neighborhood analysis can be used to trace the duplicates from pre-WGD or post-WGD and to describe in detail the tandem or repeat duplications and the gains, loss or transposition of gene copies.
Protein sequences from Saccharomyces cerevisiae (SACE), Candida glabrata (CAGL), Zygosaccharomyces rouxii (ZYRO), Kluyveromyces thermotolerans (KLTH), Saccharomyces kluyveri (SAKL), Kluyveromyces lactis (KLLA), Eremothecium gossypii (ERGO), Debaryomyces hansenii (DEHA) and Yarrowia lipolytica (YALI) were retrieved from the Génolevures database in July 2007. A bioinformatics analysis based on BLAST and Smith-Waterman alignments, clustering and Condorcet election, was applied to the Génolevures data set by NIKOLSKI and SHERMAN in order to define consensus families. All SACE Pdrp were found to belong to a single family, the GLC2004, recently relabeled GL3C0025. This family comprises 62 members (Table 1). Six peptides of 72 to 1099 amino acid residues were considered as fragments, leaving 56 full-sized members which were further analyzed. Genes are labeled according to the Génolevures nomenclature .
Principle of assignation to SONS
Within the Génolevures GLC3C0025 family, the Subsets of Orthologs by Neighborhood and Similarity (SONS) were identified as schematized in Additional files 2 and 3. Fifteen neighbor genes were retrieved on each side of the 62 PDR (called "query gene"). Each of these 30 genes is assigned to a Génolevures family. As family assignment is based on similarity criteria, neighbor genes belonging to the same GL family were considered orthologs by similarity. Two PDR query genes, sharing at least one neighbor being classified in the same GL family, were regarded as belonging to the same SONS. Indeed, in this case, the orthology based on similarity evidences is confirmed by the neighborhood information. Out of the 62 Pdrp, 40 show informative neighborhoods. Due to the complexity of the PDR evolutionary history, the SONS were manually curated to remove some weak links between chromosomal environments. The curated output of the SONS analysis is presented in Figure 2. For the sake of clarity, the representation is limited to five neighbor genes on each side of the query PDR genes. In fact, no conclusive neighborhood could be identified in the DEHA or YALI species. Therefore all neighbor analyses of Figure 2 excluded DEHA and YALI and were restricted to the seven species we examined from the Saccharomycetes complex defined by KURTZMAN, except one DEHA gene in SONS d.
The Pdrp sequences were aligned with ClustalW  on http://www.ebi.ac.uk/clustalw/ and visualized with Jalview . The N- Walker A and C-terminal ABC signature motifs were identified based on alignments with the motifs defined by DECOTTIGNIES and GOFFEAU for the Pdrp of S. cerevisiae. Identity was calculated after removing the filter for low complexity.
Multiple alignment was calculated by MUSCLE . The phylogenetic tree in Figure 1 was constructed using PROTDIST and NEIGHBOR in the PHYLIP suite . The use of Maximum likelihood (PROTML) or parsimony approach gave identical results. Only NEIGHBOR JOINING results are, therefore, presented. Representation of the phylogenetic tree made use of the Dendroscope application .
Rank GH, Bech-Hansen NT: Single nuclear gene inherited cross resistance and collateral sensitivity to 17 inhibitors of mitochondrial function in S. cerevisiae. Molecular and General Genetics MGG. 1973, 126: 93-102. 10.1007/BF00330986.
Carvajal E, vandenHazel HB, CybularzKolaczkowska A, Balzi E, Goffeau A: Molecular and phenotypic characterization of yeast PDR1 mutants that show hyperactive transcription of various ABC multidrug transporter genes. Molecular & General Genetics. 1997, 256: 406-415. 10.1007/s004380050584.
Delaveau T, Delahodde A, Carvajal E, Subik J, Jacq C: Pdr3, A New Yeast Regulatory Gene, Is Homologous to Pdr1 and Controls the Multidrug-Resistance Phenomenon. Molecular & General Genetics. 1994, 244: 501-511. 10.1007/BF00583901.
Subik J, Ulaszewski S, Goffeau A: Genetic-Mapping of Nuclear Mucidin Resistance Mutations in Saccharomyces cerevisiae. Current Genetics. 1986, 10: 665-670. 10.1007/BF00410914.
DeRisi J, Hazel van den B, Marc P, Balzi E, Brown P, Jacq C, Goffeau A: Genome microarray analysis of transcriptional activation in multidrug resistance yeast mutants. FEBS Letters. 2000, 470: 156-160. 10.1016/S0014-5793(00)01294-1.
Balzi E, Wang M, Leterme S, Vandyck L, Goffeau A: Pdr5, A Novel Yeast Multidrug-Resistance Conferring Transporter Controlled by the Transcription Regulator Pdr1. Journal of Biological Chemistry. 1994, 269: 2206-2214.
Bissinger PH, Kuchler K: Molecular-Cloning and Expression of the Saccharomyces cerevisiae Sts1 Gene-Product - A Yeast ABC Transporter Conferring Mycotoxin Resistance. Journal of Biological Chemistry. 1994, 269: 4180-4186.
Meyers S, Schauer W, Balzi E, Wagner M, Goffeau A, Golin J: Interaction of the Yeast Pleiotropic Drug-Resistance Gene-Pdr1 and Gene-Pdr5. Current Genetics. 1992, 21: 431-436. 10.1007/BF00351651.
Servos J, Haase E, Brendel M: Gene Snq2 of Saccharomyces cerevisiae, Which Confers Resistance to 4-Nitroquinoline-N-Oxide and Other Chemicals, Encodes A 169-Kda Protein Homologous to Atp-Dependent Permeases. Molecular & General Genetics. 1993, 236: 214-218. 10.1007/BF00277115.
Decottignies A, Grant AM, Nichols JW, de Wet H, McIntosh DB, Goffeau A: ATPase and multidrug transport activities of the overexpressed yeast ABC protein Yor1p. Journal of Biological Chemistry. 1998, 273: 12612-12622. 10.1074/jbc.273.20.12612.
Wolfger H, Mahe Y, ParleMcDermott A, Delahodde A, Kuchler K: The yeast ATP binding cassette (ABC) protein genes PDR10 and PDR15 are novel targets for the Pdr1 and Pdr3 transcriptional regulators. FEBS Letters. 1997, 418: 269-274. 10.1016/S0014-5793(97)01382-3.
Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, Galibert F, Hoheisel JD, Jacq C, Johnston M, Louis EJ, Mewes HW, Murakami Y, Philippsen P, Tettelin H, Oliver SG: Life with 6000 genes. Science. 1996, 274: 563-567. 10.1126/science.274.5287.546.
Decottignies A, Goffeau A: Complete inventory of the yeast ABC proteins. Nature Genetics. 1997, 15: 137-145. 10.1038/ng0297-137.
Chen CJ, Chin JE, Ueda K, Clark DP, Pastan I, Gottesman MM, Roninson IB: Internal Duplication and Homology with Bacterial Transport Proteins in the Mdr1 (P-Glycoprotein) Gene from Multidrug-Resistant Human-Cells. Cell. 1986, 47: 381-389. 10.1016/0092-8674(86)90595-7.
Gros P, Croop J, Housman D: Mammalian Multidrug Resistance Gene - Complete Cdna Sequence Indicates Strong Homology to Bacterial Transport Proteins. Cell. 1986, 47: 371-380. 10.1016/0092-8674(86)90594-5.
Kolaczkowski M, vanderRest M, CybularzKolaczkowska A, Soumillion JP, Konings WN, Goffeau A: Drugs, ionophoric peptides, and steroids as substrates of the yeast multidrug transporter Pdr5p. Journal of Biological Chemistry. 1996, 271: 31543-31548. 10.1074/jbc.271.49.31543.
Kolaczkowski M, Kolaczkowska A, Luczynski J, Witek S, Goffeau A: In vivo characterization of the drug resistance profile of the major ABC transporters and other components of the yeast pleiotropic drug resistance network. Microbial Drug Resistance-Mechanisms Epidemiology and Disease. 1998, 4: 143-158.
Leonard PJ, Rathod PK, Golin J: Loss of Function Mutation in the Yeast Multiple-Drug Resistance Gene Pdr5 Causes A Reduction in Chloramphenicol Efflux. Antimicrobial Agents and Chemotherapy. 1994, 38: 2492-2494.
Golin J, Ambudkar SV, May L: The yeast Pdr5p multidrug transporter: How does it recognize so many substrates?. Biochemical and Biophysical Research Communications. 2007, 356: 1-5. 10.1016/j.bbrc.2007.02.011.
Piper P, Mahe Y, Thompson S, Pandjaitan R, Holyoak C, Egner R, Muhlbauer M, Coote P, Kuchler K: The Pdr12 ABC transporter is required for the development of weak organic acid resistance in yeast. Embo Journal. 1998, 17: 4257-4265. 10.1093/emboj/17.15.4257.
Cannon RD, Lamping E, Holmes AR, Niimi K, Baret PV, Keniya MV, Tanabe K, Niimi M, Goffeau A, Monk BC: Efflux-Mediated Antifungal Drug Resistance. Clin Microbiol Rev. 2009, 22: 291-321. 10.1128/CMR.00051-08.
De Hertogh B, Hancy F, Goffeau A, Baret PV: Emergence of species-specific transporters during evolution of the Hemiascomycete phylum. Genetics. 2006, 172: 771-781. 10.1534/genetics.105.046813.
Gbelska Y, Krijger JJ, Breunig KD: Evolution of gene families: the multidrug resistance transporter genes in five related yeast species. Fems Yeast Research. 2006, 6: 345-355. 10.1111/j.1567-1364.2006.00058.x.
Gaur M, Choudhury D, Prasad R: Complete inventory of ABC proteins in human pathogenic yeast, Candida albicans. Journal of Molecular Microbiology and Biotechnology. 2005, 9: 3-15. 10.1159/000088141.
Wilcox LJ, Balderes DA, Wharton B, Tinkelenberg AH, Rao G, Sturley SL: Transcriptional profiling identifies two members of the ATP-binding cassette transporter superfamily required for sterol uptake in yeast. Journal of Biological Chemistry. 2002, 277: 32466-32472. 10.1074/jbc.M204707200.
Berbee ML, Taylor JW: Fungal molecular evolution: Gene trees and geologic time. Mycota Vii: Systematics and Evolution, Pt B. 2001, 7: 229-245.
Kurtzman CP: Phylogenetic circumscription of Saccharomyces, Kluyveromyces and other members of the Saccharomycetaceae, and the proposal of the new genera Lachancea, Nakaseomyces, Naumovia, Vanderwaltozyma and Zygotorulaspora. Fems Yeast Research. 2003, 4: 233-245. 10.1016/S1567-1356(03)00175-2.
Kurtzman CP, Robnett CJ: Phylogenetic relationships among yeasts of the 'Saccharomyces complex' determined from multigene sequence analyses. Fems Yeast Research. 2003, 3: 417-432. 10.1016/S1567-1356(03)00012-6.
Souciet JL, Dujon B, Gaillardin C, Johnston M, Baret PV, Cliften P, Sherman DJ, Weissenbach J, Westhof E, Wincker P, Jubin C, Poulain J, Barbe V, Segurence BG, Artiguenave F, Anthouard V, Vacherie B, Val ME, Fulton RS, Minx P, Wilson R, Durrens P, Jean G, Marck C, Martin T, Nikolski M, Rolland T, Seret ML, Casaregola S, Despons L, Fairhead C, Fischer G, Lafontaine I, Leh V, Lemaire M, de Montigny J, Neuveglise C, Thierry A, Blanc-Lenfle I, Bleykasten C, Diffels J, Fritsch E, Frangeul L, Goeffon A, Jauniaux N, Kachouri-Lafond R, Payen C, Potier S, Pribylova L, Ozanne C, Richard GF, Sacerdot C, Straub ML, Talla E: Comparative genomics of protoploid Saccharomycetaceae. Genome Res. 2009,
Byrne KP, Wolfe KH: The Yeast Gene Order Browser: Combining curated homology and syntenic context reveals gene fate in polyploid species. Genome Res. 2005, 15: 1456-1461. 10.1101/gr.3672305.
Kellis M, Patterson N, Endrizzi M, Birren B, Lander ES: Sequencing and comparison of yeast species to identify genes and regulatory elements. Nature. 2003, 423: 241-254. 10.1038/nature01644.
Thevenieau F, Le Dall MT, Nthangeni B, Mauersberger S, Marchal R, Nicaud JM: Characterization of Yarrowia lipolytica mutants affected in hydrophobic substrate utilization. Fungal Genetics and Biology. 2007, 44: 531-542. 10.1016/j.fgb.2006.09.001.
Nakayama H, Tanabe K, Bard M, Hodgson W, Wu S, Takemori D, Aoyama T, Kumaraswami NS, Metzler L, Takano Y, Chibana H, Niimi M: The Candida glabrata putative sterol transporter gene CgAUS1 protects cells against azoles in the presence of serum. Journal of Antimicrobial Chemotherapy. 2007, 60: 1264-1272. 10.1093/jac/dkm321.
Fairhead C, Dujon B: Structure of Kluyveromyces lactis subtelomeres: duplications and gene content. Fems Yeast Research. 2006, 6: 428-441. 10.1111/j.1567-1364.2006.00033.x.
Prasad R, Dewergifosse P, Goffeau A, Balzi E: Molecular-Cloning and Characterization of A Novel Gene of Candida albicans, Cdr1, Conferring Multiple Resistance to Drugs and Antifungals. Current Genetics. 1995, 27: 320-329. 10.1007/BF00352101.
DelSorbo G, Andrade AC, VanNistelrooy JGM, Vankan JAL, Balzi E, DeWaard MA: Multidrug resistance in Aspergillus nidulans involves novel ATP-binding cassette transporters. Molecular & General Genetics. 1997, 254: 417-426. 10.1007/s004380050434.
Iwaki T, Giga-Hama Y, Takegawa K: A survey of all 11 ABC transporters in fission yeast: two novel ABC transporters are required for red pigment accumulation in a Schizosaccharomyces pombe adenine biosynthetic mutant. Microbiology-Sgm. 2006, 152: 2309-2321. 10.1099/mic.0.28952-0.
Turi TG, Rose JK: Characterization of A Novel Schizosaccharomyces pombe Multidrug-Resistance Transporter Conferring Brefeldin-A Resistance. Biochemical and Biophysical Research Communications. 1995, 213: 410-418. 10.1006/bbrc.1995.2147.
Nakaune R, Adachi K, Nawata O, Tomiyama M, Akutsu K, Hibi T: A novel ATP-binding cassette transporter involved in multidrug resistance in the phytopathogenic fungus Penicillium digitatum. Applied and Environmental Microbiology. 1998, 64: 3983-3988.
Verrier PJ, Bird D, Buria B, Dassa E, Forestier C, Geisler M, Klein M, Kolukisaoglu U, Lee Y, Martinoia E, Murphy A, Rea PA, Samuels L, Schulz B, Spalding EJ, Yazaki K, Theodoulou FL: Plant ABC proteins - a unified nomenclature and updated inventory. Trends in Plant Science. 2008, 13: 151-159. 10.1016/j.tplants.2008.02.001.
Sanchez-Fernandez R, Davies TGE, Coleman JOD, Rea PA: The Arabidopsis thaliana ABC protein superfamily, a complete inventory. Journal of Biological Chemistry. 2001, 276: 30231-30244. 10.1074/jbc.M103104200.
Crouzet J, Trombik T, Fraysse AS, Boutry M: Organization and function of the plant pleiotropic drug resistance ABC transporter family. FEBS Letters. 2006, 580: 1123-1130. 10.1016/j.febslet.2005.12.043.
Ullman B: Multidrug-Resistance and P-Glycoproteins in Parasitic Protozoa. Journal of Bioenergetics and Biomembranes. 1995, 27: 77-84. 10.1007/BF02110334.
Dean M, Hamon Y, Chimini G: The human ATP-binding cassette (ABC) transporter superfamily. Journal of Lipid Research. 2001, 42: 1007-1017.
Merico A, Sulo P, Piskur J, Compagno C: Fermentative lifestyle in yeasts belonging to the Saccharomyces complex. Febs Journal. 2007, 274: 976-989. 10.1111/j.1742-4658.2007.05645.x.
Wapinski I, Pfeffer A, Friedman N, Regev A: Automatic genome-wide reconstruction of phylogenetic gene trees. Bioinformatics. 2007, 23: I549-I558. 10.1093/bioinformatics/btm193.
Nikolski M, Sherman DJ: Family relationships: should consensus reign?--consensus clustering for protein families. Bioinformatics. 2007, 23: e71-e76. 10.1093/bioinformatics/btl314.
Durrens P, Sherman DJ: A systematic nomenclature of chromosomal elements for hemiascomycete yeasts. Yeast. 2005, 22: 337-342. 10.1002/yea.1214.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucl Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.
Clamp M, Cuff J, Searle SM, Barton GJ: The Jalview Java alignment editor. Bioinformatics. 2004, 20: 426-427. 10.1093/bioinformatics/btg430.
Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. Bmc Bioinformatics. 2004, 5: 1-19. 10.1186/1471-2105-5-113.
Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.
Huson DH, Richter DC, Rausch C, Dezulian T, Franz M, Rupp R: Dendroscope: An interactive viewer for large phylogenetic trees. Bmc Bioinformatics. 2007, 8: 460-10.1186/1471-2105-8-460.
Dujon B, Sherman D, Fischer G, Durrens P, Casaregola S, Lafontaine I, de Montigny J, Marck C, Neuveglise C, Talla E, Goffard N, Frangeul L, Aigle M, Anthouard V, Babour A, Barbe V, Barnay S, Blanchin S, Beckerich JM, Beyne E, Bleykasten C, Boisrame A, Boyer J, Cattolico L, Confanioleri F, de Daruvar A, Despons L, Fabre E, Fairhead C, Ferry-Dumazet H, Groppi A, Hantraye F, Hennequin C, Jauniaux N, Joyet P, Kachouri R, Kerrest A, Koszul R, Lemaire M, Lesur I, Ma L, Muller H, Nicaud JM, Nikolski M, Oztas S, Ozier-Kalogeropoulos O, Pellenz S, Potier S, Richard GF, Straub ML, Suleau A, Swennen D, Tekaia F, Wesolowski-Louvel M, Westhof E, Wirth B, Zeniou-Meyer M, Zivanovic I, Bolotin-Fukuhara M, Thierry A, Bouchier C, Caudron B, Scarpelli C, Gaillardin C, Weissenbach J, Wincker P, Souciet JL: Genome evolution in yeasts. Nature. 2004, 430: 35-44. 10.1038/nature02579.
Dietrich FS, Voegeli S, Brachat S, Lerch A, Gates K, Steiner S, Mohr C, Pohlmann R, Luedi P, Choi S, Wing RA, Flavier A, Gaffney TD, Philippsen P: The Ashbya gossypii Genome as a Tool for Mapping the Ancient Saccharomyces cerevisiae Genome. Science. 2004, 304: 304-307. 10.1126/science.1095781.
Gattiker A, Rischatsch R, Demougin P, Voegeli S, Dietrich FS, Philippsen P, Primig M: Ashbya Genome Database 3.0: a cross-species genome and transcriptome browser for yeast biologists. Bmc Genomics. 2007, 8: 9-10.1186/1471-2164-8-9.
We would like to thank the Génolevures consortium, coordinated by Jean-Luc SOUCIET, for helpful discussions and access to the database of protein sequences from Zygosaccharomyces rouxii, Kluyveromyces thermotolerans, Saccharomyces kluyveri and Eremothecium gossypii. We acknowledge the constructive comments of two anonymous reviewers.
M-LS and JD designed and developed a computer program to identify SONS. M-LS did the manual curation of SONS and the preparation of figures. JD prepared the tables. AG identified the PDR motifs and PB did the phylogenetic study. M-LS and JD drafted manuscripts with revisions from AG and PB. All authors read and approved the final manuscript edited by M-LS.
Electronic supplementary material
Additional file 1:List of Pdrp features included in this paper. Sequences of S. cerevisiae, C. glabrata, K. lactis, D. hansenii, Y. lipolytica and E. gossypii were previously published. Z. rouxii and K. thermotolerans have both been sequenced by Génoscope, and S. kluyveri was sequenced by the Washington University Genome Sequencing Centre. The three new genomes ZYRO, KLTH, SAKL, plus KLLA have been entirely annotated and manually curated by the Génolevures consortium . ERGO annotation is by GATTIKER. (PDF 48 KB)
Additional file 2:Schema of analysis tools and databases. Rectangles represent databases and results while parallelograms represent tools used to obtain them. (PPT 103 KB)
Additional file 3:Steps of the Identification of Orthologs by Neighborhood and Similarity method. (1) The building of SONS (Subset of Orthologs defined by Neighborhood and Similarity) begins by taking a family of genes whose translation products show sequence similarity for example, the family GL3C0025. (2) For each query gene of this family we identify 15 genes, which we call neighbors, on each side and especially the family to which these neighbors belong. (3) Two query genes whose translation products belong to the same family (homologs by similarity) are members of the same SONS if they share at least one pair of neighbors that are also homologous to each other by similarity. These neighbors are highlighted in the same color. In the example, YOL075c shares some neighbors with the query gene of CAGL and has one neighbor in common with the query gene of KLTH. CAGL has no gene in common with KLTH, however YOL075c must be in the same SONS as the query genes of CAGL and KLTH. In consequence, the three query genes are clustered in the same SONS. Homologs that do not share a pair of homologous neighbors are separated into two distinct SONS. (4) See Figure 2 SONS d for the complete SONS of YOL075c. (PPT 92 KB)
Additional file 4:Relationships between the Pdrp sensu stricto based on the neighborhood analysis: comparison with YGOB and SYNERGY. Each box represents a PDR. Grey lines join PDR genes located on the same chromosome. When non-adjacent, the number of coding sequences between two PDR genes appears in white. Shared neighborhoods are represented by colored connectors. A large connector links PDR genes sharing more than one common neighbor while a thin connector links PDR genes sharing one neighbor only. The numbers shown in blue is the number of the pillar in which the PDR gene has been classified on the Yeast Gene Order Browser . The number shown in red is the number of the orthogroup in which the PDR gene has been classified by Synergy . (PDF 326 KB)
About this article
Cite this article
Seret, ML., Diffels, J.F., Goffeau, A. et al. Combined phylogeny and neighborhood analysis of the evolution of the ABC transporters conferring multiple drug resistance in hemiascomycete yeasts. BMC Genomics 10, 459 (2009). https://doi.org/10.1186/1471-2164-10-459