Molecular and functional analysis of phosphomannomutase (PMM) from higher plants and genetic evidence for the involvement of PMM in ascorbic acid biosynthesis in Arabidopsis and Nicotiana benthamiana

Phosphomannomutase (PMM) catalyzes the interconversion of mannose-6-phosphate and mannose-1-phosphate. However, systematic molecular and functional investigations on PMM from higher plants have hitherto not been reported. In this work, PMM cDNAs were isolated from Arabidopsis, Nicotiana benthamiana, soybean, tomato, rice and wheat. Amino acid sequence comparisons indicated that plant PMM proteins exhibited significant identity to their fungal and mammalian orthologs. In line with the similarity in primary structure, plant PMM complemented the sec53-6 temperature sensitive mutant of Saccharomyces cerevisiae. Histidine-tagged Arabidopsis PMM (AtPMM) purified from Escherichia coli converted mannose-1-phosphate into mannose-6-phosphate and glucose-1-phosphate into glucose-6-phosphate, with the former reaction being more efficient than the latter one. In Arabidopsis and N. benthamiana, PMM was constitutively expressed in both vegetative and reproductive organs. Reducing the PMM expression level through virus-induced gene silencing caused a substantial decrease in ascorbic acid (AsA) content in N. benthamiana leaves. Conversely, raising the PMM expression level in N. benthamiana using viral-vector-mediated ectopic expression led to a 20-50% increase in AsA content. Consistent with this finding, transgenic expression of an AtPMM-GFP fusion protein in Arabidopsis also increased AsA content by 25-33%. Collectively, this study improves our understanding on the molecular and functional properties of plant PMM and provides genetic evidence on the involvement of PMM in the biosynthesis of AsA in Arabidopsis and N. benthamiana plants.     

Evolutionary analysis for functional divergence of Jak protein kinase domains and tissue-specific genes

Jak (Janus kinase) is a nonreceptor tyrosine kinase, which plays important roles in signal transduction pathways. The unique feature of Jak is that, in addition to a fully functional tyrosine kinase domain (JH1), Jak possesses a pseudokinase domain (JH2). Although JH2 lost its catalytic function, experimental evidence has shown that this domain may have acquired some new but unknown functions. This apparent functional divergence after the (internal) domain duplication may result in dramatic changes of selective constraints at some sites. We conducted a data analysis to test this hypothesis. Our result shows that shifted selective constraints (or shifted evolutionary rates) between the JH1 and the JH2 domains are statistically significant. Predicted amino acid sites by posterior analysis can be classified into two groups: very conserved in JH1 but highly variable in JH2, and vice versa. Moreover, we have studied the evolutionary pattern of four tissue-specific genes, Jak1, Jak2, Jak3, and Tyk2, which were generated in the early stages of vertebrates. We found that after the (first) gene duplication, site-specific rate shifts between Jak2/Jak3 and Jak1/Tyk are significant, presumably as a consequence of functional divergence among these genes. The implication of our study for functional genomics is discussed.     

The normal phenotype of Pmm1-deficient mice suggests that Pmm1 is not essential for normal mouse development

Phosphomannomutases (PMMs) are crucial for the glycosylation of glycoproteins. In humans, two highly conserved PMMs exist: PMM1 and PMM2. In vitro both enzymes are able to convert mannose-6-phosphate (mannose-6-P) into mannose-1-P, the key starting compound for glycan biosynthesis. However, only mutations causing a deficiency in PMM2 cause hypoglycosylation, leading to the most frequent type of the congenital disorders of glycosylation (CDG): CDG-Ia. PMM1 is as yet not associated with any disease, and its physiological role has remained unclear. We generated a mouse deficient in Pmm1 activity and documented the expression pattern of murine Pmm1 to unravel its biological role. The expression pattern suggested an involvement of Pmm1 in (neural) development and endocrine regulation. Surprisingly, Pmm1 knockout mice were viable, developed normally, and did not reveal any obvious phenotypic alteration up to adulthood. The macroscopic and microscopic anatomy of all major organs, as well as animal behavior, appeared to be normal. Likewise, lectin histochemistry did not demonstrate an altered glycosylation pattern in tissues. It is especially striking that Pmm1, despite an almost complete overlap of its expression with Pmm2, e.g., in the developing brain, is apparently unable to compensate for deficient Pmm2 activity in CDG-Ia patients. Together, these data point to a (developmental) function independent of mannose-1-P synthesis, whereby the normal knockout phenotype, despite the stringent conservation in phylogeny, could be explained by a critical function under as-yet-unidentified challenge conditions.     

