Evolution of vitamin B2 biosynthesis. A novel class of riboflavin synthase in Archaea

The open reading frame MJ1184 of Methanococcus jannaschii with similarity to riboflavin synthase of Methanothermobacter thermoautotrophicus was cloned into an expression vector but was poorly expressed in an Escherichia coli host strain. However, a synthetic open reading frame that was optimized for expression in E.coli directed the synthesis of abundant amounts of a protein with an apparent subunit mass of 17.5 kDa. The protein was purified to apparent homogeneity. Hydrodynamic studies indicated a relative mass of 88 kDa suggesting a homopentamer structure. The enzyme was shown to catalyze the formation of riboflavin from 6,7-dimethyl-8-ribityllumazine at a rate of 24 nmol mg(-1) min(-1) at 40 degrees C. Divalent metal ions, preferably manganese or magnesium, are required for maximum activity. In contrast to pentameric archaeal type riboflavin synthases, orthologs from plants, fungi and eubacteria are trimeric proteins characterized by an internal sequence repeat with similar folding patterns. In these organisms the reaction is achieved by binding the two substrate molecules in an antiparallel orientation. With the enzyme of M.jannaschii, 13C NMR spectroscopy with 13C-labeled 6,7-dimethyl-8-ribityllumazine samples as substrates showed that the regiochemistry of the dismutation reaction is the same as observed in eubacteria and eukaryotes, however, in a non-pseudo-c2 symmetric environment. Whereas the riboflavin synthases of M.jannaschii and M.thermoautotrophicus are devoid of similarity with those of eubacteria and eukaryotes, they have significant sequence similarity with 6,7-dimethyl-8-ribityllumazine synthases catalyzing the penultimate step of riboflavin biosynthesis. 6,7-Dimethyl-8-ribityllumazine synthase and the archaeal riboflavin synthase appear to have diverged early in the evolution of Archaea from a common ancestor. Some Archaea have eubacterial type riboflavin synthases which may have been acquired by lateral gene transfer.     

CDP-2,3-Di-O-geranylgeranyl-sn-glycerol:L-serine O-archaetidyltransferase (archaetidylserine synthase) in the methanogenic archaeon Methanothermobacter thermautotrophicus

CDP-2,3-di-O-geranylgeranyl-sn-glycerol:L-serine O-archaetidyltransferase (archaetidylserine synthase) activity in cell extracts of Methanothermobacter thermautotrophicus cells was characterized. The enzyme catalyzed the formation of unsaturated archaetidylserine from CDP-unsaturated archaeol and L-serine. The identity of the reaction products was confirmed by thin-layer chromatography, fast atom bombardment-mass spectrum analysis, and chemical degradation. The enzyme showed maximal activity in the presence of 10 mM Mn2+ and 1% Triton X-100. Among various synthetic substrate analogs, both enantiomers of CDP-unsaturated archaeols with ether-linked geranylgeranyl chains and CDP-saturated archaeol with ether-linked phytanyl chains were similarly active toward the archaetidylserine synthase. The activity on the ester analog of the substrate was two to three times higher than that on the corresponding ether-type substrate. The activity of D-serine with the enzyme was 30% of that observed for L-serine. A trace amount of an acid-labile, unsaturated archaetidylserine intermediate was detected in the cells by a pulse-labeling experiment. A gene (MT1027) in M. thermautotrophicus genome annotated as the gene encoding phosphatidylserine synthase was found to be homologous to Bacillus subtilis pssA but not to Escherichia coli pssA. The substrate specificity of phosphatidylserine synthase from B. subtilis was quite similar to that observed for the M. thermautotrophicus archaetidylserine synthase, while the E. coli enzyme had a strong preference for CDP-1,2-diacyl-sn-glycerol. It was concluded that M. thermautotrophicus archaetidylserine synthase belongs to subclass II phosphatidylserine synthase (B. subtilis type) on the basis of not only homology but also substrate specificity and some enzymatic properties. The possibility that a gene encoding the subclass II phosphatidylserine synthase might be transferred from a bacterium to an ancestor of methanogens is discussed.     

A Novel Biosynthetic Pathway of Archaetidyl-myo-inositol via Archaetidyl-myo-inositol Phosphate from CDP-archaeol and d-Glucose 6-Phosphate in Methanoarchaeon Methanothermobacter thermautotrophicus Cells

