The Role of Ancestral Duplicated Genes in Adaptation to Growth on Lactate, a Non-Fermentable Carbon Source for the Yeast Saccharomyces cerevisiae.

The cell central metabolism has been shaped throughout evolutionary times when facing challenges from the availability of resources. In the budding yeast, Saccharomyces cerevisiae, a set of duplicated genes originating from an ancestral whole-genome and several coetaneous small-scale duplication events drive energy transfer through glucose metabolism as the main carbon source either by fermentation or respiration. These duplicates (~a third of the genome) have been dated back to approximately 100 MY, allowing for enough evolutionary time to diverge in both sequence and function. Gene duplication has been proposed as a molecular mechanism of biological innovation, maintaining balance between mutational robustness and evolvability of the system. However, some questions concerning the molecular mechanisms behind duplicated genes transcriptional plasticity and functional divergence remain unresolved. In this work we challenged S. cerevisiae to the use of lactic acid/lactate as the sole carbon source and performed a small adaptive laboratory evolution to this non-fermentative carbon source, determining phenotypic and transcriptomic changes. We observed growth adaptation to acidic stress, by reduction of growth rate and increase in biomass production, while the transcriptomic response was mainly driven by repression of the whole-genome duplicates, those implied in glycolysis and overexpression of ROS response. The contribution of several duplicated pairs to this carbon source switch and acidic stress is also discussed.

Transcriptional Rewiring, Adaptation, and the Role of Gene Duplication in the Metabolism of Ethanol of Saccharomyces cerevisiae.

Ethanol is the main by-product of yeast sugar fermentation that affects microbial growth parameters, being considered a dual molecule, a nutrient and a stressor. Previous works demonstrated that the budding yeast arose after an ancient hybridization process resulted in a tier of duplicated genes within its genome, many of them with implications in this ethanol "produce-accumulate-consume" strategy. The evolutionary link between ethanol production, consumption, and tolerance versus ploidy and stability of the hybrids is an ongoing debatable issue. The implication of ancestral duplicates in this metabolic rewiring, and how these duplicates differ transcriptionally, remains unsolved. Here, we study the transcriptomic adaptive signatures to ethanol as a nonfermentative carbon source to sustain clonal yeast growth by experimental evolution, emphasizing the role of duplicated genes in the adaptive process. As expected, ethanol was able to sustain growth but at a lower rate than glucose. Our results demonstrate that in asexual populations a complete transcriptomic rewiring was produced, strikingly by downregulation of duplicated genes, mainly whole-genome duplicates, whereas small-scale duplicates exhibited significant transcriptional divergence between copies. Overall, this study contributes to the understanding of evolution after gene duplication, linking transcriptional divergence with duplicates' fate in a multigene trait as ethanol tolerance.IMPORTANCE Gene duplication events have been related with increasing biological complexity through the tree of life, but also with illnesses, including cancer. Early evolutionary theories indicated that duplicated genes could explore alternative functions due to relaxation of selective constraints in one of the copies, as the other remains as ancestral-function backup. In unicellular eukaryotes like yeasts, it has been demonstrated that the fate and persistence of duplicates depend on duplication mechanism (whole-genome or small-scale events), shaping their actual genomes. Although it has been shown that small-scale duplicates tend to innovate and whole-genome duplicates specialize in ancestral functions, the implication of duplicates' transcriptional plasticity and transcriptional divergence on environmental and metabolic responses remains largely obscure. Here, by experimental adaptive evolution, we show that Saccharomyces cerevisiae is able to respond to metabolic stress (ethanol as nonfermentative carbon source) due to the persistence of duplicated genes. These duplicates respond by transcriptional rewiring, depending on their transcriptional background. Our results shed light on the mechanisms that determine the role of duplicates, and on their evolvability.

Exploring molecular evolution of Rubisco in C3 and CAM Orchidaceae and Bromeliaceae.

