Signatures of optimal codon usage in metabolic genes inform budding yeast ecology

Reverse ecology is the inference of ecological information from patterns of genomic variation. One rich, heretofore underutilized, source of ecologically relevant genomic information is codon optimality or adaptation. Bias toward codons that match the tRNA pool is robustly associated with high gene expression in diverse organisms, suggesting that codon optimization could be used in a reverse ecology framework to identify highly expressed, ecologically relevant genes. To test this hypothesis, we examined the relationship between optimal codon usage in the classic galactose metabolism (GAL) pathway and known ecological niches for 329 species of budding yeasts, a diverse subphylum of fungi. We find that optimal codon usage in the GAL pathway is positively correlated with quantitative growth on galactose, suggesting that GAL codon optimization reflects increased capacity to grow on galactose. Optimal codon usage in the GAL pathway is also positively correlated with human-associated ecological niches in yeasts of the CUG-Ser1 clade and with dairy-associated ecological niches in the family Saccharomycetaceae. For example, optimal codon usage of GAL genes is greater than 85% of all genes in the genome of the major human pathogen Candida albicans (CUG-Ser1 clade) and greater than 75% of genes in the genome of the dairy yeast Kluyveromyces lactis (family Saccharomycetaceae). We further find a correlation between optimization in the GALactose pathway genes and several genes associated with nutrient sensing and metabolism. This work suggests that codon optimization harbors information about the metabolic ecology of microbial eukaryotes. This information may be particularly useful for studying fungal dark matter—species that have yet to be cultured in the lab or have only been identified by genomic material.

Bias toward codons that match the tRNA pool is robustly associated with high gene expression in diverse organisms, suggesting that codon optimization could be used in a reverse ecology framework to identify ecologically relevant genes. This study finds that this is indeed the case for 329 species of budding yeasts, suggesting that codon optimization harbors information about the metabolic ecology of microbial eukaryotes.     

Metabolic evolution of two reducing equivalent-conserving pathways for high-yield succinate production in Escherichia coli

Reducing equivalents are an important cofactor for efficient synthesis of target products. During metabolic evolution to improve succinate production in Escherichia coli strains, two reducing equivalent-conserving pathways were activated to increase succinate yield. The sensitivity of pyruvate dehydrogenase to NADH inhibition was eliminated by three nucleotide mutations in the lpdA gene. Pyruvate dehydrogenase activity increased under anaerobic conditions, which provided additional NADH. The pentose phosphate pathway and transhydrogenase were activated by increased activities of transketolase and soluble transhydrogenase SthA. These data suggest that more carbon flux went through the pentose phosphate pathway, thus leading to production of more reducing equivalent in the form of NADPH, which was then converted to NADH through soluble transhydrogenase for succinate production. Reverse metabolic engineering was further performed in a parent strain, which was not metabolically evolved, to verify the effects of activating these two reducing equivalent-conserving pathways for improving succinate yield. Activating pyruvate dehydrogenase increased succinate yield from 1.12 to 1.31 mol/mol, whereas activating the pentose phosphate pathway and transhydrogenase increased succinate yield from 1.12 to 1.33 mol/mol. Activating these two pathways in combination led to a succinate yield of 1.5 mol/mol (88% of theoretical maximum), suggesting that they exhibited a synergistic effect for improving succinate yield.     

Metabolic characterization and transformation of the non‐dairy Lactococcus lactis strain KF147, for production of ethanol from xylose

