Highly Divergent Genes for Methanopterin-Linked C1 Transfer Reactions in Lake Washington, Assessed via Metagenomic Analysis and mRNA Detection

The origins and the evolutionary history of tetrahydromethanopterin-linked C₁ transfer reactions that are part of two environmentally important biotransformations, methylotrophy and methanogenesis, are still not well understood. In previous studies, we have expanded the known phylogenetic diversity of these reactions by identifying genes highly diverging from the ones associated with cultivated Proteobacteria, Planctomycetes, or Archaea (M. G. Kalyuzhnaya, M. E. Lidstrom, and L. Chistoserdova, Microb. Ecol. 48:463-472, 2004; M. G. Kalyuzhnaya, O. Nercessian, M. E. Lidstrom, and L. Chistoserdova, Environ. Microbiol. 7:1269-1274, 2005). Here we used a metagenomic approach to demonstrate that these divergent genes are present with high abundance in the microbial community inhabiting Lake Washington sediment. We also gained preliminary insights into the genomic composition of the organisms possessing these genes by sequencing genomic fragments from three uncultured microbes possessing the genes of interest. Phylogenetic analyses suggested that, although distantly related to each other, these organisms deeply diverge from known Bacteria and Archaea, with more relation to the former, suggesting their affiliation with a new bacterial phylum. We also demonstrate, via specific mRNA detection, that these divergent genes are expressed in the environment, pointing toward their potential role in local carbon cycling.     

Poly-beta-hydroxybutyrate biosynthesis in the facultative methylotroph methylobacterium extorquens AM1: identification and mutation of gap11, gap20, and phaR

Methylobacterium extorquens AM1, a serine cycle facultative methylotroph, accumulates poly-beta-hydroxybutyrate (PHB) as a carbon and energy reserve material during growth on both multicarbon- and single-carbon substrates. Recently, the identification and mutation of the genes involved in the biosynthesis and degradation of PHB have been described for this bacterium, demonstrating that two of the genes of the PHB cycle (phaA and phaB) are also involved in C(1) and C(2) metabolism, as part of a novel pathway for glyoxylate regeneration in the serine cycle (N. Korotkova and M. E. Lidstrom, J. Bacteriol. 183:1038-1046, 2001; N. Korotkova, L. Chistoserdova, V. Kuksa, and M. E. Lidstrom, J. Bacteriol. 184:1750-1758, 2002). In this work, three new genes involved in PHB biosynthesis in this bacterium have been investigated via mutation and phenotypic analysis: gap11, gap20, and phaR. We demonstrate that gap11 and gap20 encode two major granule-associated proteins (phasins) and that mutants with mutations in these genes are defective in PHB production and also in growth on C(2) compounds, while they show wild-type growth characteristics on C(1) or multicarbon compounds. The phaR mutant shows defects in both PHB accumulation and growth characteristics when grown on C(1) compounds and has defects in PHB accumulation but grows normally on C(3) and C(4) compounds, while both PHB accumulation and growth rate are at wild-type levels during growth on C(2) compounds. Our results suggest that this phenotype is due to altered fluxes of acetyl coenzyme A (CoA), a major intermediate in C(1), C(2), and heterotrophic metabolism in M. extorquens AM1, as well as the entry metabolite for the PHB cycle. Therefore, it seems likely that PhaR acts to control acetyl-CoA flux to PHB in this methylotrophic bacterium.     

Genetics of the serine cycle in Methylobacterium extorquens AM1: identification, sequence, and mutation of three new genes involved in C1 assimilation, orf4, mtkA, and mtkB.