BACKGROUND: The CO2-concentrating mechanism associated to Crassulacean acid metabolism (CAM) alters the catalytic context for Rubisco by increasing CO2 availability and provides an advantage in particular ecological conditions. We hypothesized about the existence of molecular changes linked to these particular adaptations in CAM Rubisco. We investigated molecular evolution of the Rubisco large (L-) subunit in 78 orchids and 144 bromeliads with C3 and CAM photosynthetic pathways. The sequence analyses were complemented with measurements of Rubisco kinetics in some species with contrasting photosynthetic mechanism and differing in the L-subunit sequence. RESULTS: We identified potential positively selected sites and residues with signatures of co-adaptation. The implementation of a decision tree model related Rubisco specific variable sites to the leaf carbon isotopic composition of the species. Differences in the Rubisco catalytic traits found among C3 orchids and between strong CAM and C3 bromeliads suggested Rubisco had evolved in response to differing CO2 concentration. CONCLUSIONS: The results revealed that the variability in the Rubisco L-subunit sequence in orchids and bromeliads is composed of coevolving sites under potential positive adaptive signal. The sequence variability was related to δ13C in orchids and bromeliads, however it could not be linked to the variability found in the kinetic properties of the studied species.

Molecular Chaperones Accelerate the Evolution of Their Protein Clients in Yeast.

Protein stability is a major constraint on protein evolution. Molecular chaperones, also known as heat-shock proteins, can relax this constraint and promote protein evolution by diminishing the deleterious effect of mutations on protein stability and folding. This effect, however, has only been stablished for a few chaperones. Here, we use a comprehensive chaperone-protein interaction network to study the effect of all yeast chaperones on the evolution of their protein substrates, that is, their clients. In particular, we analyze how yeast chaperones affect the evolutionary rates of their clients at two very different evolutionary time scales. We first study the effect of chaperone-mediated folding on protein evolution over the evolutionary divergence of Saccharomyces cerevisiae and S. paradoxus. We then test whether yeast chaperones have left a similar signature on the patterns of standing genetic variation found in modern wild and domesticated strains of S. cerevisiae. We find that genes encoding chaperone clients have diverged faster than genes encoding non-client proteins when controlling for their number of protein-protein interactions. We also find that genes encoding client proteins have accumulated more intraspecific genetic diversity than those encoding non-client proteins. In a number of multivariate analyses, controlling by other well-known factors that affect protein evolution, we find that chaperone dependence explains the largest fraction of the observed variance in the rate of evolution at both evolutionary time scales. Chaperones affecting rates of protein evolution mostly belong to two major chaperone families: Hsp70s and Hsp90s. Our analyses show that protein chaperones, by virtue of their ability to buffer destabilizing mutations and their role in modulating protein genotype-phenotype maps, have a considerable accelerating effect on protein evolution.

Chaperonin overproduction and metabolic erosion caused by mutation accumulation in Escherichia coli.

Bacterial cells adapting to a constant environment tend to accumulate mutations in portions of their genome that are not maintained by selection. This process has been observed in bacteria evolving under strong genetic drift, and especially in bacterial endosymbionts of insects. Here, we study this process in hypermutable Escherichia coli populations evolved through 250 single-cell bottlenecks on solid rich medium in a mutation accumulation experiment that emulates the evolution of bacterial endosymbionts. Using phenotype microarrays monitoring metabolic activity in 95 environments distinguished by their carbon sources, we observe how mutation accumulation has decreased the ability of cells to metabolize most carbon sources. We study if the chaperonin GroEL, which is naturally overproduced in bacterial endosymbionts, can ameliorate the process of metabolic erosion, because of its known ability to buffer destabilizing mutations in metabolic enzymes. Our results indicate that GroEL can slow down the negative phenotypic consequences of genome decay in some environments.

Inferring the regulatory network of the miRNA-mediated response to biotic and abiotic stress in melon.

