Molecular evolution of paclitaxel biosynthetic genes <Emphasis Type="Italic">TS</Emphasis> and <Emphasis Type="Italic">DBAT</Emphasis> of <Emphasis Type="Italic">Taxus</Emphasis> species

Evolutionary patterns of sequence divergence were analyzed in genes from the conifer genus Taxus (yew), encoding paclitaxel biosynthetic enzymes taxadiene synthase (TS) and 10-deacetylbaccatin III-10β-O-acetyltransferase (DBAT). N-terminal fragments of TS, full-length DBAT and internal transcribed spacer (ITS) were amplified from 15 closely related Taxus species and sequenced. Premature stop codons were not found in TS and DBAT sequences. Codon usage bias was not found, suggesting that synonymous mutations are selectively neutral. TS and DBAT gene trees are not consistent with the ITS tree, where species formed monophyletic clades. In fact, for both genes, alleles were sometimes shared across species and parallel amino acid substitutions were identified. While both TS and DBAT are, overall, under purifying selection, we identified a number of amino acids of TS under positive selection based on inference using maximum likelihood models. Positively selected amino acids in the N-terminal region of TS suggest that this region might be more important for enzyme function than previously thought. Moreover, we identify lineages with significantly elevated rates of amino acid substitution using a genetic algorithm. These findings demonstrate that the pattern of adaptive paclitaxel biosynthetic enzyme evolution can be documented between closely related Taxus species, where species-specific taxane metabolism has evolved recently.     

Do Amino Acid Biosynthetic Costs Constrain Protein Evolution in <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>?

Prokaryotic organisms preferentially utilize less energetically costly amino acids in highly expressed genes. Studies have shown that the proteome of Saccharomyces cerevisiae also exhibits this behavior, but only in broad terms. This study examines the question of metabolic efficiency as a proteome-shaping force at a finer scale, examining whether trends consistent with cost minimization as an evolutionary force are present independent of protein function and amino acid physicochemical property, and consistently with respect to amino acid biosynthetic costs. Inverse correlations between the average amino acid biosynthetic cost of the protein product and the levels of gene expression in S. cerevisiae are consistent with natural selection to minimize costs. There are, however, patterns of amino acid usage that raise questions about the strength (and possibly the universality) of this selective force in shaping S. cerevisiae’s proteome.     

Apparent longevity-related adaptation of mitochondrial amino acid content is due to nucleotide compositional shifts

Recent studies across animal phyla have suggested a possible link between amino acid compositional shifts and adaptive evolution across mitochondrial proteomes enabling longer lifespans. These studies examined associations of a gradual loss of cysteine (Cys) residues, increased usage of methionine (Met), and increased usage of threonine (Thr), with the evolution of longevity. Here, we examine all three hypotheses in a framework that considers nucleotide composition. We find that nucleotide composition is strongly correlated across codon positions, and with the above amino acid frequency patterns. We also find that the ND6 gene, which in vertebrates is the only mitochondrial gene situated on the “light-strand” shows no significant pattern for any of the amino acid associations. We also reasoned that if the mitochondrially-encoded proteins of oxidative phosphorylation (OXPHOS) were under selection for such shifts, then nuclear-encoded components should also reflect such pressure. However, we found non-correspondence of these patterns in the nuclear genes when compared to the mitochondrial genes previously associated with positive selection. These results are strongly suggestive of mutational bias, or less efficient purifying selection, as the primary driver of whole proteome shifts in amino acid composition.     

Evolution of tryptophan biosynthetic pathway in microbial genomes: a comparative genetic study

Biosynthetic pathway evolution needs to consider the evolution of a group of genes that code for enzymes catalysing the multiple chemical reaction steps leading to the final end product. Tryptophan biosynthetic pathway has five chemical reaction steps that are highly conserved in diverse microbial genomes, though the genes of the pathway enzymes show considerable variations in arrangements, operon structure (gene fusion and splitting) and regulation. We use a combined bioinformatic and statistical analyses approach to address the question if the pathway genes from different microbial genomes, belonging to a wide range of groups, show similar evolutionary relationships within and between them. Our analyses involved detailed study of gene organization (fusion/splitting events), base composition, relative synonymous codon usage pattern of the genes, gene expressivity, amino acid usage, etc. to assess inter- and intra-genic variations, between and within the pathway genes, in diverse group of microorganisms. We describe these genetic and genomic variations in the tryptophan pathway genes in different microorganisms to show the similarities across organisms, and compare the same genes across different organisms to find the possible variability arising possibly due to horizontal gene transfers. Such studies form the basis for moving from single gene evolution to pathway evolutionary studies that are important steps towards understanding the systems biology of intracellular pathways.     