Evolution of human α1-acid glycoprotein genes and surrounding Alu repeats

There is a mosaic pattern of variation between the two tandemly arranged human α₁-acid glycoprotein genes. Both the synonymous and the nonsynonymous sites of exons 3 and 4 are more divergent than the rest of the gene, suggesting that they have had a different evolutionary history. Comparisons of the two gene sequences with rat AGP indicate that exons 3 and 4 of AGP2 have been evolving without functional constraint since their divergence from AGP1. It is proposed that the conserved region of the gene has been homogenized recently by gene conversion with the homologous regions of AGP1. The Alu sequences surrounding the genes appear to have been involved in both the gene duplication and the gene conversion events.     

A zebrafish model of PMM2-CDG reveals altered neurogenesis and a substrate-accumulation mechanism for N-linked glycosylation deficiency

Congenital disorder of glycosylation (PMM2-CDG) results from mutations in pmm2, which encodes the phosphomannomutase (Pmm) that converts mannose-6-phosphate (M6P) to mannose-1-phosphate (M1P). Patients have wide-spectrum clinical abnormalities associated with impaired protein N-glycosylation. Although it has been widely proposed that Pmm2 deficiency depletes M1P, a precursor of GDP-mannose, and consequently suppresses lipid-linked oligosaccharide (LLO) levels needed for N-glycosylation, these deficiencies have not been demonstrated in patients or any animal model. Here we report a morpholino-based PMM2-CDG model in zebrafish. Morphant embryos had developmental abnormalities consistent with PMM2-CDG patients, including craniofacial defects and impaired motility associated with altered motor neurogenesis within the spinal cord. Significantly, global N-linked glycosylation and LLO levels were reduced in pmm2 morphants. Although M1P and GDP-mannose were below reliable detection/quantification limits, Pmm2 depletion unexpectedly caused accumulation of M6P, shown earlier to promote LLO cleavage in vitro. In pmm2 morphants, the free glycan by-products of LLO cleavage increased nearly twofold. Suppression of the M6P-synthesizing enzyme mannose phosphate isomerase within the pmm2 background normalized M6P levels and certain aspects of the craniofacial phenotype and abrogated pmm2-dependent LLO cleavage. In summary, we report the first zebrafish model of PMM2-CDG and uncover novel cellular insights not possible with other systems, including an M6P accumulation mechanism for underglycosylation.     

Identification of the pgmG Gene, Encoding a Bifunctional Protein with Phosphoglucomutase and Phosphomannomutase Activities, in the Gellan Gum-Producing Strain Sphingomonas paucimobilis ATCC 31461

The pgmG gene of Sphingomonas paucimobilis ATCC 31461, the industrial gellan gum-producing strain, was cloned and sequenced. It encodes a 50,059-Da polypeptide that has phosphoglucomutase (PGM) and phosphomannomutase (PMM) activities and is 37 to 59% identical to other bifunctional proteins with PGM and PMM activities from gram-negative species, including Pseudomonas aeruginosa AlgC. Purified PgmG protein showed a marked preference for glucose-1-phosphate (G1P); the catalytic efficiency was about 50-fold higher for G1P than it was for mannose-1-phosphate (M1P). The estimated apparent K_m values for G1P and M1P were high, 0.33 and 1.27 mM, respectively. The pgmG gene allowed the recovery of alginate biosynthetic ability in a P. aeruginosa mutant with a defective algC gene. This result indicates that PgmG protein can convert mannose-6-phosphate into M1P in the initial steps of alginate biosynthesis and, together with other results, suggests that PgmG may convert glucose-6-phosphate into G1P in the gellan pathway.     

Origin and Evolution of Paralogous rRNA Gene Clusters Within the Flatworm Family Dugesiidae (Platyhelminthes, Tricladida)

