Trade-Offs Predicted by Metabolic Network Structure Give Rise to Evolutionary Specialization and Phenotypic Diversification

Mitigating trade-offs between different resource-utilization functions is key to an organism’s ecological and evolutionary success. These trade-offs often reflect metabolic constraints with a complex molecular underpinning; therefore, their consequences for evolutionary processes have remained elusive. Here, we investigate how metabolic architecture induces resource-utilization constraints and how these constraints, in turn, elicit evolutionary specialization and diversification. Guided by the metabolic network structure of the bacterium Lactococcus cremoris, we selected two carbon sources (fructose and galactose) with predicted coutilization constraints. By evolving L. cremoris on either fructose, galactose, or a mix of both sugars, we imposed selection favoring divergent metabolic specializations or coutilization of both resources, respectively. Phenotypic characterization revealed the evolution of either fructose or galactose specialists in the single-sugar treatments. In the mixed-sugar regime, we observed adaptive diversification: both specialists coexisted, and no generalist evolved. Divergence from the ancestral phenotype occurred at key pathway junctions in the central carbon metabolism. Fructose specialists evolved mutations in the fbp and pfk genes that appear to balance anabolic and catabolic carbon fluxes. Galactose specialists evolved increased expression of pgmA (the primary metabolic bottleneck of galactose metabolism) and silencing of ptnABCD (the main glucose transporter) and ldh (regulator/enzyme of downstream carbon metabolism). Overall, our study shows how metabolic network architecture and historical contingency serve to predict targets of selection and inform the functional interpretation of evolved mutations. The elucidation of the relationship between molecular constraints and phenotypic trade-offs contributes to an integrative understanding of evolutionary specialization and diversification.     

Higher Gene Duplicabilities for Metabolic Proteins Than for Nonmetabolic Proteins in Yeast and <Emphasis Type="Italic">E. coli</Emphasis>

Although the evolutionary significance of gene duplication has long been appreciated, it remains unclear what factors determine gene duplicability. In this study we investigated whether metabolism is an important determinant of gene duplicability because cellular metabolism is crucial for the survival and reproduction of an organism. Using genomic data and metabolic pathway data from the yeast (Saccharomyces cerevisiae) and Escherichia coli, we found that metabolic proteins indeed tend to have higher gene duplicability than nonmetabolic proteins. Moreover, a detailed analysis of metabolic pathways in these two organisms revealed that genes in the central metabolic pathways and the catabolic pathways have, on average, higher gene duplicability than do other genes and that most genes in anabolic pathways are single-copy genes.Reviewing Editor: Dr. Rüdiger Cerff     

Relaxed Evolution in the Tyrosine Aminotransferase Gene Tat in Old World Fruit Bats (Chiroptera: Pteropodidae)

Frugivorous and nectarivorous bats fuel their metabolism mostly by using carbohydrates and allocate the restricted amounts of ingested proteins mainly for anabolic protein syntheses rather than for catabolic energy production. Thus, it is possible that genes involved in protein (amino acid) catabolism may have undergone relaxed evolution in these fruit- and nectar-eating bats. The tyrosine aminotransferase (TAT, encoded by the Tat gene) is the rate-limiting enzyme in the tyrosine catabolic pathway. To test whether the Tat gene has undergone relaxed evolution in the fruit- and nectar-eating bats, we obtained the Tat coding region from 20 bat species including four Old World fruit bats (Pteropodidae) and two New World fruit bats (Phyllostomidae). Phylogenetic reconstructions revealed a gene tree in which all echolocating bats (including the New World fruit bats) formed a monophyletic group. The phylogenetic conflict appears to stem from accelerated TAT protein sequence evolution in the Old World fruit bats. Our molecular evolutionary analyses confirmed a change in the selection pressure acting on Tat, which was likely caused by a relaxation of the evolutionary constraints on the Tat gene in the Old World fruit bats. Hepatic TAT activity assays showed that TAT activities in species of the Old World fruit bats are significantly lower than those of insectivorous bats and omnivorous mice, which was not caused by a change in TAT protein levels in the liver. Our study provides unambiguous evidence that the Tat gene has undergone relaxed evolution in the Old World fruit bats in response to changes in their metabolism due to the evolution of their special diet.     

METABOLIC FLUX IS A DETERMINANT OF THE EVOLUTIONARY RATES OF ENZYME‐ENCODING GENES