BACKGROUND: MiRNAs have emerged as key regulators of stress response in plants, suggesting their potential as candidates for knock-in/out to improve stress tolerance in agricultural crops. Although diverse assays have been performed, systematic and detailed studies of miRNA expression and function during exposure to multiple environments in crops are limited. RESULTS: Here, we present such pioneering analysis in melon plants in response to seven biotic and abiotic stress conditions. Deep-sequencing and computational approaches have identified twenty-four known miRNAs whose expression was significantly altered under at least one stress condition, observing that down-regulation was preponderant. Additionally, miRNA function was characterized by high scale degradome assays and quantitative RNA measurements over the intended target mRNAs, providing mechanistic insight. Clustering analysis provided evidence that eight miRNAs showed a broad response range under the stress conditions analyzed, whereas another eight miRNAs displayed a narrow response range. Transcription factors were predominantly targeted by stress-responsive miRNAs in melon. Furthermore, our results show that the miRNAs that are down-regulated upon stress predominantly have as targets genes that are known to participate in the stress response by the plant, whereas the miRNAs that are up-regulated control genes linked to development. CONCLUSION: Altogether, this high-resolution analysis of miRNA-target interactions, combining experimental and computational work, Illustrates the close interplay between miRNAs and the response to diverse environmental conditions, in melon.

Arabidopsis Heat Stress-Induced Proteins Are Enriched in Electrostatically Charged Amino Acids and Intrinsically Disordered Regions.

Comparison of the proteins of thermophilic, mesophilic, and psychrophilic prokaryotes has revealed several features characteristic to proteins adapted to high temperatures, which increase their thermostability. These characteristics include a profusion of disulfide bonds, salt bridges, hydrogen bonds, and hydrophobic interactions, and a depletion in intrinsically disordered regions. It is unclear, however, whether such differences can also be observed in eukaryotic proteins or when comparing proteins that are adapted to temperatures that are more subtly different. When an organism is exposed to high temperatures, a subset of its proteins is overexpressed (heat-induced proteins), whereas others are either repressed (heat-repressed proteins) or remain unaffected. Here, we determine the expression levels of all genes in the eukaryotic model system Arabidopsis thaliana at 22 and 37 °C, and compare both the amino acid compositions and levels of intrinsic disorder of heat-induced and heat-repressed proteins. We show that, compared to heat-repressed proteins, heat-induced proteins are enriched in electrostatically charged amino acids and depleted in polar amino acids, mirroring thermophile proteins. However, in contrast with thermophile proteins, heat-induced proteins are enriched in intrinsically disordered regions, and depleted in hydrophobic amino acids. Our results indicate that temperature adaptation at the level of amino acid composition and intrinsic disorder can be observed not only in proteins of thermophilic organisms, but also in eukaryotic heat-induced proteins; the underlying adaptation pathways, however, are similar but not the same.

NaCl-saturated brines are thermodynamically moderate, rather than extreme, microbial habitats.

NaCl-saturated brines such as saltern crystalliser ponds, inland salt lakes, deep-sea brines and liquids-of-deliquescence on halite are commonly regarded as a paradigm for the limit of life on Earth. There are, however, other habitats that are thermodynamically more extreme. Typically, NaCl-saturated environments contain all domains of life and perform complete biogeochemical cycling. Despite their reduced water activity, ∼0.755 at 5 M NaCl, some halophiles belonging to the Archaea and Bacteria exhibit optimum growth/metabolism in these brines. Furthermore, the recognised water-activity limit for microbial function, ∼0.585 for some strains of fungi, lies far below 0.755. Other biophysical constraints on the microbial biosphere (temperatures of >121°C; pH > 12; and high chaotropicity; e.g. ethanol at >18.9% w/v (24% v/v) and MgCl2 at >3.03 M) can prevent any cellular metabolism or ecosystem function. By contrast, NaCl-saturated environments contain biomass-dense, metabolically diverse, highly active and complex microbial ecosystems; and this underscores their moderate character. Here, we survey the evidence that NaCl-saturated brines are biologically permissive, fertile habitats that are thermodynamically mid-range rather than extreme. Indeed, were NaCl sufficiently soluble, some halophiles might grow at concentrations of up to 8 M. It may be that the finite solubility of NaCl has stabilised the genetic composition of halophile populations and limited the action of natural selection in driving halophile evolution towards greater xerophilicity. Further implications are considered for the origin(s) of life and other aspects of astrobiology.

Intrinsic adaptive value and early fate of gene duplication revealed by a bottom-up approach.