The non‐dairy lactic acid bacterium Lactococcus lactis KF147 can utilize xylose as the sole energy source. To assess whether KF147 could serve as a platform organism for converting second generation sugars into useful chemicals, the authors characterized growth and product formation for KF147 when grown on xylose. In a defined medium KF147 was found to co‐metabolize xylose and arginine, resulting in bi‐phasic growth. Especially at low xylose concentrations, arginine significantly improved growth rate. To facilitate further studies of the xylose metabolism, the authors eliminated arginine catabolism by deleting the arcA gene encoding the arginine deiminase. The fermentation product profile suggested two routes for xylose degradation, the phosphoketolase pathway and the pentose phosphate pathway. Inactivation of the phosphoketolase pathway redirected the entire flux through the pentose phosphate pathway whereas over‐expression of phosphoketolase increased the flux through the phosphoketolase pathway. In general, significant amounts of the mixed‐acid products, including lactate, formate, acetate and ethanol, were formed irrespective of xylose concentrations. To demonstrate the potential of KF147 for converting xylose into useful chemicals the authors chose to redirect metabolism towards ethanol production. A synthetic promoter library was used to drive the expression of codon‐optimized versions of the Zymomonas mobilis genes encoding pyruvate decarboxylase and alcohol dehydrogenase, and the outcome was a strain producing ethanol as the sole fermentation product with a high yield corresponding to 83% of the theoretical maximum. The results clearly indicate the great potential of using the more metabolically diverse non‐dairy L. lactis strains for bio‐production based on xylose containing feedstocks.     

Solving the riddle of codon usage preferences: a test for translational selection

Translational selection is responsible for the unequal usage of synonymous codons in protein coding genes in a wide variety of organisms. It is one of the most subtle and pervasive forces of molecular evolution, yet, establishing the underlying causes for its idiosyncratic behaviour across living kingdoms has proven elusive to researchers over the past 20 years. In this study, a statistical model for measuring translational selection in any given genome is developed, and the test is applied to 126 fully sequenced genomes, ranging from archaea to eukaryotes. It is shown that tRNA gene redundancy and genome size are interacting forces that ultimately determine the action of translational selection, and that an optimal genome size exists for which this kind of selection is maximal. Accordingly, genome size also presents upper and lower boundaries beyond which selection on codon usage is not possible. We propose a model where the coevolution of genome size and tRNA genes explains the observed patterns in translational selection in all living organisms. This model finally unifies our understanding of codon usage across prokaryotes and eukaryotes. Helicobacter pylori, Saccharomyces cerevisiae and Homo sapiens are codon usage paradigms that can be better understood under the proposed model.     

Further increased production of free fatty acids by overexpressing a predicted transketolase gene of the pentose phosphate pathway in Aspergillus oryzae faaA disruptant

Free fatty acids are useful as source materials for the production of biodiesel fuel and various chemicals such as pharmaceuticals and dietary supplements. Previously, we attained a 9.2-fold increase in free fatty acid productivity by disrupting a predicted acyl-CoA synthetase gene (faaA, AO090011000642) in Aspergillus oryzae. In this study, we achieved further increase in the productivity by overexpressing a predicted transketolase gene of the pentose phosphate pathway in the faaA disruptant. The A. oryzae genome is predicted to have three transketolase genes and overexpression of AO090023000345, one of the three genes, resulted in phenotypic change and further increase (corresponding to an increased production of 0.38 mmol/g dry cell weight) in free fatty acids at 1.4-fold compared to the faaA disruptant. Additionally, the biomass of hyphae increased at 1.2-fold by the overexpression. As a result, free fatty acid production yield per liter of liquid culture increased at 1.7-fold by the overexpression.     

The Effect of Multiple Evolutionary Selections on Synonymous Codon Usage of Genes in the Mycoplasma bovis Genome

Mycoplasma bovis is a major pathogen causing arthritis, respiratory disease and mastitis in cattle. A better understanding of its genetic features and evolution might represent evidences of surviving host environments. In this study, multiple factors influencing synonymous codon usage patterns in M. bovis (three strains’ genomes) were analyzed. The overall nucleotide content of genes in the M. bovis genome is AT-rich. Although the G and C contents at the third codon position of genes in the leading strand differ from those in the lagging strand (p<0.05), the 59 synonymous codon usage patterns of genes in the leading strand are highly similar to those in the lagging strand. The over-represented codons and the under-represented codons were identified. A comparison of the synonymous codon usage pattern of M. bovis and cattle (susceptible host) indicated the independent formation of synonymous codon usage of M. bovis. Principal component analysis revealed that (i) strand-specific mutational bias fails to affect the synonymous codon usage pattern in the leading and lagging strands, (ii) mutation pressure from nucleotide content plays a role in shaping the overall codon usage, and (iii) the major trend of synonymous codon usage has a significant correlation with the gene expression level that is estimated by the codon adaptation index. The plot of the effective number of codons against the G+C content at the third codon position also reveals that mutation pressure undoubtedly contributes to the synonymous codon usage pattern of M. bovis. Additionally, the formation of the overall codon usage is determined by certain evolutionary selections for gene function classification (30S protein, 50S protein, transposase, membrane protein, and lipoprotein) and translation elongation region of genes in M. bovis. The information could be helpful in further investigations of evolutionary mechanisms of the Mycoplasma family and heterologous expression of its functionally important proteins.     

