A genome-wide phylogenetic reconstruction of family 1 UDP-glycosyltransferases revealed the expansion of the family during the adaptation of plants to life on land

For almost a decade, our knowledge on the organisation of the family 1 UDP‐glycosyltransferases (UGTs) has been limited to the model plant A. thaliana. The availability of other plant genomes represents an opportunity to obtain a broader view of the family in terms of evolution and organisation. Family 1 UGTs are known to glycosylate several classes of plant secondary metabolites. A phylogeny reconstruction study was performed to get an insight into the evolution of this multigene family during the adaptation of plants to life on land. The organisation of the UGTs in the different organisms was also investigated. More than 1500 putative UGTs were identified in 12 fully sequenced and assembled plant genomes based on the highly conserved PSPG motif. Analyses by maximum likelihood (ML) method were performed to reconstruct the phylogenetic relationships existing between the sequences. The results of this study clearly show that the UGT family expanded during the transition from algae to vascular plants and that in higher plants the clustering of UGTs into phylogenetic groups appears to be conserved, although gene loss and gene gain events seem to have occurred in certain lineages. Interestingly, two new phylogenetic groups, named O and P, that are not present in A. thaliana were discovered.     

An evolutionary view of functional diversity in family 1 glycosyltransferases

Glycosyltransferases (GTs) (EC 2.4.x.y) catalyze the transfer of sugar moieties to a wide range of acceptor molecules, such as sugars, lipids, proteins, nucleic acids, antibiotics and other small molecules, including plant secondary metabolites. These enzymes can be classified into at least 92 families, of which family 1 glycosyltransferases (GT1), often referred to as UDP glycosyltransferases (UGTs), is the largest in the plant kingdom. To understand how UGTs expanded in both number and function during evolution of land plants, we screened genome sequences from six plants (Physcomitrella patens, Selaginella moellendorffii, Populus trichocarpa, Oryza sativa, Arabidopsis thaliana and Arabidopsis lyrata) for the presence of a conserved UGT protein domain. Phylogenetic analyses of the UGT genes revealed a significant expansion of UGTs, with lineage specificity and a higher duplication rate in vascular plants after the divergence of Physcomitrella. The UGTs from the six species fell into 24 orthologous groups that contained genes derived from the common ancestor of these six species. Some orthologous groups contained multiple UGT families with known functions, suggesting that UGTs discriminate compounds as substrates in a lineage-specific manner. Orthologous groups containing only a single UGT family tend to play a crucial role in plants, suggesting that such UGT families may have not expanded because of evolutionary constraints.     

Phylogenetic analysis of the UDP-glycosyltransferase multigene family of Arabidopsis thaliana

A class of UDP-glycosyltransferases (UGTs) defined by the presence of a C-terminal consensus sequence is found throughout the plant and animal kingdoms. Whereas mammalian enzymes use UDP-glucuronic acid, the plant enzymes typically use UDP-glucose in the transfer reactions. A diverse array of aglycones can be glucosylated by these UGTs. In plants, the aglycones include plant hormones, secondary metabolites involved in stress and defense responses, and xenobiotics such as herbicides. Glycosylation is known to regulate many properties of the aglycones such as their bioactivity, their solubility, and their transport properties within the cell and throughout the plant. As a means of providing a framework to start to understand the substrate specificities and structure-function relationships of plant UGTs, we have now applied a molecular phylogenetic analysis to the multigene family of 99 UGT sequences in Arabidopsis. We have determined the overall organization and evolutionary relationships among individual members with a surprisingly high degree of confidence. Through constructing a composite phylogenetic tree that also includes all of the additional plant UGTs with known catalytic activities, we can start to predict both the evolutionary history and substrate specificities of new sequences as they are identified. The tree already suggests that while the activities of some subgroups of the UGT family are highly conserved among different plant species, others subgroups shift substrate specificity with relative ease.     

Molecular evolutionary analyses of the Arabidopsis L7 ribosomal protein gene family