Cost-Minimization of Amino Acid Usage

The negative correlation between the frequencies of usage of amino acids and their biosynthetic cost suggests that organisms minimize costs of protein biosynthesis. Empirical results support that: (1) free-living organisms (Archaea, Bacteria, and Eucaryota) minimize the usage of heavy amino acids more than intracellular organisms (viruses, chloroplasts, and mitochondria), a result confirmed by comparing intracellular Bacteria with other Bacteria; (2) avoidance of amino acids with low impact on protein structure (Chou-Fasman indices) is greater than for those with equal molecular weight but greater structural impact: constraints on protein function limit cost-minimization; (3) amino acid weight minimization (WM) for a protein correlates positively with the protein's expression level and with its size; (4) preliminary results suggest that for different proteins, the evolutionary rate of amino acid replacements correlates negatively with WM in these proteins; (5) results suggest that WM decreases with genome-size; and (6) developmental rates correlate positively with WM (within primates and rodents), even after confounding factors were accounted for. Effects of biosynthetic cost-minimization at whole-organism levels vary with metabolic and ecological strategies. Biosynthetic cost-minimization is an adaptive hypothesis that yields a semi-mechanistic explanation for small differences in allele fitness.     

Comparative analysis of amino acid usage and protein length distribution between alternatively and non-alternatively spliced genes across six eukaryotic genomes

Alternative splicing has been discovered in nearly all metazoan organisms as a mechanism to increase the diversity of gene products. However, the origin and evolution of alternatively spliced genes are still poorly understood. To understand the mechanisms for the evolution of alternatively spliced genes, it may be important to study the differences between alternatively and non-alternatively spliced genes. The aim of this research was to compare amino acid usage and protein length distribution between alternatively and non-alternatively spliced genes across six nearly complete eukaryotic genomes, including those of human (Homo sapiens), mouse (Mus musculus), rat (Rattus norvegicus), fruit fly (Drosophila melanogaster), Caenorhabditis elegans, and bovine (Bos taurus). Our results have suggested the following: (1) across the six species, alternatively and non-alternatively spliced genes have very similar tendency for amino acids usage for not only the overall scale but also those highly expressed genes, with all of the highly expressed genes having preferred amino acids including A, E, G, K, L, P, S, V, R, T, and D. (2) For not only the overall genes but also those highly expressed ones, the average length of the protein products of alternatively spliced genes is significantly greater than that of non-alternatively spliced ones. In contrast, distributions of protein lengths for the two groups of genes are very similar among all six species. Based on these results, we propose that alternatively spliced genes may have originated from non-alternatively spliced ones through events such as DNA mutations or gene fusion.     

Amino Acid Biosynthetic Cost and Protein Conservation

Protein products of highly expressed genes tend to favor amino acids that have lower average biosynthetic costs (i.e., they exhibit metabolic efficiency). While this trend has been observed in several studies, the specific sites where cost-reducing substitutions accumulate have not been well characterized. Toward that end, weighted costs in conserved and variable positions were evaluated across a total of 9,119 homologous proteins in four mammalian orders (primate, carnivore, rodent, and artiodactyls), which together contain a total of 20,457,072 amino acids. Degree of conservation at homologous positions in these mammalian proteins and average-weighted cost across all positions within a single protein are significantly correlated. Dividing human genes into two classes (those with and those without CpG islands in their promoters) suggests that humans also preferentially utilize less costly amino acids in highly expressed genes. In contrast to the intuitive expectation that the relatively weak selective force associated with metabolic efficiency would be a selection pressure in complex multicellular organisms, the overall level of selective constraint within the variable regions of mammalian proteins allows the metabolic efficiency to derive a reduction of overall biosynthetic cost, particularly in genes with the highest levels of expression.     

