Chance and pleiotropy dominate genetic diversity in complex bacterial environments

How does environmental complexity affect the evolution of single genes? Here, we measured the effects of a set of Bacillus subtilis glutamate dehydrogenase mutants across 19 different environments-from phenotypically homogeneous single-cell populations in liquid media to heterogeneous biofilms, plant roots and soil populations. The effects of individual gene mutations on organismal fitness were highly reproducible in liquid cultures. However, 84% of the tested alleles showed opposing fitness effects under different growth conditions (sign environmental pleiotropy). In colony biofilms and soil samples, different alleles dominated in parallel replica experiments. Accordingly, we found that in these heterogeneous cell populations the fate of mutations was dictated by a combination of selection and drift. The latter relates to programmed prophage excisions that occurred during biofilm development. Overall, for each condition, a wide range of glutamate dehydrogenase mutations persisted and sometimes fixated as a result of the combined action of selection, pleiotropy and chance. However, over longer periods and in multiple environments, nearly all of this diversity would be lost-across all the environments and conditions that we tested, the wild type was the fittest allele.

Evolution of substrate specificity in a retained enzyme driven by gene loss

The connection between gene loss and the functional adaptation of retained proteins is still poorly understood. We apply phylogenomics and metabolic modeling to detect bacterial species that are evolving by gene loss, with the finding that Actinomycetaceae genomes from human cavities are undergoing sizable reductions, including loss of L-histidine and L-tryptophan biosynthesis. We observe that the dual-substrate phosphoribosyl isomerase A or priA gene, at which these pathways converge, appears to coevolve with the occurrence of trp and his genes. Characterization of a dozen PriA homologs shows that these enzymes adapt from bifunctionality in the largest genomes, to a monofunctional, yet not necessarily specialized, inefficient form in genomes undergoing reduction. These functional changes are accomplished via mutations, which result from relaxation of purifying selection, in residues structurally mapped after sequence and X-ray structural analyses. Our results show how gene loss can drive the evolution of substrate specificity from retained enzymes.

DOI:http://dx.doi.org/10.7554/eLife.22679.001

Insights into the evolution of enzyme substrate promiscuity after the discovery of (βα)₈ isomerase evolutionary intermediates from a diverse metagenome

BACKGROUND

Current sequence-based approaches to identify enzyme functional shifts, such as enzyme promiscuity, have proven to be highly dependent on a priori functional knowledge, hampering our ability to reconstruct evolutionary history behind these mechanisms. Hidden Markov Model (HMM) profiles, broadly used to classify enzyme families, can be useful to distinguish between closely related enzyme families with different specificities. The (βα)8-isomerase HisA/PriA enzyme family, involved in L-histidine (HisA, mono-substrate) biosynthesis in most bacteria and plants, but also in L-tryptophan (HisA/TrpF or PriA, dual-substrate) biosynthesis in most Actinobacteria, has been used as model system to explore evolutionary hypotheses and therefore has a considerable amount of evolutionary, functional and structural knowledge available. We searched for functional evolutionary intermediates between the HisA and PriA enzyme families in order to understand the functional divergence between these families.

RESULTS

We constructed a HMM profile that correctly classifies sequences of unknown function into the HisA and PriA enzyme sub-families. Using this HMM profile, we mined a large metagenome to identify plausible evolutionary intermediate sequences between HisA and PriA. These sequences were used to perform phylogenetic reconstructions and to identify functionally conserved amino acids. Biochemical characterization of one selected enzyme (CAM1) with a mutation within the functionally essential N-terminus phosphate-binding site, namely, an alanine instead of a glycine in HisA or a serine in PriA, showed that this evolutionary intermediate has dual-substrate specificity. Moreover, site-directed mutagenesis of this alanine residue, either backwards into a glycine or forward into a serine, revealed the robustness of this enzyme. None of these mutations, presumably upon functionally essential amino acids, significantly abolished its enzyme activities. A truncated version of this enzyme (CAM2) predicted to adopt a (βα)6-fold, and thus entirely lacking a C-terminus phosphate-binding site, was identified and shown to have HisA activity.

CONCLUSION

As expected, reconstruction of the evolution of PriA from HisA with HMM profiles suggest that functional shifts involve mutations in evolutionarily intermediate enzymes of otherwise functionally essential residues or motifs. These results are in agreement with a link between promiscuous enzymes and intragenic epistasis. HMM provides a convenient approach for gaining insights into these evolutionary processes.

Enzyme evolution beyond gene duplication: A model for incorporating horizontal gene transfer

Understanding the evolution of enzyme function after gene duplication has been a major goal of molecular biologists, biochemists and evolutionary biologists alike, for almost half a century. In contrast, the impact that horizontal gene transfer (HGT) has had on the evolution of enzyme specialization and the assembly of metabolic networks has just started to being investigated. Traditionally, evolutionary studies of enzymes have been limited to either the function of enzymes in vitro, or to sequence variability at the population level, where in almost all cases the starting conceptual framework embraces gene duplication as the mechanism responsible for the appearance of genetic redundancy. Very recently, we merged comparative phylogenomics, detection of selection signals, enzyme kinetics, X-ray crystallography and computational molecular dynamics, to characterize the sub-functionalization process of an amino acid biosynthetic enzyme prompted by an episode of HGT in bacteria. Some of the evolutionary implications of these functional studies, including a proposed model of enzyme specialization independent of gene duplication, are developed in this commentary.