Ether-type inositol phospholipids are ubiquitously distributed in Archaea membranes. The present paper describes a novel biosynthetic pathway of the archaeal inositol phospholipid. To study the biosynthesis of archaetidylinositol in vitro, we prepared two possible substrates: CDP-archaeol, which was chemically synthesized, and myo-[¹⁴C]inositol 1-phosphate, which was enzymatically prepared from [¹⁴C]glucose 6-phosphate with the inositol 1-phosphate (IP) synthase of this organism. The complete structure of the IP synthase reaction product was determined to be 1l-myo-inositol 1-phosphate, based on gas liquid chromatography with a chiral column. When the two substrates were incubated with the Methanothermobacter thermautotrophicus membrane fraction, archaetidylinositol phosphate (AIP) was formed along with a small amount of archaetidylinositol (AI). The two products were identified by fast atom bombardment-mass spectrometry and chemical analyses. AI was formed from AIP by incubation with the membrane fraction, but AIP was not formed from AI. This finding indicates that archaeal AI was synthesized from CDP-archaeol and d-glucose 6-phosphate via myo-inositol 1-phosphate and AIP. Although the relevant enzymes were not isolated, three enzymes are implied: IP synthase, AIP synthase, and AIP phosphatase. AIP synthase was homologous to yeast phosphatidylinositol synthase, and we confirmed AIP synthase activity by cloning the encoding gene (MTH1691) and expressing it in Escherichia coli. AIP synthase is a newly found member of the enzyme superfamily CDP-alcohol phosphatidyltransferase, which includes a wide range of enzymes that attach polar head groups to ester- and ether-type phospholipids of bacterial and archaeal origin. This is the first report of the biosynthesis of ether-type inositol phospholipids in Archaea.     

Ubiquitous distribution of phosphatidylinositol phosphate synthase and archaetidylinositol phosphate synthase in Bacteria and Archaea,which contain inositol phospholipid

In Eukarya, phosphatidylinositol (PI) is biosynthesized from CDP-diacylglycerol (CDP-DAG) and inositol. In Archaea and Bacteria, on the other hand, we found a novel inositol phospholipid biosynthetic pathway. The precursors, inositol 1-phosphate, CDP-archaeol (CDP-ArOH), and CDP-DAG, form archaetidylinositol phosphate (AIP) and phosphatidylinositol phosphate (PIP) as intermediates. These intermediates are dephosphorylated to synthesize archaetidylinositol (AI) and PI. To date, the activities of the key enzymes (AIP synthase, PIP synthase) have been confirmed in only three genera (two archaeal genera, Methanothermobacter and Pyrococcus, and one bacterial genus, Mycobacterium). In the present study, we demonstrated that this novel biosynthetic pathway is universal in both Archaea and Bacteria, which contain inositol phospholipid, and elucidate the specificity of PIP synthase and AIP synthase for lipid substrates. PIP and AIP synthase activity were confirmed in all recombinant cells transformed with the respective gene constructs for four bacterial species (Streptomyces avermitilis, Propionibacterium acnes, Corynebacterium glutamicum, and Rhodococcus equi) and two archaeal species (Aeropyrum pernix and Sulfolobus solfataricus). Inositol was not incorporated. CDP-ArOH was used as the substrate for PIP synthase in Bacteria, and CDP-DAG was used as the substrate for AIP synthase in Archaea, despite their fundamentally different structures. PI synthase activity was observed in two eukaryotic species, Saccharomyces cerevisiae and Homo sapiens; however, inositol 1-phosphate was not incorporated. In Eukarya, the only pathway converts free inositol and CDP-DAG directly into PI. Phylogenic analysis of PIP synthase, AIP synthase, and PI synthase revealed that they are closely related enzymes.     

Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.     

Novel archaeal phylotypes from an East African alkaline saltern

DNA has been extracted on site from the brines of the final crystallizing pond of an alkaline saltern at Lake Magadi, Kenya. Amplification of 16S rRNA genes followed by cloning, sequencing, and phylogenetic analysis has revealed the presence of two distinct new archaeal lineages. The majority of cloned sequences showed greater than 95% identity to each other, but only 88%–90% similarity to any cultivated haloalkaliphilic Archaea, and form a distinct cluster within the known Haloarchaea. Two cloned genes showed close similarity to each other but only 76% similarity to any known archaeal sequence, and therefore represent a distinct phylotype only distantly related to the euryarchaeotal branch of the Archaea. Received: August 10, 1998 / Accepted: December 10, 1998     

Ribosomal protein-sequence block structure suggests complex prokaryotic evolution with implications for the origin of eukaryotes