Cytoplasmic ribosomal protein (r-protein) genes in Arabidopsis thaliana are encoded by 80 multigene families that contain between two and seven members. Gene family members are typically similar at the protein sequence level, with the most divergent members of any gene family retaining 94% identity, on average. However, three Arabidopsis r-protein families - S15a, L7 and P2 - contain highly divergent family members. Here, we investigated the organization, structure, expression and molecular evolution of the L7 r-protein family. Phylogenetic analyses showed that L7 r-protein gene family members constitute two distinct phylogenetic groups. The first group including RPL7B, RPL7C and RPL7D has homologs in plants, animals and fungi. The second group represented by RPL7A is found in plants but has no orthologs from other fully-sequenced eukaryotic genomes. These two groups may have derived from a duplication event prior to the divergence of animals and plants. All four L7 r-protein genes are expressed and all exhibit a differential expression in inflorescence and flowers. RPL7A and RPL7B are less expressed than the other genes in all tissues analyzed. Molecular characterization of nucleic and protein sequences of L7 r-protein genes and analysis of their codon usage did not indicate any functional divergence. The probable evolution of an extra-ribosomal function of group 2 genes is discussed.     

Pathogen-responsive expression of glycosyltransferase genes UGT73B3 and UGT73B5 is necessary for resistance to Pseudomonas syringae pv tomato in Arabidopsis

The genome sequencing of Arabidopsis (Arabidopsis thaliana) has revealed that secondary metabolism plant glycosyltransferases (UGTs) are encoded by an unexpectedly large multigenic family of 120 members. Very little is known about their actual function in planta, in particular during plant pathogen interactions. Among them, members of the group D are of particular interest since they are related to UGTs involved in stress-inducible responses in other plant species. We provide here a detailed analysis of the expression profiles of this group of Arabidopsis UGTs following infection with Pseudomonas syringae pv tomato or after treatment with salicylic acid, methyljasmonate, and hydrogen peroxide. Members of the group D displayed distinct induction profiles, indicating potential roles in stress or defense responses notably for UGT73B3 and UGT73B5. Analysis of UGT expression in Arabidopsis defense-signaling mutants further revealed that their induction is methyljasmonate independent, but partially salicylic acid dependent. T-DNA tagged mutants (ugt73b3 and ugt73b5) exhibited decreased resistance to P. syringae pv tomato-AvrRpm1, indicating that expression of the corresponding UGT genes is necessary during the hypersensitive response. These results emphasize the importance of plant secondary metabolite UGTs in plant-pathogen interactions and provide foundation for future understanding of the exact role of UGTs during the hypersensitive response.     

Diversification of NRT2 and the origin of its fungal homolog

We investigated the origin and diversification of the high-affinity nitrate transporter NRT2 in fungi and other eukaryotes using Bayesian and maximum parsimony methods. To assess the higher-level relationships and origins of NRT2 in eukaryotes, we analyzed 200 amino acid sequences from the Nitrate/Nitrite Porter (NNP) Family (to which NRT2 belongs), including 55 fungal, 41 viridiplantae (green plants), 11 heterokonts (stramenopiles), and 87 bacterial sequences. To assess evolution of NRT2 within fungi and other eukaryotes, we analyzed 116 amino acid sequences of NRT2 from 58 fungi, 40 viridiplantae (green plants), 1 rhodophyte, and 5 heterokonts, rooted with 12 bacterial sequences. Our results support a single origin of eukaryotic NRT2 from 1 of several clades of mostly proteobacterial NNP transporters. The phylogeny of bacterial NNP transporters does not directly correspond with bacterial taxonomy, apparently due to ancient duplications and/or horizontal gene transfer events. The distribution of NRT2 in the eukaryotes is patchy, but the NRT2 phylogeny nonetheless supports the monophyly of major groups such as viridiplantae, flowering plants, monocots, and eudicots, as well as fungi, ascomycetes, basidiomycetes, and agaric mushrooms. At least 1 secondary origin of eukaryotic NRT2 via horizontal transfer to the fungi is suggested, possibly from a heterokont donor. Our analyses also suggest that there has been a horizontal transfer of nrt2 from a basidiomycete fungus to an ascomycete fungus and reveal a duplication of nrt2 in the ectomycorrhizal mushroom genus, Hebeloma.     

Substrate specificity of plant UDP-dependent glycosyltransferases predicted from crystal structures and homology modeling