Analysis of the 18S rDNA sequences of five species of the family Dugesiidae (phylum Platyhelminthes, suborder Tricladida, infraorder Paludicola) and eight species belonging to families Dendrocoelidae and Planaridae and to the infraorder Maricola showed that members of the family Dugesiidae have two types of 18S rDNA genes, while the rest of the species have only one. The duplication event also affected the ITS-1, 5.8S, ITS-2 region and probably the 28S gene. The mean divergence value between the type I and the type II sequences is 9% and type II 18S rDNA genes are evolving 2.3 times more rapidly than type I. The evolutionary rates of type I and type II genes were calibrated from biogeographical data, and an approximate date for the duplication event of 80–120 million years ago was calculated. The type II gene was shown, by RT-PCR, to be transcribed in adult individuals of Schmidtea polychroa, though at very low levels. This result, together with the fact that most of the functionally important positions for small-subunit rRNA in prokaryotes have been conserved, indicates that the type II gene is probably functional. Received: 24 March 1998 / Accepted: 17 March 1999     

The early vertebrate Danio rerio Mr 46000 mannose-6-phosphate receptor: biochemical and functional characterisation

Mannose-6-phosphate receptors (MPRs) have been identified in a wide range of species from humans to invertebrates such as molluscs. A characteristic of all MPRs is their common property to recognize mannose-6-phosphate residues that are labelling lysosomal enzymes and to mediate their targeting to lysosomes in mammalian cells by the corresponding receptor proteins. We present here the analysis of full-length sequences for MPR 46 from zebrafish (Danio rerio) and its functional analysis. This is the first non-mammalian MPR 46 to be characterised. The amino acid sequences of the zebrafish MPR 46 displays 70% similarity to the human MPR 46 protein. In particular, all essential cysteine residues, the transmembrane domain as well as the cytoplasmic tail residues harbouring the signals for endocytosis and Golgi-localizing, γ-ear-containing, ARF-binding protein (GGA)-mediated sorting at the trans-Golgi network, are highly conserved. The zebrafish MPR 46 has the arginine residue known to be essential for mannose-6-phosphate binding and other additional characteristic residues of the mannose-6-phosphate ligand-binding pocket. Like the mammalian MPR 46, zebrafish MPR 46 binds to the multimeric mannose-6-phosphate ligand phosphomannan and can rescue the missorting of lysosomal enzymes in mammalian MPR-deficient cells. The conserved C-terminal acidic dileucine motif (DxxLL) in the cytoplasmic domain of zebrafish MPR 46 essential for the interaction of the GGAs with the receptor domains interacts with the human GGA1-VHS domain. Interestingly, the serine residue suggested to regulate the interaction between the tail and the GGAs in a phosphorylation-dependent manner is substituted by a proline residue in fish. Electronic Supplementary Material Supplementary material is available for this article at . The zebrafish MPR 46 sequence data have been submitted to the GenBank database under accession no. DQ089037.     

Molecular evolution and functional divergence of the cytochrome P450 3 (CYP3) Family in Actinopterygii (ray-finned fish)

Background

The cytochrome P450 (CYP) superfamily is a multifunctional hemethiolate enzyme that is widely distributed from Bacteria to Eukarya. The CYP3 family contains mainly the four subfamilies CYP3A, CYP3B, CYP3C and CYP3D in vertebrates; however, only the Actinopterygii (ray-finned fish) have all four subfamilies and detailed understanding of the evolutionary relationship of Actinopterygii CYP3 family members would be valuable.

Methods and Findings

Phylogenetic relationships were constructed to trace the evolutionary history of the Actinopterygii CYP3 family genes. Selection analysis, relative rate tests and functional divergence analysis were combined to interpret the relationship of the site-specific evolution and functional divergence in the Actinopterygii CYP3 family. The results showed that the four CYP3 subfamilies in Actinopterygii might be formed by gene duplication. The first gene duplication event was responsible for divergence of the CYP3B/C clusters from ancient CYP3 before the origin of the Actinopterygii, which corresponded to the fish-specific whole genome duplication (WGD). Tandem repeat duplication in each of the homologue clusters produced stable CYP3B, CYP3C, CYP3A and CYP3D subfamilies. Acceleration of asymmetric evolutionary rates and purifying selection together were the main force for the production of new subfamilies and functional divergence in the new subset after gene duplication, whereas positive selection was detected only in the retained CYP3A subfamily. Furthermore, nearly half of the functional divergence sites appear to be related to substrate recognition, which suggests that site-specific evolution is closely related with functional divergence in the Actinopterygii CYP3 family.

Conclusions

