Laboratory evolution reveals general and specific tolerance mechanisms for commodity chemicals.

Although strain tolerance to high product concentrations is a barrier to the economically viable biomanufacturing of industrial chemicals, chemical tolerance mechanisms are often unknown. To reveal tolerance mechanisms, an automated platform was utilized to evolve Escherichia coli to grow optimally in the presence of 11 industrial chemicals (1,2-propanediol, 2,3-butanediol, glutarate, adipate, putrescine, hexamethylenediamine, butanol, isobutyrate, coumarate, octanoate, hexanoate), reaching tolerance at concentrations 60%-400% higher than initial toxic levels. Sequencing genomes of 223 isolates from 89 populations, reverse engineering, and cross-compound tolerance profiling were employed to uncover tolerance mechanisms. We show that: 1) cells are tolerized via frequent mutation of membrane transporters or cell wall-associated proteins (e.g., ProV, KgtP, SapB, NagA, NagC, MreB), transcription and translation machineries (e.g., RpoA, RpoB, RpoC, RpsA, RpsG, NusA, Rho), stress signaling proteins (e.g., RelA, SspA, SpoT, YobF), and for certain chemicals, regulators and enzymes in metabolism (e.g., MetJ, NadR, GudD, PurT); 2) osmotic stress plays a significant role in tolerance when chemical concentrations exceed a general threshold and mutated genes frequently overlap with those enabling chemical tolerance in membrane transporters and cell wall-associated proteins; 3) tolerization to a specific chemical generally improves tolerance to structurally similar compounds whereas a tradeoff can occur on dissimilar chemicals, and 4) using pre-tolerized starting isolates can hugely enhance the subsequent production of chemicals when a production pathway is inserted in many, but not all, evolved tolerized host strains, underpinning the need for evolving multiple parallel populations. Taken as a whole, this study provides a comprehensive genotype-phenotype map based on identified mutations and growth phenotypes for 223 chemical tolerant isolates.

Fermentation and purification strategies for the production of betulinic acid and its lupane-type precursors in Saccharomyces cerevisiae.

Microbial production of plant derived, biologically active compounds has the potential to provide economic and ecologic alternatives to existing low productive, plant-based processes. Current production of the pharmacologically active cyclic triterpenoid betulinic acid is realized by extraction from the bark of plane tree or birch. Here, we reengineered the reported betulinic acid pathway into Saccharomyces cerevisiae and used this novel strain to develop efficient fermentation and product purification methods. Fed-batch cultivations with ethanol excess, using either an ethanol-pulse feed or controlling a constant ethanol concentration in the fermentation medium, significantly enhanced production of betulinic acid and its triterpenoid precursors. The beneficial effect of excess ethanol was further exploited in nitrogen-limited resting cell fermentations, yielding betulinic acid concentrations of 182 mg/L, and total triterpenoid concentrations of 854 mg/L, the highest concentrations reported so far. Purification of lupane-type triterpenoids with high selectivity and yield was achieved by solid-liquid extraction without prior cell disruption using polar aprotic solvents such as acetone or ethyl acetate and subsequent precipitation with strong acids. This study highlights the potential of microbial production of plant derived triterpenoids in S. cerevisiae by combining metabolic and process engineering.

Engineering and systems-level analysis of Saccharomyces cerevisiae for production of 3-hydroxypropionic acid via malonyl-CoA reductase-dependent pathway.

