Comparative genomics of the vitamin B12 metabolism and regulation in prokaryotes

Using comparative analysis of genes, operons, and regulatory elements, we describe the cobalamin (vitamin B12) biosynthetic pathway in available prokaryotic genomes. Here we found a highly conserved RNA secondary structure, the regulatory B12 element, which is widely distributed in the upstream regions of cobalamin biosynthetic/transport genes in eubacteria. In addition, the binding signal (CBL-box) for a hypothetical B12 regulator was identified in some archaea. A search for B12 elements and CBL-boxes and positional analysis identified a large number of new candidate B12-regulated genes in various prokaryotes. Among newly assigned functions associated with the cobalamin biosynthesis, there are several new types of cobalt transporters, ChlI and ChlD subunits of the CobN-dependent cobaltochelatase complex, cobalt reductase BluB, adenosyltransferase PduO, several new proteins linked to the lower ligand assembly pathway, l-threonine kinase PduX, and a large number of other hypothetical proteins. Most missing genes detected within the cobalamin biosynthetic pathways of various bacteria were identified as nonorthologous substitutes. The variable parts of the cobalamin metabolism appear to be the cobalt transport and insertion, the CobG/CbiG- and CobF/CbiD-catalyzed reactions, and the lower ligand synthesis pathway. The most interesting result of analysis of B12 elements is that B12-independent isozymes of the methionine synthase and ribonucleotide reductase are regulated by B12 elements in bacteria that have both B12-dependent and B12-independent isozymes. Moreover, B12 regulons of various bacteria are thought to include enzymes from known B12-dependent or alternative pathways.     

Enzymatic and structural characterization of an archaeal thiamin phosphate synthase

Studies on thiamin biosynthesis have so far been achieved in eubacteria, yeast and plants, in which the thiamin structure is formed as thiamin phosphate from a thiazole and a pyrimidine moiety. This condensation reaction is catalyzed by thiamin phosphate synthase, which is encoded by the thiE gene or its orthologs. On the other hand, most archaea do not seem to have the thiE gene, but instead their thiD gene, coding for a 2-methyl-4-amino-5-hydroxymethylpyrimidine (HMP) kinase/HMP phosphate kinase, possesses an additional C-terminal domain designated thiN. These two proteins, ThiE and ThiN, do not share sequence similarity. In this study, using recombinant protein from the hyperthermophile archaea Pyrobaculum calidifontis, we demonstrated that the ThiN protein is an analog of the ThiE protein, catalyzing the formation of thiamin phosphate with the release of inorganic pyrophosphate from HMP pyrophosphate and 4-methyl-5-β-hydroxyethylthiazole phosphate (HET-P). In addition, we found that the ThiN protein can liberate an inorganic pyrophosphate from HMP pyrophosphate in the absence of HET-P. A structure model of the enzyme–product complex of P. calidifontis ThiN domain was proposed on the basis of the known three-dimensional structure of the ortholog of Pyrococcus furiosus. The significance of Arg320 and His341 residues for thiN-coded thiamin phosphate synthase activity was confirmed by site-directed mutagenesis. This is the first report of the experimental analysis of an archaeal thiamin synthesis enzyme.     

Metabolic engineering provides insight into the regulation of thiamin biosynthesis in plants