Amino acid sequence alignments of orthologous ribosomal proteins found in Bacteria, Archaea, and Eukaryota display, relative to one another, an unusual segment or block structure, with major evolutionary implications. Within each of the prokaryotic phylodomains the sequences exhibit substantial similarity, but cross-domain alignments break up into (a) universal blocks (conserved in both phylodomains), (b) bacterial blocks (unalignable with any archaeal counterparts), and (c) archaeal blocks (unalignable with any bacterial counterparts). Sequences of those eukaryotic cytoplasmic riboproteins that have orthologs in both Bacteria and Archaea, exclusively match the archaeal block structure. The distinct blocks do not correlate consistently with any identifiable functional or structural feature including RNA and protein contacts. This phylodomain-specific block pattern also exists in a number of other proteins associated with protein synthesis, but not among enzymes of intermediary metabolism. While the universal blocks imply that modern Bacteria and Archaea (as defined by their translational machinery) clearly have had a common ancestor, the phylodomain-specific blocks imply that these two groups derive from single, phylodomain-specific types that came into existence at some point long after that common ancestor. The simplest explanation for this pattern would be a major evolutionary bottleneck, or other scenario that drastically limited the progenitors of modern prokaryotic diversity at a time considerably after the evolution of a fully functional translation apparatus. The vast range of habitats and metabolisms that prokaryotes occupy today would thus reflect divergent evolution after such a restricting event. Interestingly, phylogenetic analysis places the origin of eukaryotes at about the same time and shows a closer relationship of the eukaryotic ribosome-associated proteins to crenarchaeal rather than euryarchaeal counterparts.     

Archaeal phylogeny based on ribosomal proteins

Until recently, phylogenetic analyses of Archaea have mainly been based on ribosomal RNA (rRNA) sequence comparisons, leading to the distinction of the two major archaeal phyla: the Euryarchaeota and the Crenarchaeota. Here, thanks to the recent sequencing of several archaeal genomes, we have constructed a phylogeny based on the fusion of the sequences of the 53 ribosomal proteins present in most of the archaeal species. This phylogeny was remarkably congruent with the rRNA phylogeny, suggesting that both reflected the actual phylogeny of the domain Archaea even if some nodes remained unresolved. In both cases, the branches leading to hyperthermophilic species were short, suggesting that the evolutionary rate of their genes has been slowed down by structural constraints related to environmental adaptation. In addition, to estimate the impact of lateral gene transfer (LGT) on our tree reconstruction, we used a new method that revealed that 8 genes out of the 53 ribosomal proteins used in our study were likely affected by LGT. This strongly suggested that a core of 45 nontransferred ribosomal protein genes existed in Archaea that can be tentatively used to infer the phylogeny of this domain. Interestingly, the tree obtained using only the eight ribosomal proteins likely affected by LGT was not very different from the consensus tree, indicating that LGT mainly brought random phylogenetic noise. The major difference involves organisms living in similar environments, suggesting that LGTs are mainly directed by the physical proximity of the organisms rather than by their phylogenetic proximity.     

The archaeal homolog of the Imp4 protein, a eukaryotic U3 snoRNP component

Homologs of the Imp4 protein, a component specific to the eukaryotic U3 snoRNP complex, have been found in all archaeal genomes. The archaeal and eukaryotic Imp4 proteins that are related to four other protein families, the Imp4-like, the SSF1 homologs and two sets of hypothetical proteins, are characterized by the Imp4 signature pattern. These findings, together with the presence of other snoRNPs homologs in Archaea, provide evidence for similar RNA processing and folding in Eukarya and Archaea.     

Phylogenetic diversity and ecology of environmental Archaea

On the basis of culture studies, Archaea were thought to be synonymous with extreme environments. However, the large numbers of environmental rRNA gene sequences currently flooding into databases such as GenBank show that these organisms are present in almost all environments examined to date. Large sequence databases and new fast phylogenetic software allow more precise determination of the archaeal phylogenetic tree, but also indicate that our knowledge of archaeal diversity is incomplete. Although it is apparent that Archaea can be found in all environments, the chemistry of their ecological context is mostly unknown.     

Bioenergetics of the Archaea.