Reinvestigating the codon and amino acid usage of S. cerevisiae genome: a new insight from protein secondary structure analysis

Biased usage of synonymous codons has been elucidated under the perspective of cellular tRNA abundance for quite a long time now. Taking advantage of publicly available gene expression data for Saccharomyces cerevisiae, a systematic analysis of the codon and amino acid usages in two different coding regions corresponding to the regular (helix and strand) as well as the irregular (coil) protein secondary structures, have been performed. Our analyses suggest that apart from tRNA abundance, mRNA folding stability is another major evolutionary force in shaping the codon and amino acid usage differences between the highly and lowly expressed genes in S. cerevisiae genome and surprisingly it depends on the coding regions corresponding to the secondary structures of the encoded proteins. This is obviously a new paradigm in understanding the codon usage in S. cerevisiae. Differential amino acid usage between highly and lowly expressed genes in the regions coding for the irregular protein secondary structure in S. cerevisiae is expounded by the stability of the mRNA folded structure. Irrespective of the protein secondary structural type, the highly expressed genes always tend to encode cheaper amino acids in order to reduce the overall biosynthetic cost of production of the corresponding protein. This study supports the hypothesis that the tRNA abundance is a consequence of and not a reason for the biased usage of amino acid between highly and lowly expressed genes.     

Metabolic costs of amino acid and protein production in Escherichia coli

Escherichia coli is the most popular microorganism for the production of recombinant proteins and is gaining increasing importance for the production of low-molecular weight compounds such as amino acids. The metabolic cost associated with the production of amino acids and (recombinant) proteins from glucose, glycerol and acetate was determined using three different computational techniques to identify those amino acids that put the highest burden on the biosynthetic machinery of E. coli. Comparing the costs of individual amino acids, we find that methionine is the most expensive amino acid in terms of consumed mol of ATP per molecule produced, while leucine is the most expensive amino acid when taking into account the cellular abundances of amino acids. Moreover, we show that the biosynthesis of a large number of amino acids from glucose and particularly from glycerol provides a surplus of energy, which can be used to balance the high energetic cost of amino acid polymerization.     

Hydrophobicity and Aromaticity Are Primary Factors Shaping Variation in Amino Acid Usage of Chicken Proteome

Amino acids are utilized with different frequencies both among species and among genes within the same genome. Up to date, no study on the amino acid usage pattern of chicken has been performed. In the present study, we carried out a systematic examination of the amino acid usage in the chicken proteome. Our data indicated that the relative amino acid usage is positively correlated with the tRNA gene copy number. GC contents, including GC1, GC2, GC3, GC content of CDS and GC content of the introns, were correlated with the most of the amino acid usage, especially for GC rich and GC poor amino acids, however, multiple linear regression analyses indicated that only approximately 10–40% variation of amino acid usage can be explained by GC content for GC rich and GC poor amino acids. For other intermediate GC content amino acids, only approximately 10% variation can be explained. Correspondence analyses demonstrated that the main factors responsible for the variation of amino acid usage in chicken are hydrophobicity, aromaticity and genomic GC content. Gene expression level also influenced the amino acid usage significantly. We argued that the amino acid usage of chicken proteome likely reflects a balance or near balance between the action of selection, mutation, and genetic drift.     

A strong effect of AT mutational bias on amino acid usage in Buchnera is mitigated at high-expression genes

The advent of full genome sequences provides exceptionally rich data sets to explore molecular and evolutionary mechanisms that shape divergence among and within genomes. In this study, we use multivariate analysis to determine the processes driving genome-wide patterns of amino usage in the obligate endosymbiont Buchnera and its close free-living relative Escherichia coli. In the AT-rich Buchnera genome, the primary source of variation in amino acid usage differentiates high- and low-expression genes. Amino acids of high-expression Buchnera genes are generally less aromatic and use relatively GC-rich codons, suggesting that selection against aromatic amino acids and against amino acids with AT-rich codons is stronger in high-expression genes. Selection to maintain hydrophobic amino acids in integral membrane proteins is a primary factor driving protein evolution in E. coli but is a secondary factor in Buchnera. In E. coli, gene expression is a secondary force driving amino acid usage, and a correlation with tRNA abundance suggests that translational selection contributes to this effect. Although this and previous studies demonstrate that AT mutational bias and genetic drift influence amino acid usage in Buchnera, this genome-wide analysis argues that selection is sufficient to affect the amino acid content of proteins with different expression and hydropathy levels.     

