<title>Background</title> <p>Lateral or horizontal gene transfer (HGT) is defined as the exchange of genes between different strains or species [##REF##10611671##1##]. HGT introduces new genes into a recipient genome that are either homologous to existing genes, or belong to entirely new sequence families. Large-scale genomic sequencing of prokaryotes has revealed that gene transfer is an important evolutionary mechanism for these organisms [##REF##12777514##2##,##REF##11050445##3##]. HGT has been linked to the acquisition of drug resistance by benign bacteria [##REF##13679020##4##], and also to the gain of genes that confer the ability to catabolize certain amino acids that are important virulence factors [##REF##9422600##5##]. However there is much debate as to whether lateral gene transfer is an ubiquitous influence throughout prokaryotic genome evolution [##REF##12902542##6##]. Until recently, the process of gene transfer has been assumed to be of limited significance to eukaryotes [##REF##15761667##7##]. The availability of diverse eukaryotic genome sequence data is dramatically changing our views on the important role gene transfer can play in eukaryotic evolution.</p> <p>The rapid increase in fungal sequence data has promoted this kingdom to the forefront of comparative genomics [##REF##16188435##8##]. Whereas there is some documented evidence for HGT between fungal species [##REF##16832356##9##, ####REF##16055562##10##, ##REF##16870684##11##, ##REF##18218086##12##, ##REF##16367844##13##, ##REF##17513882##14##, ##REF##17971860##15##, ##REF##16675503##16##, ##REF##17652520##17####17652520##17##] or from bacteria to fungi [##REF##12546782##18##, ####REF##15947202##19##, ##REF##18073433##20##, ##REF##17277171##21##, ##REF##15229592##22##, ##REF##15014982##23##, ##REF##11447209##24##, ##REF##15118835##25##, ##REF##15483318##26##, ##REF##10723736##27##, ##REF##9287428##28####9287428##28##] [see additional file ##SUPPL##0##1##], overall very few incidences have been identified. There are two possible explanations: either gene transfer is indeed extremely rare amongst fungi, or it has not yet been thoroughly studied. To address this question we investigated the frequency of successful recent interdomain HGT events between prokaryotes and yeast species belonging to the CTG clade. We chose this course of action as we expect recent interdomain HGT events to be more readily identified and supported than more ancient transfers.</p> <p>For the purposes of this study, we define CTG species as the immediate relatives of <italic>C. albicans</italic>, including <italic>C. tropicalis, C. parapsilosis, Clavispora lusitaniae, Pichia guilliermondii</italic>, and <italic>Lodderomyces elongisporus</italic>. These species have been completely sequenced, share a relatively recent common ancestor [##REF##17121679##29##], and the codon CUG is translated as serine rather than leucine [##REF##10188281##30##].</p> <p>We used syntenic, phylogenetic and sequence based analyses to identify two cases of interdomain HGT between prokaryotes and <italic>C. parapsilosis</italic>, most likely involving the proteobacteria phylum. Our results suggest that extant CTG species do not readily take up exogenous DNA.</p>
<title>Methods</title> <title>Sequence data</title> <p>The complete <italic>C. albicans </italic>(SC5314) genome (Assembly 19) was obtained from the <italic>Candida </italic>genome database [##UREF##4##65##]. The Broad institute have sequenced and annotated five CTG species (<italic>C. albicans </italic>(WO-1), <italic>C. tropicalis</italic>, <italic>L. elongisporus</italic>, <italic>P. guilliermondii</italic>, and <italic>Cl. lusitaniae</italic>). These genomes were obtained directly from the Broad Institute [##UREF##5##66##]. Gene sets for the <italic>C. dubliniensis </italic>were downloaded from GeneDB [##UREF##1##44##].</p> <p>The incomplete <italic>C. parapsilosis </italic>geneome was downloaded from the Sanger Institute [##UREF##6##67##]. Gene annotations were performed using two separate approaches. The first involved a reciprocal best BLAST [##REF##9254694##31##] search with a cutoff E- value of 10<sup>-7 </sup>of <italic>Candida albicans </italic>SC5314 protein coding genes against the unannotated <italic>C. parapsilosis </italic>genome. Top BLAST hits longer than 300 nucleotides were retained as putative open reading frames. The second approach involved a pipeline of analysis that combined several different gene prediction programs, including <italic>ab initio </italic>programs SNAP [##REF##15144565##68##], Genezilla [##REF##15145805##69##], and AUGUSTUS [##REF##15145805##69##], with gene models from Exonerate [##REF##15713233##70##] and Genewise [##REF##15123596##71##] based on alignments of proteins and Expressed Sequence Tags. Putative gene sets from both approaches were imported into Artemis [##REF##11120685##72##] and cross corroborated manually. The resultant gene sets contained 5,823 protein-coding genes. The <italic>C. parapsilosis </italic>genome was also annotated by the Broad Institute, and where possible we have used the gene names they assigned.</p> <p>The UniProt database (v11.1) was downloaded [##UREF##7##73##]. Database searches against GenBank refer to release 164.0.</p> <title>Blast based approach to detect potential horizontally transferred genes</title> <p>Taking one CTG species at a time, we located gene families of interest by comparing individual protein coding genes against the UniProt database (v11.1) using the BlastP algorithm [##REF##9254694##31##] with a cutoff expectation (E) value of 10<sup>-20</sup>. To use all available sequence data, CTG proteins with a top database hit to a bacterial protein in UniProt were extracted for a second round of database searching against GenBank (E value of 10<sup>-20</sup>). Proteins which also had a top database hit to a bacterial protein in GenBank were considered as possible incidences of horizontal gene transfer. All putative homologs were extracted from GenBank and searched against the relevant CTG genome to ensure a reciprocal best Blast hit. For completeness, CTG proteins not yet deposited in GenBank were added to gene families of interest where appropriate.</p> <p>Accession numbers for all sequences used in this analysis can be found in additional material [see additional file ##SUPPL##1##2##].</p> <title>Phylogenetic methods</title> <p>Gene families were aligned using MUSCLE (v3.6) [##REF##15034147##74##] using the default settings. Obvious alignment ambiguities were corrected manually.</p> <p>Phylogenetic relationships were inferred using maximum likelihood methods. Appropriate protein models of substitution were selected for each gene family using ModelGenerator [##REF##16563161##75##]. One hundred bootstrap replicates were then carried out with the appropriate protein model using the software program PHYML [##REF##14530136##76##] and summarized using the majority-rule consensus method.</p> <p>We performed the approximately unbiased test of phylogenetic tree selection [##REF##12079646##53##], to assess whether differences in topology between constrained and unconstrained gene trees are no greater than expected by chance.</p> <title>Codon usage analysis and spectral analysis</title> <p>To determine if the putative HGT genes had a different codon usage pattern to the host genome an analysis of variation in synonymous codon usage was undertaken using the GCUA software [##REF##9632833##77##]. Individual correspondence analyses of raw codon counts for the <italic>Candida parapsilosis</italic>, <italic>Ustilago maydis, Malassezia globosa, Aspergillus flavus, Aspergillus niger, Gibberella zeae, Aspergillus oryzae, Phaeosphaeria nodorum</italic>, and <italic>Schizosaccharomyces pombe </italic>genomes were performed, with the first four principal axes being used to evaluate synonymous codon usage patterns. Similar analyses were also carried out on members of the proline racemase and phenazine F gene families displayed in Figures ##FIG##1##2## and ##FIG##3##4##. We used spectrum [##UREF##8##78##] to perform a spectral analysis on a subset of the phenazine data.</p>
<title>Results and discussion</title> <title>Identification of horizontal gene transfer candidates through Blast database search</title> <p>We compared all available CTG gene sets against UniProt using BlastP [##REF##9254694##31##]. CTG genes with top database hits to bacterial species were identified as putative horizontally transferred genes and the resultant Blast files were inspected manually. A <italic>D. hansenii </italic>gene (protein ydhR precursor) with a top database hit to a bacterial sequence was not considered for further analsyes as it has previously been described [##REF##15229592##22##]. After this process two genes from <italic>C. parapsilosis </italic>were considered for further analysis; one encodes a putative proline racemase, and the second encodes a member of the phenazine F superfamily. Related family members were identified by a second round of database searching against GenBank to ensure all available genomic data was utilized.</p> <title>Proline racemase phylogeny and characterization</title> <p>The <italic>C. parapsilosis </italic>gene (designated CPAG_02038) is most similar to a proline racemase homolog from <italic>Burkholderia cenocepacia </italic>AU 1054 protein (66% pairwise identity; Figure ##FIG##0##1A##). Amino acid racemases catalyze the interconversion of L- and D-amino acids by abstraction of the α-amino proton of the enzyme bound substrate [##REF##5722267##32##]. CPAG_02038 lies within a large contig and is also present in a previously published genome survey of <italic>C. parapsilosis </italic>[##REF##15947193##33##], suggesting its presence does not the result from contamination. We could not locate any related genes in any other CTG genome (using BlastP or TBlastN). Family members are widely distributed throughout the prokaryotes however, and are also located within the Pezizomycotina.</p> <p>We extracted 321 putative proline racemases from 207 organisms, including members of the α, β, γ, and δ-proteobacteria, Actinobacteria, Fungi, Protozoa and Metazoa. Numerous species were found to have several family members [see additional file ##SUPPL##1##2##]; all were included for complete comparative purposes. A maximum likelihood (ML) phylogeny was reconstructed from an alignment of all the PR proteins (Figure ##FIG##1##2##).</p> <p>There are a large number of polytomies displayed in Figure ##FIG##1##2##. These probably result from duplication of PR genes followed by diversifying selection, leading to a high degree of sequence heterogeneity. For example, <italic>Agrobacterium tumefaciens </italic>str. C58 contains three PR homologs [see additional file ##SUPPL##1##2##], with an average amino acid pairwise percentage identity of ~31%. <italic>Burkholderia cenocepacia </italic>AU 1054 contains 2 proline racemase homologs [see additional file ##SUPPL##1##2##], which are only 28% identical. To help resolve the evolutionary history amongst PR homologs we reconstructed an additional ML phylogeny based on a reduced dataset (Figure ##FIG##2##3##). We also reconstructed a Bayesian phylogeny using the heterogeneous CAT site model. The CAT model can account for site-specific features of sequence evolution and has been found to be more robust than other methods against phylogenetic artifacts such as long branch attraction [##REF##17288577##34##]. The resultant Bayesian phylogeny is highly congruent with the ML phylogeny (not shown).</p> <p>The putative <italic>C. parapsilosis </italic>PR homolog lies in a strongly supported (100% Bootstrap support (BP)) clade with <italic>Burkholderia </italic>species (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-A). <italic>Burkholderia </italic>are β-proteobacteria. However, no other β-proteobacteria, or indeed any other bacterial genus were found within clade-A (Figures ##FIG##1##2## &amp;##FIG##2##3##).</p> <p>Although no PR homologs were identified in other CTG species, or indeed in any other of the Saccharomycotina, there are homologs in family members of the Pezizomycotina. A Pezizomycotina specific subclade is evident in our phylogeny containing <italic>Phaeosphaeria nodorum</italic>, <italic>Aspergillus niger </italic>and <italic>Gibberella zeae </italic>(Figures ##FIG##1##2## &amp;##FIG##2##3## clade-B 100% BP). This subclade is found in a strongly supported clade with members of the Actinobacteria (Figure ##FIG##1##2## 100% BP), containing <italic>Brevibacterium linens </italic>and an unclassified marine actinobacterium and excluding <italic>Rubrobacter xylanophilus </italic>(Figure ##FIG##1##2## 87% BP). This suggests that these Pezizomycotina species obtained their PR gene from the Actinobacteridae subclass rather than the Rubrobacteridae subclass. This transfer event is another independent HGT event of a PR gene into fungi, and we hypothesize it occurred early in the Pezizomycotina lineage, as it is shared by three distantly related species. Its patchy phyletic distribution suggests it has been subsequently lost in other Pezizomycotina species.</p> <p>There are also PR homologs in the Metazoans. These are found in a eukaryote clade that also contains a number of Pezizomycotina representatives (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-C 93% BP). Several scenarios can explain this phylogenetic positioning. Firstly, the PR gene may have been present in the last universal common ancestor of all eukaryotes but has been differentially lost in all lineages except those leading to modern day Metazoa and Pezizomycotina. Alternatively, an ancient gene transfer from bacteria to the last common ancestor (LCA) of Metazoa and Fungi could have occurred, with subsequent gene loss amongst different Metazoan and Fungal lineages. A third hypothesis is that two independent gene transfers have occurred into the Metazoan and Pezizomycotina lineages from unsampled bacterial donors. Finally, a transfer from unsampled bacteria into one of the eukaryote clades (either Metazoa or Pezizomycotina) may have occurred with subsequent transfer from one eukaryotic group to the other.</p> <p><italic>A. niger</italic>, <italic>A. oryzae </italic>and <italic>G. zeae </italic>all contain multiple PR homologs [see additional file ##SUPPL##1##2##]. One <italic>A. niger</italic>, one <italic>G. zeae </italic>and the three <italic>A. oryzae </italic>PR homologs are nested in a strongly supported Pezizomycotina specific subclade (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-D 100% BP). This subclade if found within a larger predominately proteobacterial clade (Figure ##FIG##1##2## 74% BP). This infers that there was an independent gene transfer event of a bacterial PR homolog into an ancestral Pezizomycotina species.</p> <p>The phylogenetic position of the <italic>C. parapsilosis </italic>PR homolog (Figures ##FIG##1##2## &amp;##FIG##2##3##) resemble that described for the adenosine deaminase (ADA) gene in the <italic>Dekkera bruxellensis </italic>genome [##REF##17277171##21##]. In that analysis, the authors suggest that <italic>D. bruxellensis </italic>and <italic>Burkholderia </italic>species received the ADA gene from a species not yet represented in the public sequence databases. Our PR phylogeny suggests a similar event may have occurred within clade-A, which contains only <italic>C. parapsilosis </italic>and <italic>Burkholderia </italic>species (Figures ##FIG##1##2## &amp;##FIG##2##3##). <italic>Burkholderia </italic>species are known to have a genomic repertoire that allows the transfer and receipt of exogenous DNA [##REF##12748267##35##] and a number of studies have reported successful gene transfers into <italic>Burkholderia </italic>species [##REF##15493822##36##,##UREF##0##37##]. It is possible therefore that there have been other successful gene transfers into this bacterial lineage.</p> <p>The vast majority of amino acids found in living cells correspond to the L-stereoisomer [##REF##16446443##38##]. However, D-amino acids are long known to be found in the cell walls of Gram positive and negative bacteria, where they are essential components of peptidoglycan [##REF##8749373##39##]. Apart from low levels of D-amino acids derived from spontaneous racemization as a result of aging [##REF##9949871##40##], it was assumed that only L-amino acid enantiomers were present in eukaryotes [##REF##12735293##41##]. However, recent studies have reported the presence of numerous D-amino acids in an array of organisms, including mammals [##REF##10557334##42##]. The first eukaryotic (proline) amino acid racemase has recently been described from the human pathogen <italic>Trypanosoma cruzi </italic>[##REF##10932226##43##]. A high degree of sequence similarity was observed between the <italic>T. cruzi </italic>and bacterial homologs [##REF##10932226##43##]. Our phylogeny infers that <italic>T. cruzi </italic>obtained its PR homolog through interdomain HGT from a member of the Firmicutes (subclass Clostridia), as it is grouped beside members of this group with a high degree of support (Figure ##FIG##1##2## clade-E 96% BP). We performed database searches [##UREF##1##44##], against other Protozoan genomes including <italic>Trypanosoma brucei</italic>, <italic>Trypanosoma congolense </italic>and <italic>Trypanosoma annulata</italic>. We failed to locate a homolog in all species except for <italic>T. vivax</italic>.</p> <p>Previous analysis has shown that <italic>T. cruzi </italic>and <italic>T. vivax </italic>are not each others closest phylogenetic neighbors, relative to the other species sampled [##REF##10633190##45##]. This suggests an ancestral Trypanosoma gained the PR gene and multiple losses in different Trypanosoma lineages has subsequently occurred.</p> <title>Gene order around PR homologs</title> <p>The <italic>C. parapsilosis </italic>PR homolog lies close to an ortholog (CPAG_02041) of orf19.1135 from <italic>C. albicans </italic>(Figure ##FIG##3##4##). The gene order to the left of this ORF is conserved in all CTG species, the order to the right is conserved in most CTG species apart from <italic>C. parapsilosis </italic>and <italic>L. elongisporus</italic>. <italic>C. parapsilosis </italic>and <italic>L. elongisporus </italic>are closely related [##REF##17121679##29##], and an examination of synteny suggests that the PR gene (together with a second ORF, cpar5437) were inserted between CPAG_2041 and CPAG_2037 (Figure ##FIG##3##4##). cpar5437 encodes a neutral amino acid (AA) transporter. The presence of an AA transporter beside the PR homolog is interesting. If the putative proline racemase has a role in amino acid metabolism, then the presence of the transporter may be the result of an adaptive translocation to enhance the activity of the PR gene. Unlike the PR ORF the AA transporter is fungal in origin. Most CTG species contain a single neutral AA transporter; however <italic>C. parapsilosis </italic>and <italic>D. hansenii </italic>have four.</p> <p>We located tRNA genes for nearly all CTG species beside the large conserved syntenic block (Figure ##FIG##3##4##). It has been shown that tRNA genes are associated with genomic breakpoints [##REF##10839539##46##]. We hypothesize that a genomic rearrangement has occurred at this site in the LCA of <italic>C. parapsilosis </italic>and <italic>L. elongisporus</italic>. We cannot determine if the bacterial PR homolog was inserted into the LCA of <italic>L. elongisporus</italic>/<italic>C. parapsilosis </italic>and subsequently lost in <italic>L. elongisporus</italic>, or gained by <italic>C. parapsilosis </italic>after speciation.</p> <p>We also investigated the gene order around the Pezizomycotina PR homologs [see additional file ##SUPPL##2##3##]. Gene synteny around the PR homologs found in clade-D (Figures ##FIG##1##2## &amp;##FIG##2##3##) is not conserved (not shown). Interestingly however, both <italic>A. niger </italic>and <italic>G. zeae </italic>in clade-D (Figures ##FIG##1##2## &amp;##FIG##2##3##) have genes containing a FAD dependent oxidoreductase domain in close proximity to their PR homologs (not shown). According to Pfam [##UREF##2##47##], FAD dependent oxidases include D-amino acid oxidases, that catalyze the oxidation of neutral and basic D-amino acids into their corresponding keto acids. The presence of these oxidases may be another example of an adaptive translocation to enhance the activity of the PR gene in these Pezizomycotina species.</p> <p><italic>A. oryzae </italic>has three PR homologs (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-C). All of these have orthologs in its close relative <italic>A. flavus </italic>(Figure ##FIG##2##3## clade-C), and synteny around these is conserved [see additional file ##SUPPL##2##3## clade-D]. The remaining two species in clade-C are <italic>A. niger </italic>and <italic>G. zeae</italic>. There is no evidence of conserved gene order within these species, or with <italic>A. oryzae </italic>or <italic>A. flavus</italic>. Gene order around the <italic>A. flavus </italic>and <italic>A. terreus </italic>PR homologs found in the Metazoan/Pezizomycotina clade (Figures ##FIG##1##2## &amp;##FIG##2##3##) is also conserved [see additional file ##SUPPL##2##3##], as is the order between <italic>A. fumigatus </italic>and <italic>N. fishceri </italic>[see additional file ##SUPPL##2##3##]. We could not locate amino acid transporters or FAD dependent oxidases beside any of the PR homologs found in clades B or C (Figure ##FIG##1##2##).</p> <title>Proline racemase codon usage</title> <p>It has been shown that recently acquired genes often display an atypical codon preference when compared to other genes in the genome [##REF##1762151##48##,##REF##10830951##49##]. However, the transferred PR homologs have a codon usage consistent with the rest of their genomes [see additional file ##SUPPL##3##4##]. We undertook an analysis of variation in synonymous codon usage on all PR genes shown in Figure ##FIG##1##2##. Homologs from related species cluster together [see additional file ##SUPPL##4##5##]. For example, the Actinobactria, the Firmicutes and the <italic>Burkholderia </italic>species all inhabit unique areas in two dimensional correspondence analysis space [see additional file ##SUPPL##4##5##].</p> <p>The majority of fungal and Metazoan PRs are clustered together [see additional file ##SUPPL##4##5##]. The <italic>C. parapsilosis </italic>PR homolog has a codon usage distinct from the other Pezizomycotina fungal PR homologs [see additional file ##SUPPL##4##5##], which is unsurprising as <italic>C. parapsilosis </italic>belongs to the Saccharomycotina subphylum. The <italic>C. parapsilosis </italic>homolog is also separate from the <italic>Burkholderia </italic>(β-proteobacteria) genes with which it forms a closely related phylogenetic group (Figures ##FIG##1##2## &amp;##FIG##2##3##). This suggests that the gene may have originated from a genome with no other close relatives among the species analyzed here.</p> <title>Proline racemase activity</title> <p>The PR active site from <italic>Trypanosoma cruzi, Clostridium sticklandii, Agrobacterium tumefaciens</italic>, <italic>Brucella melitensis </italic>and <italic>Pseudomonas aeruginosa </italic>all contain cysteine at amino acid position 330 [##REF##10932226##43##,##REF##1174513##50##]. This amino acid is essential for enzymatic function, because substitution with serine abolishes activity [##REF##12735293##41##]. However, PR homologs from human, mouse, <italic>Rhizobium </italic>and <italic>Brucella </italic>contain a threonine instead of a cysteine at position 330 [##REF##12735293##41##]. We observed that cysteine is found in the equivalent position in many of the bacterial proteins. The Pezizomycotina PR genes found in clade-B and clade-D contain a cysteine at the active site (Figure ##FIG##2##3##). The PR homologs found in the Metazoan/Pezizomycotina clade (clade-B) have a threonine at position 330. Similarly, the <italic>C. parapsilosis </italic>PR homolog, together with its relatives from <italic>Burkholderia </italic>all contain a threonine (Figure ##FIG##2##3##). However, <italic>Burkholderia </italic>species have multiple PR homologs [see additional file ##SUPPL##1##2##] with a cysteine as the active site (not shown). It is not clear what effect the substitution has on enzyme activity. It has been suggested that homologs containing threonine at the active site are not true PRs [##REF##12735293##41##], but may instead belong to a superfamily. We cannot detect any difference in the ability of <italic>C. parapsilosis</italic>, the other CTG species or any of the Pezizomycotina species to utilize D-proline as growth media (data not shown). We therefore cannot confidently infer the function of the PR homologs in the fungi analyzed here.</p> <title>Phenazine F phylogeny and characterisation</title> <p>The <italic>C. parapsilosis </italic>gene (designated CPAG_03462) is most similar to a <italic>Photorhabdus luminescens </italic>phenazine F (PhzF) protein with 61% pairwise identity (Figure ##FIG##0##1B##). Phenazines are biologically active compounds, all of which have a characteristic tricyclic ring system and have been shown to confer a selective growth advantage to organisms which secrete them, as they possess broad-spectrum antibiotic activity towards bacteria, fungi and higher eukaryotes [##REF##15545603##51##]. In <italic>Pseudomonas</italic>, the best studied phenazine producer, PhzF is part of an operon required for the conversion of chorismic acid to phenazine-1-carboxylate (PCA) [##REF##15449932##52##]. PhzF homologs were identified in most of the CTG species tested as well as several other fungal species. However, we could not identify a PhzF homolog in the <italic>L. elongisporus </italic>genome, even when multiple TBlastN and BlastN searches were used.</p> <p>PhzF homologs were extracted from GenBank for subsequent phylogenetic analysis. In total 181 representative protein coding sequences distributed amongst 154 organisms were used. These taxa were distributed amongst α, β, γ and δ-proteobacteria, Actinobacteria, Fungi, Firmicutes a well as other bacterial groups.</p> <p>We aligned all sequences and reconstructed a PhzF ML phylogeny (Figure ##FIG##4##5##). The <italic>C. parapsilosis </italic>PhzF homolog is found in a clade with members of the β-proteobacteria <italic>(Burkholderia multiovorans, Burkholderia cepacia, Burkholderia ambifaria</italic>), α-proteobacteria (<italic>Roseovarius</italic>) and the γ-proteobacteria (<italic>Azotobacter vinelandii, Acinetobacter baumannii, Shewanella baltica and Photorhabdus luminescens</italic>) (81% BP). In contrast, all other PhzF homologs from CTG species are in a completely separate clade (Figure ##FIG##4##5##). These form a sister group (63% BP) to PhzF homologs from other Saccharomycotina species (<italic>C. glabrata, Saccharomyces cerevisiae, Kluyveromyces lactis </italic>and <italic>Vanderwaltozyma polyspora</italic>). All three clades are grouped together in a larger clade with high support (75% BP).</p> <p>The sister group relationship between the PhzF homologs from the Ascomycota and the proteobacteria clade is intriguing (Figure ##FIG##4##5##), as it suggests that an ancestral Saccharomycotina species gained the PhzF homolog from a proteobacteria. The bacterial PhzF gene has subsequently been retained after multiple speciation events, but lost in <italic>C. parapsilosis</italic>. We hypothesize that <italic>C. parapsilosis </italic>has recently reacquired a bacterial PhzF homolog from a proteobacterial source, as it is grouped (81% BP) within a proteobacterial subclade. To test this hypothesis we reconstructed constrained trees that placed <italic>C. parapsilosis </italic>together with the remaining Ascomycota species [see additional file ##SUPPL##5##6## C-H]. The AU test of phylogenetic tree selection [##REF##12079646##53##], showed that the original unconstrained tree (groups <italic>C. parapsilosis </italic>with proteobacteria) receives the optimal likelihood tree score, and the differences in likelihood scores when compared to the constrained trees [see additional file ##SUPPL##5##6##], are significant (P &lt; 0.05). This is also supported by spectral analysis [see additional file ##SUPPL##6##7##].</p> <p>Our phylogeny shows that the <italic>Schizosaccharomyces pombe </italic>PhzF homolog is found in a clade containing all CTG PhzF homologs (Figure ##FIG##4##5## 99% BP). Furthermore it is grouped beside <italic>D. hansenii </italic>(66% BP). <italic>S. pombe </italic>is not a member of the Saccharomycotina, it belongs to the Taphrinomycotina subphylum. The genome sequences of <italic>Schizosaccharomyces japonicus </italic>and <italic>Schizosaccharomyces octosporus </italic>have recently been completed [##UREF##3##54##]. We could not locate a PhzF homolog in <italic>S. japonicus </italic>but did locate a homolog in <italic>S. octosporus </italic>using a TBlastN search strategy. Phylogenetic analysis has shown that <italic>S. pombe </italic>and <italic>S. octosporus </italic>are more closely related to one another than to <italic>S. japonicus </italic>[##REF##12527786##55##]. Therefore we hypothesize that the LCA ancestor of <italic>S. pombe </italic>and <italic>S. octosporus </italic>gained the PhzF gene from an ancestral <italic>D. hansenii</italic>-like species after speciation from <italic>S. japonicus</italic>. We reconstructed a constrained tree that placed <italic>S. pombe </italic>outside the Saccharomycotina clade [see additional file ##SUPPL##5##6B##]. The approximately unbiased test of phylogenetic tree selection (AU test) [##REF##12079646##53##], showed that the phylogenetic inferences of the unconstrained tree are significantly better (P &lt; 0.05) than the constrained tree [see additional file ##SUPPL##5##6##]. This infers that <italic>S. pombe </italic>has obtained a PhzF homolog from a member of the CTG clade.</p> <p>A small basidiomycete clade is evident amongst prokaryote species (Figure ##FIG##4##5##). Both <italic>Ustilago maydis </italic>and <italic>Malassezia globosa </italic>belong to the Ustilaginomycotina subphylum. Therefore our phylogeny infers that an ancestral Ustilaginomycotina species gained a PhzF gene from an unknown bacterial source, and both species have retained this after speciation.</p> <p>A correspondence analysis of synonymous codon usage for all PhzF homologs was also performed and is shown in additional information [see additional file ##SUPPL##7##8##]. The <italic>S. pombe </italic>PhzF homolog has a codon usage pattern very similar to the <italic>D. hansenii </italic>protein.</p> <title>Gene order around PhzF</title> <p>Analysis of the genes adjacent to the PhzF homolog in <italic>C. parapsilosis </italic>shows that there is a high conservation of gene synteny and supports our hypothesis that PhzF was recently acquired in this species (Figure ##FIG##5##6##). Homologs in the other CTG species are located in completely different regions of the genome relative to <italic>C. parapsilosis </italic>(not shown). For example, the <italic>C. albicans </italic>PhzF homolog is located between orf19.5619 and orf19.5621, whereas the <italic>C. parapsilosis </italic>homolog is found between orf19.6689 &amp; orf19.6687 relative to <italic>C. albicans </italic>SC5314 (Figure ##FIG##5##6##). However, the <italic>L. elongisporus </italic>genome contains no PhzF homolog, either at a position equivalent to the <italic>C. parapsilosis </italic>copy or elsewhere in the genome.</p> <p>We propose that the LCA of <italic>L. elongisporus </italic>and <italic>C. parapsilosis </italic>lost the PhzF gene present in the other CTG species, and a second (new) copy was subsequently gained by <italic>C. parapsilosis </italic>after speciation. We have partial sequence data (unpublished) from <italic>Candida orthopsilosis</italic>, a species so closely related to <italic>C. parapsilosis </italic>that it was once designated <italic>C. parapsilosis </italic>group II [##REF##7665651##56##]. We located a <italic>C. orthopsilosis </italic>PR homolog that is 83% identical (at the amino acid level) to the <italic>C. parapsilosis </italic>copy. This implies that the common ancestor of <italic>C. parapsilosis </italic>and <italic>C. orthopsilosis </italic>acquired the bacterial PhzF homolog after speciation from <italic>L. elongisporus</italic>.</p> <p>Mechanisms of gene transfer into fungi are poorly understood. To date no DNA uptake mechanism has been identified in CTG species. Interkingdom conjugation between bacteria and yeast has been observed however [##REF##2666856##57##, ####REF##8045911##58##, ##REF##8720926##59####8720926##59##]. Similarly, <italic>Saccharomyces cerevisiae </italic>has been shown to be transformant competent under certain conditions [##REF##10953082##60##]. CTG species are known interact with bacteria <italic>in vivo </italic>[##REF##12077418##61##], and it is therefore possible that interkingdom conjugation and transformation may facilitate DNA transfer in <italic>C. parapsilosis</italic>. These mechanisms may also be applicable to the Pezizomycotina species examined in this analysis.</p>
<title>Results and discussion</title> <title>Identification of horizontal gene transfer candidates through Blast database search</title> <p>We compared all available CTG gene sets against UniProt using BlastP [##REF##9254694##31##]. CTG genes with top database hits to bacterial species were identified as putative horizontally transferred genes and the resultant Blast files were inspected manually. A <italic>D. hansenii </italic>gene (protein ydhR precursor) with a top database hit to a bacterial sequence was not considered for further analsyes as it has previously been described [##REF##15229592##22##]. After this process two genes from <italic>C. parapsilosis </italic>were considered for further analysis; one encodes a putative proline racemase, and the second encodes a member of the phenazine F superfamily. Related family members were identified by a second round of database searching against GenBank to ensure all available genomic data was utilized.</p> <title>Proline racemase phylogeny and characterization</title> <p>The <italic>C. parapsilosis </italic>gene (designated CPAG_02038) is most similar to a proline racemase homolog from <italic>Burkholderia cenocepacia </italic>AU 1054 protein (66% pairwise identity; Figure ##FIG##0##1A##). Amino acid racemases catalyze the interconversion of L- and D-amino acids by abstraction of the α-amino proton of the enzyme bound substrate [##REF##5722267##32##]. CPAG_02038 lies within a large contig and is also present in a previously published genome survey of <italic>C. parapsilosis </italic>[##REF##15947193##33##], suggesting its presence does not the result from contamination. We could not locate any related genes in any other CTG genome (using BlastP or TBlastN). Family members are widely distributed throughout the prokaryotes however, and are also located within the Pezizomycotina.</p> <p>We extracted 321 putative proline racemases from 207 organisms, including members of the α, β, γ, and δ-proteobacteria, Actinobacteria, Fungi, Protozoa and Metazoa. Numerous species were found to have several family members [see additional file ##SUPPL##1##2##]; all were included for complete comparative purposes. A maximum likelihood (ML) phylogeny was reconstructed from an alignment of all the PR proteins (Figure ##FIG##1##2##).</p> <p>There are a large number of polytomies displayed in Figure ##FIG##1##2##. These probably result from duplication of PR genes followed by diversifying selection, leading to a high degree of sequence heterogeneity. For example, <italic>Agrobacterium tumefaciens </italic>str. C58 contains three PR homologs [see additional file ##SUPPL##1##2##], with an average amino acid pairwise percentage identity of ~31%. <italic>Burkholderia cenocepacia </italic>AU 1054 contains 2 proline racemase homologs [see additional file ##SUPPL##1##2##], which are only 28% identical. To help resolve the evolutionary history amongst PR homologs we reconstructed an additional ML phylogeny based on a reduced dataset (Figure ##FIG##2##3##). We also reconstructed a Bayesian phylogeny using the heterogeneous CAT site model. The CAT model can account for site-specific features of sequence evolution and has been found to be more robust than other methods against phylogenetic artifacts such as long branch attraction [##REF##17288577##34##]. The resultant Bayesian phylogeny is highly congruent with the ML phylogeny (not shown).</p> <p>The putative <italic>C. parapsilosis </italic>PR homolog lies in a strongly supported (100% Bootstrap support (BP)) clade with <italic>Burkholderia </italic>species (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-A). <italic>Burkholderia </italic>are β-proteobacteria. However, no other β-proteobacteria, or indeed any other bacterial genus were found within clade-A (Figures ##FIG##1##2## &amp;##FIG##2##3##).</p> <p>Although no PR homologs were identified in other CTG species, or indeed in any other of the Saccharomycotina, there are homologs in family members of the Pezizomycotina. A Pezizomycotina specific subclade is evident in our phylogeny containing <italic>Phaeosphaeria nodorum</italic>, <italic>Aspergillus niger </italic>and <italic>Gibberella zeae </italic>(Figures ##FIG##1##2## &amp;##FIG##2##3## clade-B 100% BP). This subclade is found in a strongly supported clade with members of the Actinobacteria (Figure ##FIG##1##2## 100% BP), containing <italic>Brevibacterium linens </italic>and an unclassified marine actinobacterium and excluding <italic>Rubrobacter xylanophilus </italic>(Figure ##FIG##1##2## 87% BP). This suggests that these Pezizomycotina species obtained their PR gene from the Actinobacteridae subclass rather than the Rubrobacteridae subclass. This transfer event is another independent HGT event of a PR gene into fungi, and we hypothesize it occurred early in the Pezizomycotina lineage, as it is shared by three distantly related species. Its patchy phyletic distribution suggests it has been subsequently lost in other Pezizomycotina species.</p> <p>There are also PR homologs in the Metazoans. These are found in a eukaryote clade that also contains a number of Pezizomycotina representatives (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-C 93% BP). Several scenarios can explain this phylogenetic positioning. Firstly, the PR gene may have been present in the last universal common ancestor of all eukaryotes but has been differentially lost in all lineages except those leading to modern day Metazoa and Pezizomycotina. Alternatively, an ancient gene transfer from bacteria to the last common ancestor (LCA) of Metazoa and Fungi could have occurred, with subsequent gene loss amongst different Metazoan and Fungal lineages. A third hypothesis is that two independent gene transfers have occurred into the Metazoan and Pezizomycotina lineages from unsampled bacterial donors. Finally, a transfer from unsampled bacteria into one of the eukaryote clades (either Metazoa or Pezizomycotina) may have occurred with subsequent transfer from one eukaryotic group to the other.</p> <p><italic>A. niger</italic>, <italic>A. oryzae </italic>and <italic>G. zeae </italic>all contain multiple PR homologs [see additional file ##SUPPL##1##2##]. One <italic>A. niger</italic>, one <italic>G. zeae </italic>and the three <italic>A. oryzae </italic>PR homologs are nested in a strongly supported Pezizomycotina specific subclade (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-D 100% BP). This subclade if found within a larger predominately proteobacterial clade (Figure ##FIG##1##2## 74% BP). This infers that there was an independent gene transfer event of a bacterial PR homolog into an ancestral Pezizomycotina species.</p> <p>The phylogenetic position of the <italic>C. parapsilosis </italic>PR homolog (Figures ##FIG##1##2## &amp;##FIG##2##3##) resemble that described for the adenosine deaminase (ADA) gene in the <italic>Dekkera bruxellensis </italic>genome [##REF##17277171##21##]. In that analysis, the authors suggest that <italic>D. bruxellensis </italic>and <italic>Burkholderia </italic>species received the ADA gene from a species not yet represented in the public sequence databases. Our PR phylogeny suggests a similar event may have occurred within clade-A, which contains only <italic>C. parapsilosis </italic>and <italic>Burkholderia </italic>species (Figures ##FIG##1##2## &amp;##FIG##2##3##). <italic>Burkholderia </italic>species are known to have a genomic repertoire that allows the transfer and receipt of exogenous DNA [##REF##12748267##35##] and a number of studies have reported successful gene transfers into <italic>Burkholderia </italic>species [##REF##15493822##36##,##UREF##0##37##]. It is possible therefore that there have been other successful gene transfers into this bacterial lineage.</p> <p>The vast majority of amino acids found in living cells correspond to the L-stereoisomer [##REF##16446443##38##]. However, D-amino acids are long known to be found in the cell walls of Gram positive and negative bacteria, where they are essential components of peptidoglycan [##REF##8749373##39##]. Apart from low levels of D-amino acids derived from spontaneous racemization as a result of aging [##REF##9949871##40##], it was assumed that only L-amino acid enantiomers were present in eukaryotes [##REF##12735293##41##]. However, recent studies have reported the presence of numerous D-amino acids in an array of organisms, including mammals [##REF##10557334##42##]. The first eukaryotic (proline) amino acid racemase has recently been described from the human pathogen <italic>Trypanosoma cruzi </italic>[##REF##10932226##43##]. A high degree of sequence similarity was observed between the <italic>T. cruzi </italic>and bacterial homologs [##REF##10932226##43##]. Our phylogeny infers that <italic>T. cruzi </italic>obtained its PR homolog through interdomain HGT from a member of the Firmicutes (subclass Clostridia), as it is grouped beside members of this group with a high degree of support (Figure ##FIG##1##2## clade-E 96% BP). We performed database searches [##UREF##1##44##], against other Protozoan genomes including <italic>Trypanosoma brucei</italic>, <italic>Trypanosoma congolense </italic>and <italic>Trypanosoma annulata</italic>. We failed to locate a homolog in all species except for <italic>T. vivax</italic>.</p> <p>Previous analysis has shown that <italic>T. cruzi </italic>and <italic>T. vivax </italic>are not each others closest phylogenetic neighbors, relative to the other species sampled [##REF##10633190##45##]. This suggests an ancestral Trypanosoma gained the PR gene and multiple losses in different Trypanosoma lineages has subsequently occurred.</p> <title>Gene order around PR homologs</title> <p>The <italic>C. parapsilosis </italic>PR homolog lies close to an ortholog (CPAG_02041) of orf19.1135 from <italic>C. albicans </italic>(Figure ##FIG##3##4##). The gene order to the left of this ORF is conserved in all CTG species, the order to the right is conserved in most CTG species apart from <italic>C. parapsilosis </italic>and <italic>L. elongisporus</italic>. <italic>C. parapsilosis </italic>and <italic>L. elongisporus </italic>are closely related [##REF##17121679##29##], and an examination of synteny suggests that the PR gene (together with a second ORF, cpar5437) were inserted between CPAG_2041 and CPAG_2037 (Figure ##FIG##3##4##). cpar5437 encodes a neutral amino acid (AA) transporter. The presence of an AA transporter beside the PR homolog is interesting. If the putative proline racemase has a role in amino acid metabolism, then the presence of the transporter may be the result of an adaptive translocation to enhance the activity of the PR gene. Unlike the PR ORF the AA transporter is fungal in origin. Most CTG species contain a single neutral AA transporter; however <italic>C. parapsilosis </italic>and <italic>D. hansenii </italic>have four.</p> <p>We located tRNA genes for nearly all CTG species beside the large conserved syntenic block (Figure ##FIG##3##4##). It has been shown that tRNA genes are associated with genomic breakpoints [##REF##10839539##46##]. We hypothesize that a genomic rearrangement has occurred at this site in the LCA of <italic>C. parapsilosis </italic>and <italic>L. elongisporus</italic>. We cannot determine if the bacterial PR homolog was inserted into the LCA of <italic>L. elongisporus</italic>/<italic>C. parapsilosis </italic>and subsequently lost in <italic>L. elongisporus</italic>, or gained by <italic>C. parapsilosis </italic>after speciation.</p> <p>We also investigated the gene order around the Pezizomycotina PR homologs [see additional file ##SUPPL##2##3##]. Gene synteny around the PR homologs found in clade-D (Figures ##FIG##1##2## &amp;##FIG##2##3##) is not conserved (not shown). Interestingly however, both <italic>A. niger </italic>and <italic>G. zeae </italic>in clade-D (Figures ##FIG##1##2## &amp;##FIG##2##3##) have genes containing a FAD dependent oxidoreductase domain in close proximity to their PR homologs (not shown). According to Pfam [##UREF##2##47##], FAD dependent oxidases include D-amino acid oxidases, that catalyze the oxidation of neutral and basic D-amino acids into their corresponding keto acids. The presence of these oxidases may be another example of an adaptive translocation to enhance the activity of the PR gene in these Pezizomycotina species.</p> <p><italic>A. oryzae </italic>has three PR homologs (Figures ##FIG##1##2## &amp;##FIG##2##3## clade-C). All of these have orthologs in its close relative <italic>A. flavus </italic>(Figure ##FIG##2##3## clade-C), and synteny around these is conserved [see additional file ##SUPPL##2##3## clade-D]. The remaining two species in clade-C are <italic>A. niger </italic>and <italic>G. zeae</italic>. There is no evidence of conserved gene order within these species, or with <italic>A. oryzae </italic>or <italic>A. flavus</italic>. Gene order around the <italic>A. flavus </italic>and <italic>A. terreus </italic>PR homologs found in the Metazoan/Pezizomycotina clade (Figures ##FIG##1##2## &amp;##FIG##2##3##) is also conserved [see additional file ##SUPPL##2##3##], as is the order between <italic>A. fumigatus </italic>and <italic>N. fishceri </italic>[see additional file ##SUPPL##2##3##]. We could not locate amino acid transporters or FAD dependent oxidases beside any of the PR homologs found in clades B or C (Figure ##FIG##1##2##).</p> <title>Proline racemase codon usage</title> <p>It has been shown that recently acquired genes often display an atypical codon preference when compared to other genes in the genome [##REF##1762151##48##,##REF##10830951##49##]. However, the transferred PR homologs have a codon usage consistent with the rest of their genomes [see additional file ##SUPPL##3##4##]. We undertook an analysis of variation in synonymous codon usage on all PR genes shown in Figure ##FIG##1##2##. Homologs from related species cluster together [see additional file ##SUPPL##4##5##]. For example, the Actinobactria, the Firmicutes and the <italic>Burkholderia </italic>species all inhabit unique areas in two dimensional correspondence analysis space [see additional file ##SUPPL##4##5##].</p> <p>The majority of fungal and Metazoan PRs are clustered together [see additional file ##SUPPL##4##5##]. The <italic>C. parapsilosis </italic>PR homolog has a codon usage distinct from the other Pezizomycotina fungal PR homologs [see additional file ##SUPPL##4##5##], which is unsurprising as <italic>C. parapsilosis </italic>belongs to the Saccharomycotina subphylum. The <italic>C. parapsilosis </italic>homolog is also separate from the <italic>Burkholderia </italic>(β-proteobacteria) genes with which it forms a closely related phylogenetic group (Figures ##FIG##1##2## &amp;##FIG##2##3##). This suggests that the gene may have originated from a genome with no other close relatives among the species analyzed here.</p> <title>Proline racemase activity</title> <p>The PR active site from <italic>Trypanosoma cruzi, Clostridium sticklandii, Agrobacterium tumefaciens</italic>, <italic>Brucella melitensis </italic>and <italic>Pseudomonas aeruginosa </italic>all contain cysteine at amino acid position 330 [##REF##10932226##43##,##REF##1174513##50##]. This amino acid is essential for enzymatic function, because substitution with serine abolishes activity [##REF##12735293##41##]. However, PR homologs from human, mouse, <italic>Rhizobium </italic>and <italic>Brucella </italic>contain a threonine instead of a cysteine at position 330 [##REF##12735293##41##]. We observed that cysteine is found in the equivalent position in many of the bacterial proteins. The Pezizomycotina PR genes found in clade-B and clade-D contain a cysteine at the active site (Figure ##FIG##2##3##). The PR homologs found in the Metazoan/Pezizomycotina clade (clade-B) have a threonine at position 330. Similarly, the <italic>C. parapsilosis </italic>PR homolog, together with its relatives from <italic>Burkholderia </italic>all contain a threonine (Figure ##FIG##2##3##). However, <italic>Burkholderia </italic>species have multiple PR homologs [see additional file ##SUPPL##1##2##] with a cysteine as the active site (not shown). It is not clear what effect the substitution has on enzyme activity. It has been suggested that homologs containing threonine at the active site are not true PRs [##REF##12735293##41##], but may instead belong to a superfamily. We cannot detect any difference in the ability of <italic>C. parapsilosis</italic>, the other CTG species or any of the Pezizomycotina species to utilize D-proline as growth media (data not shown). We therefore cannot confidently infer the function of the PR homologs in the fungi analyzed here.</p> <title>Phenazine F phylogeny and characterisation</title> <p>The <italic>C. parapsilosis </italic>gene (designated CPAG_03462) is most similar to a <italic>Photorhabdus luminescens </italic>phenazine F (PhzF) protein with 61% pairwise identity (Figure ##FIG##0##1B##). Phenazines are biologically active compounds, all of which have a characteristic tricyclic ring system and have been shown to confer a selective growth advantage to organisms which secrete them, as they possess broad-spectrum antibiotic activity towards bacteria, fungi and higher eukaryotes [##REF##15545603##51##]. In <italic>Pseudomonas</italic>, the best studied phenazine producer, PhzF is part of an operon required for the conversion of chorismic acid to phenazine-1-carboxylate (PCA) [##REF##15449932##52##]. PhzF homologs were identified in most of the CTG species tested as well as several other fungal species. However, we could not identify a PhzF homolog in the <italic>L. elongisporus </italic>genome, even when multiple TBlastN and BlastN searches were used.</p> <p>PhzF homologs were extracted from GenBank for subsequent phylogenetic analysis. In total 181 representative protein coding sequences distributed amongst 154 organisms were used. These taxa were distributed amongst α, β, γ and δ-proteobacteria, Actinobacteria, Fungi, Firmicutes a well as other bacterial groups.</p> <p>We aligned all sequences and reconstructed a PhzF ML phylogeny (Figure ##FIG##4##5##). The <italic>C. parapsilosis </italic>PhzF homolog is found in a clade with members of the β-proteobacteria <italic>(Burkholderia multiovorans, Burkholderia cepacia, Burkholderia ambifaria</italic>), α-proteobacteria (<italic>Roseovarius</italic>) and the γ-proteobacteria (<italic>Azotobacter vinelandii, Acinetobacter baumannii, Shewanella baltica and Photorhabdus luminescens</italic>) (81% BP). In contrast, all other PhzF homologs from CTG species are in a completely separate clade (Figure ##FIG##4##5##). These form a sister group (63% BP) to PhzF homologs from other Saccharomycotina species (<italic>C. glabrata, Saccharomyces cerevisiae, Kluyveromyces lactis </italic>and <italic>Vanderwaltozyma polyspora</italic>). All three clades are grouped together in a larger clade with high support (75% BP).</p> <p>The sister group relationship between the PhzF homologs from the Ascomycota and the proteobacteria clade is intriguing (Figure ##FIG##4##5##), as it suggests that an ancestral Saccharomycotina species gained the PhzF homolog from a proteobacteria. The bacterial PhzF gene has subsequently been retained after multiple speciation events, but lost in <italic>C. parapsilosis</italic>. We hypothesize that <italic>C. parapsilosis </italic>has recently reacquired a bacterial PhzF homolog from a proteobacterial source, as it is grouped (81% BP) within a proteobacterial subclade. To test this hypothesis we reconstructed constrained trees that placed <italic>C. parapsilosis </italic>together with the remaining Ascomycota species [see additional file ##SUPPL##5##6## C-H]. The AU test of phylogenetic tree selection [##REF##12079646##53##], showed that the original unconstrained tree (groups <italic>C. parapsilosis </italic>with proteobacteria) receives the optimal likelihood tree score, and the differences in likelihood scores when compared to the constrained trees [see additional file ##SUPPL##5##6##], are significant (P &lt; 0.05). This is also supported by spectral analysis [see additional file ##SUPPL##6##7##].</p> <p>Our phylogeny shows that the <italic>Schizosaccharomyces pombe </italic>PhzF homolog is found in a clade containing all CTG PhzF homologs (Figure ##FIG##4##5## 99% BP). Furthermore it is grouped beside <italic>D. hansenii </italic>(66% BP). <italic>S. pombe </italic>is not a member of the Saccharomycotina, it belongs to the Taphrinomycotina subphylum. The genome sequences of <italic>Schizosaccharomyces japonicus </italic>and <italic>Schizosaccharomyces octosporus </italic>have recently been completed [##UREF##3##54##]. We could not locate a PhzF homolog in <italic>S. japonicus </italic>but did locate a homolog in <italic>S. octosporus </italic>using a TBlastN search strategy. Phylogenetic analysis has shown that <italic>S. pombe </italic>and <italic>S. octosporus </italic>are more closely related to one another than to <italic>S. japonicus </italic>[##REF##12527786##55##]. Therefore we hypothesize that the LCA ancestor of <italic>S. pombe </italic>and <italic>S. octosporus </italic>gained the PhzF gene from an ancestral <italic>D. hansenii</italic>-like species after speciation from <italic>S. japonicus</italic>. We reconstructed a constrained tree that placed <italic>S. pombe </italic>outside the Saccharomycotina clade [see additional file ##SUPPL##5##6B##]. The approximately unbiased test of phylogenetic tree selection (AU test) [##REF##12079646##53##], showed that the phylogenetic inferences of the unconstrained tree are significantly better (P &lt; 0.05) than the constrained tree [see additional file ##SUPPL##5##6##]. This infers that <italic>S. pombe </italic>has obtained a PhzF homolog from a member of the CTG clade.</p> <p>A small basidiomycete clade is evident amongst prokaryote species (Figure ##FIG##4##5##). Both <italic>Ustilago maydis </italic>and <italic>Malassezia globosa </italic>belong to the Ustilaginomycotina subphylum. Therefore our phylogeny infers that an ancestral Ustilaginomycotina species gained a PhzF gene from an unknown bacterial source, and both species have retained this after speciation.</p> <p>A correspondence analysis of synonymous codon usage for all PhzF homologs was also performed and is shown in additional information [see additional file ##SUPPL##7##8##]. The <italic>S. pombe </italic>PhzF homolog has a codon usage pattern very similar to the <italic>D. hansenii </italic>protein.</p> <title>Gene order around PhzF</title> <p>Analysis of the genes adjacent to the PhzF homolog in <italic>C. parapsilosis </italic>shows that there is a high conservation of gene synteny and supports our hypothesis that PhzF was recently acquired in this species (Figure ##FIG##5##6##). Homologs in the other CTG species are located in completely different regions of the genome relative to <italic>C. parapsilosis </italic>(not shown). For example, the <italic>C. albicans </italic>PhzF homolog is located between orf19.5619 and orf19.5621, whereas the <italic>C. parapsilosis </italic>homolog is found between orf19.6689 &amp; orf19.6687 relative to <italic>C. albicans </italic>SC5314 (Figure ##FIG##5##6##). However, the <italic>L. elongisporus </italic>genome contains no PhzF homolog, either at a position equivalent to the <italic>C. parapsilosis </italic>copy or elsewhere in the genome.</p> <p>We propose that the LCA of <italic>L. elongisporus </italic>and <italic>C. parapsilosis </italic>lost the PhzF gene present in the other CTG species, and a second (new) copy was subsequently gained by <italic>C. parapsilosis </italic>after speciation. We have partial sequence data (unpublished) from <italic>Candida orthopsilosis</italic>, a species so closely related to <italic>C. parapsilosis </italic>that it was once designated <italic>C. parapsilosis </italic>group II [##REF##7665651##56##]. We located a <italic>C. orthopsilosis </italic>PR homolog that is 83% identical (at the amino acid level) to the <italic>C. parapsilosis </italic>copy. This implies that the common ancestor of <italic>C. parapsilosis </italic>and <italic>C. orthopsilosis </italic>acquired the bacterial PhzF homolog after speciation from <italic>L. elongisporus</italic>.</p> <p>Mechanisms of gene transfer into fungi are poorly understood. To date no DNA uptake mechanism has been identified in CTG species. Interkingdom conjugation between bacteria and yeast has been observed however [##REF##2666856##57##, ####REF##8045911##58##, ##REF##8720926##59####8720926##59##]. Similarly, <italic>Saccharomyces cerevisiae </italic>has been shown to be transformant competent under certain conditions [##REF##10953082##60##]. CTG species are known interact with bacteria <italic>in vivo </italic>[##REF##12077418##61##], and it is therefore possible that interkingdom conjugation and transformation may facilitate DNA transfer in <italic>C. parapsilosis</italic>. These mechanisms may also be applicable to the Pezizomycotina species examined in this analysis.</p>
<title>Conclusion</title> <p>We investigated the frequency of recent interkingdom gene transfer between CTG and bacterial species. We located two strongly supported incidences of HGT, both within the <italic>C. parapsilosis </italic>genome. We also located independent transfers into the Pezizomycotina, Basidiomycotina and Protozoan lineages.</p> <p>We cannot determine the exact origin of the PR homolog (CPAG_02038) found in the <italic>C. parapsilosis </italic>genome. However, based on its phylogenetic position it either originated from a <italic>Burkholderia </italic>source, or more likely an organism not yet represented in the sequence databases. Our PR phylogenetic analysis also suggests there were two independent transfers into Pezizomycotina species, one from an Actinobacterial source, and the second is from an unknown proteobacterial source. There is also evidence that <italic>T. cruzi </italic>has obtained its PR homolog from a Firmicutes species. The transferred PR genes analyzed here belong to a superfamily of proline racemases, although we cannot determine their exact function in the fungal species examined. Their proximity to an amino acid transporter (in <italic>C. parapsilosis</italic>) and a FAD dependent oxidoreductase (in <italic>A. niger </italic>and <italic>G. zeae</italic>) suggests they do have a role in amino acid metabolism. Furthermore, evidence of multiple independent transfers into fungi suggests the protein does confer a biological advantage, although we cannot determine what is. The bacteria-derived PR gene has the potential to be a novel antifungal drug target as there would be no undesired host protein-drug interactions.</p> <p>The bacterial PhzF homolog (CPAG_03462) found in <italic>C. parapsilosis </italic>most likely originated from a proteobacterial source. Most CTG species examined contained PhzF homologs, with the exception of <italic>L. elongisporus</italic>. The crystal structure the PhzF homolog in <italic>S. cerevisiae </italic>has been determined and while its function remains unknown, it is not thought to be involved in phenazine production [##REF##16021630##62##]. We postulate that the PhzF homolog present in other CTG species was initially lost by the ancestor of <italic>C. parapsilosis </italic>and <italic>L. elongisporus</italic>, but subsequently regained by <italic>C. parapsilosis </italic>through HGT. The loss of eukaryote genes and subsequent reacquisition of a prokaryotic copy has previously been described in yeast, and can confer specific metabolic capabilities. An analysis of the biotin biosynthesis pathway discovered that the ancestor of <italic>Candida</italic>, <italic>Debaryomyces</italic>, <italic>Kluyveromyces </italic>and <italic>Saccharomyces </italic>lost the majority of the pathway after the divergence from the ancestor of <italic>Y. lipolytica</italic>. However, <italic>Saccharomyces </italic>species have rebuilt the biotin pathway through gene duplication/neofunctionalization after horizontal gene transfer from α and γ proteobacterial sources [##REF##18073433##20##]. The acquisition of the <italic>URA1 </italic>gene (encoding dihydroorotate dehydrogenase) from <italic>Lactobacillus </italic>and replacement of the endogenous gene in <italic>S. cerevisiae</italic>, allowed growth under anaerobic conditions [##REF##15947202##19##]. Similarly, acquisition of <italic>BDS1 </italic>(alkyl-aryl-sulfatase) from proteobacteria may have enabled the survival of <italic>S. cerevisiae </italic>in a harsh soil environment [##REF##15947202##19##]. Our PhzF phylogeny suggests that the PhzF homolog found in most CTG species originated from an ancient HGT event, from a member of the proteobacteria. Our analysis also shows that <italic>S. pombe </italic>has obtained a PhzF homolog from a CTG species, most likely one closely related to <italic>D. hansenii</italic>. There is also phylogenetic evidence showing that an ancestral Ustilaginomycotina species gained a PhzF gene from an unknown bacterial source. We cannot however, determine the biological advantage to the organisms.</p> <p>Although it was not the major goal of this study, we did locate HGT from bacteria into fungal genomes outside the CTG clade, and also inter-fungal transfers. In a previous analysis of HGT in diplomonads, fifteen genes were found to have undergone HGT [##REF##12546782##18##]. There is phylogenetic evidence that these genes have undergone independent transfers into other eukaryotic lineages including Fungi. Therfore, in eukaryotes just as HGT has affected some species more than others [##REF##16701456##63##], there may be groups of genes that are more likely to be taken up through HGT than others. We cannot test this directly however, as we have not identified all cases of HGT from bacteria to fungi outside the CTG clade.</p> <p>Our analysis indicates that recent interkingdom gene transfer into extant CTG species is negligible. This supports a previous hypothesis that genetic code alterations blocks horizontal gene transfer [##REF##17932489##64##]. It should be noted however that we searched for recent bacterial gene transfers into individual CTG species, and not for more ancient transfers. We took this approach because the presence of recently gained bacterial genes in a eukaryote genome should be readily detected compared to older transfers. Similarly, we have not investigated eukaryote-to-eukaryote transfers. It is therefore possible that we have underestimated the overall rate of HGT into the CTG lineage. The discovery of HGT in other fungal lineages implies that HGT plays an important role in fungal evolution and deserves further analysis. In particular a strategy which can detect ancient gene transfers would be meaningful.</p>
<title>Background</title> <p>To date very few incidences of interdomain gene transfer into fungi have been identified. Here, we used the emerging genome sequences of <italic>Candida albicans </italic>WO-1, <italic>Candida tropicalis, Candida parapsilosis, Clavispora lusitaniae, Pichia guilliermondii</italic>, and <italic>Lodderomyces elongisporus </italic>to identify recent interdomain HGT events. We refer to these as CTG species because they translate the CTG codon as serine rather than leucine, and share a recent common ancestor.</p> <title>Results</title> <p>Phylogenetic and syntenic information infer that two <italic>C. parapsilosis </italic>genes originate from bacterial sources. One encodes a putative proline racemase (PR). Phylogenetic analysis also infers that there were independent transfers of bacterial PR enzymes into members of the Pezizomycotina, and protists. The second HGT gene in <italic>C. parapsilosis </italic>belongs to the phenazine F (PhzF) superfamily. Most CTG species also contain a fungal PhzF homolog. Our phylogeny suggests that the CTG homolog originated from an ancient HGT event, from a member of the proteobacteria. An analysis of synteny suggests that <italic>C. parapsilosis </italic>has lost the endogenous fungal form of PhzF, and subsequently reacquired it from a proteobacterial source. There is evidence that <italic>Schizosaccharomyces pombe </italic>and Basidiomycotina also obtained a PhzF homolog through HGT.</p> <title>Conclusion</title> <p>Our search revealed two instances of well-supported HGT from bacteria into the CTG clade, both specific to <italic>C. parapsilosis</italic>. Therefore, while recent interkingdom gene transfer has taken place in the CTG lineage, its occurrence is rare. However, our analysis will not detect ancient gene transfers, and we may have underestimated the global extent of HGT into CTG species.</p>
<title>Authors' contributions</title> <p>DAF, MEL and GB were involved in the design phase. MEL predicted genes in unannotated genomes. DAF sourced homologs, examined synteny and performed phylogenetic analyses. DAF and GB drafted the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>The authors wish to acknowledge the Wellcome Trust Sanger Institute and Broad institute of MIT &amp; Harvard for releasing data ahead of publication. We would like to acknowledge the financial support of the Irish Research Council for Science, Engineering and Technology (IRCSET), the Irish Health Research Board (HRB) and Science Foundation Ireland (SFI). We wish to acknowledge the SFI/HEA Irish Centre for High-End Computing (ICHEC) for the provision of computational facilities and support. We thank Mike Lorenz and Paul Dyer for fungal strains. We also thank Jason Stajich for help with gene annotations.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>A) </bold>An alignment of PR proteins from <italic>C. parapsilosis </italic>(CpPR CPAG_02038) and <italic>Burkholderia cenocepacia </italic>(BcPR) was generated with MUSCLE. These proteins are 66% identical <bold>B) </bold>An alignment of PhzF proteins from <italic>Candida parapsilosis </italic>(CpPhzF CPAG_03462) and <italic>Photorhabdus luminescens </italic>(PlPhzF), these are 61% identical.</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Proline racemase maximum likelihood phylogeny</bold>. The optimum model of protein substitution was found to be WAG+G. The number of gamma rate categories was 4 (alpha = 1.163). Bootstrap resampling (100 iterations) was undertaken and are displayed. For display purposes branches with less than 50% support were collapsed. Letters (A-E) in parentheses are used to distinguish clades and are discussed in the text. Branches are colored according to their taxonomy. Fungal branches and species names are colored green.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Reduced Proline racemase maximum likelihood phylogeny with active site alignment</bold>. Bootstrap resampling (100 iterations) was undertaken and percentages are displayed. Fungal branches are shown in green. An alignment around the active site is also displayed. Clade letters in parentheses correspond to those in Figure 2. The phylogeny is rooted around the Metazoan/Pezizomycotina specific clade (clade-C), all members of this clade have a threonine at the active site. <italic>C. parapsilosis </italic>and its phylogenetic neighbors have a threonine instead of a cysteine at the active site (clade-A). <italic>A. oryzae, A. niger </italic>and <italic>G. zeae </italic>all contain cysteine at the active site (clade-D). <italic>A. flavus</italic>, <italic>A. oryzae</italic>, <italic>A. niger</italic>, <italic>A. nidulans </italic>and <italic>G. zeae </italic>also have cysteine at the active site (clade-C).</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Gene order around <italic>C. parapsilosis </italic>proline racemase gene</bold>. Species names and identifiers are shown in each box. Gene identifiers relate to annotations from the Broad Institute [66]. On the left hand side orthologous genes are stacked under one another in pillars. Relative positions of t-RNA genes are shown and may indicate a breakpoint. After the breakpoint, synteny is conserved between <italic>C. albicans</italic>, <italic>C. dubliniensis</italic>, <italic>C. tropicalis</italic>, <italic>D. hansenii </italic>and <italic>Cl. lusitaniae</italic>. Synteny between <italic>C. parapsilosis </italic>and <italic>L. elongisporus </italic>is conserved but differs to the other CTG species. <italic>C. parapsilosis </italic>has a proline racemase (PR CPAG_02038) and a neutral amino acid transporter (AA cpar5437) insertion in this region. cpar5437 is absent from the Broad gene list but present in our manual gene call.</p></caption></fig> <fig id="F5" position="float"><label>Figure 5</label><caption><p><bold>PhzF maximum likelihood phylogeny</bold>. The optimum model of protein substitution was found to be WAG+G. The number of gamma rate categories was 4 (alpha = 0.873). Bootstrap resampling (100 iterations) was undertaken and are displayed. For display purposes branches with less than 50% support were collapsed. Branches are colored according to their taxonomy. Fungal branches are shown in green. The <italic>S. pombe </italic>PhzF homolog is highlighted with a red rectangle.</p></caption></fig> <fig id="F6" position="float"><label>Figure 6</label><caption><p><bold>Gene order around <italic>C. parapsilosis </italic>PhzF gene</bold>. Species names and gene identifiers are shown in each box. Orthologous genes are stacked under one another in pillars. The <italic>C. parapsilosis </italic>PhzF homolog (CPAG_03462) is highlighted with a red box. Synteny relative to the <italic>C. parapsilosis </italic>PhzF homolog is conserved in all species. Other CTG PhzF homologs are found in completely different regions of the genome relative to <italic>C. parapsilosis</italic>. <italic>L. elongisporus </italic>is the only CTG species missing a PhzF gene, and there is no evidence for a pseudogene in the genome.</p></caption></fig>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Examples of reported incidences of interkingdom gene transfer between prokaryotes and fungi</bold>. One <italic>Kluyveromyces lactis </italic>gene (KLLA0D19949g) previously highlighted [##REF##15229592##22##], been omitted as it is no longer recognized as an ORF. <italic>Y. lipolytica </italic>genes denoted with a * and ^ indicate possible gene duplications after HGT.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p><bold>GenBank accession numbers for PR (A) and PhzF (B) sequences used in this analysis</bold>. Species identified with an * use the accession numbers created by the Broad Institute [##UREF##5##66##] or the Wellcome Trust Sanger Institute [##UREF##6##67##].</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p><bold>Gene order around Pezizomycotina proline racemase genes</bold>. Species names and identifiers are shown in each box. PR genes are labeled. Gene identifiers relate to annotations from the Broad Institute. Clade letters in parentheses correspond to those in Figure ##FIG##1##2##. There is evidence for conserved gene synteny between some species such as <italic>A. oryzae </italic>and <italic>A. flavus </italic>(clade-C). <italic>A. flavus</italic>/<italic>A. terreus </italic>and <italic>N. fischeri</italic>/<italic>A. fumigatus </italic>in the Metazoan/Pezizomycotina clade (B). The <italic>A. flavus </italic>gene denoted by a * is absent from the Broad gene set but we were able to locate it with a BlastX search.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p><bold>Correspondence analysis of codon usage</bold>. Correspondence analysis of codon usage in the <italic>C. parapsilosis </italic>(1), <italic>U. maydis </italic>(2), <italic>M. globosa </italic>(3), <italic>A. flavus </italic>(4), <italic>A. niger </italic>(5), <italic>G. zeae </italic>(6), <italic>A. oryzae </italic>(7), <italic>P. nodorum </italic>(8), and <italic>S. pombe </italic>(9) genomes. Transferred genes are highlighted. All have a codon usage similar to the rest of their genomes which is unsurprising as transferred genes have been shown to ameliorate their codon usage to their hosts [##REF##9089078##79##].</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p><bold>Correspondence analysis of codon usage in the proline racemase gene family analyzed in this study</bold>. Major groups are color-coded. The <italic>C. parapsilosis </italic>PR gene has a codon usage pattern distinct from other fungal species in this analysis. It is also quite distinct from the <italic>Burkholderia </italic>(β-proteobacteria) species, which were found to be its phylogenetic neighbors.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p><bold>Trees for approximately unbiased test for PhzF homologs</bold>. Tree A is the original unconstrained topology, which groups <italic>C. parapsilosis </italic>with proteobacteria. Topology B is a constrained tree that places S. pombe outside the Saccharomycotina clade. Topologies C-H are constrained and place <italic>C. parapsilosis </italic>amongst the other Saccharomycotina species. Log likelihood scores for each tree are given. To assess the likelihood that any differences in topology between the inferred trees is no more significant that that expected by chance, we performed the approximately unbiased test. The AU test shows that the unconstrained tree receives the optimal likelihood tree score. Furthermore, the differences in likelihood scores when compared to the constrained trees are significant (P &lt; 0.05). Therefore based on these results the placement of the <italic>C. parapsilosis </italic>homolog in the proteobacterial clade to the exclusion of the Saccharomycotina and S. pombe within the Saccharomycotina clade is significant.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p><bold>PhzF spectral analysis</bold>. Analysis was performed on the Saccharomycotina and selected proteobacterial clade. Bars above the x-axis represent frequency of support for each split. Bars below the x-axis represent the sum of all corresponding conflicts. Clad grams above columns represent the corresponding splits in the data. There is no support for the placement of <italic>C. parapsilosis </italic>with the other Saccharomycotina species.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p><bold>Correspondence analysis of codon usage in the PhzF gene family analyzed in this study</bold>. Major groups are color-coded. The <italic>C. parapsilosis </italic>PhzF gene has a codon usage pattern similar to other CTG species analyzed. It is quite distinct from the proteobacterial species that were found to be its phylogenetic neighbors.</p></caption></supplementary-material>
<graphic xlink:href="1471-2148-8-181-1"/> <graphic xlink:href="1471-2148-8-181-2"/> <graphic xlink:href="1471-2148-8-181-3"/> <graphic xlink:href="1471-2148-8-181-4"/> <graphic xlink:href="1471-2148-8-181-5"/> <graphic xlink:href="1471-2148-8-181-6"/>
<media xlink:href="1471-2148-8-181-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-181-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-181-S3.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-181-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-181-S5.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-181-S6.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-181-S7.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-181-S8.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Tuanyok", "Auerbach", "Brettin", "Bruce", "Munk", "Detter", "Pearson", "Hornstra", "Sermswan", "Wuthiekanun", "Peacock", "Currie", "Keim", "Wagner"], "given-names": ["A", "RK", "TS", "DC", "AC", "JC", "T", "H", "RW", "V", "SJ", "BJ", "P", "DM"], "article-title": ["A horizontal gene transfer event defines two distinct groups within Burkholderia pseudomallei that have dissimilar geographic distributions"], "source": ["J Bacteriol"], "year": ["2007"]}, {"article-title": ["GeneDB"], "ext-link": ["http://www.genedb.org"]}, {"article-title": ["Pfam"], "ext-link": ["http://pfam.sanger.ac.uk"]}, {"article-title": ["The Schizosaccharomyces group at the Broad Institute"], "ext-link": ["http://www.broad.mit.edu/annotation/genome/schizosaccharomyces_group/MultiHome.html"]}, {"article-title": ["The Candida Genome Database"], "ext-link": ["http://www.candidagenome.org"]}, {"article-title": ["The Candida group at the Broad Institute "], "ext-link": ["http://www.broad.mit.edu/annotation/genome/candida_group/MultiHome.html"]}, {"article-title": ["The Wellcome Trust Sanger Institute "], "ext-link": ["http://www.sanger.ac.uk/"]}, {"article-title": ["The UniProt database"], "ext-link": ["ftp://ftp.uniprot.org/pub/databases/"]}, {"surname": ["Charleston"], "given-names": ["MA"], "article-title": ["Spectrum: spectral analysis of phylogenetic data"], "source": ["Bioinformatics (Oxford, England)"], "year": ["1998"], "volume": ["14"], "issue": ["1"], "fpage": ["98"], "lpage": ["99"], "pub-id": ["10.1093/bioinformatics/14.1.98"]}]
BMC Evol Biol. 2008 Jun 24; 8:181