Thiamin (or thiamine) is a water-soluble B-vitamin (B1), which is required, in the form of thiamin pyrophosphate, as an essential cofactor in crucial carbon metabolism reactions in all forms of life. To ensure adequate metabolic functioning, humans rely on a sufficient dietary supply of thiamin. Increasing thiamin levels in plants via metabolic engineering is a powerful strategy to alleviate vitamin B1 malnutrition and thus improve global human health. These engineering strategies rely on comprehensive knowledge of plant thiamin metabolism and its regulation. Here, multiple metabolic engineering strategies were examined in the model plant Arabidopsis thaliana. This was achieved by constitutive overexpression of the three biosynthesis genes responsible for B1 synthesis, HMP-P synthase (THIC), HET-P synthase (THI1), and HMP-P kinase/TMP pyrophosphorylase (TH1), either separate or in combination. By monitoring the levels of thiamin, its phosphorylated entities, and its biosynthetic intermediates, we gained insight into the effect of either strategy on thiamin biosynthesis. Moreover, expression analysis of thiamin biosynthesis genes showed the plant’s intriguing ability to respond to alterations in the pathway. Overall, we revealed the necessity to balance the pyrimidine and thiazole branches of thiamin biosynthesis and assessed its biosynthetic intermediates. Furthermore, the accumulation of nonphosphorylated intermediates demonstrated the inefficiency of endogenous thiamin salvage mechanisms. These results serve as guidelines in the development of novel thiamin metabolic engineering strategies.     

Involvement of thiaminase II encoded by the THI20 gene in thiamin salvage of Saccharomyces cerevisiae

The physiological significance of thiaminase II, which catalyzes the hydrolysis of thiamin, has remained elusive for several decades. The C-terminal domains of THI20 family proteins (THI20/21/22) and the whole region of PET18 gene product of Saccharomyces cerevisiae are homologous to bacterial thiaminase II. On the other hand, the N-terminal domains of THI20 and THI21 encode 2-methyl-4-amino-5-hydroxymethylpyrimidine kinase and 2-methyl-4-amino-5-hydroxymethylpyrimidine phosphate kinase involved in the thiamin synthetic pathway. In this study, it was first indicated that the C-terminal domains of the THI20 family and PET18 are not required for de novo thiamin synthesis in S. cerevisiae, using a quadruple deletion strain expressing the N-terminal domain of THI20. Biochemical analysis using cell-free extracts and recombinant proteins demonstrated that yeast thiaminase II activity is exclusively encoded by THI20. It appeared that Thi20p has an affinity for the pyrimidine moiety of thiamin, and HMP produced by the thiaminase II activity is immediately phosphorylated. Thi20p was found to participate in the formation of thiamin from two synthetic antagonists, pyrithiamin and oxythiamin, by hydrolyzing both antagonists and phosphorylating HMP to give HMP pyrophosphate. Furthermore, 2-methyl-4-amino-5-aminomethylpyrimidine, a presumed naturally occurring thiamin precursor, was effectively converted to HMP by incubation with Thi20p. It is proposed that the thiaminase II activity of Thi20p is involved in the thiamin salvage pathway by catalyzing the hydrolysis of HMP precursors in S. cerevisiae.     

A new thiamin salvage pathway

The physiological function for thiaminase II, a thiamin-degrading enzyme, has eluded investigators for more than 50 years. Here, we demonstrate that this enzyme is involved in the regeneration of the thiamin pyrimidine rather than in thiamin degradation, and we identify a new pathway involved in the salvage of base-degraded forms of thiamin. This pathway is widely distributed among bacteria, archaea and eukaryotes. In this pathway, thiamin hydrolysis products such as N-formyl-4-amino-5-aminomethyl-2-methylpyrimidine (formylaminopyrimidine; 15) are transported into the cell using the ThiXYZ transport system, deformylated by the ylmB-encoded amidohydrolase and hydrolyzed to 4-amino-5-hydroxymethyl-2-methylpyrimidine (HMP; 6)-an intermediate on the de novo thiamin biosynthetic pathway. To our knowledge this is the first example of a thiamin salvage pathway involving thiamin analogs generated by degradation of one of the heterocyclic rings of the cofactor.     

Discovery of a SAR11 growth requirement for thiamin's pyrimidine precursor and its distribution in the Sargasso Sea