The split of fish-specific CYP3 subfamilies was related to the fish-specific WGD, and site-specific acceleration of asymmetric evolutionary rates and purifying selection was the main force for the origin of the new subfamilies and functional divergence in the new subset after gene duplication. Site-specific evolution in substrate recognition was related to functional divergence in the Actinopterygii CYP3 family.     

Characterization of a thermostable enzyme with phosphomannomutase/phosphoglucomutase activities from the hyperthermophilic archaeon Pyrococcus horikoshii OT3

The phosphomannomutase/phosphoglucomutase (PMM/PGM) enzyme catalyzes reversibly the intra-molecular phosphoryl interconverting reaction of mannose-6-phosphate and mannose-1-phosphate or glucose-6-phosphate and glucose-1-phosphate. Glucose-6-phosphate and glucose-1-phosphate are known to be utilized for energy metabolism and cell surface construction, respectively. PMM/PGM has been isolated from many microorganisms. By performing similarity searches using existing PMM/PGM sequences, the homologous ORFs PH0923 and PH1210 were identified from the genomic data of Pyrococcus horikoshii OT3. Since PH0923 appears to be part of an operon consisting of four carbohydrate metabolic enzymes, PH0923 was selected as the first target for the investigation of PMM/PGM activity in P. horikoshii OT3. The coding region of PH0923 was cloned and the purified recombinant protein was utilized for an examination of its biochemical properties. The enzyme retained half its initial activity after treatment at 95 degrees C for 90 min. Detailed analyses of activities showed that this protein is capable of utilizing a variety of metal ions that are not utilized by previously characterized PMM/PGM proteins. A mutated protein with an alanine residue replacing the active site serine residue indicated that this residue plays an important but non-essential role in PMM/PGM activity.     

Structure of Leishmania mexicana phosphomannomutase highlights similarities with human isoforms

Phosphomannomutase (PMM) catalyses the conversion of mannose-6-phosphate to mannose-1-phosphate, an essential step in mannose activation and the biosynthesis of glycoconjugates in all eukaryotes. Deletion of PMM from Leishmania mexicana results in loss of virulence, suggesting that PMM is a promising drug target for the development of anti-leishmanial inhibitors. We report the crystallization and structure determination to 2.1 A of L. mexicana PMM alone and in complex with glucose-1,6-bisphosphate to 2.9 A. PMM is a member of the haloacid dehalogenase (HAD) family, but has a novel dimeric structure and a distinct cap domain of unique topology. Although the structure is novel within the HAD family, the leishmanial enzyme shows a high degree of similarity with its human isoforms. We have generated L. major PMM knockouts, which are avirulent. We expressed the human pmm2 gene in the Leishmania PMM knockout, but despite the similarity between Leishmania and human PMM, expression of the human gene did not restore virulence. Similarities in the structure of the parasite enzyme and its human isoforms suggest that the development of parasite-selective inhibitors will not be an easy task.     

HYPERMUTABILITY OF HOXA13A AND FUNCTIONAL DIVERGENCE FROM ITS PARALOG ARE ASSOCIATED WITH THE ORIGIN OF A NOVEL DEVELOPMENTAL FEATURE IN ZEBRAFISH AND RELATED TAXA (CYPRINIFORMES)

Gene duplication is widely regarded as the predominant mechanism by which genes with new functions and associated phenotypic novelties arise. A whole genome duplication occurred shortly before the most recent common ancestor of teleosts, the most diverse chordate group, resulting in duplication and retention of many Hox cluster genes. Because they play a key role in determination of body plan morphology, it has been widely assumed that Hox genes play a key role in the evolution of diverse metazoan body plans. However, it is not clear whether certain aspects of molecular evolution, such as asymmetric divergence and neofunctionalization, contribute to the initial retention of paralogs. We investigate the molecular evolution and functional divergence of the duplicated HoxA13 paralogs in zebrafish to determine when asymmetric divergence and functional divergence occurred after the duplication event. Our findings demonstrate the contribution of gene duplication to the evolution of novel features through evolutionary mechanisms other than those traditionally investigated, such as positive selection occurring immediately after gene duplication. Rather, we find a latent build up of molecular changes in a gene associated with the development of a novel feature in a very diverse group of fishes.     

Cloning and characterization of phosphoglucomutase and phosphomannomutase derived from Sphingomonas chungbukensis DJ77