In a recent paper we reported the sequence of the beginning of a serine cycle gene cluster on the Methylobacterium extorquens AM1 chromosome, containing the genes encoding serine glyoxylate aminotransferase (sgaA), hydroxypyruvate reductase (hprA), and 5,10-methylenetetrahydrofolate dehydrogenase (mtdA) (L. V. Chistoserdova and M. E. Lidstrom J. Bacteriol. 176:1957-1968, 1994). Here we present the sequence of the adjacent downstream region containing three full and one partial open reading frames. The first of the full open reading frames (orf4) remains unidentified, while the other two (mtkA and mtkB) code for the two subunits of malate thiokinase, and the fourth, a partial open reading frame (ppcA), apparently encodes phosphoenolpyruvate carboxylase. Mutants containing insertion mutations in orf4, mtdA, and mtdB all were unable to grow on C1 compounds, showing that these three newly identified genes are indispensable for the operation of the serine cycle. Mutants in orf4 were also unable to grow on C2 compounds, but growth was restored by glyoxylate, suggesting that orf4 might be required for the conversion of acetyl coenzyme A to glyoxylate.     

Genetics of the serine cycle in Methylobacterium extorquens AM1: identification of sgaA and mtdA and sequences of sgaA, hprA, and mtdA.

In a previous paper, we reported identification of the 5' part of hprA of Methylobacterium extorquens AM1, which encodes the serine cycle enzyme hydroxypyruvate reductase (L. V. Chistoserdova and M. E. Lidstrom, J. Bacteriol. 174:71-77, 1992). Here we present the complete sequence of hprA and partial sequence of genes adjacent to hprA. Upstream of hprA, the 3' part of an open reading frame was discovered, separated from hprA by 263 bp. This open reading frame was identified as the gene encoding another serine cycle enzyme, serine glyoxylate aminotransferase (sgaA). Cells containing an insertion mutation into sgaA were unable to grow on C1 compounds, demonstrating that the gene is required for C1 metabolism. Sequencing downstream of hprA has revealed the presence of another open reading frame (mtdA), which is probably cotranscribed with hprA. This open reading frame was identified as the gene required for the synthesis of 5,10-methylenetetrahydrofolate dehydrogenase. Our data suggest that this enzyme plays an integral role in methylotrophic metabolism in M. extorquens AM1, either in formaldehyde oxidation or as part of the serine cycle.     

Liverwort genomes display extensive structural variations

PIKE, L. M., HU, A., RENZAGLIA, K. S. & MUSICH, P. R., 1992. Liverwort genomes display extensive structural variations. Analyses of the total genomic DNA of eight species of liverworts and two species of green algae by thermal denaturation and CsCl buoyant density gradient centrifugation reveal a high degree of structural complexity and interspecific heterogeneity. The hepatic taxa exhibit two or more DNA components of varying base composition. Average G4-C contents of total cellular DNA calculated from melting profiles are similarly variable, ranging from 38% to 53% G + C. The green alga Chara , a member of the ancestral line to land plants, shows similarities with liverworts in possessing multiple DNA components of comparable complexity, whereas Hydrodiciyon DNA displays a single component. Detailed hybridization analyses of individual density gradient fractions using α-tubulin, rRNA and ribulose 1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL) gene probes were performed to locate the low-copy number and moderately repetitive nuclear genes, and the chloroplast chromosome, respectively. The location of each gene within the density gradient is highly variable among the organisms examined; a-tubulin occurs in fractions ranging from 44–64% G + C, rDNA in 50–64% G + C fractions, and the RbcL gene is located in fractions from 30–59% G + C. For a given species, the two nuclear genes normally overlap in their distributions within the gradient. In most instances, neither gene occurs in the major DNA components, indicating that these components may contain repetitive DNAs. The observed variation in the density of the rbcL gene implies substantial reorganization of the chloroplast genome. The overall differences in the genomic components within and between taxa provide insight into the dynamics of DNA structure that have occurred during the extended evolutionary history of these organisms.     

MeaB is a component of the methylmalonyl-CoA mutase complex required for protection of the enzyme from inactivation