Relationships between evolutionary rates and gene properties on a genomic, functional, pathway, or system level are being explored to unravel the principles of the evolutionary process. In particular, functional network properties have been analyzed to recognize the constraints they may impose on the evolutionary fate of genes. Here we took as a case study the core metabolic network in human erythrocytes and we analyzed the relationship between the evolutionary rates of its genes and the metabolic flux distribution throughout it. We found that metabolic flux correlates with the ratio of nonsynonymous to synonymous substitution rates. Genes encoding enzymes that carry high fluxes have been more constrained in their evolution, while purifying selection is more relaxed in genes encoding enzymes carrying low metabolic fluxes. These results demonstrate the importance of considering the dynamical functioning of gene networks when assessing the action of selection on system‐level properties.     

Evolutionary rate patterns of the Gibberellin pathway genes

Background  

Analysis of molecular evolutionary patterns of different genes within metabolic pathways allows us to determine whether these genes are subject to equivalent evolutionary forces and how natural selection shapes the evolution of proteins in an interacting system. Although previous studies found that upstream genes in the pathway evolved more slowly than downstream genes, the correlation between evolutionary rate and position of the genes in metabolic pathways as well as its implications in molecular evolution are still less understood.     

Evolution of arginine deiminase (ADI) pathway genes

We have analyzed the evolution of the three genes encoding structural enzymes of the arginine deiminase (ADI) pathway, arginine deiminase (ADI), ornithine transcarbamoylase (OTC), and carbamate kinase (CK) in a wide range of organisms, including Archaea, Bacteria, and Eukarya. This catabolic route was probably present in the last common ancestor to all the domains of life. The results obtained indicate that these genes have undergone a complex evolutionary history, including horizontal transfer events, duplications, and losses. Therefore, these genes are not adequate to infer organismal relationships at deep branching levels, but they provide an insight into how catabolic genes evolved and were assembled into metabolic pathways. Our results suggest that the three genes evolved independently and were later assembled into a single cluster with functional interdependence, thus, providing support for the gene recruitment hypothesis. Furthermore, the molecular phylogenetic analysis of OTC suggests a new classification of these genes into three subfamilies.     

Rapid evolution through gene duplication and subfunctionalization of the testes-specific alpha4 proteasome subunits in Drosophila

Gene duplication is an important mechanism for acquiring new genes and creating genetic novelty in organisms. Evidence suggests that duplicated genes are retained at a much higher rate than originally thought and that functional divergence of gene copies is a major factor promoting their retention in the genome. We find that two Drosophila testes-specific alpha4 proteasome subunit genes (alpha4-t1 and alpha4-t2) have a higher polymorphism within species and are significantly more diverged between species than the somatic alpha4 gene. Our data suggest that following gene duplication, the alpha4-t1 gene experienced relaxed selective constraints, whereas the alpha4-t2 gene experienced positive selection acting on several codons. We report significant heterogeneity in evolutionary rates among all three paralogs at homologous codons, indicating that functional divergence has coincided with genic divergence. Reproductive subfunctionalization may allow for a more rapid evolution of reproductive traits and a greater specialization of testes function. Our data add to the increasing evidence that duplicated genes experience lower selective constraints and in some cases positive selection following duplication. Newly duplicated genes that are freer from selective constraints may provide a mechanism for developing new interactions and a pathway for the evolution of new genes.     

Experiments in microbial evolution: new enzymes, new metabolic activities

Biological evolution has resulted in a richness and diversity of species. Among microorganisms this is most evident in the wealth and diversity of biochemical transformations. Evidence for evolutionary relationships may be obtained from comparative studies, but with microorganisms it is also possible to follow evolution in action. Microbial populations adapt rapidly to changes in the environment and the evolution of new metabolic activities can be observed in laboratory experiments. The enzymes of many catabolic pathways are synthesized in response to the presence of inducing substrates. New catabolic activities may be acquired by mutations in regulatory genes resulting in alterations in the specificity of induction, or in enzyme synthesis in the absence of inducer. Mutations in structural genes may given rise to enzymes with altered substrate specificities. In bacteria, catabolic genes may be carried on plasmids and the exchange of plasmids among bacterial populations increases the evolutionary potential. Experiments in microbial evolution have produced strains with novel catabolic activities involving regulatory or structural gene mutations, gene duplications and plasmid exchange. Enzymes studied in this way include amidase, ribitol dehydrogenase, evolved beta-galactosidase, and enzymes of the catabolic pathways for pentoses and pentitols and haloaromatic compounds.     

Extensive Functional Diversification of the Populus Glutathione S-Transferase Supergene Family