In the late 1970s, on the basis of rRNA phylogeny, Archaea (archaebacteria) was identified as a distinct domain of life besides Bacteria (eubacteria) and Eucarya. Though forming a separate domain, Archaea display an enormous diversity of lifestyles and metabolic capabilities. Many archaeal species are adapted to extreme environments with respect to salinity, temperatures around the boiling point of water, and/or extremely alkaline or acidic pH. This has posed the challenge of studying the molecular and mechanistic bases on which these organisms can cope with such adverse conditions. This review considers our cumulative knowledge on archaeal mechanisms of primary energy conservation, in relationship to those of bacteria and eucarya. Although the universal principle of chemiosmotic energy conservation also holds for Archaea, distinct features have been discovered with respect to novel ion-transducing, membrane-residing protein complexes and the use of novel cofactors in bioenergetics of methanogenesis. From aerobically respiring Archaea, unusual electron-transporting supercomplexes could be isolated and functionally resolved, and a proposal on the organization of archaeal electron transport chains has been presented. The unique functions of archaeal rhodopsins as sensory systems and as proton or chloride pumps have been elucidated on the basis of recent structural information on the atomic scale. Whereas components of methanogenesis and of phototrophic energy transduction in halobacteria appear to be unique to Archaea, respiratory complexes and the ATP synthase exhibit some chimeric features with respect to their evolutionary origin. Nevertheless, archaeal ATP synthases are to be considered distinct members of this family of secondary energy transducers. A major challenge to future investigations is the development of archaeal genetic transformation systems, in order to gain access to the regulation of bioenergetic systems and to overproducers of archaeal membrane proteins as a prerequisite for their crystallization.     

Crystal structure of an archaeal pentameric riboflavin synthase in complex with a substrate analog inhibitor: stereochemical implications

Whereas eubacterial and eukaryotic riboflavin synthases form homotrimers, archaeal riboflavin synthases from Methanocaldococcus jannaschii and Methanothermobacter thermoautrophicus are homopentamers with sequence similarity to the 6,7-dimethyl-8-ribityllumazine synthase catalyzing the penultimate step in riboflavin biosynthesis. Recently it could be shown that the complex dismutation reaction catalyzed by the pentameric M. jannaschii riboflavin synthase generates riboflavin with the same regiochemistry as observed for trimeric riboflavin synthases. Here we present crystal structures of the pentameric riboflavin synthase from M. jannaschii and its complex with the substrate analog inhibitor, 6,7-dioxo-8-ribityllumazine. The complex structure shows five active sites located between adjacent monomers of the pentamer. Each active site can accommodate two substrate analog molecules in anti-parallel orientation. The topology of the two bound ligands at the active site is well in line with the known stereochemistry of a pentacyclic adduct of 6,7-dimethyl-8-ribityllumazine that has been shown to serve as a kinetically competent intermediate. The pentacyclic intermediates of trimeric and pentameric riboflavin synthases are diastereomers.     

Potential coexistence of both bacterial and eukaryotic small RNA biogenesis and functional related protein homologs in Archaea

RNA silencing plays crucial roles in both bacteria and eukaryotes, yet its machinery appears to differ in these two kingdoms. A couple of Argonaute protein homologs have been reported in some archaeal species in recent years. As Argonaute protein is the key component of eukaryotic RNA silencing pathways, such findings suggested the possibility of existence of eukaryotic RNA silencing like pathways in Archaea, which present the life forms between prokaryotes and eukaryotes. To further explore such hypothesis, we systematically screened 71 fully sequenced archaeal genomes, and identified some proteins containing homologous regions to the functional domains of eukaryotie RNA silencing pathway key proteins. The phylogenetic relationships of these proteins were analyzed. The conserved time-tional amino acids between archaeal and eukaryotic Piwi domains suggested their functional similarity. Our results provide new clues to the evolution of RNA silencing pathways.     

Origin of eukaryotic cell nuclei by symbiosis of Archaea in Bacteria is revealed by homology-hit analysis

The origin of eukaryotic cell nuclei by symbiosis of Archaea in Bacteria was proposed on the basis of the phylogenetic topologies of genes. However, it was not possible to conclude whether or not the genes involved were authentic representative genes. Furthermore, using the BLAST and FASTA programs, the similarity of open reading frame (ORF) groups between three domains (Eukarya, Archaea and Bacteria) was estimated at one threshold. Therefore, their similarities at other thresholds could not be clarified. Here we use our newly developed 'homology-hit analysis' method, which uses multiple thresholds, to determine the origin of the nucleus. We removed mitochondria-related ORFs from yeast ORFs, and determined the number of yeast orthologous ORFs in each functional category to the ORFs in six Archaea and nine Bacteria at several thresholds (E-values) using the BLAST. Our results indicate that yeast ORFs related to the nucleus may share their origins with archaeal ORFs, whereas ORFs that are related to the cytoplasm may share their origins with bacterial ORFs. Our results thus strongly support the idea of nucleus symbiosis.     

Structural analysis of the plasmid pTA1 isolated from the thermoacidophilic archaeon Thermoplasma acidophilum