Evolutionary Patterns of Codon Usage in the Chloroplast Gene rbcL

In this study we reconstruct the evolution of codon usage bias in the chloroplast gene rbcL using a phylogeny of 92 green-plant taxa. We employ a measure of codon usage bias that accounts for chloroplast genomic nucleotide content, as an attempt to limit plausible explanations for patterns of codon bias evolution to selection- or drift-based processes. This measure uses maximum likelihood-ratio tests to compare the performance of two models, one in which a single codon is overrepresented and one in which two codons are overrepresented. The measure allowed us to analyze both the extent of bias in each lineage and the evolution of codon choice across the phylogeny. Despite predictions based primarily on the low G+C content of the chloroplast and the high functional importance of rbcL, we found large differences in the extent of bias, suggesting differential molecular selection that is clade specific. The seed plants and simple leafy liverworts each independently derived a low level of bias in rbcL, perhaps indicating relaxed selectional constraint on molecular changes in the gene. Overrepresentation of a single codon was typically plesiomorphic, and transitions to overrepresentation of two codons occurred commonly across the phylogeny, possibly indicating biochemical selection. The total codon bias in each taxon, when regressed against the total bias of each amino acid, suggested that twofold amino acids play a strong role in inflating the level of codon usage bias in rbcL, despite the fact that twofolds compose a minority of residues in this gene. Those amino acids that contributed most to the total codon usage bias of each taxon are known through amino acid knockout and replacement to be of high functional importance. This suggests that codon usage bias may be constrained by particular amino acids and, thus, may serve as a good predictor of what residues are most important for protein fitness. Present address (Joshua T. Herbeck): JBP Center for Comparative Molecular Biology and Evolution, Marine Biological Laboratory, Woods Hole, MA 02543, USA     

Evolutionary Patterning of Hemagglutinin Gene Sequence of 2009 H1N1 Pandemic

The 2009 H1N1 swine flu is the first pandemic in decades. Infectivity of the influenza virus for human host depends largely on its ability to evade antibodies specific for viral protein called hemagglutinin (HA) that mediates attachment to the host. In the present study we analysed large number of HA gene sequences available in Flu Database maintained at NCBI. Our sequence based analysis clearly demonstrates that the amino acid usage pattern may dramatically change during the course of evolution, and there exists a clear link between a particular pattern of amino acid usage of HA genes and its potential to become infectious. Structural studies revealed how binding efficiency between the HA and sialic acid may alter the pandemic potential of infection. Our work highlights the evolutionary significance and biochemical basis of the selective advantage of certain amino acids of HA in 2009 and provides a link between the characteristics changes in HA protein and their potential to pronounce a global menace to public health.     

Comparative analysis of essential genes and nonessential genes in Escherichia coli K12

Genes can be classified as essential or nonessential based on their indispensability for a living organism. Previous researches have suggested that essential genes evolve more slowly than nonessential genes and the impact of gene dispensability on a gene’s evolutionary rate is not as strong as expected. However, findings have not been consistent and evidence is controversial regarding the relationship between the gene indispensability and the rate of gene evolution. Understanding how different classes of genes evolve is essential for a full understanding of evolutionary biology, and may have medical relevance in the design of new antibacterial agents. We therefore performed an investigation into the properties of essential and nonessential genes. Analysis of evolutionary conservation, protein length distribution and amino acid usage between essential and nonessential genes in Escherichia coli K12 demonstrated that essential genes are relatively preserved throughout the bacterial kingdom when compared to nonessential genes. Furthermore, results show that essential genes, compared to nonessential genes, have a significantly higher proportion of large (>534 amino acids) and small proteins (<139 amino acids) relative to medium-sized proteins. The pattern of amino acids usage shows a similar trend for essential and nonessential genes, although some notable exceptions are observed. These findings help to clarify our understanding of the evolutionary mechanisms of essential and nonessential genes, relevant to the study of mutagenesis and possibly allowing prediction of gene properties in other poorly understood organisms.     