Adenosylcobalamin-dependent methylmalonyl-CoA mutase catalyzes the interconversion of methylmalonyl-CoA and succinyl-CoA. In humans, deficiencies in the mutase lead to methylmalonic aciduria, a rare disease that is fatal in the first year of life. Such inherited deficiencies can result from mutations in the mutase structural gene or from mutations that impair the acquisition of cobalamins. Recently, a human gene of unknown function, MMAA, has been implicated in methylmalonic aciduria (Dobson, C. M., Wai, T., Leclerc, D., Wilson, A., Wu, X., Dore, C., Hudson, T., Rosenblatt, D. S., and Gravel, R. A. (2002) Proc. Natl. Acad. Sci. U. S. A. 99, 15554-15559). MMAA orthologs are widespread in bacteria, archaea, and eukaryotes. In Methylobacterium extorquens AM1, a mutant defective in the MMAA homolog meaB was unable to grow on C(1) and C(2) compounds because of the inability to convert methylmalonyl-CoA to succinyl-CoA (Korotkova N., Chistoserdova, L., Kuksa, V., and Lidstrom, M. E. (2002) J. Bacteriol. 184, 1750-1758). Here we demonstrate that this defect is not due to the absence of adenosylcobalamin but due to an inactive form of methylmalonyl-CoA mutase. The presence of active mutase in double mutants defective in MeaB and in the synthesis of either R-methylmalonyl-CoA or adenosylcobalamin indicates that MeaB is necessary for protection of mutase from inactivation during catalysis. MeaB and methylmalonyl-CoA mutase from M. extorquens were cloned and purified in their active forms. We demonstrated that MeaB forms a complex with methylmalonyl-CoA mutase and stimulates in vitro mutase activity. These results support the hypothesis that MeaB functions to protect methylmalonyl-CoA mutase from irreversible inactivation.     

Identification of a fourth formate dehydrogenase in Methylobacterium extorquens AM1 and confirmation of the essential role of formate oxidation in methylotrophy

A mutant of Methylobacterium extorquens AM1 with lesions in genes for three formate dehydrogenase (FDH) enzymes was previously described by us (L. Chistoserdova, M. Laukel, J.-C. Portais, J. A. Vorholt, and M. E. Lidstrom, J. Bacteriol. 186:22-28, 2004). This mutant had lost its ability to grow on formate but still maintained the ability to grow on methanol. In this work, we further investigated the phenotype of this mutant. Nuclear magnetic resonance experiments with [¹³C]formate, as well as ¹⁴C-labeling experiments, demonstrated production of labeled CO₂ in the mutant, pointing to the presence of an additional enzyme or a pathway for formate oxidation. The tungsten-sensitive phenotype of the mutant suggested the involvement of a molybdenum-dependent enzyme. Whole-genome array experiments were conducted to test for genes overexpressed in the triple-FDH mutant compared to the wild type, and a gene (fdh4A) was identified whose translated product carried similarity to an uncharacterized putative molybdopterin-binding oxidoreductase-like protein sharing relatively low similarity with known formate dehydrogenase alpha subunits. Mutation of this gene in the triple-FDH mutant background resulted in a methanol-negative phenotype. When the gene was deleted in the wild-type background, the mutant revealed diminished growth on methanol with accumulation of high levels of formate in the medium, pointing to an important role of FDH4 in methanol metabolism. The identity of FDH4 as a novel FDH was also confirmed by labeling experiments that revealed strongly reduced CO₂ formation in growing cultures. Mutation of a small open reading frame (fdh4B) downstream of fdh4A resulted in mutant phenotypes similar to the phenotypes of fdh4A mutants, suggesting that fdh4B is also involved in formate oxidation.     

Preferred amino acids and thermostability