Vitamin traffic, the production of organic growth factors by some microbial community members and their use by other taxa, is being scrutinized as a potential explanation for the variation and highly connected behavior observed in ocean plankton by community network analysis. Thiamin (vitamin B₁), a cofactor in many essential biochemical reactions that modify carbon–carbon bonds of organic compounds, is distributed in complex patterns at subpicomolar concentrations in the marine surface layer (0–300 m). Sequenced genomes from organisms belonging to the abundant and ubiquitous SAR11 clade of marine chemoheterotrophic bacteria contain genes coding for a complete thiamin biosynthetic pathway, except for thiC, encoding the 4-amino-5-hydroxymethyl-2-methylpyrimidine (HMP) synthase, which is required for de novo synthesis of thiamin''s pyrimidine moiety. Here we demonstrate that the SAR11 isolate ‘Candidatus Pelagibacter ubique'', strain HTCC1062, is auxotrophic for the thiamin precursor HMP, and cannot use exogenous thiamin for growth. In culture, strain HTCC1062 required 0.7 zeptomoles per cell (ca. 400 HMP molecules per cell). Measurements of dissolved HMP in the Sargasso Sea surface layer showed that HMP ranged from undetectable (detection limit: 2.4 pM) to 35.7 pM, with maximum concentrations coincident with the deep chlorophyll maximum. In culture, some marine cyanobacteria, microalgae and bacteria exuded HMP, and in the Western Sargasso Sea, HMP profiles changed between the morning and evening, suggesting a dynamic biological flux from producers to consumers.     

Molecular characterization of the thi3 gene involved in thiamine biosynthesis in Zea mays: cDNA sequence and enzymatic and structural properties of the recombinant bifunctional protein with 4-amino-5-hydroxymethyl-2-methylpyrimidine (phosphate) kinase and thiamine monophosphate synthase activities

A thiamine biosynthesis gene, thi3, from maize Zea mays has been identified through cloning and sequencing of cDNA and heterologous overexpression of the encoded protein, THI3, in Escherichia coli. The recombinant THI3 protein was purified to homogeneity and shown to possess two essentially different enzymatic activities of HMP(-P) [4-amino-5-hydroxymethyl-2-methylpyrimidine (phosphate)] kinase and TMP (thiamine monophosphate) synthase. Both activities were characterized in terms of basic kinetic constants, with interesting findings that TMP synthase is uncompetitively inhibited by excess of one of the substrates [HMP-PP (HMP diphosphate)] and ATP. A bioinformatic analysis of the THI3 sequence suggested that these activities were located in two distinct, N-terminal kinase and C-terminal synthase, domains. Models of the overall folds of THI3 domains and the arrangements of active centre residues were obtained with the SWISS-MODEL protein modelling server, on the basis of the known three-dimensional structures of Salmonella enterica serotype Typhimurium HMP(-P) kinase and Bacillus subtilis TMP synthase. The essential roles of Gln98 and Met134 residues for HMP kinase activity and of Ser444 for TMP synthase activity were experimentally confirmed by site-directed mutagenesis.     

Carotenoid biosynthetic pathway: molecular phylogenies and evolutionary behavior of crt genes in eubacteria

Phylogenetic analysis of carotenoid biosynthetic pathway genes and their evolutionary rate variations were studied among eubacterial taxa. The gene sequences for the enzymes involved in this pathway were obtained for major phylogenetic groups of eubacteria (green sulfur bacteria, green nonsulphur bacteria, Gram-positive bacteria, proteobacteria, flavobacteria, cyanobacteria) and archeabacteria. These gene datasets were distributed under five major steps of carotenoid biosynthesis in eubacteria; isoprenoid precursor biosynthesis, phytoene synthesis, dehydrogenation of phytoene, lycopene cyclization, formation of acyclic xanthophylls, formation of cyclic xanthophylls and carotenoid biosynthesis regulation. The NJ algorithm was used on protein coding DNA sequences to deduce the evolutionary relationship for the respective crt genes among different eubacterial lineages. The rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)) and synonymous nucleotide substitutions per synonymous site (d(S)) were calculated for different clades of the respective phylogenetic tree for specific crt genes. The phylogenetic analysis suggests that evolutionary pattern of crt genes in eubacteria is characterized by lateral gene transfer and gene duplication events. The d(N) values indicate that carotenoid biosynthetic genes are more conserved in proteobacteria than in any other eubacterial phyla. Furthermore, of the genes involved in carotenoid biosynthesis pathway, structural genes evolve slowly than the regulatory genes in eubacteria.     