Thermoplasma acidophilum is a thermoacidophilic archaeon that grows optimally at pH1.8 and 56°C and has no cell wall. Plasmid pTA1 was found in some strains of the species. We sequenced plasmid pTA1 and analyzed the open reading frames (ORFs). pTA1 was found to be a circular DNA molecule of 15,723 bp. Eighteen ORFs were found; none of the gene products except ORF1 had sequence similarity to known proteins. ORF1 showed similarity to Cdc6, which is involved in genome-replication initiation in Eukarya and Archaea. T. acidophilum has two Cdc6 homologues in the genome. The homologue found in pTA1 is most similar to Tvo3, one of the three Cdc6 homologues found in the genome of Thermoplasma volcanium, among all of the Cdc6 family proteins. The phylogenetic analysis suggested that plasmid pTA1 is possibly originated from the chromosomal DNA of Thermoplasma.     

In vitro biosynthesis of ether-type glycolipids in the methanoarchaeon Methanothermobacter thermautotrophicus

The biosynthesis of archaeal ether-type glycolipids was investigated in vitro using Methanothermobacter thermautotrophicus cell-free homogenates. The sole sugar moiety of glycolipids and phosphoglycolipids of the organism is the beta-D-glucosyl-(1-->6)-D-glucosyl (gentiobiosyl) unit. The enzyme activities of archaeol:UDP-glucose beta-glucosyltransferase (monoglucosylarchaeol [MGA] synthase) and MGA:UDP-glucose beta-1,6-glucosyltransferase (diglucosylarchaeol [DGA] synthase) were found in the methanoarchaeon. The synthesis of DGA is probably a two-step glucosylation: (i) archaeol + UDP-glucose --> MGA + UDP, and (ii) MGA + UDP-glucose --> DGA + UDP. Both enzymes required the addition of K(+) ions and archaetidylinositol for their activities. DGA synthase was stimulated by 10 mM MgCl(2), in contrast to MGA synthase, which did not require Mg(2+). It was likely that the activities of MGA synthesis and DGA synthesis were carried out by different proteins because of the Mg(2+) requirement and their cellular localization. MGA synthase and DGA synthase can be distinguished in cell extracts greatly enriched for each activity by demonstrating the differing Mg(2+) requirements of each enzyme. MGA synthase preferred a lipid substrate with the sn-2,3 stereostructure of the glycerol backbone on which two saturated isoprenoid chains are bound at the sn-2 and sn-3 positions. A lipid substrate with unsaturated isoprenoid chains or sn-1,2-dialkylglycerol configuration exhibited low activity. Tetraether-type caldarchaetidylinositol was also actively glucosylated by the homogenates to form monoglucosyl caldarchaetidylinositol and a small amount of diglucosyl caldarchaetidylinositol. The addition of Mg(2+) increased the formation of diglucosyl caldarchaetidylinositol. This suggested that the same enzyme set synthesized the sole sugar moiety of diether-type glycolipids and tetraether-type phosphoglycolipids.     

Phylogenetic tree of bacterial and eucaryotic FtsZ-proteins based on the homology of their primary sequences

In this report we present phylogenetic tree based on the sequences of FtsZ-proteins of Archaea, Bacteria and Eukaryota. We have analyzed the relative positions of separate sequences and their clusters and their distance to the hypothetical root of phylogenetic tree as well. It has been demonstrated that most of procaryotic and eucaryotic FtsZ-protein sequences are located at the distance of 0.3-0.8 substitutions per site from the hypothetical root of phylogenetic tree. The tubulins are outlined from the root at 7-8 fold larger distances in comparison with the separate sequences. This confirms the slight aminoacid homology between these two groups of cytoskeletal proteins. At the same time the significant similarity between FtsZ-proteins of cyanobacteria and of plant plastids has been shown.     

Novel coding regions in four complete archaeal genomes.

In the process of analysing the four available complete archaeal genomes, we have noted that certain regions characterised as 'non-coding' exhibit significant sequence similarity to other protein sequences from Archaea and other species. Using established technology, we have identified a number of potential protein coding regions in these putative 'non-coding' regions. We have detected 524 such cases, of which 113 regions appear to code for proteins present in archaeal or other species, while the remaining 411 regions are mostly start/stop definition conflicts. Of the 113 protein coding regions, only 21 code for proteins with homologues of known function. The number of novel coding sequences identified herein amounts to 1. 5% of the total genome entries, while the conflicting cases represent an additional 5%. The observed differences between the four complete archaeal genomes seem to reflect disparate approaches to genome annotation. Genome sequence collections should be regularly checked to improve gene prediction by sequence similarity and greater effort is required to make gene definitions consistent across related species.