BACKGROUND: In the future, oil- and gas-derived polymers may be replaced with bio-based polymers, produced from renewable feedstocks using engineered cell factories. Acrylic acid and acrylic esters with an estimated world annual production of approximately 6 million tons by 2017 can be derived from 3-hydroxypropionic acid (3HP), which can be produced by microbial fermentation. For an economically viable process 3HP must be produced at high titer, rate and yield and preferably at low pH to minimize downstream processing costs. RESULTS: Here we describe the metabolic engineering of baker's yeast Saccharomyces cerevisiae for biosynthesis of 3HP via a malonyl-CoA reductase (MCR)-dependent pathway. Integration of multiple copies of MCR from Chloroflexus aurantiacus and of phosphorylation-deficient acetyl-CoA carboxylase ACC1 genes into the genome of yeast increased 3HP titer fivefold in comparison with single integration. Furthermore we optimized the supply of acetyl-CoA by overexpressing native pyruvate decarboxylase PDC1, aldehyde dehydrogenase ALD6, and acetyl-CoA synthase from Salmonella enterica SEacs (L641P). Finally we engineered the cofactor specificity of the glyceraldehyde-3-phosphate dehydrogenase to increase the intracellular production of NADPH at the expense of NADH and thus improve 3HP production and reduce formation of glycerol as by-product. The final strain produced 9.8 ± 0.4 g L(-1) 3HP with a yield of 13% C-mol C-mol(-1) glucose after 100 h in carbon-limited fed-batch cultivation at pH 5. The 3HP-producing strain was characterized by (13)C metabolic flux analysis and by transcriptome analysis, which revealed some unexpected consequences of the undertaken metabolic engineering strategy, and based on this data, future metabolic engineering directions are proposed. CONCLUSIONS: In this study, S. cerevisiae was engineered for high-level production of 3HP by increasing the copy numbers of biosynthetic genes and improving flux towards precursors and redox cofactors. This strain represents a good platform for further optimization of 3HP production and hence an important step towards potential commercial bio-based production of 3HP.

Assembly of a novel biosynthetic pathway for production of the plant flavonoid fisetin in Escherichia coli.

Plant secondary metabolites are an underutilized pool of bioactive molecules for applications in the food, pharma and nutritional industries. One such molecule is fisetin, which is present in many fruits and vegetables and has several potential health benefits, including anti-cancer, anti-viral and anti-aging activity. Moreover, fisetin has recently been shown to prevent Alzheimer's disease in mice and to prevent complications associated with diabetes type I. Thus far the biosynthetic pathway of fisetin in plants remains elusive. Here, we present the heterologous assembly of a novel fisetin pathway in Escherichia coli. We propose a novel biosynthetic pathway from the amino acid, tyrosine, utilizing nine heterologous enzymes. The pathway proceeds via the synthesis of two flavanones never produced in microorganisms before--garbanzol and resokaempferol. We show for the first time a functional biosynthetic pathway and establish E. coli as a microbial platform strain for the production of fisetin and related flavonols.

Establishing a synthetic pathway for high-level production of 3-hydroxypropionic acid in Saccharomyces cerevisiae via ß-alanine.

Microbial fermentation of renewable feedstocks into plastic monomers can decrease our fossil dependence and reduce global CO2 emissions. 3-Hydroxypropionic acid (3HP) is a potential chemical building block for sustainable production of superabsorbent polymers and acrylic plastics. With the objective of developing Saccharomyces cerevisiae as an efficient cell factory for high-level production of 3HP, we identified the ß-alanine biosynthetic route as the most economically attractive according to the metabolic modeling. We engineered and optimized a synthetic pathway for de novo biosynthesis of ß-alanine and its subsequent conversion into 3HP using a novel ß-alanine-pyruvate aminotransferase discovered in Bacillus cereus. The final strain produced 3HP at a titer of 13.7±0.3gL(-1) with a 0.14±0.0C-molC-mol(-1) yield on glucose in 80h in controlled fed-batch fermentation in mineral medium at pH 5, and this work therefore lays the basis for developing a process for biological 3HP production.

Highly Active and Specific Tyrosine Ammonia-Lyases from Diverse Origins Enable Enhanced Production of Aromatic Compounds in Bacteria and Saccharomyces cerevisiae.

Phenylalanine and tyrosine ammonia-lyases form cinnamic acid and p-coumaric acid, which are precursors of a wide range of aromatic compounds of biotechnological interest. Lack of highly active and specific tyrosine ammonia-lyases has previously been a limitation in metabolic engineering approaches. We therefore identified 22 sequences in silico using synteny information and aiming for sequence divergence. We performed a comparative in vivo study, expressing the genes intracellularly in bacteria and yeast. When produced heterologously, some enzymes resulted in significantly higher production of p-coumaric acid in several different industrially important production organisms. Three novel enzymes were found to have activity exclusively for phenylalanine, including an enzyme from the low-GC Gram-positive bacterium Brevibacillus laterosporus, a bacterial-type enzyme from the amoeba Dictyostelium discoideum, and a phenylalanine ammonia-lyase from the moss Physcomitrella patens (producing 230 µM cinnamic acid per unit of optical density at 600 nm [OD600]) in the medium using Escherichia coli as the heterologous host). Novel tyrosine ammonia-lyases having higher reported substrate specificity than previously characterized enzymes were also identified. Enzymes from Herpetosiphon aurantiacus and Flavobacterium johnsoniae resulted in high production of p-coumaric acid in Escherichia coli (producing 440 µM p-coumaric acid OD600 unit(-1) in the medium) and in Lactococcus lactis. The enzymes were also efficient in Saccharomyces cerevisiae, where p-coumaric acid accumulation was improved 5-fold over that in strains expressing previously characterized tyrosine ammonia-lyases.