The enzymes phosphoglucomutase (PGM) and phosphomannomutase (PMM) play an important role in the synthesis of extracellular polysaccharide. By colony hybridization of the fosmid library of Sphingomonas chungbukensis DJ77, an open reading frame (ORF-1) of 1,626 nucleotides, whose predicted product is highly homologous with other PGM proteins from several bacterial species, was identified. An additional open reading frame (ORF-2) of 1,437 nucleotides was identified, and its encoded protein shows a high level of similarity with the PGM/PMM protein family. The two genes were cloned into a bacterial expression vector pET-15b (+) and expressed in Escherichia coli as fusion proteins with (His)(6)-tag. Both recombinant proteins (designated as SP-1 and SP-2 for ORF-1 and ORF-2, respectively) exhibited PGM and PMM activities. The molecular masses of subunits of SP-1 and SP-2 were estimated to be around 58 and 51 kDa from SDS-PAGE, respectively. However, molecular masses of SP-1 and SP-2 in their native condition were determined to be approximately 59.5 and 105.4 kDa, according to non-denaturing PAGE, respectively. The SP-1 protein has a preference for glucose-1-phosphate rather than mannose-1-phosphate, while the preferred substrate of SP-2 is mannose-1-phosphate. Thus, the existence of two proteins with bifunctional PGM/PMM activities was first found S. chungbukensis DJ77.     

Molecular evolution of theSaccharomyces cerevisiae histone gene loci

Summary The core histone genes ofSaccharomyces cerevisiae are arranged as duplicate nonallelic sets of specifically paired genes. The identity of structural organization between the duplicated gene pairs would have its simplest evolutionary origin in the duplication of a complete locus in a single event. In such a case, the time since the duplication of one of the genes should be identical to that since duplication of the gene adjacent to it on the chromosome. A calculation of the evolutionary distances between the coding DNA sequences of the histone genes leads to a duplication paradox: The extents of sequence divergence in the silent component of third-base positions for adjacent pairs of genes are not identical. Estimates of the evolutionary distance between the two H3-H4 noncoding intergene DNA sequences are large; the divergence between the two separate sequences is indistinguishable from the divergence between either of the regions and a randomly generated permutation of itself. These results suggest that the duplication event may have occurred much earlier than previously estimated. The potential age of the duplication, and the attractive simplicity of the duplication of both the H3-H4 and the H2A-H2B gene pairs having taken place in a single event, leads to the hypothesis that modern haploidS. cerevisiae may have evolved by diploidization or fusion of two ancient fungi.     

Molecular evolution of the 14-3-3 protein family

Members of the highly conserved and ubiquitous 14-3-3 protein family modulate a wide variety of cellular processes. To determine the evolutionary relationships among specific 14-3-3 proteins in different plant, animal, and fungal species and to initiate a predictive analysis of isoform-specific differences in light of the latest functional and structural studies of 14-3-3, multiple alignments were constructed from forty-six 14-3-3 sequences retrieved from the GenBank and SwissProt databases and a newly identified second 14-3-3 gene fromCaenorhabditis elegans. The alignment revealed five highly conserved sequence blocks. Blocks 2–5 correlate well with the alpha helices 3, 5, 7, and 9 which form the proposed internal binding domain in the three-dimensional structure model of the functioning dimer. Amino acid differences within the functional and structural domains of plant and animal 14-3-3 proteins were identified which may account for functional diversity amongst isoforms. Protein phylogenic trees were constructed using both the maximum parsimony and neighbor joining methods of the PHYLIP(3.5c) package; 14-3-3 proteins fromEntamoeba histolytica, an amitochondrial protozoa, were employed as an outgroup in our analysis. Epsilon isoforms from the animal lineage form a distinct grouping in both trees, which suggests an early divergence from the other animal isoforms. Epsilons were found to be more similar to yeast and plant isoforms than other animal isoforms at numerous amino acid positions, and thus epsilon may have retained functional characteristics of the ancestral protein. The known invertebrate proteins group with the nonepsilon mammalian isoforms. Most of the current 14-3-3 isoform diversity probably arose through independent duplication events after the divergence of the major eukaryotic kingdoms. Divergence of the seven mammalian isoforms beta, zeta, gamma, eta, epsilon, tau, and sigma (stratifin/ HME1) occurred before the divergence of mammalian and perhaps before the divergence of vertebrate species. A possible ancestral 14-3-3 sequence is proposed. Correspondence to: D.C. Shakes     