The population genetic mechanisms governing the preservation of gene duplicates, especially in the critical very initial phase, have remained largely unknown. Here, we demonstrate that gene duplication confers per se a weak selective advantage in scenarios of fitness trade-offs. Through a precise quantitative description of a model system, we show that a second gene copy serves to reduce gene expression inaccuracies derived from pervasive molecular noise and suboptimal gene regulation. We then reveal that such an accuracy in the phenotype yields a selective advantage in the order of 0.1% on average, which would allow the positive selection of gene duplication in populations with moderate/large sizes. This advantage is greater at higher noise levels and intermediate concentrations of the environmental molecule, when fitness trade-offs become more evident. Moreover, we discuss how the genome rearrangement rates greatly condition the eventual fixation of duplicates. Overall, our theoretical results highlight an original adaptive value for cells carrying new-born duplicates, broadly analyze the selective conditions that determine their early fates in different organisms, and reconcile population genetics with evolution by gene duplication.

Correction: Positively selected amino acid replacements within the RuBisCO enzyme of oak trees are associated with ecological adaptations.

[This corrects the article DOI: 10.1371/journal.pone.0183970.].

Positively selected amino acid replacements within the RuBisCO enzyme of oak trees are associated with ecological adaptations.

Phylogenetic analysis by maximum likelihood (PAML) has become the standard approach to study positive selection at the molecular level, but other methods may provide complementary ways to identify amino acid replacements associated with particular conditions. Here, we compare results of the decision tree (DT) model method with ones of PAML using the key photosynthetic enzyme RuBisCO as a model system to study molecular adaptation to particular ecological conditions in oaks (Quercus). We sequenced the chloroplast rbcL gene encoding RuBisCO large subunit in 158 Quercus species, covering about a third of the global genus diversity. It has been hypothesized that RuBisCO has evolved differentially depending on the environmental conditions and leaf traits governing internal gas diffusion patterns. Here, we show, using PAML, that amino acid replacements at the residue positions 95, 145, 251, 262 and 328 of the RuBisCO large subunit have been the subject of positive selection along particular Quercus lineages associated with the leaf traits and climate characteristics. In parallel, the DT model identified amino acid replacements at sites 95, 219, 262 and 328 being associated with the leaf traits and climate characteristics, exhibiting partial overlap with the results obtained using PAML.

Deconstruction of the beaten Path-Sidestep interaction network provides insights into neuromuscular system development.

An 'interactome' screen of all Drosophila cell-surface and secreted proteins containing immunoglobulin superfamily (IgSF) domains discovered a network formed by paralogs of Beaten Path (Beat) and Sidestep (Side), a ligand-receptor pair that is central to motor axon guidance. Here we describe a new method for interactome screening, the Bio-Plex Interactome Assay (BPIA), which allows identification of many interactions in a single sample. Using the BPIA, we 'deorphanized' four more members of the Beat-Side network. We confirmed interactions using surface plasmon resonance. The expression patterns of beat and side genes suggest that Beats are neuronal receptors for Sides expressed on peripheral tissues. side-VI is expressed in muscle fibers targeted by the ISNb nerve, as well as at growth cone choice points and synaptic targets for the ISN and TN nerves. beat-V genes, encoding Side-VI receptors, are expressed in ISNb and ISN motor neurons.

Expression properties exhibit correlated patterns with the fate of duplicated genes, their divergence, and transcriptional plasticity in Saccharomycotina.

Gene duplication is an important source of novelties and genome complexity. What genes are preserved as duplicated through long evolutionary times can shape the evolution of innovations. Identifying factors that influence gene duplicability is therefore an important aim in evolutionary biology. Here, we show that in the yeast Saccharomyces cerevisiae the levels of gene expression correlate with gene duplicability, its divergence, and transcriptional plasticity. Genes that were highly expressed before duplication are more likely to be preserved as duplicates for longer evolutionary times and wider phylogenetic ranges than genes that were lowly expressed. Duplicates with higher expression levels exhibit greater divergence between their gene copies. Duplicates that exhibit higher expression divergence are those enriched for TATA-containing promoters. These duplicates also show transcriptional plasticity, which seems to be involved in the origin of adaptations to environmental stresses in yeast. While the expression properties of genes strongly affect their duplicability, divergence and transcriptional plasticity are enhanced after gene duplication. We conclude that highly expressed genes are more likely to be preserved as duplicates due to their promoter architectures, their greater tolerance to expression noise, and their ability to reduce the noise-plasticity conflict.