CrEdit: CRISPR mediated multi-loci gene integration in Saccharomyces cerevisiae.

BACKGROUND: One of the bottlenecks in production of biochemicals and pharmaceuticals in Saccharomyces cerevisiae is stable and homogeneous expression of pathway genes. Integration of genes into the genome of the production organism is often a preferred option when compared to expression from episomal vectors. Existing approaches for achieving stable simultaneous genome integrations of multiple DNA fragments often result in relatively low integration efficiencies and furthermore rely on the use of selection markers. RESULTS: Here, we have developed a novel method, CrEdit (CRISPR/Cas9 mediated genome Editing), which utilizes targeted double strand breaks caused by CRISPR/Cas9 to significantly increase the efficiency of homologous integration in order to edit and manipulate genomic DNA. Using CrEdit, the efficiency and locus specificity of targeted genome integrations reach close to 100% for single gene integration using short homology arms down to 60 base pairs both with and without selection. This enables direct and cost efficient inclusion of homology arms in PCR primers. As a proof of concept, a non-native ß-carotene pathway was reconstructed in S. cerevisiae by simultaneous integration of three pathway genes into individual intergenic genomic sites. Using longer homology arms, we demonstrate highly efficient and locus-specific genome integration even without selection with up to 84% correct clones for simultaneous integration of three gene expression cassettes. CONCLUSIONS: The CrEdit approach enables fast and cost effective genome integration for engineering of S. cerevisiae. Since the choice of the targeting sites is flexible, CrEdit is a powerful tool for diverse genome engineering applications.

Novel biosensors based on flavonoid-responsive transcriptional regulators introduced into Escherichia coli.

This study describes the construction of two flavonoid biosensors, which can be applied for metabolic engineering of Escherichia coli strains. The biosensors are based on transcriptional regulators combined with autofluorescent proteins. The transcriptional activator FdeR from Herbaspirillum seropedicae SmR1 responds to naringenin, while the repressor QdoR from Bacillus subtilis is inactivated by quercetin and kaempferol. Both biosensors showed over a 7-fold increase of the fluorescent signal after addition of their specific effectors, and a linear correlation between the fluorescence intensity and externally added flavonoid concentration. The QdoR-biosensor was successfully applied for detection of kaempferol production in vivo at the single cell level by fluorescence-activated cell sorting. Furthermore, the amount of kaempferol produced highly correlated with the specific fluorescence of E. coli cells containing a flavonol synthase from Arabidopsis thaliana (fls1). We expect the designed biosensors to be applied for isolation of genes involved in flavonoid biosynthetic pathways.

Evolution reveals a glutathione-dependent mechanism of 3-hydroxypropionic acid tolerance.

Biologically produced 3-hydroxypropionic acid (3 HP) is a potential source for sustainable acrylates and can also find direct use as monomer in the production of biodegradable polymers. For industrial-scale production there is a need for robust cell factories tolerant to high concentration of 3 HP, preferably at low pH. Through adaptive laboratory evolution we selected S. cerevisiae strains with improved tolerance to 3 HP at pH 3.5. Genome sequencing followed by functional analysis identified the causal mutation in SFA1 gene encoding S-(hydroxymethyl)glutathione dehydrogenase. Based on our findings, we propose that 3 HP toxicity is mediated by 3-hydroxypropionic aldehyde (reuterin) and that glutathione-dependent reactions are used for reuterin detoxification. The identified molecular response to 3 HP and reuterin may well be a general mechanism for handling resistance to organic acid and aldehydes by living cells.

EasyClone: method for iterative chromosomal integration of multiple genes in Saccharomyces cerevisiae.