Rapid Diversification of FoxP2 in Teleosts through Gene Duplication in the Teleost-Specific Whole Genome Duplication Event

As one of the most conserved genes in vertebrates, FoxP2 is widely involved in a number of important physiological and developmental processes. We systematically studied the evolutionary history and functional adaptations of FoxP2 in teleosts. The duplicated FoxP2 genes (FoxP2a and FoxP2b), which were identified in teleosts using synteny and paralogon analysis on genome databases of eight organisms, were probably generated in the teleost-specific whole genome duplication event. A credible classification with FoxP2, FoxP2a and FoxP2b in phylogenetic reconstructions confirmed the teleost-specific FoxP2 duplication. The unavailability of FoxP2b in Danio rerio suggests that the gene was deleted through nonfunctionalization of the redundant copy after the Otocephala-Euteleostei split. Heterogeneity in evolutionary rates among clusters consisting of FoxP2 in Sarcopterygii (Cluster 1), FoxP2a in Teleostei (Cluster 2) and FoxP2b in Teleostei (Cluster 3), particularly between Clusters 2 and 3, reveals asymmetric functional divergence after the gene duplication. Hierarchical cluster analyses of hydrophobicity profiles demonstrated significant structural divergence among the three clusters with verification of subsequent stepwise discriminant analysis, in which FoxP2 of Leucoraja erinacea and Lepisosteus oculatus were classified into Cluster 1, whereas FoxP2b of Salmo salar was grouped into Cluster 2 rather than Cluster 3. The simulated thermodynamic stability variations of the forkhead box domain (monomer and homodimer) showed remarkable divergence in FoxP2, FoxP2a and FoxP2b clusters. Relaxed purifying selection and positive Darwinian selection probably were complementary driving forces for the accelerated evolution of FoxP2 in ray-finned fishes, especially for the adaptive evolution of FoxP2a and FoxP2b in teleosts subsequent to the teleost-specific gene duplication.     

Two independent duplications forming the Cyp307a genes in Drosophila

The conserved relationship between orthologs of many cytochrome P450 genes involved in ecdysone synthesis is not reflected in the evolution of the Drosophila Cyp307a genes. In Drosophila melanogaster Cyp307a1 (spook) and Cyp307a2 (spookier) both play essential roles in ecdysone synthesis and may possess biochemically redundant functions. Using phylogenetic analyses we show that the Drosophila Cyp307a genes were formed from two independent duplication events depicting a complicated evolutionary scenario. An initial duplication, from a Cyp307a2 ancestral gene produced the Cyp307a1 gene that has been maintained only in the Sophophoran subgenus. A second duplication in the Drosophila subgenus formed an additional paralog, Cyp307a3. Microsynteny is conserved for Cyp307a2 throughout the Drosophila species, but is not conserved between Cyp307a1 and Cyp307a3. These are located in different genomic positions in the Sophophora and Drosophila subgenera, respectively. Cyp307a3 appears to encode a functional gene product and is expressed in a different spatial and temporal manner to Cyp307a1. This suggests some level of functional divergence between the Cyp307a paralogs in different Drosophila species.     

Functional Evolution of Phosphatidylethanolamine Binding Proteins in Soybean and Arabidopsis

Gene duplication provides resources for novel gene functions. Identification of the amino acids responsible for functional conservation and divergence of duplicated genes will strengthen our understanding of their evolutionary course. Here, we conducted a systemic functional investigation of phosphatidylethanolamine binding proteins (PEBPs) in soybean (Glycine max) and Arabidopsis thaliana. Our results demonstrated that after the ancestral duplication, the lineage of the common ancestor of the FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) subfamilies functionally diverged from the MOTHER OF FT AND TFL1 (MFT) subfamily to activate flowering and repress flowering, respectively. They also underwent further specialization after subsequent duplications. Although the functional divergence increased with duplication age, we observed rapid functional divergence for a few pairs of young duplicates in soybean. Association analysis between amino acids and functional variations identified critical amino acid residues that led to functional differences in PEBP members. Using transgenic analysis, we validated a subset of these differences. We report clear experimental evidence for the functional evolution of the PEBPs in the MFT, FT, and TFL1 subfamilies, which predate the origin of angiosperms. Our results highlight the role of amino acid divergence in driving evolutionary novelty after duplication.