Most organisms grow at temperatures from 20 to 50 degrees C, but some prokaryotes, including Archaea and Bacteria, are capable of withstanding higher temperatures, from 60 to >100 degrees C. Their biomolecules, especially proteins, must be sufficiently stable to function under these extreme conditions; however, the basis for thermostability remains elusive. We investigated the preferential usage of certain groupings of amino acids and codons in thermally adapted organisms, by comparative proteome analysis, using 28 complete genomes from 18 mesophiles (M), 4 thermophiles (T), and 6 hyperthermophiles (HT). Whenever the percent of glutamate (E) and lysine (K) increased in the HT proteomes, the percent of glutamine (Q) and histidine (H) decreased, so that the E + K/Q + H ratio was >4.5; it was <2.5 in the M proteomes, and 3.2 to 4.6 in T. The E + K/Q + H ratios for chaperonins, potentially thermostable proteins, were higher than their proteome ratios, whereas for DNA ligases, which are not necessarily thermostable, they followed the proteome ratios. Analysis of codon usage revealed that HT had more AGR codons for Arg than they did CGN codons, which were more common in mesophiles. The E + K/Q + H ratio may provide a useful marker for distinguishing HT, T and M prokaryotes, and the high percentage of the amino acid couple E + K, consistently associated with a low percentage of the pair Q + H, could contribute to protein thermostability. The preponderance of AGR codons for Arg is a signature of all HT so far analyzed. The E + K/Q + H ratio and the codon bias for Arg are apparently not related to phylogeny. HT members of the Bacteria show the same values as the HT members of the Archaea; the values for T organisms are related to their lifestyle (intermediate temperature) and not to their domain (Archaea) and the values for M are similar in Eukarya, Bacteria and Archaea.     

Seven GC-rich microbial genomes adopt similar codon usage patterns regardless of their phylogenetic lineages

Seven GC-rich (group I) and three AT-rich (group II) microbial genomes are analyzed in this paper. The seven microbes in group I belong to different phylogenetic lineages, even different domains of life. The common feature is that they are highly GC-rich organisms, with more than 60% genomic GC content. Group II includes three bacteria, which belong to the same subdivision as Pseudomonas aeruginosa in group I. The genomic GC content of the three bacteria is in the range of 26-50%. It is shown that although the phylogenetic lineages of the organisms in group I are remote, the common feature of highly genomic GC content forces them to adopt similar codon usage patterns, which constitutes the basis of an algorithm using a set of universal parameters to recognize known genes in the seven genomes. The common codon usage pattern of function known genes in the seven genomes is GGS type, where G, G, and S are the bases of G, non-G, and G/C, respectively. On the contrary, although the phylogenetic lineages of the three bacteria in group II are quite close, the codon usage patterns of function known genes in these genomes are obviously distinct. There are no universal parameters to identify known genes in the three genomes in group II. It can be deduced that the genomic GC content is more important than phylogenetic lineage in gene recognition programs. We hope that the work might be useful for understanding the common characteristics in the organization of microbial genomes.     

Introns of Entamoeba histolytica and Entamoeba dispar.

The genome of Entamoeba histolytica is considered to possess very few intervening sequences (introns), as only 5 intron-containing genes from this protozoan parasite have been reported so far. However, while sequencing a number of genomic contigs as well as three independent genes coding for ribosomal protein L27a, we have identified 9 additional intron-containing genes of E. histolytica and the closely related species Entamoeba dispar, indicating that introns are more common in these organisms than previously suggested. The various amoeba introns are relatively short comprising between 46 and 115 nucleotides only and have a higher AT-content compared to the corresponding exon sequences. In contrast to higher eukaryotes, amoeba introns do not contain a well-conserved branch point consensus, and have extended donor and acceptor splice sites of the sequences G     

The structure and evolution of the ribosomal proteins encoded in the spc operon of the archaeon (Crenarchaeota) Sulfolobus acidocaldarius.