Development of strains for efficient production of chemicals and pharmaceuticals requires multiple rounds of genetic engineering. In this study, we describe construction and characterization of EasyClone vector set for baker's yeast Saccharomyces cerevisiae, which enables simultaneous expression of multiple genes with an option of recycling selection markers. The vectors combine the advantage of efficient uracil excision reaction-based cloning and Cre-LoxP-mediated marker recycling system. The episomal and integrative vector sets were tested by inserting genes encoding cyan, yellow, and red fluorescent proteins into separate vectors and analyzing for co-expression of proteins by flow cytometry. Cells expressing genes encoding for the three fluorescent proteins from three integrations exhibited a much higher level of simultaneous expression than cells producing fluorescent proteins encoded on episomal plasmids, where correspondingly 95% and 6% of the cells were within a fluorescence interval of Log10 mean ± 15% for all three colors. We demonstrate that selective markers can be simultaneously removed using Cre-mediated recombination and all the integrated heterologous genes remain in the chromosome and show unchanged expression levels. Hence, this system is suitable for metabolic engineering in yeast where multiple rounds of gene introduction and marker recycling can be carried out.

An automated workflow for multi-omics screening of microbial model organisms.

Multi-omics datasets are becoming of key importance to drive discovery in fundamental research as much as generating knowledge for applied biotechnology. However, the construction of such large datasets is usually time-consuming and expensive. Automation might enable to overcome these issues by streamlining workflows from sample generation to data analysis. Here, we describe the construction of a complex workflow for the generation of high-throughput microbial multi-omics datasets. The workflow comprises a custom-built platform for automated cultivation and sampling of microbes, sample preparation protocols, analytical methods for sample analysis and automated scripts for raw data processing. We demonstrate possibilities and limitations of such workflow in generating data for three biotechnologically relevant model organisms, namely Escherichia coli, Saccharomyces cerevisiae, and Pseudomonas putida.

Dynamic control of gene expression in Saccharomyces cerevisiae engineered for the production of plant sesquitepene α-santalene in a fed-batch mode.

Microbial cells engineered for efficient production of plant sesquiterpenes may allow for sustainable and scalable production of these compounds that can be used as e.g. perfumes and pharmaceuticals. Here, for the first time a Saccharomyces cerevisiae strain capable of producing high levels of α-santalene, the precursor of a commercially interesting compound, was constructed through a rationally designed metabolic engineering approach. Optimal sesquiterpene production was obtained by modulating the expression of one of the key metabolic steps of the mevalonate (MVA) pathway, squalene synthase (Erg9). To couple ERG9 expression to glucose concentration its promoter was replaced by the HXT1 promoter. In a second approach, the HXT2 promoter was used to express an ERG9 antisense construct. Using the HXT1 promoter to control ERG9 expression, it was possible to divert the carbon flux from sterol synthesis towards α-santalene improving the productivity by 3.4 fold. Combining this approach together with the overexpression of a truncated form of 3-hydroxyl-3-methyl-glutaryl-CoA reductase (HMGR) and deletion of lipid phosphate phosphatase encoded by LPP1 led to a strain with a productivity of 0.18mg/gDCWh. The titer was further increased by deleting DPP1 encoding a second FPP consuming pyrophosphate phosphatase yielding a final productivity and titer, respectively, of 0.21mg/gDCWh and 92mg/l of α-santalene.

Diversion of flux toward sesquiterpene production in Saccharomyces cerevisiae by fusion of host and heterologous enzymes.

The ability to transfer metabolic pathways from the natural producer organisms to the well-characterized cell factory Saccharomyces cerevisiae is well documented. However, as many secondary metabolites are produced by collaborating enzymes assembled in complexes, metabolite production in yeast may be limited by the inability of the heterologous enzymes to collaborate with the native yeast enzymes. This may cause loss of intermediates by diffusion or degradation or due to conversion of the intermediate through competitive pathways. To bypass this problem, we have pursued a strategy in which key enzymes in the pathway are expressed as a physical fusion. As a model system, we have constructed several fusion protein variants in which farnesyl diphosphate synthase (FPPS) of yeast has been coupled to patchoulol synthase (PTS) of plant origin (Pogostemon cablin). Expression of the fusion proteins in S. cerevisiae increased the production of patchoulol, the main sesquiterpene produced by PTS, up to 2-fold. Moreover, we have demonstrated that the fusion strategy can be used in combination with traditional metabolic engineering to further increase the production of patchoulol. This simple test case of synthetic biology demonstrates that engineering the spatial organization of metabolic enzymes around a branch point has great potential for diverting flux toward a desired product.