Identifying how genes and their functions evolve after duplication is central to understanding gene family radiation. In this study, we systematically examined the functional diversification of the glutathione S-transferase (GST) gene family in Populus trichocarpa by integrating phylogeny, expression, substrate specificity, and enzyme kinetic data. GSTs are ubiquitous proteins in plants that play important roles in stress tolerance and detoxification metabolism. Genome annotation identified 81 GST genes in Populus that were divided into eight classes with distinct divergence in their evolutionary rate, gene structure, expression responses to abiotic stressors, and enzymatic properties of encoded proteins. In addition, when all the functional parameters were examined, clear divergence was observed within tandem clusters and between paralogous gene pairs, suggesting that subfunctionalization has taken place among duplicate genes. The two domains of GST proteins appear to have evolved under differential selective pressures. The C-terminal domain seems to have been subject to more relaxed functional constraints or divergent directional selection, which may have allowed rapid changes in substrate specificity, affinity, and activity, while maintaining the primary function of the enzyme. Our findings shed light on mechanisms that facilitate the retention of duplicate genes, which can result in a large gene family with a broad substrate spectrum and a wide range of reactivity toward different substrates.     

Identification of neurotransmitter receptor genes under significantly relaxed selective constraint by orthologous gene comparisons between humans and rodents

Neurotransmitter receptors (neuroreceptors) are classified into two types, G protein-coupled receptors (GPCRs) and ligand-gated ion channels. The former occupies a small part of the large GPCR superfamily, whereas the latter consists of three superfamilies. In these superfamilies, humans and rodents share almost the same set of neuroreceptor genes. This neuroreceptor gene set is good material to examine the degree of selective constraint exerted on each member gene of a given superfamily. If there are any neuroreceptor genes under the degree of selective constraint that is very different from that of the other member genes, they may have influenced the functional features characteristic of human neural activities. With the aim of identifying such neuroreceptor genes, we collected sequence data of orthologous neuroreceptor genes for humans, mice, and rats by database searches. This data set included ortholog pairs for 141 kinds of neuroreceptor genes, which covered almost the whole set of neuroreceptor genes known to be expressed in the human brain. The degree of selective constraint was estimated by computing the ratio (d(N)/d(S)) of the number of nonsynonymous substitutions to that of synonymous substitutions. We found that the d(N)/d(S) ratio ranged widely and its distribution fitted a gamma distribution. In particular, we found that four neuroreceptor genes are under the significantly relaxed selective constraint. They are an ionotropic glutamate receptor subunit NMDA-2C, two GABA(A) receptor subunits, i.e., GABA(A)-epsilon and GABA(A)-theta, and a dopamine receptor D4. Interestingly, these neuroreceptors have been reported to be associated with cognitive abilities such as memory and learning, and responsiveness to novel stimuli. These cognitive abilities can influence the behavioral features of an individual. Thus, it suggests that the relaxed-constraint neuroreceptor genes have evolved, assuring that the nervous system responds to a variety of stimuli with proper flexibility.     

Primary and quaternary structure of the catabolic ornithine carbamoyltransferase from Pseudomonas aeruginosa. Extensive sequence homology with the anabolic ornithine carbamoyltransferases of Escherichia coli

We have determined the complete nucleotide sequence of the arcB gene from Pseudomonas aeruginosa strain PAO and we have purified the arcB product, the catabolic ornithine carbamoyltransferase (EC 2.1.3.3), to apparent homogeneity from the same strain. The N-terminal amino acid sequence, the total amino acid composition and the subunit size of the purified enzyme were in agreement with nucleotide sequencing results, which predict a polypeptide of 336 amino acids (Mr 38,108). Crosslinking experiments suggest that the native enzyme (apparent Mr approx. 420,000) basically consists of a trimer aggregating to form nonamers or dodecamers. The arcB gene of P. aeruginosa had strong homology with the argF and argI genes which code for the anabolic ornithine carbamoyltransferase isoenzymes in Escherichia coli; 63% of the nucleotides and 57% of the amino acids were absolutely conserved in arcB and argF. This indicates a close evolutionary relationship between these genes although their products have different physiological functions in the cell. Under conditions of induction (energy depletion) the catabolic ornithine carbamoyltransferase represented greater than or equal to 10% of the total cellular protein. Like other highly expressed Pseudomonas genes, the arcB gene was found not to use seven codons which correspond to minor or weakly interacting tRNA species in E. coli.     

Faster-X Evolution of Gene Expression in Drosophila