Effect of transketolase modifications on carbon flow to the purine-nucleotide pathway in Corynebacterium ammoniagenes

Transketolase, one of the enzymes in the nonoxidative branch of the pentose phosphate pathway, operates to shuttle ribose 5-phosphate and glycolytic intermediates together with transaldolase, and might be involved in the availability of ribose 5-phosphate, a precursor of nucleotide biosynthesis. The tkt and tal genes encoding transketolase and transaldolase, respectively, were cloned from the typical nucleotide- and nucleoside-producing organism Corynebacterium ammoniagenes by a PCR approach using oligonucleotide primers derived from conserved regions of each amino acid sequence from other organisms. Enzymatic and molecular analyses revealed that the two genes were clustered on the genome together with the glucose 6-phosphate dehydrogenase gene (zwf). The effect of transketolase modifications on the production of inosine and 5'-xanthylic acid was investigated in industrial strains of C. ammoniagenes. Multiple copies of plasmid-borne tkt caused about tenfold increases in transketolase activity and resulted in 10-20% decreased yields of products relative to the parents. In contrast, site-specific disruption of tkt enabled both producers to accumulate 10-30% more products concurrently with a complete loss of transketolase activity and the expected phenotype of shikimate auxotrophy. These results indicate that transketolase normally shunts ribose 5-phosphate back into glycolysis in these biosynthetic processes and interception of this shunt allows cells to redirect carbon flux through the oxidative pentose pathway from the intermediate towards the purine-nucleotide pathway.     

Metabolic engineering of Corynebacterium glutamicum for production of sunscreen shinorine

Ultraviolet-absorbing chemicals are useful in cosmetics and skin care to prevent UV-induced skin damage. We demonstrate here that heterologous production of shinorine, which shows broad absorption maxima in the UV-A and UV-B region. A shinorine producing Corynebacterium glutamicum strain was constructed by expressing four genes from Actinosynnema mirum DSM 43827, which are responsible for the biosynthesis of shinorine from sedoheptulose-7-phosphate in the pentose phosphate pathway. Deletion of transaldolase encoding gene improved shinorine production by 5.2-fold. Among the other genes in pentose phosphate pathway, overexpression of 6-phosphogluconate dehydrogenase encoding gene further increased shinorine production by 60% (19.1 mg/L). The genetic engineering of the pentose phosphate pathway in C. glutamicum improved shinorine production by 8.3-fold in total, and could be applied to produce the other chemicals derived from sedoheptulose-7-phosphate.     

Dissecting the function of Atg1 complex in Dictyostelium autophagy reveals a connection with the pentose phosphate pathway enzyme transketolase

The network of protein–protein interactions of the Dictyostelium discoideum autophagy pathway was investigated by yeast two-hybrid screening of the conserved autophagic proteins Atg1 and Atg8. These analyses confirmed expected interactions described in other organisms and also identified novel interactors that highlight the complexity of autophagy regulation. The Atg1 kinase complex, an essential regulator of autophagy, was investigated in detail here. The composition of the Atg1 complex in D. discoideum is more similar to mammalian cells than to Saccharomyces cerevisiae as, besides Atg13, it contains Atg101, a protein not conserved in this yeast. We found that Atg101 interacts with Atg13 and genetic disruption of these proteins in Dictyostelium leads to an early block in autophagy, although the severity of the developmental phenotype and the degree of autophagic block is higher in Atg13-deficient cells. We have also identified a protein containing zinc-finger B-box and FNIP motifs that interacts with Atg101. Disruption of this protein increases autophagic flux, suggesting that it functions as a negative regulator of Atg101. We also describe the interaction of Atg1 kinase with the pentose phosphate pathway enzyme transketolase (TKT). We found changes in the activity of endogenous TKT activity in strains lacking or overexpressing Atg1, suggesting the presence of an unsuspected regulatory pathway between autophagy and the pentose phosphate pathway in Dictyostelium that seems to be conserved in mammalian cells.     