Whole genome sequencing of Saccharomyces cerevisiae: from genotype to phenotype for improved metabolic engineering applications.

BACKGROUND: The need for rapid and efficient microbial cell factory design and construction are possible through the enabling technology, metabolic engineering, which is now being facilitated by systems biology approaches. Metabolic engineering is often complimented by directed evolution, where selective pressure is applied to a partially genetically engineered strain to confer a desirable phenotype. The exact genetic modification or resulting genotype that leads to the improved phenotype is often not identified or understood to enable further metabolic engineering. RESULTS: In this work we performed whole genome high-throughput sequencing and annotation can be used to identify single nucleotide polymorphisms (SNPs) between Saccharomyces cerevisiae strains S288c and CEN.PK113-7D. The yeast strain S288c was the first eukaryote sequenced, serving as the reference genome for the Saccharomyces Genome Database, while CEN.PK113-7D is a preferred laboratory strain for industrial biotechnology research. A total of 13,787 high-quality SNPs were detected between both strains (reference strain: S288c). Considering only metabolic genes (782 of 5,596 annotated genes), a total of 219 metabolism specific SNPs are distributed across 158 metabolic genes, with 85 of the SNPs being nonsynonymous (e.g., encoding amino acid modifications). Amongst metabolic SNPs detected, there was pathway enrichment in the galactose uptake pathway (GAL1, GAL10) and ergosterol biosynthetic pathway (ERG8, ERG9). Physiological characterization confirmed a strong deficiency in galactose uptake and metabolism in S288c compared to CEN.PK113-7D, and similarly, ergosterol content in CEN.PK113-7D was significantly higher in both glucose and galactose supplemented cultivations compared to S288c. Furthermore, DNA microarray profiling of S288c and CEN.PK113-7D in both glucose and galactose batch cultures did not provide a clear hypothesis for major phenotypes observed, suggesting that genotype to phenotype correlations are manifested post-transcriptionally or post-translationally either through protein concentration and/or function. CONCLUSIONS: With an intensifying need for microbial cell factories that produce a wide array of target compounds, whole genome high-throughput sequencing and annotation for SNP detection can aid in better reducing and defining the metabolic landscape. This work demonstrates direct correlations between genotype and phenotype that provides clear and high-probability of success metabolic engineering targets. The genome sequence, annotation, and a SNP viewer of CEN.PK113-7D are deposited at http://www.sysbio.se/cenpk.

Characterization of different promoters for designing a new expression vector in Saccharomyces cerevisiae.

The widely used pESC vector series (Stratagene, La Jolla, CA, USA) with the bidirectional GAL1/GAL10 promoter provides the possibility of simultaneously expressing two different genes from a single vector in Saccharomyces cerevisiae. This system can be induced by galactose and is repressed by glucose. Since S. cerevisiae prefers glucose as a carbon source, and since its growth rate is higher in glucose than in galactose-containing media, we compared and evaluated seven different promoters expressed during growth on glucose (pTEF1, pADH1, pTPI1, pHXT7, pTDH3, pPGK1 and pPYK1) with two strong galactose-induced promoters (pGAL1 and pGAL10), using lacZ as a reporter gene and measuring LacZ activity in batch and continuous cultivation. TEF1 and PGK1 promoters showed the most constant activity pattern at different glucose concentrations. Based on these results, we designed and constructed two new expression vectors which contain the two constitutive promoters, TEF1 and PGK1, in opposite orientation to each other. These new vectors retain all the features from the pESC-URA plasmid except that gene expression is mediated by constitutive promoters.

Enhancement of farnesyl diphosphate pool as direct precursor of sesquiterpenes through metabolic engineering of the mevalonate pathway in Saccharomyces cerevisiae.

The mevalonate pathway in the yeast Saccharomyces cerevisiae was deregulated in order to enhance the intracellular pool of farnesyl diphosphate (FPP), the direct precursor for the biosynthesis of sesquiterpenes. Over-expression of the catalytic domain of HMG1, both from the genome and plasmid, resulted in higher production of cubebol, a plant originating sesquiterpene, and increased squalene accumulation. Down-regulation of ERG9 by replacing its native promoter with the regulatable MET3 promoter, enhanced cubebol titers but simultaneous over-expression of tHMG1 and repression of ERG9 did not further improve cubebol production. Furtheremore, the concentrations of squalene and ergosterol were measured in the engineered strains. Unexpectedly, significant accumulation of squalene and restoring the ergosterol biosynthesis were observed in the ERG9 repressed strains transformed with the plasmids harboring cubebol synthase gene. This could be explained by a toxicity effect of cubebol, possibly resulting in higher transcription levels for the genes under control of MET3 promoter, which could lead to accumulation of squalene and ergosterol.