DNA sequences on X chromosomes often have a faster rate of evolution when compared to similar loci on the autosomes, and well articulated models provide reasons why the X-linked mode of inheritance may be responsible for the faster evolution of X-linked genes. We analyzed microarray and RNA–seq data collected from females and males of six Drosophila species and found that the expression levels of X-linked genes also diverge faster than autosomal gene expression, similar to the “faster-X” effect often observed in DNA sequence evolution. Faster-X evolution of gene expression was recently described in mammals, but it was limited to the evolutionary lineages shortly following the creation of the therian X chromosome. In contrast, we detect a faster-X effect along both deep lineages and those on the tips of the Drosophila phylogeny. In Drosophila males, the dosage compensation complex (DCC) binds the X chromosome, creating a unique chromatin environment that promotes the hyper-expression of X-linked genes. We find that DCC binding, chromatin environment, and breadth of expression are all predictive of the rate of gene expression evolution. In addition, estimates of the intraspecific genetic polymorphism underlying gene expression variation suggest that X-linked expression levels are not under relaxed selective constraints. We therefore hypothesize that the faster-X evolution of gene expression is the result of the adaptive fixation of beneficial mutations at X-linked loci that change expression level in cis. This adaptive faster-X evolution of gene expression is limited to genes that are narrowly expressed in a single tissue, suggesting that relaxed pleiotropic constraints permit a faster response to selection. Finally, we present a conceptional framework to explain faster-X expression evolution, and we use this framework to examine differences in the faster-X effect between Drosophila and mammals.     

Genetic and structural organization of the aminophenol catabolic operon and its implication for evolutionary process

The aminophenol (AP) catabolic operon in Pseudomonas putida HS12 mineralizing nitrobenzene was found to contain all the enzymes responsible for the conversion of AP to pyruvate and acetyl coenzyme A via extradiol meta cleavage of 2-aminophenol. The sequence and functional analyses of the corresponding genes of the operon revealed that the AP catabolic operon consists of one regulatory gene, nbzR, and the following nine structural genes, nbzJCaCbDGFEIH, which encode catabolic enzymes. The NbzR protein, which is divergently transcribed with respect to the structural genes, possesses a leucine zipper motif and a MarR homologous domain. It was also found that NbzR functions as a repressor for the AP catabolic operon through binding to the promoter region of the gene cluster in its dimeric form. A comparative study of the AP catabolic operon with other meta cleavage operons led us to suggest that the regulatory unit (nbzR) was derived from the MarR family and that the structural unit (nbzJCaCbDGFEIH) has evolved from the ancestral meta cleavage gene cluster. It is also proposed that these two functional units assembled through a modular type gene transfer and then have evolved divergently to acquire specialized substrate specificities (NbzCaCb and NbzD) and catalytic function (NbzE), resulting in the creation of the AP catabolic operon. The evolutionary process of the AP operon suggests how bacteria have efficiently acquired genetic diversity and expanded their metabolic capabilities by modular type gene transfer.     

Human testis-specific genes are under relaxed negative selection

Recent studies have suggested that selective forces and constraints acting on genes varied during human evolution depending on the organ in which they are expressed. To gain insight into the evolution of organ determined negative selection forces, we compared the non-synonymous SNP diversity of genes expressed in different organs. Based on a HAPMAP dataset, we determined for each SNP its frequency in 11 human populations and, in each case, predicted whether or not the change it produces is deleterious. We have shown that, for all organs under study, SNPs predicted to be deleterious are present at a significantly lower frequency than SNPs predicted to be tolerated. However, testis-specific genes contain a higher proportion of deleterious SNPs than other organs. This study shows that negative selection is acting on the whole human genome, but that the action of negative selection is relaxed on testis-specific genes. This result adds to and expands the hypothesis of a recent evolutionary change in the human male reproductive system and its behavior.     

Patterns of variation during adaptation in functionally linked loci

An understanding of the distribution of natural patterns of genetic variation is relevant to such fundamental biological fields as evolution and development. One recent approach to understanding such patterns has been to focus on the constraints that may arise as a function of the network or pathway context in which genes are embedded. Despite theoretical expectations of higher evolutionary constraint for genes encoding upstream versus downstream enzymes in metabolic pathways, empirical results have varied. Here we combine two complementary models from population genetics and enzyme kinetics to explore genetic variation as a function of pathway position when selection acts on whole-pathway flux. We are able to qualitatively reproduce empirically observed patterns of polymorphism and divergence and suggest that expectations should vary depending on the evolutionary trajectory of a population. Upstream genes are initially more polymorphic and diverge faster after an environmental change, while we see the opposite trend as the population approaches its fitness optimum.