Plant family 1 UDP-dependent glycosyltransferases (UGTs) catalyze the glycosylation of a plethora of bioactive natural products. In Arabidopsis thaliana, 120 UGT encoding genes have been identified. The crystal-based 3D structures of four plant UGTs have recently been published. Despite low sequence conservation, the UGTs show a highly conserved secondary and tertiary structure. The sugar acceptor and sugar donor substrates of UGTs are accommodated in the cleft formed between the N- and C-terminal domains. Several regions of the primary sequence contribute to the formation of the substrate binding pocket including structurally conserved domains as well as loop regions differing both with respect to their amino acid sequence and sequence length. In this review we provide a detailed analysis of the available plant UGT crystal structures to reveal structural features determining substrate specificity. The high 3D structural conservation of the plant UGTs render homology modeling an attractive tool for structure elucidation. The accuracy and utility of UGT structures obtained by homology modeling are discussed and quantitative assessments of model quality are performed by modeling of a plant UGT for which the 3D crystal structure is known. We conclude that homology modeling offers a high degree of accuracy. Shortcomings in homology modeling are also apparent with modeling of loop regions remaining as a particularly difficult task.     

A multigene family of glycosyltransferases in a model plant,Arabidopsis thaliana

Glycosyltransferases transfer sugars from NDP-sugar donors to acceptors. The multigene family of transferases described in this paper typically transfer glucose from UDP-glucose to low-molecular-mass acceptors in the cytosol of plant cells. There are 107 sequences in the genome of Arabidopsis thaliana that contain a consensus, suggesting they belong to this Group 1 multigene family. The family has been analysed phylogenetically, and a functional genomics approach has been applied to explore the relatedness of sequence similarity to catalytic specificity and stereoselectivity. Enzymes belonging to this class of transferases glycosylate a vast array of acceptors, including natural products such as secondary metabolites and hormones, as well as xenobiotics absorbed by the plant, such as herbicides and pesticides. Conjugation to glucose potentially changes the activity of the acceptor molecule and invariably changes its location within the plant cell. Using the genomics approach described, a platform of knowledge has been constructed that will enable an understanding to be gained on the role of these enzymes in cellular homoeostasis, as well as their activity in biotransformations in vitro that require strict regioselectivity of glycosylation.     

Was the ANITA rooting of the angiosperm phylogeny affected by long-branch attraction? Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya

Five groups of basal angiosperms, Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya (ANITA), were identified in several recent studies as representing a series of the earliest-diverging lineages of the angiosperm phylogeny. All of these studies except one employed a multigene analysis approach and used gymnosperms as the outgroup to determine the ingroup topology. The high level of divergence between gymnosperms and angiosperms, however, has long been implicated in the difficulty of reconstructing relationships at the base of angiosperm phylogeny using DNA sequences, for fear of long-branch attraction (LBA). In this study, we replaced the gymnosperm sequences from the five-gene matrix (mitochondrial atp1 and matR, plastid atpB and rbcL, and nuclear 18S rDNA) used in our earlier study with four categories of divergent sequences--random sequences with equal base frequencies or equally AT- and GC-rich contents, homopolymers and heteropolymers, misaligned gymnosperm sequences, and aligned lycopod and bryophyte sequences--to evaluate whether the gymnosperms were an appropriate outgroup to angiosperms in our earlier study that identified the ANITA rooting. All 24 analyses performed rooted the angiosperm phylogeny at either Acorus or Alisma (or Alisma-Triglochin-Potamogeton in one case due to use of a slightly different alignment) and placed the monocots as a basal grade, producing genuine LBA results. These analyses demonstrate that the identification of ANITA as the basalmost extant angiosperms was based on historical signals preserved in the gymnosperm sequences and that the gymnosperms were an appropriate outgroup with which to root the angiosperm phylogeny in the multigene sequence analysis. This strategy of evaluating the appropriateness of an outgroup using artificial sequences and a series of outgroups with increments of divergence levels can be applied to investigations of phylogenetic patterns at the bases of other major clades, such as land plants, animals, and eukaryotes.     

Evolutionary analysis of glycosyl hydrolase family 28 (GH28) suggests lineage-specific expansions in necrotrophic fungal pathogens