The genes for nine ribosomal proteins, L24, L5, S14, S8, L6, L18, S5, L30, and L15, have been isolated and sequenced from the spc operon in the archaeon (Crenarchaeota) Sulfolobus acidocaldarius, and the putative amino acid sequence of the proteins coded by these genes has been determined. In addition, three other genes in the spc operon, coding for ribosomal proteins S4E, L32E, and L19E (equivalent to rat ribosomal proteins S4, L32, and L19), were sequenced and the structure of the putative proteins was determined. The order of the ribosomal protein genes in the spc operon of the Crenarchaeota kingdom of Archaea is identical to that present in the Euryarchaeota kingdom of Archaea and also identical to that found in bacteria, except for the genes for r-proteins S4E, L32E, and L19E, which are absent in bacteria. Although AUG is the initiation codon in most of the spc genes, GUG (val) and UUG (leu) are also used as initiation codons in S. acidocaldarius. Over 70% of the codons in the Sulfolobus spc operon have A or U in the third position, reflecting the low GC content of Sulfolobus DNA. Phylogenetic analysis indicated that the archaeal r-proteins are a sister group of their eucaryotic counterparts but did not resolve the question of whether the Archaea is monophyletic, as suggested by the L6P, L15P, and L18P trees, or the question of whether the Crenarchaeota is separate from the Euryarchaeota and closer to the Eucarya, as suggested by the S8P, S5P, and L24P trees. In the case of the three Sulfolobus r-proteins that do not have a counterpart in the bacterial ribosome (S4E, L32E, and L19E), the archaeal r-proteins showed substantial identity to their eucaryotic equivalents, but in all cases the archaeal proteins formed a separate group from the eucaryotic proteins.     

A complex gene superfamily encodes actin in petunia.

We have shown by several independent criteria that actin is encoded by a very large and complex superfamily of genes in Petunia. Several cDNA and genomic probes encoding actins from diverse organisms (Dictyostelium, Drosophila, chicken and soybean) hybridize to hundreds of restriction fragments in the petunia genome. Actin-hybridizing sequences were isolated from a petunia genomic library at a rate of at least 200 per genome equivalent. Twenty randomly selected actin-hybridizing clones were characterized in more detail. DNA sequence data from four representative and highly divergent clones, PAc2, PAc3, PAc4 and PAc7, demonstrate that these actin-like sequences are related to functional actin genes. Intron positions typical of other known plant actin genes are conserved in these clones. Four of six clones analyzed (PAc1, PAc2, PAc3, PAc4) hybridize to leaf mRNA of the same size (1.7 kb) as that reported for other plant actin mRNAs and to a slightly smaller mRNA species (1.5 kb). Five distinct subfamilies of actin-related genes were characterized which varied in size from a few members to several dozen members. It is clear from our data that other actin gene subfamilies must also exist within the genome. Possible mechanisms of actin gene amplification and genome turnover are discussed.     

Analysis of gene islands involved in methanopterin-linked C1 transfer reactions reveals new functions and provides evolutionary insights

In this study, the occurrence and chromosomal clustering of genes encoding C(1) transfer reactions linked to tetrahydromethanopterin (H(4)MPT) were analyzed in a variety of proteobacteria and in representatives of the Planctomycetes via genomic analysis or via partial sequencing by cosmid walking. Although a tendency for clustering was found common for the genes of interest, significant variations in gene order and the degree of clustering were uncovered both between and within different groups of Proteobacteria and between Proteobacteria and Planctomycetes. Phylogenetic analyses suggested that the evolution of genes encoding H(4)MPT-linked reactions in Proteobacteria involved lateral transfers within Proteobacteria and possibly between Proteobacteria and other phyla. Gene cluster comparisons revealed a number of novel genes potentially involved in the C(1) transfer reactions, and these were analyzed by mutation and expression analyses. Four genes, a homolog of pabB, and three genes conserved between methanogenic Archaea and Bacteria possessing H(4)MPT-linked functions, orfY, orf1, and afpA were shown to be involved in formaldehyde oxidation/detoxification, as judged by specific mutant phenotypes. In particular, pabB contributes to the biosynthesis of para-aminobenzoic acid, a precursor of both tetrahydrofolate and H(4)MPT, and afpA apparently encodes a novel dihydromethanopterin reductase, based on mutant complementation experiments.