Compositional dynamics and codon usage pattern of BRCA1 gene across nine mammalian species

The BRCA1 gene is located on the human chromosome 17q21.31 and plays important role in biological processes. The aminoacyl-tRNA synthetases (AARS) are a family of heterogenous enzymes responsible protein synthesis and whose secondary functions include a role in autoimmune myositis. Our findings reveal that the compositional constraint and the preference of more A/T –ending codons determine the codon usage patterns in BRCA1 gene while more G/C-ending codons influence the codon usage pattern of AARS gene among mammals. The codon usage bias in BRCA1 and AARS genes is low. The codon CGC encoding arginine amino acid and the codon TTA encoding leucine were uniformly distributed in BRCA1 and AARS genes, respectively in mammals including human. Natural selection might have played a major role while mutation pressure might have played a minor role in shaping the codon usage pattern of BRCA1 and AARS genes.     

Synonymous codon usage in different protein secondary structural classes of human genes: Implication for increased non-randomness of GC<Subscript>3</Subscript> rich genes towards protein stability

The relationship between the synonymous codon usage and different protein secondary structural classes were investigated using 401 Homo sapiens proteins extracted from Protein Data Bank (PDB). A simple Chi-square test was used to assess the significance of deviation of the observed and expected frequencies of 59 codons at the level of individual synonymous families in the four different protein secondary structural classes. It was observed that synonymous codon families show non-randomness in codon usage in four different secondary structural classes. However, when the genes were classified according to their GC₃ levels there was an increase in non-randomness in high GC₃ group of genes. The non-randomness in codon usage was further tested among the same protein secondary structures belonging to four different protein folding classes of high GC₃ group of genes. The results show that in each of the protein secondary structural unit there exist some synonymous family that shows class specific codon-usage pattern. Moreover, there is an increased non-random behaviour of synonymous codons in sheet structure of all secondary structural classes in high GC₃ group of genes. Biological implications of these results have been discussed.     

High control coefficient of transketolase in the nonoxidative pentose phosphate pathway of human erythrocytes: NMR, antibody, and computer simulation studies.

The degree of control exerted by transketolase over metabolite flux in the nonoxidative pentose phosphate pathway in human erythrocytes was investigated using transketolase antiserum to modulate the activity of that enzyme. 31P NMR enabled the simultaneous measurement of the levels of pentose phosphate pathway metabolites following incubation of hemolysates with ribose 5-phosphate. The variations in metabolic flux which occurred as the transketolase activity of hemolysate samples was altered indicated that a high degree of control was exerted by transketolase. Investigations using transaldolase-depleted hemolysates showed that transaldolase exhibits a lesser degree of control over pathway flux. Experimental data were compared with simulations generated by a computer model encompassing the reactions of the classical nonoxidative pentose phosphate pathway. The sensitivity coefficients (also called "control strengths" or "flux-control coefficients") calculated from the computer simulations were 0.74 and 0.03 for transketolase and transaldolase, respectively.     

Biochemistry and physiology of nitrogen fixation with particular emphasis on nitrogen-fixing phototrophs