Glycosyl hydrolase family 28 (GH28) is a set of structurally related enzymes that hydrolyze glycosidic bonds in pectin, and are important extracellular enzymes for both pathogenic and saprotrophic fungi. Yet, very little is understood about the evolutionary forces driving the diversification of GH28s in fungal genomes. We reconstructed the evolutionary history of family GH28 in fungi by examining the distribution of GH28 copy number across the phylogeny of fungi, and by reconstructing the phylogeny of GH28 genes. We also examined the relationship between lineage-specific GH28 expansions and fungal ecological strategy, testing the hypothesis that GH28 evolution in fungi is driven by ecological strategy (pathogenic vs. non-pathogenic) and pathogenic niche (necrotrophic vs. biotrophic). Our results showed that GH28 phylogeny of Ascomycota and Basidiomycota sequences was structured by specific biochemical function, with endo-polygalacturonases and endo-rhamnogalacturonases forming distinct, apparently ancient clades, while exo-polygalacturonases are more widely distributed. In contrast, Mucoromycotina and Stramenopile sequences formed taxonomically-distinct clades. Large, lineage-specific variation in GH28 copy number indicates that the evolution of this gene family is consistent with the birth-and-death model of gene family evolution, where diversity of GH28 loci within genomes was generated through multiple rounds of gene duplication followed by functional diversification and loss of some gene family members. Although GH28 copy number was correlated with genome size, our findings suggest that ecological strategy also plays an important role in determining the GH28 repertoire of fungi. Both necrotrophic and biotrophic fungi have larger genomes than non-pathogens, yet only necrotrophs possess more GH28 enzymes than non-pathogens. Hence, lineage-specific GH28 expansion is the result of both variation in genome size across fungal species and diversifying selection within the necrotrophic plant pathogen ecological niche. GH28 evolution among necrotrophs has likely been driven by a co-evolutionary arms race with plants, whereas the need to avoid plant immune responses has resulted in purifying selection within biotrophic fungi.     

Monophyly of beta-tubulin and H+-ATPase gene variants in Glomus intraradices: consequences for molecular evolutionary studies of AM fungal genes

Arbuscular mycorrhizal fungi (AMF) are an ecologically important group of fungi. Previous studies showed the presence of divergent copies of beta-tubulin and V-type vacuolar H+-ATPase genes in AMF genomes and suggested horizontal gene transfer from host plants or mycoparasites to AMF. We sequenced these genes from DNA isolated from an in vitro cultured isolate of Glomus intraradices that was free of any obvious contaminants. We found two highly variable beta-tubulin sequences and variable H+-ATPase sequences. Despite this high variation, comparison of the sequences with those in gene banks supported a glomeromycotan origin of G. intraradices beta-tubulin and H+-ATPase sequences. Thus, our results are in sharp contrast with the previously reported polyphyletic origin of those genes. We present evidence that some highly divergent sequences of beta-tubulin and H+-ATPase deposited in the databases are likely to be contaminants. We therefore reject the prediction of horizontal transfer to AMF genomes. High differences in GC content between glomeromycotan sequences and sequences grouping in other lineages are shown and we suggest they can be used as an indicator to detect such contaminants. H+-ATPase phylogeny gave unexpected results and failed to resolve fungi as a natural group. beta-Tubulin phylogeny supported Glomeromeromycota as sister group of the Chytridiomycota. Contrasts between our results and trees previously generated using rDNA sequences are discussed.     

Phylogenetic relationships in class I of the superfamily of bacterial, fungal, and plant peroxidases.

Molecular phylogeny among catalase-peroxidases, cytochrome c peroxidases, and ascorbate peroxidases was analysed. Sixty representative sequences covering all known subgroups of class I of the superfamily of bacterial, fungal, and plant heme peroxidases were selected. Each sequence analysed contained the typical peroxidase motifs evolved to bind effectively the prosthetic heme group, enabling peroxidatic activity. The N-terminal and C-terminal domains of catalase-peroxidases matching the ancestral tandem gene duplication event were treated separately in the phylogenetic analysis to reveal their specific evolutionary history. The inferred unrooted phylogenetic tree obtained by three different methods revealed the existence of four clearly separated clades (C-terminal and N-terminal domains of catalase-peroxidases, ascorbate peroxidases, and cytochrome c peroxidases) which were segregated early in the evolution of this superfamily. From the results, it is obvious that the duplication event in the gene for catalase-peroxidase occurred in the later phase of evolution, in which the individual specificities of the peroxidase families distinguished were already formed. Evidence is presented that class I of the heme peroxidase superfamily is spread among prokaryotes and eukaryotes, obeying the birth-and-death process of multigene family evolution.     