High-Resolution Scanning of Optimal Biosensor Reporter Promoters in Yeast.

Small-molecule binding allosteric transcription factors (aTFs) derived from bacteria enable real-time monitoring of metabolite abundances, high-throughput screening of genetic designs, and dynamic control of metabolism. Yet, engineering of reporter promoter designs of prokaryotic aTF biosensors in eukaryotic cells is complex. Here we investigate the impact of aTF binding site positions at single-nucleotide resolution in >300 reporter promoter designs in Saccharomyces cerevisiae. From this we identify biosensor output landscapes with transient and distinct aTF binding site position effects for aTF repressors and activators, respectively. Next, we present positions for tunable reporter promoter outputs enabling metabolite-responsive designs for a total of four repressor-type and three activator-type aTF biosensors with dynamic output ranges up to 8- and 26-fold, respectively. This study highlights aTF binding site positions in reporter promoters as key for successful biosensor engineering and that repressor-type aTF biosensors allows for more flexibility in terms of choice of binding site positioning compared to activator-type aTF biosensors.

Enhancing sesquiterpene production in Saccharomyces cerevisiae through in silico driven metabolic engineering.

A genome-scale metabolic model was used to identify new target genes for enhanced biosynthesis of sesquiterpenes in the yeast Saccharomyces cerevisiae. The effect of gene deletions on the flux distributions in the metabolic model of S. cerevisiae was assessed using OptGene as the modeling framework and minimization of metabolic adjustments (MOMA) as objective function. Deletion of NADPH-dependent glutamate dehydrogenase encoded by GDH1 was identified as the best target gene for the improvement of sesquiterpene biosynthesis in yeast. Deletion of this gene enhances the available NADPH in the cytosol for other NADPH requiring enzymes, including HMG-CoA reductase. However, since disruption of GDH1 impairs the ammonia utilization, simultaneous over-expression of the NADH-dependent glutamate dehydrogenase encoded by GDH2 was also considered in this study. Deletion of GDH1 led to an approximately 85% increase in the final cubebol titer. However, deletion of this gene also caused a significant decrease in the maximum specific growth rate. Over-expression of GDH2 did not show a further effect on the final cubebol titer but this alteration significantly improved the growth rate compared to the GDH1 deleted strain.

In silico strategy to rationally engineer metabolite production: A case study for threonine in Escherichia coli.

Genetic engineering of metabolic pathways is a standard strategy to increase the production of metabolites of economic interest. However, such flux increases could very likely lead to undesirable changes in metabolite concentrations, producing deleterious perturbations on other cellular processes. These negative effects could be avoided by implementing a balanced increase of enzyme concentrations according to the Universal Method [Kacser and Acerenza (1993) Eur J Biochem 216:361-367]. Exact application of the method usually requires modification of many reactions, which is difficult to achieve in practice. Here, improvement of threonine production via pyruvate kinase deletion in Escherichia coli is used as a case study to demonstrate a partial application of the Universal Method, which includes performing sensitivity analysis. Our analysis predicts that manipulating a few reactions is sufficient to obtain an important increase in threonine production without major perturbations of metabolite concentrations.

Genome-wide systematic identification of methyltransferase recognition and modification patterns.

Genome-wide analysis of DNA methylation patterns using single molecule real-time DNA sequencing has boosted the number of publicly available methylomes. However, there is a lack of tools coupling methylation patterns and the corresponding methyltransferase genes. Here we demonstrate a high-throughput method for coupling methyltransferases with their respective motifs, using automated cloning and analysing the methyltransferases in vectors carrying a strain-specific cassette containing all potential target sites. To validate the method, we analyse the genomes of the thermophile Moorella thermoacetica and the mesophile Acetobacterium woodii, two acetogenic bacteria having substantially modified genomes with 12 methylation motifs and a total of 23 methyltransferase genes. Using our method, we characterize the 23 methyltransferases, assign motifs to the respective enzymes and verify activity for 11 of the 12 motifs.