Summary This paper presents an overview of aspects of N₂-fixation in phototrophic N₂-fixers. Nitrogenase is little different in phototrophs from other organisms. Evidence suggests that fixed carbon dissimilation rather than direct photoreduction from oxidised inorganic compounds or exogenous photosynthetic electron donors is the major route of reductant supply to nitrogenase in phototrophs; inRhodospirillum rubrum pyruvate is a possible electron donor to nitrogenase; in cyanobacteria the oxidative pentose phosphate pathway is important, although some recent evidence implicates glycolysis and the tricarboxylic acid cycle in reductant supply in heterocystous cyanobacteria. In photosynthetic organisms light modulation of various enzymes occurs-some Calvin cycle enzymes are light activated, some oxidative pentose phosphate pathway and glycolytic enzymes are deactivated and some tricarboxylic acid cycle enzymes are activated. Reduced levels of thioredoxin in heterocysts may contribute to the sustained functioning of the oxidative pentose phosphate pathway in heterocysts in the light and dark. In photosynthetic bacteria such asRhodospirillum rubrum an activating enzyme which removes a modifying group from inactive Fe protein can activate nitrogenase. O₂ and NH ₄ ⁺ both inhibit N₂-fixation and there is some evidence in cyanobacteria that O₂ stability of whole cell nitrogenase can be achieved by prolonged incubation of cultures at high O₂.     

The Distribution of Synonymous Codon Choice in the Translation Initiation Region of Dengue Virus

Dengue is the most common arthropod-borne viral (Arboviral) illness in humans. The genetic features concerning the codon usage of dengue virus (DENV) were analyzed by the relative synonymous codon usage, the effective number of codons and the codon adaptation index. The evolutionary distance between DENV and the natural hosts (Homo sapiens, Pan troglodytes, Aedes albopictus and Aedes aegypti) was estimated by a novel formula. Finally, the synonymous codon usage preference for the translation initiation region of this virus was also analyzed. The result indicates that the general trend of the 59 synonymous codon usage of the four genotypes of DENV are similar to each other, and this pattern has no link with the geographic distribution of the virus. The effect of codon usage pattern of Aedes albopictus and Aedes aegypti on the formation of codon usage of DENV is stronger than that of the two primates. Turning to the codon usage preference of the translation initiation region of this virus, some codons pairing to low tRNA copy numbers in the two primates have a stronger tendency to exist in the translation initiation region than those in the open reading frame of DENV. Although DENV, like other RNA viruses, has a high mutation to adapt its hosts, the regulatory features about the synonymous codon usage have been ‘branded’ on the translation initiation region of this virus in order to hijack the translational mechanisms of the hosts.     

Biosynthesis of ribose-5-phosphate and erythrose-4-phosphate in archaea: a phylogenetic analysis of archaeal genomes

A phylogenetic analysis of the genes encoding enzymes in the pentose phosphate pathway (PPP), the ribulose monophosphate (RuMP) pathway, and the chorismate pathway of aromatic amino acid biosynthesis, employing data from 13 complete archaeal genomes, provides a potential explanation for the enigmatic phylogenetic patterns of the PPP genes in archaea. Genomic and biochemical evidence suggests that three archaeal species (Methanocaldococcus jannaschii, Thermoplasma acidophilum and Thermoplasma volcanium) produce ribose-5-phosphate via the nonoxidative PPP (NOPPP), whereas nine species apparently lack an NOPPP but may employ a reverse RuMP pathway for pentose synthesis. One species (Halobacterium sp. NRC-1) lacks both the NOPPP and the RuMP pathway but may possess a modified oxidative PPP (OPPP), the details of which are not yet known. The presence of transketolase in several archaeal species that are missing the other two NOPPP genes can be explained by the existence of differing requirements for erythrose-4-phosphate (E4P) among archaea: six species use transketolase to make E4P as a precursor to aromatic amino acids, six species apparently have an alternate biosynthetic pathway and may not require the ability to make E4P, and one species (Pyrococcus horikoshii) probably does not synthesize aromatic amino acids at all.     

转酮酶在癌症中的作用及其机制研究进展

肿瘤细胞的一大重要特征是代谢水平的改变。戊糖磷酸途径作为细胞产生NADPH和五碳糖的主要通路,在肿瘤发生发展过程中发挥着重要功能。转酮酶是戊糖磷酸途径中的关键酶之一,越来越多的研究表明其与癌症病人预后具有显著相关性。作为肿瘤诊断、治疗的潜在靶标,转酮酶具有重要的研究价值。我们就目前癌症研究中对转酮酶的研究进展做简要综述。