Methionine adenosyltransferase as a useful molecular systematics tool revealed by phylogenetic and structural analyses

Structural and phylogenetic relationships among Bacteria and Eukaryota were analyzed by examining 292 methionine adenosyltransferase (MAT) amino acid sequences with respect to the crystal structure of this enzyme established for Escherichia coli and rat liver. Approximately 30% of MAT residues were found to be identical in all species. Five highly conserved amino acid sequence blocks did not vary in the MAT family. We detected specific structural features that correlated with sequence signatures for several clades, allowing taxonomical identification by sequence analysis. In addition, the number of amino acid residues in the loop connecting beta-strands A2 and A3 served to clearly distinguish sequences between eukaryotes and eubacteria. The molecular phylogeny of MAT genes in eukaryotes can be explained in terms of functional diversification coupled to gene duplication or alternative splicing and adaptation through strong structural constraints. Sequence analyses and intron/exon junction positions among nematodes, arthropods and vertebrates support the traditional Coelomata hypothesis. In vertebrates, the liver MAT I isoenzyme has gradually adapted its sequence towards one providing a more specific liver function. MAT phylogeny also served to cluster the major bacterial groups, demonstrating the superior phylogenetic performance of this ubiquitous, housekeeping gene in reconstructing the evolutionary history of distant relatives.     

Functional genomic analysis of Arabidopsis thaliana glycoside hydrolase family 35

Catalysing the hydrolysis of terminal beta-galactosyl residues from carbohydrates, galactolipids, and glycoproteins, glycoside hydrolase family 35 (beta-galactosidases; BGALs) are widely distributed in plants and believed to play many key roles, including modification of cell wall components. Completion of the Arabidopsis thaliana genome sequencing project has, for the first time, allowed an examination of the total number, gene structure, and evolutionary patterns of all Family 35 members in a representative (model) angiosperm. Reiterative database searches established a multigene family of 17 members (designated BGAL1-BGAL17). Using these genes as query sequences, BLAST and Hidden Markov Model searches identified BGAL genes among 22 other eukaryotes, whose genomic sequences are known. The Arabidopsis (n=17) and rice (n=15) BGAL families were much larger than those of Chlamydomonas, fungi, and animals (n=0-4), and a lineage-specific expansion of BGAL genes apparently occurred after divergence of the Arabidopsis and rice lineages. All plant BGAL genes, with the exception of Arabidopsis BGAL17 and rice Os 9633.m04334, form a monophyletic group. Arabidopsis BGAL expression levels are much higher in mature leaves, roots, flowers, and siliques but are lower in young seedlings. BGAL8, BGAL11, BGAL13, BGAL14, and BGAL16 are expressed only in flowers. Catalytically active BGAL4 was produced in the E. coli and baculoviral expression systems, purified to electrophoretic homogeneity, and partially characterized. The purified enzyme hydrolyzed p- and o-nitrophenyl-beta-d-galactosides. It also cleaved beta-(1,3)-, beta-(1,4)-, and beta-(1,6)-linked galactobiosides and galactotriosides, showing a marked preference for beta-(1,3)- and beta-(1,4)-linkages.     

Comparative analysis of HD2 type histone deacetylases in higher plants

Zea mays (L.) histone deacetylase HD2 was identified as a new type of histone deacetylase (HDAC) unrelated to the well-known Rpd3p and Hdalp families but with sequence homology to peptidyl-prolyl cis-trans isomerases (PPIases). Here we show that HD2 is a multigene family with highly related members in various plant species. Gene analysis revealed a similar exon/intron structure in Arabidopsis thaliana (L.) Heynh. and Z. mays, and most of the sequences analyzed were demonstrated to possess an intron of the very rare AT-AC type.