Selection on Synthesis Cost Affects Interprotein Amino Acid Usage in All Three Domains of Life

Most investigations of the forces shaping protein evolution have focussed on protein function. However, cells are typically 50%–75% protein by dry weight, with protein expression levels distributed over five orders of magnitude. Cells may, therefore, be under considerable selection pressure to incorporate amino acids that are cheap to synthesize into proteins that are highly expressed. Such selection pressure has been demonstrated to alter amino acid usage in a few organisms, but whether “cost selection” is a general phenomenon remains unknown. One reason for this is that reliable protein expression level data is not available for most organisms. Accordingly, I have developed a new method for detecting cost selection. This method depends solely on interprotein gradients in amino acid usage. Applying it to an analysis of 43 whole genomes from all three domains of life, I show that selection on the synthesis cost of amino acids is a pervasive force in shaping the composition of proteins. Moreover, some amino acids have different price tags for different organisms—the cost of amino acids is changed for organisms living in hydrothermal vents compared with those living at the sea surface or for organisms that have difficulty acquiring elements such as nitrogen compared with those that do not—so I also investigated whether differences between organisms in amino acid usage might reflect differences in synthesis or acquisition costs. The results suggest that organisms evolve to alter amino acid usage in response to environmental conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Hector Musto]     

Codon usage bias in prokaryotic pyrimidine-ending codons is associated with the degeneracy of the encoded amino acids

Synonymous codons are unevenly distributed among genes, a phenomenon termed codon usage bias. Understanding the patterns of codon bias and the forces shaping them is a major step towards elucidating the adaptive advantage codon choice can confer at the level of individual genes and organisms. Here, we perform a large-scale analysis to assess codon usage bias pattern of pyrimidine-ending codons in highly expressed genes in prokaryotes. We find a bias pattern linked to the degeneracy of the encoded amino acid. Specifically, we show that codon-pairs that encode two- and three-fold degenerate amino acids are biased towards the C-ending codon while codons encoding four-fold degenerate amino acids are biased towards the U-ending codon. This codon usage pattern is widespread in prokaryotes, and its strength is correlated with translational selection both within and between organisms. We show that this bias is associated with an improved correspondence with the tRNA pool, avoidance of mis-incorporation errors during translation and moderate stability of codon-anticodon interaction, all consistent with more efficient translation.     

Yeast engineering to express sex pheromone gland genes of the oriental fruit moth,Grapholita molesta

Mating disruption by using sex pheromone is an ecofriendly alternative way to control insect pests. To be effective, large amounts of sex pheromone are needed, leading to a relatively high production cost. To reduce the cost for chemical synthesis of sex pheromone, yeast engineering technology has been devised. This study used a baker's yeast, Saccharomyces cerevisiae, to express genes associated with sex pheromone biosynthesis of the Oriental fruit moth, Grapholita molesta. Compared to other fatty acid biosynthetic pathways, two steps that are unique to pheromone gland of G. molesta are proposed: desaturation at even number catalyzed by desaturase (Gm-DES) and terminal reduction catalyzed by fatty acyl reductase (Gm-FAR). Gm-DES and Gm-FAR were cloned into a yeast expression vector, pYES2.1. They were used to transform S. cerevisiae by a double transfection method. The transformed yeast was induced with 2% galactose to over-express these two exogenous genes. Their expression was confirmed by RT-PCR and western blotting. To facilitate pheromone production, transformed yeasts were supplied with myristic acid during over-expression. Resulting fatty acid composition was analyzed by GC-MS after fatty acid methyl ester derivatization. Control yeast produced mostly saturated fatty acids. However, a single gene (Gm-DES)-transformed yeast produced unsaturated fatty acids at ∆9 such as Z9-tetradecenoic acid (Z9-14:1), palmitoleic acid (Z9-16:1), and oleic acid (Z9-18:1) in addition to saturated fatty acids. The double-transformed yeast produced an additional component, alcohol form of oleic acid (Z9-18:OH). These results suggest that Gm-DES can catalyze desaturation of fatty acids at ∆9 and Gm-FAR can reduce terminal carboxylic acid into alcohol.