Evolution of substrate specificity in a recipient's enzyme following horizontal gene transfer

Despite the prominent role of horizontal gene transfer (HGT) in shaping bacterial metabolism, little is known about the impact of HGT on the evolution of enzyme function. Specifically, what is the influence of a recently acquired gene on the function of an existing gene? For example, certain members of the genus Corynebacterium have horizontally acquired a whole l-tryptophan biosynthetic operon, whereas in certain closely related actinobacteria, for example, Mycobacterium, the trpF gene is missing. In Mycobacterium, the function of the trpF gene is performed by a dual-substrate (βα)8 phosphoribosyl isomerase (priA gene) also involved in l-histidine (hisA gene) biosynthesis. We investigated the effect of a HGT-acquired TrpF enzyme upon PriA's substrate specificity in Corynebacterium through comparative genomics and phylogenetic reconstructions. After comprehensive in vivo and enzyme kinetic analyses of selected PriA homologs, a novel (βα)8 isomerase subfamily with a specialized function in l-histidine biosynthesis, termed subHisA, was confirmed. X-ray crystallography was used to reveal active-site mutations in subHisA important for narrowing of substrate specificity, which when mutated to the naturally occurring amino acid in PriA led to gain of function. Moreover, in silico molecular dynamic analyses demonstrated that the narrowing of substrate specificity of subHisA is concomitant with loss of ancestral protein conformational states. Our results show the importance of HGT in shaping enzyme evolution and metabolism.

What can genome-scale metabolic network reconstructions do for prokaryotic systematics?

It has recently been proposed that in addition to Nomenclature, Classification and Identification, Comprehending Microbial Diversity may be considered as the fourth tenet of microbial systematics [Staley JT (2010) The Bulletin of BISMiS, 1(1): 1-5]. As this fourth goal implies a fundamental understanding of microbial speciation, this perspective article argues that translation of bacterial genome sequences into metabolic features may contribute to the development of modern polyphasic taxonomic approaches. Genome-scale metabolic network reconstructions (GSMRs), which are the result of computationally predicted and experimentally confirmed stoichiometric matrices incorporating all enzyme and metabolite components encoded by a genome sequence, provide a platform that can illustrate bacterial speciation. As the topology and the composition of GSMRs are expected to be the result of adaptive evolution, the features of these networks may provide the prokaryotic taxonomist with novel tools for reaching the fourth tenet of microbial systematics. Through selected examples from the Actinobacteria, which have been inferred from GSMRs and experimentally confirmed after phenotypic characterisation, it will be shown that this level of information can be incorporated into modern polyphasic taxonomic approaches. In conclusion, three specific examples are illustrated to show how GSMRs will revolutionize prokaryotic systematics, as has previously occurred in many other fields of microbiology.

Identification and analysis of residues contained on beta --> alpha loops of the dual-substrate (beta alpha)8 phosphoribosyl isomerase A specific for its phosphoribosyl anthranilate isomerase activity

A good model to experimentally explore evolutionary hypothesis related to enzyme function is the ancient-like dual-substrate (beta alpha)(8) phosphoribosyl isomerase A (PriA), which takes part in both histidine and tryptophan biosynthesis in Streptomyces coelicolor and related organisms. In this study, we determined the Michaelis-Menten enzyme kinetics for both isomerase activities in wild-type PriA from S. coelicolor and in selected single-residue monofunctional mutants, identified after Escherichia coli in vivo complementation experiments. Structural and functional analyses of a hitherto unnoticed residue contained on the functionally important beta --> alpha loop 5, namely, Arg(139), which was postulated on structural grounds to be important for the dual-substrate specificity of PriA, is presented for the first time. Indeed, enzyme kinetics analyses done on the mutant variants PriA_Ser(81)Thr and PriA_Arg(139)Asn showed that these residues, which are contained on beta --> alpha loops and in close proximity to the N-terminal phosphate-binding site, are essential solely for the phosphoribosyl anthranilate isomerase activity of PriA. Moreover, analysis of the X-ray crystallographic structure of PriA_Arg(139)Asn elucidated at 1.95 A herein strongly implicates the occurrence of conformational changes in this beta --> alpha loop as a major structural feature related to the evolution of the dual-substrate specificity of PriA. It is suggested that PriA has evolved by tuning a fine energetic balance that allows the sufficient degree of structural flexibility needed for accommodating two topologically dissimilar substrates--within a bifunctional and thus highly constrained active site--without compromising its structural stability.

The structure/function relationship of a dual-substrate (betaalpha)8-isomerase

Two structures of phosphoribosyl isomerase A (PriA) from Streptomyces coelicolor, involved in both histidine and tryptophan biosynthesis, were solved at 1.8A resolution. A closed conformer was obtained, which represents the first complete structure of PriA, revealing hitherto unnoticed molecular interactions and the occurrence of conformational changes. Inspection of these conformers, including ligand-docking simulations, allowed identification of residues involved in substrate recognition, chemical catalysis and conformational changes. These predictions were validated by mutagenesis and functional analysis. Arg19 and Ser81 were shown to play critical roles within the carboxyl and amino phosphate-binding sites, respectively; the catalytic residues Asp11 and Asp130 are responsible for both activities; and Thr166 and Asp171, which make an unusual contact, are likely to elicit the conformational changes needed for adopting the active site architectures. This represents the first report of the structure/function relationship of this (betaalpha)8-isomerase.