Functional assignment of gene AAC16202.1 from Rhodobacter capsulatus SB1003: New insights into the bacterial SDR sorbitol dehydrogenases family

Short-chain dehydrogenases/reductases (SDR) constitute one of the largest enzyme superfamilies with over 60,000 non-redundant sequences in the database, many of which need a correct functional assignment. Among them, the gene AAC16202.1 (NCBI) from Rhodobacter capsulatus SB1003 has been assigned in Uniprot both as a sorbitol dehydrogenase (#D5AUY1) and, as an N-acetyl-d-mannosamine dehydrogenase (#O66112), both enzymes being of biotechnological interest. When the gene was overexpressed in Escherichia coli Rosetta (DE3)pLys, the purified enzyme was not active toward N-acetyl-d-mannosamine, whereas it was active toward d-sorbitol and d-fructose. However, the relative activities toward xylitol and l-iditol (0.45 and 6.9%, respectively) were low compared with that toward d-sorbitol. Thus, the enzyme could be considered sorbitol dehydrogenase (SDH) with very low activity toward xylitol, which could increase its biotechnological interest for determining sorbitol without the unspecific cross-determination of added xylitol in food and pharma compositions. The tetrameric enzyme (120 kDa) showed similar catalytic efficiency (2.2 × 10³ M⁻¹ s⁻¹) to other sorbitol dehydrogenases for d-sorbitol, with an optimum pH of 9.0 and an optimum temperature of 37 °C. The enzyme was also more thermostable than other reported SDH, ammonium sulfate being the best stabilizer in this respect, increasing the melting temperature (T_m) up to 52.9 °C. The enzyme can also be considered as a new member of the Zn²⁺ independent SDH family since no effect on activity was detected in the presence of divalent cations or chelating agents. Finally, its in silico analysis enabled the specific conserved sequence blocks that are the fingerprints of bacterial sorbitol dehydrogenases and mainly located at C-terminal of the protein, to be determined for the first time. This knowledge will facilitate future data curation of present databases and a better functional assignment of newly described sequences.     

Freshwater selenium-methylating bacterial thiopurine methyltransferases: diversity and molecular phylogeny

The diversity of bacterial thiopurine methyltransferases (bTPMT) among five natural Se-methylating freshwaters was investigated by polymerase chain reaction (PCR) screenings and sequencings. DNA sequence analyses confirmed the cloned products' identity and revealed a broad diversity of freshwater TPMTs. Neighbour-joining (NJ) phylogenetic analyses combining these sequences, all GenBank entries closely related to these sequences and deduced TPMTs obtained in this work from selected gamma-proteobacteria showed TPMTs to form a distinct radiation, closely related to UbiG methyltransferases. Inside the TPMT phylogenetic cluster, eukaryote sequences diverged early from the bacterial ones, and all the bacterial database entries belonged to a subgroup of gamma-proteobacteria, with an apparent lateral transfer of a particular allele to beta-proteobacteria of Bordetella. The NJ phylogenetic tree revealed 22 bTPMT lineages, 10 of which harboured freshwater sequences. All lineages showed deep and long branches indicative of major genetic drifts outside regions encoding highly conserved domains. Selected residues among these highly variable domains could reflect adaptations for particular ecological niches. PCR lineage-specific primers differentiated Se-methylating freshwaters according to their 'tpm lineage' signatures. Most freshwater tpm alleles were found to be distinct from those available in the databases, but a group of tpm was found encoding TPMTs identical to an Aeromonas veronii TPMT characterized in this work.     

Phylogeny of γ‐proteobacteria: resolution of one branch of the universal tree?

The reconstruction of bacterial evolutionary relationships has proven to be a daunting task because variable mutation rates and horizontal gene transfer (HGT) among species can cause grave incongruities between phylogenetic trees based on single genes. Recently, a highly robust phylogenetic tree was constructed for 13 gamma-proteobacteria using the combined alignments of 205 conserved orthologous proteins.1 Only two proteins had incongruent tree topologies, which were attributed to HGT between Pseudomonas species and Vibrio cholerae or enterics. While the evolutionary relationships among these species appears to be resolved, further analysis suggests that HGT events with other bacterial partners likely occurred; this alters the implicit assumption of gamma-proteobacteria monophyly. Thus, any thorough reconstruction of bacterial evolution must not only choose a suitable set of molecular markers but also strive to reduce potential bias in the selection of species.     

Evolutionary placement of Xanthomonadales based on conserved protein signature sequences

Xanthomonadales comprises one of the largest phytopathogenic bacterial groups, and is currently classified within the gamma-proteobacteria. However, the phylogenetic placement of this group is not clearly resolved, and the results of different studies contradict one another. In this work, the evolutionary position of Xanthomonadales was determined by analyzing the presence of shared insertions and deletions (INDELs) in highly conserved proteins. Several distinctive insertions found in most of the members of the gamma-proteobacteria are absent in Xanthomonadales and groups such as Legionelalles, Chromatiales, Methylococcales, Thiotrichales and Cardiobacteriales. These INDELs were most likely introduced after the branching of Xanthomonadales from most of the gamma-proteobacteria and provide evidence for the phylogenetic placement of the early gamma-proteobacteria. Moreover, other proteins contain insertions exclusive to the Xanthomonadales order, confirming that this is a monophyletic group and provide important specific genetic markers. Thus, the data presented clearly support the Xanthomonadales group as an independent subdivision, and constitute one of the deepest branching lineage within the gamma-proteobacteria clade.     

Structural role of conserved Asn179 in the short-chain dehydrogenase/reductase scaffold.

Short-chain dehydrogenases/reductases (SDR) constitute a large family of enzymes found in all forms of life. Despite a low level of sequence identity, the three-dimensional structures determined display a nearly superimposable alpha/beta folding pattern. We identified a conserved asparagine residue located within strand betaF and analyzed its role in the short-chain dehydrogenase/reductase architecture. Mutagenetic replacement of Asn179 by Ala in bacterial 3beta/17beta-hydroxysteroid dehydrogenase yields a folded, but enzymatically inactive enzyme, which is significantly more resistant to denaturation by guanidinium hydrochloride. Crystallographic analysis of the wild-type enzyme at 1.2-A resolution reveals a hydrogen bonding network, including a buried and well-ordered water molecule connecting strands betaE to betaF, a common feature found in 16 of 21 known three-dimensional structures of the family. Based on these results, we hypothesize that in mammalian 11beta-hydroxysteroid dehydrogenase the essential Asn-linked glycosylation site, which corresponds to the conserved segment, displays similar structural features and has a central role to maintain the SDR scaffold.     

Crystal structure of sorbitol dehydrogenase

Sorbitol dehydrogenase (SDH) is a distant relative to the alcohol dehydrogenases (ADHs) with sequence identities around 20%. SDH is a tetramer with one zinc ion per subunit. We have crystallized rat SDH and determined the structure by molecular replacement using a tetrameric bacterial ADH as search object. The conformation of the bound coenzyme is extended and similar to NADH bound to mammalian ADH but the interactions with the NMN-part have several differences with those of ADH. The active site zinc coordination in SDH is significantly different than in mammalian ADH but similar to the one found in the bacterial tetrameric NADP(H)-dependent ADH of Clostridiim beijerinckii. The substrate cleft is significantly more polar than for mammalian ADH and a number of residues are ideally located to position the sorbitol molecule in the active site. The SDH molecule can be considered to be a dimer of dimers, with subunits A-B and C-D, where the dimer interactions are similar to those in mammalian ADH. The tetramers are composed of two of these dimers, which interact with their surfaces opposite the active site clefts, which are accessible on the opposite side. In contrast to the dimer interactions, the tetramer-forming interactions are small with only few hydrogen bonds between side-chains.     

Characterization of a novel airway epithelial cell-specific short chain alcohol dehydrogenase/reductase gene whose expression is up-regulated by retinoids and is involved in the metabolism of retinol

Multiple retinoic acid responsive cDNAs were isolated from a high density cDNA microarray membrane, which was developed from a cDNA library of human tracheobronchial epithelial cells. Five selected cDNA clones encoded the sequence of the same novel gene. The predicted open reading frame of the novel gene encoded a protein of 319 amino acids. The deduced amino acid sequence contains four motifs that are conserved in the short-chain alcohol dehydrogenase/reductase (SDR) family of proteins. The novel gene shows the greatest homology to a group of dehydrogenases that can oxidize retinol (retinol dehydrogenases). The mRNA of the novel gene was found in trachea, colon, tongue, and esophagus. In situ hybridization of airway tissue sections demonstrated epithelial cell-specific gene expression, especially in the ciliated cell type. Both all-trans-retinoic acid and 9-cis-retinoic acid were able to elevate the expression of the novel gene in primary human tracheobronchial epithelial cells in vitro. This elevation coincided with an enhanced retinol metabolism in these cultures. COS cells transfected with an expression construct of the novel gene were also elevated in the metabolism of retinol. The results suggested that the novel gene represents a new member of the SDR family that may play a critical role in retinol metabolism in airway epithelia as well as in other epithelia of colon, tongue, and esophagus.     

Remarkable ancient divergences amongst neglected lorisiform primates

Lorisiform primates (Primates: Strepsirrhini: Lorisiformes) represent almost 10% of the living primate species and are widely distributed in sub‐Saharan Africa and South/South‐East Asia; however, their taxonomy, evolutionary history, and biogeography are still poorly understood. In this study we report the largest molecular phylogeny in terms of the number of represented taxa. We sequenced the complete mitochondrial cytochrome b gene for 86 lorisiform specimens, including ～80% of all the species currently recognized. Our results support the monophyly of the Galagidae, but a common ancestry of the Lorisinae and Perodicticinae (family Lorisidae) was not recovered. These three lineages have early origins, with the Galagidae and the Lorisinae diverging in the Oligocene at about 30 Mya and the Perodicticinae emerging in the early Miocene. Our mitochondrial phylogeny agrees with recent studies based on nuclear data, and supports Euoticus as the oldest galagid lineage and the polyphyletic status of Galagoides. Moreover, we have elucidated phylogenetic relationships for several species never included before in a molecular phylogeny. The results obtained in this study suggest that lorisiform diversity remains substantially underestimated and that previously unnoticed cryptic diversity might be present within many lineages, thus urgently requiring a comprehensive taxonomic revision of this primate group. © 2015 The Linnean Society of London     

Crystal structure of sorbitol dehydrogenase

Sorbitol dehydrogenase (SDH) is a distant relative to the alcohol dehydrogenases (ADHs) with sequence identities around 20%. SDH is a tetramer with one zinc ion per subunit. We have crystallized rat SDH and determined the structure by molecular replacement using a tetrameric bacterial ADH as search object. The conformation of the bound coenzyme is extended and similar to NADH bound to mammalian ADH but the interactions with the NMN-part have several differences with those of ADH. The active site zinc coordination in SDH is significantly different than in mammalian ADH but similar to the one found in the bacterial tetrameric NADP(H)-dependent ADH of Clostridiim beijerinckii. The substrate cleft is significantly more polar than for mammalian ADH and a number of residues are ideally located to position the sorbitol molecule in the active site. The SDH molecule can be considered to be a dimer of dimers, with subunits A–B and C–D, where the dimer interactions are similar to those in mammalian ADH. The tetramers are composed of two of these dimers, which interact with their surfaces opposite the active site clefts, which are accessible on the opposite side. In contrast to the dimer interactions, the tetramer-forming interactions are small with only few hydrogen bonds between side-chains.     

Molecular evidence of sorbitol dehydrogenase in tomato, a non-Rosaceae plant

The enzyme NAD-dependent sorbitol dehydrogenase (SDH) is well characterized in the Rosaceae family of fruit trees, which synthesizes sorbitol as a translocatable photosynthate. Expressed sequence tags of SDH-like sequences have also been generated from various non-Rosaceae species that do not synthesize sorbitol as a primary photosynthetic product, but the physiological roles of the encoded proteins in non-Rosaceae plants are unknown. Therefore, we isolated an SDH-like cDNA (SDL) from tomato (Lycopersicon esculentum Mill.). Genomic Southern blot analysis suggested that SDL exists in the tomato genome as a single-copy gene. Northern blot analysis showed that SDL is ubiquitously expressed in tomato plants. Recombinant SDL protein was produced and purified for enzymatic characterization. SDL catalyzed the interconversion of sorbitol and fructose with NAD (H). SDL showed highest activity for sorbitol among the several substrates tested. SDL showed no activity with NADP+. Thus, SDL was identified as a SDH, although the Km values and substrate specificity of SDL were significantly different from those of SDH purified from the Japanese pear (Pyrus pyrifolia), a Rosaceae fruit tree. In addition, tomato was transformed with antisense SDL to evaluate the contribution of SDL to SDH activity in tomato. The transformation decreased SDH activity to approximately 50% on average. Taken together, these results provide molecular evidence of SDH in tomato, and SDL was renamed LeSDH.     

Cloning and characterization of a novel mouse short-chain dehydrogenase/reductases cDNA mHsdl2#, encoding a protein with a SDR domaid and a SCP2 domain

The short-chain dehydrogenases/reductases (SDRs) play important roles in body's metabolism. We cloned a novel mouse SDR cDNA which encodes a deduced HSD-like protein with a conserved SDR domain and a SCP2 domain. The 1.8 kb cDNA consists of 11 exons and is mapped to mouse chromosome 4B3. The corresponding gene is widely expressed in normal mouse tissues and its expression level in liver increases after inducement with cholesterol food. The predicted mouse HSDL2 protein, which has a peroxisomal target signal, is localized in the cytoplasm of NIH 3T3 cells.     

Widespread occurrence of non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase among gram-positive bacteria.

The non-phosphorylating glyceraldehyde 3-phosphate dehydrogenase (GAPDHN, NADP+-specific, EC 1.2.1.9) is present in green eukaryotes and some Streptococcus strains. The present report describes the results of activity and immunoblot analyses, which were used to generate the first survey of bacterial GAPDHN distribution in a number of Bacillus, Streptococcus and Clostridium strains. Putative gapN genes were identified after PCR amplification of partial 700-bp sequences using degenerate primers constructed from highly conserved protein regions. Alignment of the amino acid sequences of these fragments with those of known sequences from other eukaryotic and prokaryotic GAPDHNs, demonstrated the presence of conserved residues involved in catalytic activity that are not conserved in aldehyde dehydrogenases, a protein family closely linked to GAPDHNs. The results confirm that the basic structural features of the members of the GAPDHN family have been conserved throughout evolution and that no identity exists with phosphorylating GAPDHs. Furthermore, phylogenetic trees generated from multiple sequence alignments suggested a close relationship between plant and bacterial GAPDHN families.     

Phylogenetic analysis of gamma-proteobacteria inferred from nucleotide sequence comparisons of the house-keeping genes adk, aroE and gdh: comparisons with phylogeny inferred from 16S rRNA gene sequences

Nucleotide sequence comparisons of three house-keeping genes, adenylate kinase (adk), shikimate dehydrogenase (aroE), and glucose-6-phosphate dehydrogenase (gdh), were used to infer the phylogeny of 33 gamma-proteobacteria. Phylogenetic trees inferred from each gene, and from the concatenated sequences of all three genes, are, in general, similar to a 16S rRNA gene-inferred tree. Similar grouping of bacteria are revealed at the family, genus, species and strain levels in all five trees. The house-keeping genes, however, show a higher rate of nucleotide sequence substitutions. Consequently, they can possibly probe deeper branches of a phylogenetic tree than the 16S rRNA gene. However, because their nucleotide sequences are not as highly conserved among gamma-proteobacteria, family- or genus-specific primers would need to be designed for the amplification of any of these three house-keeping genes. Since these genes are used in multilocus sequence typing, it is expected that the number of sequences publicly available for many taxa will increase over time proving them very useful either at complementing 16S rRNA-inferred phylogenies or for specific, targeted, phylogenetic analysis.     

Halohydrin dehalogenases are structurally and mechanistically related to short-chain dehydrogenases/reductases

Halohydrin dehalogenases, also known as haloalcohol dehalogenases or halohydrin hydrogen-halide lyases, catalyze the nucleophilic displacement of a halogen by a vicinal hydroxyl function in halohydrins to yield epoxides. Three novel bacterial genes encoding halohydrin dehalogenases were cloned and expressed in Escherichia coli, and the enzymes were shown to display remarkable differences in substrate specificity. The halohydrin dehalogenase of Agrobacterium radiobacter strain AD1, designated HheC, was purified to homogeneity. The k(cat) and K(m) values of this 28-kDa protein with 1,3-dichloro-2-propanol were 37 s(-1) and 0.010 mM, respectively. A sequence homology search as well as secondary and tertiary structure predictions indicated that the halohydrin dehalogenases are structurally similar to proteins belonging to the family of short-chain dehydrogenases/reductases (SDRs). Moreover, catalytically important serine and tyrosine residues that are highly conserved in the SDR family are also present in HheC and other halohydrin dehalogenases. The third essential catalytic residue in the SDR family, a lysine, is replaced by an arginine in halohydrin dehalogenases. A site-directed mutagenesis study, with HheC as a model enzyme, supports a mechanism for halohydrin dehalogenases in which the conserved Tyr145 acts as a catalytic base and Ser132 is involved in substrate binding. The primary role of Arg149 may be lowering of the pK(a) of Tyr145, which abstracts a proton from the substrate hydroxyl group to increase its nucleophilicity for displacement of the neighboring halide. The proposed mechanism is fundamentally different from that of the well-studied hydrolytic dehalogenases, since it does not involve a covalent enzyme-substrate intermediate.     

Cloning, nucleotide sequence, and overexpression of smoS, a component of a novel operon encoding an ABC transporter and polyol dehydrogenases of Rhodobacter sphaeroides Si4.

The gene coding for sorbitol dehydrogenase (SDH) of Rhodobacter sphaeroides Si4 was located 55 nucleotides upstream of the mannitol dehydrogenase gene (mtlK) within a previously unrecognized polyol operon. This operon probably consists of all the proteins necessary for transport and metabolization of various polyols. The gene encoding SDH (smoS) was cloned and sequenced. Analysis of the deduced amino acid sequence revealed homology to enzymes of the short-chain dehydrogenase/reductase protein family. For structure analysis of this unique bacterial enzyme, smoS was subcloned into the overexpression vector pET-24a(+) and then overproduced in Escherichia coli BL21(DE3), which yielded a specific activity of 24.8 U/mg of protein and a volumetric yield of 38,000 U/liter. Compared to values derived with the native host, R. sphaeroides, these values reflected a 270-fold increase in expression of SDH and a 971-fold increase in the volumetric yield. SDH was purified to homogeneity, with a recovery of 49%, on the basis of a three-step procedure. Upstream from smoS, another gene (smoK), which encoded a putative ATP-binding protein of an ABC transporter, was identified.     

A phylogenomic analysis of the shikimate dehydrogenases reveals broadscale functional diversification and identifies one functionally distinct subclass

The shikimate dehydrogenases (SDH) represent a widely distributed enzyme family with an essential role in secondary metabolism. This superfamily had been previously subdivided into 4 enzyme groups (AroE, YdiB, SdhL, and RifI), which show clear biochemical and functional differences ranging from amino acid biosynthesis to antibiotic production. Despite the importance of this group, little is known about how such essential enzymatic functions can evolve and diversify. We dissected the enzyme superfamily with a phylogenomic analysis of approximately 250 fully sequenced genomes, making use of previously characterized representatives from each enzyme class, and the key substrate-binding residues known to distinguish substrate specificity. We identified 5 major evolutionary and functional SDH subgroups and several other potentially unique functional classes within this complex enzyme family and then validated the functional distinctiveness of each group by characterizing the 5 SDH homologs found in Pseudomonas putida KT2440 biochemically. We identified an entirely novel functionally distinct subgroup, which we designated Ael1 (AroE-like1) and also delineated a new group of shikimate/quinate dehydrogenases (YdiB2), which is phylogenetically distinct from the previously described Escherichia coli YdiB. The combination of biochemical, phylogenetic, and genomic approaches has revealed the broad extent to which the SDH enzyme superfamily has diversified. Five functional groups were validated with the potential for at least 5 additional subgroups. Our analysis also identified a new SDH functional group, which appears to have evolved recently from an ancestral AroE, illustrating a very prominent role of horizontal transmission and neofunctionalizaton in the evolutionary and functional diversification of this enzyme family.     

Phylogenetic analyses of proton-translocating transhydrogenases

The proton-translocating nicotinamide nucleotide transhydrogenases (TH) provide a simple model for understanding chemically coupled transmembrane proton translocation. To further our understanding of TH structure-function relationships, we have identified all sequenced homologous of these vectorial enzymes and have conducted sequence comparison studies. The NAD-binding domains of TH are homologous to bacterial alanine dehydrogenases (ADH) and eukaryotic saccharopine dehydrogenases (SDH) as well as N5(carboxyethyl)-L-ornithine synthase of Lactococcus lactis and dipicolinate synthase of Bacillus subtilis. A multiple alignment, a phylogenetic tree, and two signature sequences for this family, designated the TH-ADH-SDH or TAS superfamily, have been derived. Additionally, the TH family has been characterized. Phylogenetic analyses suggested that these proteins have evolved without inter-system shuffling. However, interdomain splicing-fusion events have occurred during the evolution of several of these systems. Analyses of the multiple alignment for the TH family revealed that domain conservation occurs in the order: NADP-binding domain (domain III) > NAD-binding domain (domain I) > proton-translocating transmembrane domain (domain II). A topologic model for the proton-translocating transmembrane domain consistent with published data is presented, and a possible involvement of specific transmembrane alpha-helical segments in channel formation is suggested.     

Cloning and Characterization of a Novel Mouse Short-Chain Dehydrogenases/Reductases cDNA, mHsdl2, Encoding a Protein with an SDR Domain and an SCP2 Domain

The short-chain dehydrogenases/reductases (SDRs) play an important role in the body's metabolism. We have cloned a novel mouse SDR cDNA, which encodes a deduced HSD-like protein with a conserved SDR domain and an SCP2 domain. The 1.8 kb cDNA consists of 11 exons and is mapped to mouse chromosome 4B3. The corresponding gene is widely expressed in normal mouse tissues and its expression level in the liver increases after inducement with cholesterol food. The predicted mouse HSDL2 protein, which has a peroxisomal target signal, is localized in the cytoplasm of NIH 3T3 cells.     

Molecular basis for thermoprotection in Bemisia: structural differences between whitefly ketose reductase and other medium-chain dehydrogenases/reductases

The silverleaf whitefly (Bemisia argentifolii, Bellows and Perring) accumulates sorbitol as a thermoprotectant in response to elevated temperature. Sorbitol synthesis in this insect is catalyzed by an unconventional ketose reductase (KR) that uses NADPH to reduce fructose. A cDNA encoding the NADPH-KR from adult B. argentifolii was cloned and sequenced to determine the primary structure of this enzyme. The cDNA encoded a protein of 352 amino acids with a calculated molecular mass of 38.2 kDa. The deduced amino acid sequence of the cDNA shared 60% identity with sheep NAD(+)-dependent sorbitol dehydrogenase (SDH). Residues in SDH involved in substrate binding were conserved in the whitefly NADPH-KR. An important structural difference between the whitefly NADPH-KR and NAD(+)-SDHs occurred in the nucleotide-binding site. The Asp residue that coordinates the adenosyl ribose hydroxyls in NAD(+)-dependent dehydrogenases (including NAD(+)-SDH), was replaced by an Ala in the whitefly NADPH-KR. The whitefly NADPH-KR also contained two neutral to Arg substitutions within four residues of the Asp to Ala substitution. Molecular modeling indicated that addition of the Arg residues and loss of the Asp decreased the electric potential of the adenosine ribose-binding pocket, creating an environment favorable for NADPH-binding. Because of the ability to use NADPH, the whitefly NADPH-KR synthesizes sorbitol under physiological conditions, unlike NAD(+)-SDHs, which function in sorbitol catabolism.     

The Vein Patterning 1 (VEP1) gene family laterally spread through an ecological network

Lateral gene transfer (LGT) is a major evolutionary mechanism in prokaryotes. Knowledge about LGT--particularly, multicellular--eukaryotes has only recently started to accumulate. A widespread assumption sees the gene as the unit of LGT, largely because little is yet known about how LGT chances are affected by structural/functional features at the subgenic level. Here we trace the evolutionary trajectory of VEin Patterning 1, a novel gene family known to be essential for plant development and defense. At the subgenic level VEP1 encodes a dinucleotide-binding Rossmann-fold domain, in common with members of the short-chain dehydrogenase/reductase (SDR) protein family. We found: i) VEP1 likely originated in an aerobic, mesophilic and chemoorganotrophic α-proteobacterium, and was laterally propagated through nets of ecological interactions, including multiple LGTs between phylogenetically distant green plant/fungi-associated bacteria, and five independent LGTs to eukaryotes. Of these latest five transfers, three are ancient LGTs, implicating an ancestral fungus, the last common ancestor of land plants and an ancestral trebouxiophyte green alga, and two are recent LGTs to modern embryophytes. ii) VEP1's rampant LGT behavior was enabled by the robustness and broad utility of the dinucleotide-binding Rossmann-fold, which provided a platform for the evolution of two unprecedented departures from the canonical SDR catalytic triad. iii) The fate of VEP1 in eukaryotes has been different in different lineages, being ubiquitous and highly conserved in land plants, whereas fungi underwent multiple losses. And iv) VEP1-harboring bacteria include non-phytopathogenic and phytopathogenic symbionts which are non-randomly distributed with respect to the type of harbored VEP1 gene. Our findings suggest that VEP1 may have been instrumental for the evolutionary transition of green plants to land, and point to a LGT-mediated 'Trojan Horse' mechanism for the evolution of bacterial pathogenesis against plants. VEP1 may serve as tool for revealing microbial interactions in plant/fungi-associated environments.