Genome Mutational and Transcriptional Hotspots Are Traps for Duplicated Genes and Sources of Adaptations.

Gene duplication generates new genetic material, which has been shown to lead to major innovations in unicellular and multicellular organisms. A whole-genome duplication occurred in the ancestor of Saccharomyces yeast species but 92% of duplicates returned to single-copy genes shortly after duplication. The persisting duplicated genes in Saccharomyces led to the origin of major metabolic innovations, which have been the source of the unique biotechnological capabilities in the Baker's yeast Saccharomyces cerevisiae. What factors have determined the fate of duplicated genes remains unknown. Here, we report the first demonstration that the local genome mutation and transcription rates determine the fate of duplicates. We show, for the first time, a preferential location of duplicated genes in the mutational and transcriptional hotspots of S. cerevisiae genome. The mechanism of duplication matters, with whole-genome duplicates exhibiting different preservation trends compared to small-scale duplicates. Genome mutational and transcriptional hotspots are rich in duplicates with large repetitive promoter elements. Saccharomyces cerevisiae shows more tolerance to deleterious mutations in duplicates with repetitive promoter elements, which in turn exhibit higher transcriptional plasticity against environmental perturbations. Our data demonstrate that the genome traps duplicates through the accelerated regulatory and functional divergence of their gene copies providing a source of novel adaptations in yeast.

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects.

An open question in evolutionary biology is how does the selection-drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host-symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shift.

Glycerol stress in Saccharomyces cerevisiae: Cellular responses and evolved adaptations.

Glycerol synthesis is key to central metabolism and stress biology in Saccharomyces cerevisiae, yet the cellular adjustments needed to respond and adapt to glycerol stress are little understood. Here, we determined impacts of acute and chronic exposures to glycerol stress in S. cerevisiae. Glycerol stress can result from an increase of glycerol concentration in the medium due to the S. cerevisiae fermenting activity or other metabolic activities. Acute glycerol-stress led to a 50% decline in growth rate and altered transcription of more than 40% of genes. The increased genetic diversity in S. cerevisiae population, which had evolved in the standard nutrient medium for hundreds of generations, led to an increase in growth rate and altered transcriptome when such population was transferred to stressful media containing a high concentration of glycerol; 0.41 M (0.990 water activity). Evolution of S. cerevisiae populations during a 10-day period in the glycerol-containing medium led to transcriptome changes and readjustments to improve control of glycerol flux across the membrane, regulation of cell cycle, and more robust stress response; and a remarkable increase of growth rate under glycerol stress. Most of the observed regulatory changes arose in duplicated genes. These findings elucidate the physiological mechanisms, which underlie glycerol-stress response, and longer-term adaptations, in S. cerevisiae; they also have implications for enigmatic aspects of the ecology of this otherwise well-characterized yeast.

The Phenotypic Plasticity of Duplicated Genes in Saccharomyces cerevisiae and the Origin of Adaptations.

Gene and genome duplication are the major sources of biological innovations in plants and animals. Functional and transcriptional divergence between the copies after gene duplication has been considered the main driver of innovations . However, here we show that increased phenotypic plasticity after duplication plays a more major role than thought before in the origin of adaptations. We perform an exhaustive analysis of the transcriptional alterations of duplicated genes in the unicellular eukaryote Saccharomyces cerevisiae when challenged with five different environmental stresses. Analysis of the transcriptomes of yeast shows that gene duplication increases the transcriptional response to environmental changes, with duplicated genes exhibiting signatures of adaptive transcriptional patterns in response to stress. The mechanism of duplication matters, with whole-genome duplicates being more transcriptionally altered than small-scale duplicates. The predominant transcriptional pattern follows the classic theory of evolution by gene duplication; with one gene copy remaining unaltered under stress, while its sister copy presents large transcriptional plasticity and a prominent role in adaptation. Moreover, we find additional transcriptional profiles that are suggestive of neo- and subfunctionalization of duplicate gene copies. These patterns are strongly correlated with the functional dependencies and sequence divergence profiles of gene copies. We show that, unlike singletons, duplicates respond more specifically to stress, supporting the role of natural selection in the transcriptional plasticity of duplicates. Our results reveal the underlying transcriptional complexity of duplicated genes and its role in the origin of adaptations.