The genome sequence of Methanosphaera stadtmanae reveals why this human intestinal archaeon is restricted to methanol and H2 for methane formation and ATP synthesis

Methanosphaera stadtmanae has the most restricted energy metabolism of all methanogenic archaea. This human intestinal inhabitant can generate methane only by reduction of methanol with H2 and is dependent on acetate as a carbon source. We report here the genome sequence of M. stadtmanae, which was found to be composed of 1,767,403 bp with an average G+C content of 28% and to harbor only 1,534 protein-encoding sequences (CDS). The genome lacks 37 CDS present in the genomes of all other methanogens. Among these are the CDS for synthesis of molybdopterin and for synthesis of the carbon monoxide dehydrogenase/acetyl-coenzyme A synthase complex, which explains why M. stadtmanae cannot reduce CO2 to methane or oxidize methanol to CO2 and why this archaeon is dependent on acetate for biosynthesis of cell components. Four sets of mtaABC genes coding for methanol:coenzyme M methyltransferases were found in the genome of M. stadtmanae. These genes exhibit homology to mta genes previously identified in Methanosarcina species. The M. stadtmanae genome also contains at least 323 CDS not present in the genomes of all other archaea. Seventy-three of these CDS exhibit high levels of homology to CDS in genomes of bacteria and eukaryotes. These 73 CDS include 12 CDS which are unusually long (>2,400 bp) with conspicuous repetitive sequence elements, 13 CDS which exhibit sequence similarity on the protein level to CDS encoding enzymes involved in the biosynthesis of cell surface antigens in bacteria, and 5 CDS which exhibit sequence similarity to the subunits of bacterial type I and III restriction-modification systems.     

Common ancestry and novel genetic traits of Francisella novicida-like isolates from North America and Australia as revealed by comparative genomic analyses

Francisella novicida is a close relative of Francisella tularensis, the causative agent of tularemia. The genomes of F. novicida-like clinical isolates 3523 (Australian strain) and Fx1 (Texas strain) were sequenced and compared to F. novicida strain U112 and F. tularensis strain Schu S4. The strain 3523 chromosome is 1,945,310 bp and contains 1,854 protein-coding genes. The strain Fx1 chromosome is 1,913,619 bp and contains 1,819 protein-coding genes. NUCmer analyses revealed that the genomes of strains Fx1 and U112 are mostly colinear, whereas the genome of strain 3523 has gaps, translocations, and/or inversions compared to genomes of strains Fx1 and U112. Using the genome sequence data and comparative analyses with other members of the genus Francisella, several strain-specific genes that encode putative proteins involved in RTX toxin production, polysaccharide biosynthesis/modification, thiamine biosynthesis, glucuronate utilization, and polyamine biosynthesis were identified. The RTX toxin synthesis and secretion operon of strain 3523 contains four open reading frames (ORFs) and was named rtxCABD. Based on the alignment of conserved sequences upstream of operons involved in thiamine biosynthesis from various bacteria, a putative THI box was identified in strain 3523. The glucuronate catabolism loci of strains 3523 and Fx1 contain a cluster of nine ORFs oriented in the same direction that appear to constitute an operon. Strains U112 and Schu S4 appeared to have lost the loci for RTX toxin production, thiamine biosynthesis, and glucuronate utilization as a consequence of host adaptation and reductive evolution. In conclusion, comparative analyses provided insights into the common ancestry and novel genetic traits of these strains.