Essentiality Is a Strong Determinant of Protein Rates of Evolution during Mutation Accumulation Experiments in Escherichia coli.

The Neutral Theory of Molecular Evolution is considered the most powerful theory to understand the evolutionary behavior of proteins. One of the main predictions of this theory is that essential proteins should evolve slower than dispensable ones owing to increased selective constraints. Comparison of genomes of different species, however, has revealed only small differences between the rates of evolution of essential and nonessential proteins. In some analyses, these differences vanish once confounding factors are controlled for, whereas in other cases essentiality seems to have an independent, albeit small, effect. It has been argued that comparing relatively distant genomes may entail a number of limitations. For instance, many of the genes that are dispensable in controlled lab conditions may be essential in some of the conditions faced in nature. Moreover, essentiality can change during evolution, and rates of protein evolution are simultaneously shaped by a variety of factors, whose individual effects are difficult to isolate. Here, we conducted two parallel mutation accumulation experiments in Escherichia coli, during 5,500-5,750 generations, and compared the genomes at different points of the experiments. Our approach (a short-term experiment, under highly controlled conditions) enabled us to overcome many of the limitations of previous studies. We observed that essential proteins evolved substantially slower than nonessential ones during our experiments. Strikingly, rates of protein evolution were only moderately affected by expression level and protein length.

The Molecular Chaperone DnaK Is a Source of Mutational Robustness.

Molecular chaperones, also known as heat-shock proteins, refold misfolded proteins and help other proteins reach their native conformation. Thanks to these abilities, some chaperones, such as the Hsp90 protein or the chaperonin GroEL, can buffer the deleterious phenotypic effects of mutations that alter protein structure and function. Hsp70 chaperones use a chaperoning mechanism different from that of Hsp90 and GroEL, and it is not known whether they can also buffer mutations. Here, we show that they can. To this end, we performed a mutation accumulation experiment in Escherichia coli, followed by whole-genome resequencing. Overexpression of the Hsp70 chaperone DnaK helps cells cope with mutational load and completely avoid the extinctions we observe in lineages evolving without chaperone overproduction. Additionally, our sequence data show that DnaK overexpression increases mutational robustness, the tolerance of its clients to nonsynonymous nucleotide substitutions. We also show that this elevated mutational buffering translates into differences in evolutionary rates on intermediate and long evolutionary time scales. Specifically, we studied the evolutionary rates of DnaK clients using the genomes of E. coli, Salmonella enterica, and 83 other gamma-proteobacteria. We find that clients that interact strongly with DnaK evolve faster than weakly interacting clients. Our results imply that all three major chaperone classes can buffer mutations and affect protein evolution. They illustrate how an individual protein like a chaperone can have a disproportionate effect on the evolution of a proteome.

Evolution by gene duplication of Medicago truncatula PISTILLATA-like transcription factors.

PISTILLATA (PI) is a member of the B-function MADS-box gene family, which controls the identity of both petals and stamens in Arabidopsis thaliana. In Medicago truncatula (Mt), there are two PI-like paralogs, known as MtPI and MtNGL9. These genes differ in their expression patterns, but it is not known whether their functions have also diverged. Describing the evolution of certain duplicated genes, such as transcription factors, remains a challenge owing to the complex expression patterns and functional divergence between the gene copies. Here, we report a number of functional studies, including analyses of gene expression, protein-protein interactions, and reverse genetic approaches designed to demonstrate the respective contributions of each M. truncatula PI-like paralog to the B-function in this species. Also, we have integrated molecular evolution approaches to determine the mode of evolution of Mt PI-like genes after duplication. Our results demonstrate that MtPI functions as a master regulator of B-function in M. truncatula, maintaining the overall ancestral function, while MtNGL9 does not seem to have a role in this regard, suggesting that the pseudogenization could be the functional evolutionary fate for this gene. However, we provide evidence that purifying selection is the primary evolutionary force acting on this paralog, pinpointing the conservation of its biochemical function and, alternatively, the acquisition of a new role for this gene.