Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotheraphy
    • Applied and Environmental Mircobiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • Data Policy
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotheraphy
    • Applied and Environmental Mircobiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • Data Policy
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Molecular Biology and Physiology

Oxidative Pathways of Deoxyribose and Deoxyribonate Catabolism

Morgan N. Price, Jayashree Ray, Anthony T. Iavarone, Hans K. Carlson, Elizabeth M. Ryan, Rex R. Malmstrom, Adam P. Arkin, Adam M. Deutschbauer
Korneel Rabaey, Editor
Morgan N. Price
Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Morgan N. Price
Jayashree Ray
Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anthony T. Iavarone
QB3/Chemistry Mass Spectrometry Facility, University of California, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hans K. Carlson
Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elizabeth M. Ryan
DOE Joint Genome Institute, Walnut Creek, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rex R. Malmstrom
DOE Joint Genome Institute, Walnut Creek, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Adam P. Arkin
Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, California, USADepartment of Bioengineering, University of California, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Adam M. Deutschbauer
Environmental Genomics and Systems Biology, Lawrence Berkeley National Laboratory, Berkeley, California, USADepartment of Plant and Microbial Biology, University of California, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Korneel Rabaey
Ghent University
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00297-18
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Using genome-wide mutant fitness assays in diverse bacteria, we identified novel oxidative pathways for the catabolism of 2-deoxy-d-ribose and 2-deoxy-d-ribonate. We propose that deoxyribose is oxidized to deoxyribonate, oxidized to ketodeoxyribonate, and cleaved to acetyl coenzyme A (acetyl-CoA) and glyceryl-CoA. We have genetic evidence for this pathway in three genera of bacteria, and we confirmed the oxidation of deoxyribose to ketodeoxyribonate in vitro. In Pseudomonas simiae, the expression of enzymes in the pathway is induced by deoxyribose or deoxyribonate, while in Paraburkholderia bryophila and in Burkholderia phytofirmans, the pathway proceeds in parallel with the known deoxyribose 5-phosphate aldolase pathway. We identified another oxidative pathway for the catabolism of deoxyribonate, with acyl-CoA intermediates, in Klebsiella michiganensis. Of these four bacteria, only P. simiae relies entirely on an oxidative pathway to consume deoxyribose. The deoxyribose dehydrogenase of P. simiae is either nonspecific or evolved recently, as this enzyme is very similar to a novel vanillin dehydrogenase from Pseudomonas putida that we identified. So, we propose that these oxidative pathways evolved primarily to consume deoxyribonate, which is a waste product of metabolism.

IMPORTANCE Deoxyribose is one of the building blocks of DNA and is released when cells die and their DNA degrades. We identified a bacterium that can grow with deoxyribose as its sole source of carbon even though its genome does not contain any of the known genes for breaking down deoxyribose. By growing many mutants of this bacterium together on deoxyribose and using DNA sequencing to measure the change in the mutants’ abundance, we identified multiple protein-coding genes that are required for growth on deoxyribose. Based on the similarity of these proteins to enzymes of known function, we propose a 6-step pathway in which deoxyribose is oxidized and then cleaved. Diverse bacteria use a portion of this pathway to break down a related compound, deoxyribonate, which is a waste product of metabolism. Our study illustrates the utility of large-scale bacterial genetics to identify previously unknown metabolic pathways.

INTRODUCTION

DNA is present in every cell, and once released into the environment by cell lysis, it is hydrolyzed by a variety of secreted or extracellular enzymes (1). DNA is hydrolyzed to oligonucleotides and then to individual deoxynucleotides or deoxynucleosides (1). Deoxynucleotides can be hydrolyzed extracellularly to deoxynucleosides (2), and deoxynucleosides can be cleaved to release 2-deoxy-d-ribose and nucleobases. Some nucleosidases are specific for 2′-deoxyribonucleosides over ribonucleosides (3), which suggests that the cleavage of deoxynucleosides is their physiological role.

We are studying the catabolism of the end products of DNA hydrolysis, including deoxyribose. The only previously known pathway for catabolizing deoxyribose (of which we are aware) involves deoxyribose kinase (deoK) (4) and deoxyribose 5-phosphate aldolase (deoC) (5), which yield acetaldehyde and d-glyceraldehyde 3-phosphate (Fig. 1A). However we recently screened 28 bacteria for growth on deoxyribose and found that the plant-root-associated bacterium Pseudomonas simiae WCS417 grows with deoxyribose as the sole source of carbon (6), even though its genome does not contain either deoK or deoC. We also tested which genes are important for the growth of P. simiae under various conditions, including growth with deoxyribose as the sole source of carbon (6), by using a pool of over 100,000 randomly barcoded transposon mutants of P. simiae. If a gene is important for growth under a condition, then the abundance of its mutants will decrease during growth, and the abundance of each mutant in the pool can be measured by amplifying and sequencing the DNA barcodes (7).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

An oxidative pathway for deoxyribose utilization in Pseudomonas simiae WCS417. (A) The previously described kinase/aldolase pathway. (B) Gene fitness in P. simiae. Each point shows the gene’s fitness value (x axis) in a different experiment. Fitness values under −1 indicate a defect in growth. Values outside of the plot's range are shown at −4 or +1. Within each gene's section, the y axis is random. Experiments with deoxyribose or deoxyribonate as the carbon source are indicated by symbols. (C) The proposed pathway for deoxyribose catabolism in P. simiae. Putative enzymes are color coded to match the gene descriptions in panel B.

Here, we use these genetic data to identify a novel oxidative pathway for deoxyribose catabolism. An intermediate in the proposed pathway is 2-deoxy-d-ribonate, which may also be widely available in the environment, as deoxyribonate has been detected in human urine (8). Deoxyribonate might form by the oxidation of deoxyribose, by the repair of oxidatively damaged DNA (9), or from a spontaneous reaction of hydrogen peroxide with 2-keto-3-deoxy-d-gluconate (which is an intermediate in hexuronate metabolism and in variants of the Entner-Doudoroff pathway of glycolysis) (10).

We used the pool of mutants to study the growth of P. simiae on deoxyribonate and found that P. simiae uses many of the same genes to catabolize deoxyribonate as deoxyribose. Furthermore, these genes are induced during growth on either deoxyribose or deoxyribonate. We also reconstituted the enzymatic conversion of deoxyribose to deoxyribonate to ketodeoxyribonate in vitro.

We then used pooled mutant fitness assays to study the catabolism of deoxyribose and deoxyribonate in other bacteria. In Paraburkholderia bryophila 376MFSha3.1 and in Burkholderia phytofirmans PsJN, similar genes to those in P. simiae are involved in the catabolism of both deoxyribose and deoxyribonate. Klebsiella michiganensis M5al uses an analogous pathway with acyl coenzyme A (acyl-CoA) intermediates to consume deoxyribonate (but not deoxyribose). Thus, diverse bacteria use oxidative pathways to consume deoxyribose and deoxyribonate.

RESULTS

The genetic basis of growth on deoxyribose and deoxyribonate in Pseudomonas simiae.To supplement the previously described mutant fitness data for P. simiae WCS417 (6), we tested deoxyribose as the carbon source on 2 additional days, and we also tested deoxyribonate as the carbon source. For each experiment, we grew a library of over 100,000 barcoded mutant strains together and used DNA sequencing to measure the change in each strain’s abundance. Overall, our data set for P. simiae included 6 genome-wide assays of fitness in deoxyribose, 3 assays of fitness in deoxyribonate, and 153 other fitness assays from diverse conditions, including 49 other carbon sources.

From these data, we identified 13 genes that were consistently important for utilizing deoxyribose and were not important for fitness under most other conditions (see Materials and Methods). The fitness data for these genes are shown in Fig. 1B: each fitness value is the log2 change of the relative abundance of mutants in that gene during the experiment. In the genome of P. simiae, 11 of the 13 genes cluster into three groups of nearby genes that are on the same strand and are probably operons.

Of the 13 genes, eight are likely to code for catabolic enzymes:

  • two subunits of a molybdenum-dependent dehydrogenase (PS417_10890 and PS417_10885) that are distantly related to isoquinoline 1-oxidoreductase (iorAB [49% or 29% amino acid identity]) (11)

  • a cytochrome c that is in the same putative operon (PS417_10880) and might accept electrons from the iorAB-like dehydrogenase

  • a dehydrogenase of the short-chain dehydrogenase/reductase (SDR) superfamily (PS417_07245)

  • a β-keto acid cleavage enzyme (PS417_07250) (12)

  • a lactonase (PS417_07255)

  • glycerate kinase (PS417_13970)

  • and acetyl-CoA C-acetyltransferase (PS417_10515).

The remaining genes that are specifically important for deoxyribose utilization are two putative transporters, two regulatory genes, and mobA (PS417_21490), which is involved in the biosynthesis of molybdenum cofactor and is probably required for the activity of the molybdenum-dependent dehydrogenase (PS417_10890 to PS417_10885). Most of these genes are also important for growth on deoxyribonate (Fig. 1B). The exceptions are the genes coding for the two subunits of the molybdenum-dependent dehydrogenase, the associated cytochrome c gene, mobA, and perhaps the lactonase gene (Fig. 1B); we therefore propose that these genes act upstream of deoxyribonate formation (Fig. 1C).

To confirm the genome-wide fitness data, we obtained individual mutant strains for some of the key genes and grew them with either d-glucose, deoxyribose, or deoxyribonate as the sole source of carbon. We found that a transposon mutant with mutation of the iorB-like subunit of the molybdenum-dependent dehydrogenase (PS417_10885) grows on glucose or deoxyribonate but not on deoxyribose (Fig. 2). Mutants of the SDR dehydrogenase (PS417_07245) or the β-keto acid cleavage enzyme (PS417_07250) grow on glucose but not on deoxyribose or deoxyribonate (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Growth curves for individual transposon mutants. We compared the growth of wild-type P. simiae WCS417 and of individual transposon mutant strains in defined media with three carbon sources: d-glucose, deoxyribose, and deoxyribonate. The mutant of PS417_04200 (2-ketoglutaric semialdehyde dehydrogenase) is included as a control. Each curve is the median of six replicates. Each carbon source was provided at 10 mM.

The proposed pathway for deoxyribose catabolism in P. simiae.The genetic data suggest that P. simiae WCS417 catabolizes deoxyribose by an oxidative pathway (Fig. 1C). Because the molybdenum-dependent dehydrogenase was not required for growth on deoxyribonate, and because the lactonase was more important for growth on deoxyribose than on deoxyribonate (Fig. 1B), it appears that the molybdenum-dependent dehydrogenase oxidizes deoxyribose to a lactone, which is hydrolyzed to deoxyribonate by the lactonase. We tentatively predict that the lactone intermediate is 2-deoxy-d-ribono-1,5-lactone instead of a 1,4-lactone because the pyranose form of deoxyribose is more stable in solution at 30°C (13).

Oxidation of the deoxyribonate by the SDR enzyme could yield 2-deoxy-3-keto-d-ribonate, which is a β-keto acid and could be cleaved by the β-keto acid cleavage enzyme (PS417_07250). β-Keto acid cleavage enzymes remove an acetyl group from a β-keto acid and transfer the acetyl group to acetyl-CoA (12). In the proposed pathway, the β-keto acid cleavage enzyme is expected to yield d-glyceryl-CoA and acetoacetate. The requirement for the glycerate kinase suggests that the d-glyceryl-CoA is converted to glycerate. We looked for a thiolase or CoA-transferase that might act on d-glyceryl-CoA and that is important for deoxyribose utilization in the mutant fitness data. There are many candidate genes in the genome, but none with that phenotype. We suspect that the conversion of d-glyceryl-CoA to d-glycerate is genetically redundant and therefore cannot be identified by assaying mutant strains with only one transposon insertion. The glycerate kinase of P. simiae (PS417_13970) is 51% identical to glycerate kinase 1 from Escherichia coli (garK), which forms 2-phosphoglycerate (14, 15). 2-Phosphoglycerate is an intermediate in glycolysis and thus links the proposed pathway to central metabolism.

In the preceding analysis, we proposed reactions for all of the putative enzymes with strong negative phenotypes on deoxyribose (relative to other conditions), except for the acetyl-CoA C-acetyltransferase (PS417_10515). That enzyme converts acetoacetate (which is produced by the β-keto acid cleavage enzyme) back to acetyl-CoA, but acetoacetate must first be activated to acetoacetyl-CoA. Acetoacetate might be activated by a CoA-transferase (PS417_10525 to PS417_10520) or a CoA-synthetase (there are several candidates in the genome). None of these genes is important for the utilization of deoxyribose, so we propose that this activity is genetically redundant. Overall, the proposed pathway accounts for all of the putative enzymes that are specifically important for deoxyribose utilization in P. simiae.

Oxidation of deoxyribose to ketodeoxyribonate in vitro.We tested the activity of the deoxyribose dehydrogenase and the deoxyribonate dehydrogenase in vitro. (For details, see the Appendix.) First, we overexpressed the iorA- and iorB-like subunits of deoxyribose dehydrogenase (PS417_10890 and PS417_10885) in E. coli and tested the resulting cell lysate. We found that this preparation oxidized deoxyribose to deoxyribonate with phenazine methosulfate as the electron acceptor. The lactonase (PS417_07255) was not necessary: deoxyribonolactone may hydrolyze spontaneously in water, or it might be hydrolyzed by enzymes in the cell lysate, or the deoxyribose dehydrogenase may form deoxyribonate instead of deoxyribonolactone. Second, we purified the deoxyribonate dehydrogenase (PS417_07245) and showed that it oxidized deoxyribonate with NADH as the electron acceptor. Tandem mass spectrometry (MS/MS) showed that the oxidized product was a ketodeoxyribonate, but we were not able to determine which position was oxidized. We also tested the two dehydrogenases on each others’ substrates (deoxyribose dehydrogenase on deoxyribonate or deoxyribonate dehydrogenase on deoxyribose) and did not observe any activity. Thus, biochemical assays confirmed the proposed activities of deoxyribose dehydrogenase and deoxyribonate dehydrogenase.

Accessory genes for deoxyribose catabolism in P. simiae.Based on the proposed pathway, we can explain the roles of the other genes with specific phenotypes on deoxyribose (Fig. 1B). First, the putative operon that includes the lactonase, the SDR deoxyribonate dehydrogenase, and the β-keto acid cleavage enzyme also includes two genes from the major facilitator superfamily (MFS) of transporters (PS417_07260 and PS417_07265). We propose that these two genes are required for the uptake of deoxyribonate or deoxyribonolactone. We suspect that deoxyribose is oxidized in the periplasm, because the iorB-like subunit has a potential signal for Tat-dependent export (RRRFLA) near its N terminus, and the associated cytochrome c protein has a membrane anchor and may be exported. In contrast, the lactonase and the deoxyribonate dehydrogenase are predicted to be in the cytoplasm (16). If deoxyribose is oxidized in the periplasm, then deoxyribonate could be the substrate of the transporter when either deoxyribose or deoxyribonate is the carbon source, which would explain why these transporter genes were important under both conditions. However we do not understand why two different transporter genes appear to be required. These transporter genes are at the end of the putative operon, so polar effects are not likely.

Second, a gntR-like transcription factor (PS417_07240) is present at the beginning of the deoxyribonate oxidation cluster. By comparing the sequences upstream of this cluster and upstream of similar clusters from other strains of Pseudomonas or from Burkholderia graminis, we identified the conserved palindrome GTGATCAC. This sequence occurs at positions −47 to −40 relative to the predicted start codon of PS417_07240 and at similar positions relative to its homologs. The same motif was identified among the binding sites of the transcriptional activator AkgR from Rhodobacter sphaeroides (17), which is 32% identical to PS417_07240. We propose that PS417_07240 regulates the expression of the cluster by binding to this motif.

Third, the transcription factor PS417_13975 (sdaR) is probably an activator for the expression of the glycerate kinase (PS417_13970) gene, which is downstream of sdaR. Similarly, a close homolog of PS417_13975 from Pseudomonas fluorescens Pf-5 (PFL_3379 [88% identical]) is predicted in RegPrecise (18) to regulate a downstream glycerate kinase.

Finally, mobA (PS417_21490) is expected to be necessary for the attachment of a nucleotide to molybdenum cofactor to form a molybdenum dinucleotide cofactor. This cofactor is probably required for the activity of deoxyribose dehydrogenase, which is homologous to the coxLS subunits of a CO dehydrogenase that has a molybdenum-cytosine dinucleotide cofactor (19). While mobA is required for the activity of deoxyribose dehydrogenase, it is not important for oxidation of deoxyribonate (Fig. 1B), which is expected because none of the other enzymes in the pathway uses a molybdenum cofactor. Similarly, genes for the biosynthesis of molybdopterin (a precursor to molybdenum cofactor) are important for utilizing deoxyribose but not deoxyribonate (moaABCE: PS417_21395, PS417_08670, PS417_04865, and PS417_04875). These genes’ phenotypes are not specific to deoxyribose (and hence are not shown in Fig. 1B) because they are also important for fitness under other conditions such as inosine utilization and nitrate stress (6).

Induction of the oxidative pathway in P. simiae.Next we asked if the enzymes in the oxidative pathway are induced when P. simiae grows on deoxyribose or deoxyribonate. As shown in Fig. 3, all of these proteins seem to be expressed more highly during growth on deoxyribose or deoxyribonate than during growth on glucose. (The iorA-like subunit of deoxyribose dehydrogenase and the associated cytochrome were not detected; they may be membrane associated.) The differences in expression levels were statistically significant for the lactonase, the deoxyribonate dehydrogenase, and acetyl-CoA C-acetyltransferase (false-discovery rate of <5% by analysis of variance [ANOVA]).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Levels of protein expression in Pseudomonas simiae WCS417. We show the normalized abundance of enzymes in the deoxyribose utilization pathway, as estimated from the integrated ion abundances of the detected peptide ions. There are three biological replicates of each condition, and the y axis within each panel is random. Note the log scale for the x axis.

Comparative genetics of the oxidative pathway.To determine if other bacteria use the oxidative pathway, we examined the genomic organization of the putative enzymes and their homologs. We focused on homologs from bacteria for which mutant libraries are available, and we obtained genome-wide mutant fitness data on deoxyribose, deoxyribonate, or deoxynucleosides for three additional bacteria: Burkholderia phytofirmans PsJN, Paraburkholderia bryophila 376MFSha3.1, and Klebsiella michiganensis M5al. This section gives an overview of the genes that were important for fitness under these conditions, the chromosomal organization of these genes, and the catabolic pathways that we propose.

Although these three additional bacteria, as well as P. simiae, grow on deoxyribonate by oxidative pathways, their other catabolic capabilities vary. Figure 4A shows the relationships between the proposed pathways for the catabolism of deoxyribose, deoxyribonate, or deoxynucleosides in the four bacteria based on our analysis of the genetic data. The full rationale for the assignment of each bacterium to a particular pathway is contained in the following sections.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Genome organization of the deoxyribose utilization pathways. (A) Overview of the catabolism of deoxyribose and related compounds by 4 bacteria. For each transformation, we show which of the bacteria performs it. (B) Deoxyribonate oxidation in four bacteria. (C) The deoxyribose-phosphate aldolase pathway in three bacteria. (D) Deoxyribose dehydrogenase from P. simiae and vanillin dehydrogenase from P. putida. In panels B to D, adjacent genes are connected by horizontal lines, genes that are potential orthologs are the same color, and the scale is the same. All of the putative enzymes shown are important for the utilization of deoxyribose or related compounds, except for the iorAB-like cluster in P. putida and the cleavage enzyme and deoxyribonyl-CoA dehydrogenase in B. phytofirmans. We lack fitness data for deoC from P. bryophila.

First, the deoxyribonate oxidation genes are clustered together, with some variation, in diverse bacteria (Fig. 4B). The clusters from two related bacteria, Paraburkholderia bryophila 376MFSha3.1 and Burkholderia phytofirmans PsJN, contain the deoxyribonate dehydrogenase and β-keto acid cleavage enzyme genes, and these genes were important for deoxyribonate utilization. The genomes of both of these bacteria also encode the kinase/aldolase pathway for deoxyribose utilization (Fig. 4C). Genes from both pathways were important for deoxyribose utilization, so we propose that these bacteria consume deoxyribose by the kinase/aldolase and oxidative pathways in parallel.

Second, Klebsiella michiganensis M5al has a gene cluster that shares some genes with the deoxyribonate oxidation cluster, but K. michiganensis does not encode deoxyribonate dehydrogenase or the β-keto acid cleavage enzyme (Fig. 4B). Nevertheless, this gene cluster was important for deoxyribonate utilization by K. michiganensis. We propose that K. michiganensis consumes deoxyribonate by another oxidative pathway with acyl-CoA intermediates. Also, K. michiganensis does not grow on deoxyribose, but it can grow on deoxynucleosides via DeoABD, which can convert them to deoxyribose-5-phosphate and deoxyribose 5-phosphate aldolase (Fig. 4C).

Third, the deoxyribose dehydrogenase gene of P. simiae is part of a conserved gene cluster, but its biological role does not seem to be conserved. The conserved cluster includes iorA- and iorB-like subunits and a putative cytochrome c (Fig. 4D). (In contrast, the iorAB genes of isoquinoline 1-oxidoreductase are more distantly related, and no gene for a cytochrome c was found near iorB [11] or near its closest homologs in complete genomes [WP_083944220.1 or WP_113935277.1].) The deoxyribose dehydrogenase-like gene cluster is found in some other Pseudomonas species and in some genera of Proteobacteria such as Bordetella and Agrobacterium. Although the cluster is found in Pseudomonas putida KT2440, P. putida does not grow on deoxyribose or deoxyribonate, and its genome does not contain the deoxyribonate oxidation genes. While studying the seemingly unrelated process of vanillin (4-hydroxy-3-methoxybenzaldehyde) catabolism by P. putida, we identified mutant phenotypes for this cluster. We propose that the cluster from P. putida (PP_3621 to PP_3623) encodes a novel vanillin dehydrogenase.

We will now describe the genetic data for each of the bacteria in more detail.

Oxidation of deoxyribose and deoxyribonate in Paraburkholderia bryophila and Burkholderia phytofirmans.Paraburkholderia bryophila 376MFSha3.1 and Burkholderia phytofirmans PsJN contain close homologs of the deoxyribonate dehydrogenase and the β-keto acid cleavage enzyme (50 to 60% identical to the proteins in P. simiae), but they do not contain close homologs of the deoxyribose dehydrogenase or the lactonase. In P. bryophila, the standard kinase/aldolase pathway was important for growth on deoxyribose, but not on deoxyribonate or glucose (Fig. 5). Enzymes for the oxidative catabolism of deoxyribonate—the deoxyribonate dehydrogenase (H281DRAFT_00644) and the β-keto acid cleavage enzyme (H218DRAFT_00641)—were important for growth on both deoxyribose and deoxyribonate. Also important for growth on both substrates were glycerate kinase (H281DRAFT_00894) and genes for the conversion of acetoacetate to acetoacetyl-CoA to acetyl-CoA (H281DRAFT_04495 and H281DRAFT_00852). These genes were not important for glucose utilization. These data strongly suggest that the oxidative pathway in P. bryophila is the same as in P. simiae.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Mutant fitness data for the kinase/aldolase and oxidative pathways in Paraburkholderia bryophila 376MFSha3.1. Fitness values for the deoxyribose 5-phosphate aldolase (deoC) and for subunit B of acetoacetyl-CoA transferase were not estimated due to inadequate coverage. The data are plotted as in Fig. 1B.

For most of the genes in the oxidative pathway, the mutant phenotypes were stronger (the fitness values were more negative) during growth on deoxyribonate than on deoxyribose (Fig. 5). It appears that P. bryophila uses the oxidative pathway to grow on deoxyribonate, but deoxyribose is metabolized by the oxidative pathway in parallel with the kinase/aldolase pathway.

We did not identify any candidates for deoxyribose dehydrogenase in the fitness data from P. bryophila. This activity may be genetically redundant or nonspecific, or the kinase/aldolase pathway may suffice for efficient growth on deoxyribose.

In the genome of P. bryophila, the oxidative enzymes are clustered with a transporter (H281DRAFT_00642), a regulator of the araC family (H281DRAFT_00640), and a putative β-hydroxyacyl-CoA dehydrogenase (H281DRAFT_00643) (Fig. 4B). These genes are also important for the utilization of both deoxyribose and deoxyribonate, although the β-hydroxyacyl-CoA dehydrogenase mutants have a milder phenotype (Fig. 5). Based on the phenotype of a homolog of this protein in K. michiganensis (see below), we suspect that the β-hydroxyacyl-CoA dehydrogenase acts on deoxyribonyl-CoA. So, we speculate that P. bryophila may use a third parallel pathway for deoxyribose catabolism.

When we performed fitness assays with B. phytofirmans, we found similar results to P. bryophila for deoxyribonate utilization (Fig. 6). In particular, deoxyribonate dehydrogenase (BPHYT_RS04775), the β-keto acid cleavage enzyme (BPHYT_RS04760), and glycerate kinase (BPHYT_RS09440) are all important during growth on deoxyribonate, although β-keto acid cleavage enzyme mutants have a milder phenotype. One unexplained result is that the acetyl-CoA C-acetyltransferase (BPHYT_RS09150) is detrimental to fitness on deoxyribonate (mutants in this gene increased more than 4-fold in abundance relative to other mutants in the pooled assay). During growth of B. phytofirmans on deoxyribose, we observe strong phenotypes for the enzymes in the kinase/aldolase pathway (Fig. 6). Genes in the oxidative pathway are somewhat important for fitness during growth on deoxyribose at 10 mM (but not at 5 mM). Although these defects are modest, the defects are larger during growth on 10 mM deoxyribose than in most of 80 other conditions that we previously studied (data from reference 6) (Fig. 6).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Mutant fitness data for the aldolase and oxidative pathways in Burkholderia phytofirmans PsJN. Fitness values outside of the plot's range are shown at −5 or +2. The data are plotted as in Fig. 1B.

Overall, we conclude that both Paraburkholderia bryophila and Burkholderia phytofirmans catabolize deoxyribonate by the same oxidative pathway that P. simiae uses. Furthermore, the oxidative pathway contributes to deoxyribose utilization in P. bryophila and probably in B. phytofirmans as well. However, in contrast to P. simiae, both P. bryophila and B. phytofirmans also use the aldolase-kinase pathway for deoxyribose catabolism.

Oxidation of deoxyribonate in Klebsiella michiganensis.As mentioned above, the genome of Klebsiella michiganensis contains a gene cluster that shares some genes with the deoxyribonate oxidation cluster of P. simiae (Fig. 4B). In K. michiganensis, the cluster is important for the utilization of deoxyribonate but not for the utilization of deoxynucleosides or deoxynucleotides (Fig. 7) or under any of 92 other conditions that we previously profiled (6). The cluster includes a putative CoA-synthetase (BWI76_RS23695), a β-hydroxyacyl-CoA dehydrogenase (BWI76_RS23705), a thiolase (BWI76_RS23710), a CoA-transferase (BWI76_23700), and glycerate kinase (BWI76_RS23730). We propose that the CoA-synthetase converts deoxyribonate to deoxyribonyl-CoA, the dehydrogenase oxidizes this to 3-ketodeoxyribonyl-CoA, and the thiolase cleaves this to glyceryl-CoA and acetyl-CoA (see Fig. S1 in the supplemental material). These transformations are analogous to the oxidative pathway of the other three bacteria, but with CoA bound intermediates. In all four bacteria, deoxyribonate is oxidized and cleaved to yield glyceryl-CoA and acetyl-CoA, but the point at which the CoA thioesterification occurs varies. In K. michiganensis, the cluster also includes a CoA-transferase that may convert the glyceryl-CoA to glycerate and a glycerate kinase that forms 2-phosphoglycerate, which feeds into glycolysis (as in the other bacteria).

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

Mutant fitness data for the utilization of deoxyribonate and various deoxynucleotides/deoxynucleosides in Klebsiella michiganensis M5al. dAMP is 2′-deoxyadenosine 5′-phosphate, dA is 2′-deoxyadenosine, and dC is 2′-deoxycytidine.

FIG S1

The proposed pathway of deoxyribonate oxidation in Klebsiella michiganensis M5al. We identified the closest characterized homolog of each putative enzyme using PaperBLAST (30). The CoA-synthetase is 30% identical to the long-chain fatty-acid CoA-synthetase of B. subtilis (lcfB, or O07610). The dehydrogenase is 45% identical to a 3-oxoacyl-(acyl-carrier protein) reductase (fabG, or P0A0H9), which carries out a reaction in the opposite (reductive) direction from that shown above and with the CoA attached to acyl carrier protein. However an ACP substrate is normally part of biosynthesis, which is difficult to reconcile with this gene's mutant phenotype or its presence in the deoxyribonate utilization cluster. The thiolase is 47 to 50% identical to thiolases that act on acetoacetyl-CoA (see P45359 in Swiss-Prot or see PfGW456L13_2411 in the Fitness Browser). The CoA-transferase is 42% identical to a propionate CoA-transferase (Q9L3F7), which uses propionyl-CoA and acetate as substrates, so we show acetate as the CoA acceptor, but it might use another acceptor, such as succinate. Download FIG S1, PDF file, 0.2 MB.
Copyright © 2019 Price et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The deoxyribonate utilization cluster of K. michiganensis also includes a close homolog (70% amino acid identity) of the putative deoxyribonolactonase from P. simiae. This gene (BWI76_RS23720) was also somewhat important for deoxyribonate utilization (fitnesses of −1.3 and −1.1). It is possible that the latter phenotype is due to polar effects that inhibit the expression of downstream genes (which include glycerate kinase). Or, because deoxyribonolactone and deoxyribonate interconvert in solution, the lactonase’s activity may be useful. We also observed a similar mild phenotype for the lactonase during deoxyribonate utilization in P. simiae (Fig. 1B).

As in P. simiae, the cluster in K. michiganensis contains two MFS transporter genes (Fig. 4B). These are probably orthologs of the two transporters in P. simiae’s cluster. (Each pair of proteins is over 50% identical between the two organisms.) As in P. simiae, both of the transporters in K. michiganensis are important for growth on deoxyribonate, although BWI76_RS23715 mutants have a stronger phenotype than BWI76_RS23725 mutants (average fitness of −4.0 versus −1.4).

Although the genome of K. michiganesis contains a deoxyribose-phosphate aldolase (BWI76_RS04130), it does not seem to encode deoxyribose kinase. This might explain why it did not grow with deoxyribose as the carbon source. We did obtain fitness data with three deoxynucleotides or deoxynucleosides as carbon sources: 2′-deoxyadenosine 5′-phosphate, 2′-deoxyadenosine, and 2′-deoxycytidine (Fig. 7). These data confirmed that K. michiganensis uses the aldolase pathway to break down deoxynucleosides. Deoxynucleosides are probably converted to deoxyribose 5-phosphate by the combination of a phosphorylase (either DeoA or DeoD; BWI76_RS04135 or BWI76_RS04145) which forms deoxyribose 1-phosphate and a phosphopentomutase (DeoB, BWI76_RS04140) which converts it to deoxyribose 5-phosphate. deoABD were important for growth on deoxynucleotides and deoxynucleosides (all fitness values were under −1, except for deoD on 2′-deoxycytidine). Based on these data, we propose that K. michiganensis uses an oxidative pathway with acyl-CoA intermediates to metabolize deoxyribonate, it uses the aldolase pathway to metabolize deoxynucleosides, and it cannot metabolize deoxyribose.

Deoxyribose dehydrogenase is closely related to a novel vanillin dehydrogenase from Pseudomonas putida KT2440.We found that all three genes from the deoxyribose dehydrogenase-like cluster of Pseudomonas putida KT2440 (PP_3621 to PP_3623) were specifically important for growth on 5 mM vanillin (all fitness values of ≤−2 in two replicate experiments). These genes were not important for growth with other aromatic compounds as carbon sources (vanillic acid, benzoic acid, ferulic acid, or p-coumaric acid) (see Fig. S2 in the supplemental material).

FIG S2

Fitness data for aromatic compound degradation by Pseudomonas putida KT2440. Vanillin (4-hydroxy-3-methoxybenzaldehyde) is oxidized to vanillate, demethylated to protocatechuate (3,4-dihydroxybenzoate), and catabolized via the ortho-cleavage pathway and 3-oxoadipate to yield acetyl-CoA and succinyl-CoA. p-Coumarate [3-(4-hydroxyphenyl)-2-propenoic acid; also known as 4-hydroxycinnamic acid] appears to be activated to p-coumaryl-CoA, hydrated and cleaved to p-hydroxybenzaldehyde, oxidized to p-hydroxybenzoate, and hydroxylated to protocatechuate. Since mutants of vdh (PP_3357) are deficient in growth on p-coumaric acid, its primary substrate is probably p-hydroxybenzaldehyde, not vanillin. (The right-most glucose experiment, shown here as a control, is from M. G. Thompson, J. M. Blake-Hedges, P. Cruz-Morales, J. F. Barajas, S. C. Curran, et al., bioRxiv, https://doi.org/10.1101/450254.) Download FIG S2, PDF file, 0.3 MB.
Copyright © 2019 Price et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Because PP_3621 to PP_3623 are important for growth on vanillin but not on vanillic acid (which is the oxidized form), we propose that it is the major vanillin dehydrogenase in P. putida. Although another vanillin dehydrogenase from P. putida (vdh, or PP_3357) has been studied biochemically, vdh is not required for vanillin utilization (20). (Our data instead suggest a role for vdh as a p-hydroxybenzaldehyde dehydrogenase during the catabolism of p-coumaric acid [Fig. S2].) Vanillin utilization by P. putida requires the uptake of molybdate (21), which we can now explain because molybdate is a precursor for the molybdenum cofactor of PP_3621 to PP_3623. The vanillin dehydrogenase is expected to produce vanillate, which is demethylated (by PP_3736 to PP_3737) to protocatechuic acid, which feeds into the ortho-cleavage pathway that begins with protocatechuate 3,4-dioxygenase (PP_4655 to PP_4656). The demethylase and dioxygenase genes are also specifically important during growth on vanillin (all fitness values under −2 [Fig. S2]).

The genome of P. putida does not contain the deoxyribonate oxidation genes, and P. putida does not grow with either deoxyribose or deoxyribonate as the sole carbon source. Conversely, none of the four other bacteria that we studied grew with vanillin as the sole carbon source, and their genomes do not encode orthologs of the vanillin O-demethylase. However all three components of P. simiae’s deoxyribose dehydrogenase are very similar to the components of P. putida’s vanillin dehydrogenase (84 to 92% amino acid identity). We speculate that both enzymes are active on a broader range of substrates.

DISCUSSION

Oxidation and cleavage of deoxyribose and deoxyribonate.Using a genetic approach, we identified pathways for the oxidation of deoxyribose and/or deoxyribonate in four genera of bacteria. Although we have biochemical evidence for the conversion of deoxyribose to a ketodeoxyribonate in Pseudomonas simiae, we do not have direct evidence for the cleavage to d-glyceryl-CoA and acetyl-CoA. However the importance of glycerate kinase for deoxyribonate utilization by four different bacteria is difficult to explain otherwise. Also, for the three bacteria that use a β-keto cleavage enzyme during growth on deoxyribonate, experimental data from a homologous enzyme is consistent with our proposal. Specifically, in vitro assays of the β-keto cleavage enzyme BKACE_178 (UniProt accession no. A6X2V8), which is 61% identical to PS417_07250, found that it had some activity with 5-hydroxy-β-ketohexanoate, β-ketopentanoate, β-ketoisocaproate, and β-ketohexanoate as substrates (12). These substrates are similar to 2-deoxy-3-keto-ribonate (which is 4,5-dihydroxy-β-ketopentanoate). Furthermore, the oxidative cluster in Klebsiella michiganensis contains the glycerate kinase gene along with close homologs of three genes from the oxidative cluster of P. simiae (the lactonase gene and both MFS transporter genes [Fig. 4B]). This suggests that these gene clusters evolved to catabolize deoxyribonate via glycerate and 2-phosphoglycerate.

What are the natural substrates of the oxidative pathways?We identified deoxyribose dehydrogenase in just one of the four bacteria that grow on deoxyribonate. Furthermore, the natural role of the deoxyribose dehydrogenase of P. simiae remains uncertain. Its close similarity to the major vanillin dehydrogenase of Pseudomonas putida suggests that both enzymes might act on a broader range of substrates. In addition, P. simiae grows much more slowly on deoxyribose than on deoxyribonate, with a doubling time of around 12 h instead of around 2 h (Fig. 2), which suggests that deoxyribose might not be the natural substrate. Overall, it appears that deoxyribose utilization is not the main reason that these bacteria maintain a cluster for deoxyribonate oxidation.

We also wondered if 2-deoxyglucose might be a substrate for the oxidative pathway. A strain of Pseudomonas was reported to grow on 2-deoxyglucose by oxidizing it to 2-deoxy-3-ketogluconate (22). This transformation is chemically analogous to the conversion of 2-deoxyribose to 2-deoxy-3-ketoribonate. We tested the five bacteria considered here with 2.5 to 20 mM 2-deoxyglucose as a carbon source and found that none of them grew. Also, 2-deoxyglucose is not, as far as we know, a naturally occurring compound.

Instead, we propose that deoxyribonate is the natural substrate for the oxidative pathways. There are several plausible pathways by which deoxyribonate may form. First, nonspecific oxidation of deoxyribose, which may occur in Paraburkholderia bryophila and in Burkholderia phytofirmans, could yield deoxyribonate. Although we did not identify the genes responsible for this nonspecific activity, other sugar 1-dehydrogenase enzymes are known to be promiscuous, such as bifunctional l-arabinose/d-galactose 1-dehydrogenases (6, 23, 24). Also, because the deoxyribose dehydrogenase from P. simiae probably acts on substrates in the periplasm, some bacteria might obtain energy (but not carbon) by oxidizing deoxyribose in the periplasm without metabolizing it further. Deoxyribose in DNA is abundant, so even a low rate of nonspecific oxidation could select for deoxyribonate catabolism. All four of the deoxyribonate-catabolizing bacteria in this study are associated with plant roots. Also, the deoxyribonate utilization cluster is present in about 20% of the Pseudomonas genomes in PATRIC (25) and in 34% of plant-associated Pseudomonas, but it is rarely found in Pseudomonas isolates from humans (under 1%). DNA is estimated to be ∼3% of the dry weight of a typical plant cell (26), which implies that deoxyribose in DNA is around 1% of biomass. Second, deoxyribonate or deoxyribonolactone could be released after oxidative damage of DNA (9). Third, deoxyribonate is proposed to form spontaneously from the metabolic intermediate 2-keto-3-deoxy-d-gluconate and hydrogen peroxide (10). Although we do not know of any measurement of the concentration of deoxyribonate in plants or in soils, 2-deoxyribonate is present in human urine at about 40-fold less than the concentration of creatinine (8), or around 25 mg/day excreted by the typical adult.

MATERIALS AND METHODS

Strains and growth media.Pseudomonas simiae WCS417 and Paraburkholderia bryophila 376MFSha3.1 were provided by Jeff Dangl (University of North Carolina). Pools of barcoded transposon mutants of Pseudomonas simiae WCS417, Pseudomonas putida KT2440, Burkholderia phytofirmans PsJN, and Klebsiella michiganensis M5al were described previously (6, 27). Each mutant has a transposon insertion that includes kanR (providing kanamycin resistance), as well as a 20-nucleotide barcode flanked by common priming sites.

Individual strains from the library of mutants for P. simiae were obtained from an arrayed collection (28). Briefly, the library was spread on agar plates containing LB (Luria-Bertani) plus kanamycin, and colonies were picked and arrayed into 384-well plates containing LB with 7.5% glycerol and kanamycin and allowed to grow overnight. Next, 25 nl of each well was collected as part of a multiplexing strategy involving pooling of rows, columns, and plates, and these pools were subjected to PCR amplification and sequencing of the insertion barcodes. The location of each mutant was determined using an in-house script to identify the row, column, and plate pools sharing the same barcode. The insertion barcodes of the specific mutants used in this study were confirmed by Sanger sequencing of colonies recovered from glycerol stocks, using PCR with BarSeq_P1 and BarSeq_P2_noindex primers to amplify the region and BarSeq_P1_SangerSeq as the sequencing primer (Table 1). Furthermore, the presence of an insertion within the expected gene was confirmed by PCR amplification of the transposon junction using a mariner transposon-specific primer and a gene-specific primer (see the “confirm” primers in Table 1).

View this table:
  • View inline
  • View popup
TABLE 1

Primers used in this study

A transposon mutant library of Paraburkholderia bryophila 376MFSha3.1 (Burk376_ML3) was constructed via conjugation with an E. coli donor strain (APA752) carrying the randomly barcoded mariner transposon delivery vector pKMW3, as described previously (7). Briefly, we conjugated a 1:1 ratio of mid-log-phase P. byrophila to mid-log-phase APA752 for 12 h at 30°C on LB plates supplemented with diaminopimelic acid (DAP). The E. coli donor strain is derived from WM3064, a DAP auxotroph. After conjugation, we selected for mutants by plating on LB plates supplemented with 100 μg/ml kanamycin at 30°C. After pooling about 100,000 colonies, we subinoculated the cells into 50 ml of LB with 100 μg/ml kanamycin at a starting optical density at 600 nm (OD600) of 0.2, grew the library to saturation at 30°C, and made multiple freezer stocks in glycerol. The mapping of transposon insertion locations and the identification of their associated DNA barcodes were performed as described previously (7). Burk376_ML3 contains 106,195 mapped strains with unique DNA barcodes and contains at least one central insertion (in 10 to 90% of the gene’s length) for 5,624 of 6,555 predicted protein-coding genes. Of 63,242 strains with insertions in the central parts of genes, 59,045 were at sufficient abundance in the recovered library and were used to compute gene fitness.

All bacteria were cultured in LB medium during recovery from glycerol stocks. For the recovery of mutant libraries or individual mutant strains, 50 μg/ml kanamycin was added to the LB. All mutant libraries except K. michiganensis were recovered until the cells reached mid-log phase (OD600 of around 1.0). For K. michiganensis, we recovered the mutant library until the culture reached saturation (OD600 of around 4.0). All bacterial growth described in this study was conducted at 30°C. For cloning work, LB was used as the rich medium. All chemicals for medium preparations, mutant fitness assays, growth assays, and biochemical reactions were purchased from Sigma unless otherwise mentioned. For deoxyribonate, we used the lithium salt.

To assay fitness or growth on a specific carbon source, we generally used a defined medium with 5 to 20 mM of the carbon source as well as 30 mM PIPES [piperazine-N,N′-bis(2-ethanesulfonic acid)] sesquisodium salt, 0.25 g/liter ammonium chloride, 0.1 g/liter potassium chloride, 0.6 g/liter sodium phosphate monobasic monohydrate, Wolfe’s vitamins, and Wolfe’s minerals (http://fit.genomics.lbl.gov). For some fitness assays of P. putida, we instead used a minimal MOPS medium [40 mM 3-(N-morpholino)propanesulfonic acid, 4 mM Tricine, 1.32 mM potassium phosphate dibasic, 10 μM iron(II) sulfate heptahydrate, 9.5 mM ammonium chloride, 0.276 mM aluminum potassium sulfate dodecahydrate, 0.5 μM calcium chloride, 0.525 mM magnesium chloride hexahydrate, 50 mM sodium chloride, 3 nM ammonium heptamolybdate tetrahydrate, 0.4 μM boric acid, 30 nM cobalt chloride hexahydrate, 10 nM copper(II) sulfate pentahydrate, 80 nM manganese(II) chloride tetrahydrate, and 10 nM zinc sulfate heptahydrate].

For the growth assays with individual mutants of P. simiae (Fig. 2), we washed the cells 3 times with defined media before starting the growth experiments. (The cells were pregrown in LB.) We grew these individual strains in a 96-well microplate within a Tecan Sunrise instrument, with readings every 15 min, as described previously (6).

Genome-wide mutant fitness assays.For pooled fitness assays (RB-TnSeq), cells from the recovered mutant library were washed 3 times with defined media without a carbon source and then inoculated into the medium of choice at an OD600 of 0.02 and grown until saturation in a transparent 24-well microplate plate in a Multitron shaker. Total genomic DNA was extracted and barcodes were PCR amplified as described previously (7), but with some modifications to the primers. For most of the experiments, the P1 primer included variable spacing between the Illumina adapter and the sequence that flanks the barcode (2 to 5 N’s instead of always 5 N’s [Table 1]). Barcodes were sequenced using an Illumina HiSeq 4000, which can lead to some bleed-through between samples (in our experience, 0.1% to 0.3% of reads from BarSeq; for the likely mechanism behind bleed-through, see reference 29). To eliminate this bleed-through, some experiments used a P1 primer that contained an additional sequence to verify that it came from the expected sample. (The P2 primer contains the tag that is used for demultiplexing by Illumina software; in our earlier design, the P1 primer was the same mix for all the samples.) Specifically, the new P1 primers contain the same sequence as BarSeq_P1 before the N’s, 1 to 4 N’s (which varies with the index), the reverse (not the reverse complement) of the 6-nucleotide (nt) index sequence, and the sequence from BarSeq_P1 after the N’s.

Analysis of the mutant fitness data.The fitness data were analyzed as described previously (7) using publicly available scripts (https://bitbucket.org/berkeleylab/feba). Briefly, the fitness of a strain is the normalized log2 ratio of the number of reads for its barcode in the sample after growth versus the control sample (taken after the cells recover from the freezer). The fitness of a gene is the weighted average of the fitness of strains with insertions in the central 10 to 90% of the gene. For example, fitness values were computed for 4,414 of the 5,033 predicted nonessential proteins in P. simiae, and for the typical (median) protein, the fitness values of five independent mutant strains were averaged. For each experiment, the gene fitness values are normalized so that the mode of the distribution is at zero.

To assess the statistical significance of each fitness value, we used a t-like test statistic of the form fitness/sqrt(estimated variance) (7). For most of the bacteria, a gene was considered to have a specific phenotype in an experiment if fitness < −1, t < −5, 95% of fitness values are between −1 and 1, and |fitness| > 95th percentile(|fitness|) + 0.5, as described previously (6). However, for P. simiae, 9/162 experiments (6%) are for deoxyribose or deoxyribonate, so the 95th percentile threshold for a specific phenotype was not appropriate. We used the 90th percentile instead.

In P. simiae, we found 19 genes with a specific phenotype in at least one of six deoxyribose experiments. We manually excluded 6 of these genes from Fig. 1B because they did not have a consistent fitness defect. These genes were as follows: PS417_08675, or moaA, which is involved in molybdenum cofactor biosynthesis; PS417_23250, or tsaA, which is involved in tRNA modification; and four regulatory or signaling genes, PS417_00445, PS417_07615, PS417_13295, and PS417_22670.

Sequence analysis.Protein sequences were analyzed by using PaperBLAST (30) to search for characterized homologs, with the Fitness Browser (http://fit.genomics.lbl.gov [6]), which incorporates results from Pfam (31), TIGRFAMs (32), KEGG (33), and SEED/RAST (34), by using the Conserved Domain Database (35) to search for similarity to gene families and by using PSORTb (16) to predict protein localization.

To identify a motif for the putative binding sites of the transcription factor PS417_07240, we first used the MicrobesOnline tree-browser (36) to identify orthologs with conserved synteny. We identified similar proteins from P. fluorescens SBW25, P. fluorescens Pf-5, Pseudomonas syringae B728a, and B. graminis C4D1M (44 to 99% amino acid identity) that were also upstream of an SDR dehydrogenase and a β-keto acid cleavage enzyme. We ran MEME (37) on the 200 nucleotides upstream of these five genes.

Cloning and protein purification.We used the expression vector pET32a for cloning and protein purification as described previously (38). We PCR amplified four genes from P. simiae WCS417 (PS417_10885, PS417_10890, PS417_07255, and PS417_07245) using the amplify primers listed in Table 1 and cloned the PCR products into pET32a using Gibson assembly (39). For PS417_07255, PS417_07245, and PS417_10890, we transformed the Gibson assembly products directly into the E. coli expression strain BL21(DE3). However, for PS417_10885, we first transformed into E. coli TOP10 to identify the correct plasmid. We then transformed this plasmid into the E. coli BL21(DE3) expression strain. Protein was purified using a cobalt affinity column (Clonetech Talon kit) and quantified using a Nanodrop spectrophotometer (Thermo Fisher). Purified PS417_07245 and PS417_07255 had the expected molecular weights (29 kDa and 39 kDa, respectively). Although we were not able to purify PS417_10885 or PS417_10890, protein gels of cell lysates from E. coli strains that expressed these proteins showed predominant bands at the expected molecular weights (80 and 17 kDa, respectively).

Enzymatic assays.Enzymatic assays were performed with purified enzymes (2 to 5 μM), except that clarified cell lysates of E. coli expressing PS417_10885 or PS417_10890 were mixed together to obtain deoxyribose dehydrogenase. Each enzymatic assay was done in 50 μl of 100 mM ammonium bicarbonate (NH4HCO3) buffer with three independent replicates, run for 1 h at room temperature, and then diluted 1,000-fold in ammonium bicarbonate buffer and frozen at −20°C. NADH oxidase was purchased from Nacalai USA (San Diego, CA).

Liquid chromatography and mass spectrometry of end products.Samples were analyzed using a Thermo-Dionex UltiMate3000 RSLCnano liquid chromatograph that was connected in-line with an LTQ-Orbitrap-XL mass spectrometer equipped with a nanoelectrospray ionization (nano-ESI) source (Thermo Fisher Scientific, Waltham, MA). The liquid chromatograph was equipped with a C18 analytical column (length, 150 mm; inner diameter, 0.075 mm; particle size, 3 µm; pore size, 100 Å; Thermo Acclaim) and a 1-µl sample loop. Acetonitrile, formic acid (Optima grade, 99.9%; Fisher), and water purified to a resistivity of 18.2 MΩ·cm (at 25°C) using a Milli-Q Gradient ultrapure water purification system (Millipore, Billerica, MA) were used to prepare mobile-phase solvents. Solvent A was 99.9% water–0.1% formic acid, and solvent B was 99.9% acetonitrile–0.1% formic acid (vol/vol). The elution program consisted of isocratic flow at 1% B for 3 min, a linear gradient to 95% B over 1 min, isocratic flow at 95% B for 3 min, and isocratic flow at 1% B for 13 min, at a flow rate of 300 nl/min. Full-scan mass spectra were acquired in the positive-ion mode over the range m/z = 70 to 1,500 using the Orbitrap mass analyzer, in profile format, with a mass resolution setting of 60,000 (at m/z = 400, measured at full width at half-maximum peak height). For tandem mass spectrometry (MS/MS) analysis, precursor ions were fragmented using collision-induced dissociation (CID) under the following conditions: MS/MS spectra acquired using the linear ion trap, centroid format, isolation width of 2 m/z units, normalized collision energy of 35%, default charge state of 1+, activation Q of 0.25, and activation time of 30 ms. Data acquisition and analysis were performed using Xcalibur software (version 2.0.7; Thermo).

Proteomics.P. simiae was inoculated at an initial OD600 of 0.04 into 1 ml of defined medium that contained 10 mM either glucose, deoxyribose, or lithium deoxyribonate. At late log phase (OD600 of 0.4 to 0.7), cells were pelleted by centrifugation at 4,000 × g and resuspended in 100 mM NH4HCO3 (pH 7.5) buffer (Sigma-Aldrich, St. Louis, MO). Fifty microliters of a 1-μg/μl sample of each cell suspension was mixed with 10 μl 100 mM NH4HCO3 (pH 7.5) and 25 μl of a 0.2% aqueous solution of Rapigest SF (Waters, Milford, MA), and the sample was placed in a block heater at 80°C for 15 min and vortexed briefly to lyse cells and solubilize proteins. Five microliters of a 1-μg/μl solution of Trypsin Gold porcine protease (Promega, Madison, WI) was added, and the sample was incubated overnight at 37°C for the tryptic digest of proteins. The next morning, 10 μl of trifluoroacetic acid (Sigma-Aldrich) was added, and the sample was incubated for 90 min at 37°C to hydrolyze the Rapigest SF. The sample was centrifuged at 21,130 × g for 30 min at 6°C, and the supernatant was transferred to an autosampler vial for mass spectrometry analysis.

Trypsin-digested proteins were analyzed using a Synapt G2-Si high-definition ion mobility mass spectrometer that was equipped with a nanoelectrospray ionization source and connected in line with an Acquity M-class ultraperformance liquid chromatography (UPLC) system (Waters, Milford, MA). Data-independent, ion mobility-enabled mass spectra and tandem mass spectra (40–42) were acquired in the positive-ion mode. Data acquisition was controlled using MassLynx software (version 4.1). Tryptic peptide identification, relative protein quantification using a label-free approach (43, 44), and calculation of statistical analysis of variance (ANOVA) P values were performed using Progenesis QI for Proteomics software (version 4.0; Waters). The abundance of 1,277 proteins was quantified.

Data availability.The fitness data are available from the Fitness Browser (http://fit.genomics.lbl.gov). The fitness data and the proteomics data are archived at figshare (https://doi.org/10.6084/m9.figshare.7304159.v1). The 96 primer sequences described in the section “Genome-wide mutant fitness assays” are available with the analysis code (https://bitbucket.org/berkeleylab/feba) in the source file named primers/barseq3.index2.

APPENDIX

Oxidation of deoxyribose to a ketodeoxyribonate in vitro.To test the activity of deoxyribose dehydrogenase from P. simiae in vitro, we recombinantly expressed the iorA- and iorB-like components (PS417_10890 and PS147_10885) with polyhistidine (6× His) tags in E. coli. We tried to use a cobalt affinity column to purify each subunit, but did not succeed. Instead, we combined cell lysates from two strains of E. coli that overexpressed each subunit.

We incubated these cell lysates with 1 mM deoxyribose and with 100 μM phenazine methosulfate as the electron acceptor. The iorAB-like dehydrogenase converted deoxyribose to a compound with the same mass-to-charge ratio (m/z) as deoxyribonate (see Fig. S3 in the supplemental material). Tandem mass spectrometry confirmed that this product had a fragmentation pattern very similar to that of a 2-deoxy-d-ribonate standard (see Fig. S4 in the supplemental material). In contrast, addition of an equivalent amount of cell lysate from E. coli with an empty vector yielded only a small amount of deoxyribonate (less than a tenth as much [Fig. S3]). This shows that the deoxyribose dehydrogenase activity was primarily due to the iorAB-like dehydrogenase from P. simiae and not to an enzyme that is native to E. coli. No deoxyribonate formed if no lysate was added, so the conversion does not occur spontaneously in solution. Using polyhistidine tags, we were able to purify the lactonase (PS417_07255). The conversion of deoxyribose to deoxyribonate occurred with or without the added lactonase, perhaps because deoxyribonolactone in water can spontaneously hydrolyze to deoxyribonate or because lactonases are present in the cell lysate. In both reactions, and in the starting sample of deoxyribose, we detected a small amount of a compound with an m/z of 155.0 that could be deoxyribonolactone, although the abundance of this ion was too low for tandem mass spectrometry.

FIG S3

Reconstitution of the oxidation of deoxyribose to deoxyribonate in vitro. Deoxyribose, a potential deoxyribonolactone, and deoxyribonate were detected as sodium adducts (m/z = 157.05, 155.0, and 173.04, respectively). Each point is the average of 3 replicates, and the error bar is the standard deviation. Download FIG S3, PDF file, 0.01 MB.
Copyright © 2019 Price et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S4

The tandem mass spectra of the product of deoxyribose dehydrogenase and of deoxyribonate. Download FIG S4, PDF file, 0.1 MB.
Copyright © 2019 Price et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

We also purified the deoxyribonate dehydrogenase (PS417_07245) using a polyhistidine tag and a cobalt affinity column. We incubated the purified enzyme with 1 mM deoxyribonate and with 1 mM NAD+ as the electron acceptor. We observed a product with the molecular weight of 2-deoxy-3-keto-d-ribonate (a putative sodium adduct with an m/z of 171.02). The putative ketodeoxyribonate was not detected in the absence of the deoxyribonate dehydrogenase or when either deoxyribose or deoxyribonate was incubated with the iorAB-like dehydrogenase. Tandem mass spectrometry showed the loss of CO2 or H2O, as expected for a sugar acid (see Fig. S5 in the supplemental material). The amount of product was small and difficult to quantify (perhaps 3 to 4% of the ion intensity of the input deoxyribonate), but in the absence of the enzyme, the product was below the limit of detection. We also tested the deoxyribonate dehydrogenase for activity with deoxyribose as the substrate, but we did not observe any formation of deoxyribonate or ketodeoxyribonate.

FIG S5

Tandem mass spectrum (MS/MS) of the putative 2-deoxy-3-keto-d-ribonate from deoxyribonate incubated with purified deoxyribonate dehydrogenase. Download FIG S5, PDF file, 0.1 MB.
Copyright © 2019 Price et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Because the deoxyribonate dehydrogenase reaction is similar to the 2-deoxygluconate 3-dehydrogenase reaction, which is thermodynamically unfavorable (22), we tried to reoxidize the NADH in order to drive the reaction forward. Specifically, we incorporated 0.002 U/μl of NADH oxidase, which converts NADH + O2 to NAD+ + H2O2. This change increased the yield of the putative ketodeoxyribonate product to 12 to 14%. The retention time of the product by high-performance liquid chromatography (HPLC) was identical to that of deoxyribonate, so we did not try to separate them.

ACKNOWLEDGMENTS

We thank Mitchell Thompson for growing the P. putida mutant library on several carbon sources.

This material by ENIGMA (Ecosystems and Networks Integrated with Genes and Molecular Assemblies; http://enigma.lbl.gov), a Scientific Focus Area Program at Lawrence Berkeley National Laboratory, is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research, under contract no. DE-AC02-05CH11231. The QB3/Chemistry Mass Spectrometry Facility at UC Berkeley receives support from NIH (grant no. 1S10OD020062-01). This work used the Vincent J. Coates Genomics Sequencing Laboratory at UC Berkeley, supported by NIH Instrumentation Grant S10 OD018174.

FOOTNOTES

    • Received November 19, 2018.
    • Accepted January 12, 2019.
  • Copyright © 2019 Price et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Levy-Booth DJ,
    2. Campbell RG,
    3. Gulden RH,
    4. Hart MM,
    5. Powell JR,
    6. Klironomos JN,
    7. Pauls KP,
    8. Swanton CJ,
    9. Trevors JT,
    10. Dunfield KE
    . 2007. Cycling of extracellular DNA in the soil environment. Soil Biol Biochem 39:2977–2991. doi:10.1016/j.soilbio.2007.06.020.
    OpenUrlCrossRef
  2. 2.↵
    1. Rittmann D,
    2. Sorger-Herrmann U,
    3. Wendisch VF
    . 2005. Phosphate starvation-inducible gene ushA encodes a 5′ nucleotidase required for growth of Corynebacterium glutamicum on media with nucleotides as the phosphorus source. Appl Environ Microbiol 71:4339–4344. doi:10.1128/AEM.71.8.4339-4344.2005.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. Koszalka GW,
    2. Krenitsky TA
    . 1979. Nucleosidases from Leishmania donovani. Pyrimidine ribonucleosidase, purine ribonucleosidase, and a novel purine 2′-deoxyribonucleosidase. J Biol Chem 254:8185–8193.
    OpenUrlFREE Full Text
  4. 4.↵
    1. Tourneux L,
    2. Bucurenci N,
    3. Saveanu C,
    4. Kaminski PA,
    5. Bouzon M,
    6. Pistotnik E,
    7. Namane A,
    8. Marlière P,
    9. Bârzu O,
    10. Li De La Sierra I,
    11. Neuhard J,
    12. Gilles AM
    . 2000. Genetic and biochemical characterization of Salmonella enterica serovar Typhi deoxyribokinase. J Bacteriol 182:869–873.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Valentin-Hansen P,
    2. Boëtius F,
    3. Hammer-Jespersen K,
    4. Svendsen I
    . 1982. The primary structure of Escherichia coli K12 2-deoxyribose 5-phosphate aldolase. Nucleotide sequence of the deoC gene and the amino acid sequence of the enzyme. Eur J Biochem 125:561–566.
    OpenUrlPubMedWeb of Science
  6. 6.↵
    1. Price MN,
    2. Wetmore KM,
    3. Waters RJ,
    4. Callaghan M,
    5. Ray J,
    6. Liu H,
    7. Kuehl JV,
    8. Melnyk RA,
    9. Lamson JS,
    10. Suh Y,
    11. Carlson HK,
    12. Esquivel Z,
    13. Sadeeshkumar H,
    14. Chakraborty R,
    15. Zane GM,
    16. Rubin BE,
    17. Wall JD,
    18. Visel A,
    19. Bristow J,
    20. Blow MJ,
    21. Arkin AP,
    22. Deutschbauer AM
    . 2018. Mutant phenotypes for thousands of bacterial genes of unknown function. Nature 557:503–509. doi:10.1038/s41586-018-0124-0.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Wetmore KM,
    2. Price MN,
    3. Waters RJ,
    4. Lamson JS,
    5. He J,
    6. Hoover CA,
    7. Blow MJ,
    8. Bristow J,
    9. Butland G,
    10. Arkin AP,
    11. Deutschbauer A
    . 2015. Rapid quantification of mutant fitness in diverse bacteria by sequencing randomly bar-coded transposons. mBio 6:e00306-15. doi:10.1128/mBio.00306-15.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Chamberlin BA,
    2. Sweeley CC
    . 1987. Metabolic profiles of urinary organic acids recovered from absorbent filter paper. Clin Chem 33:572–576.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Kappen LS,
    2. Goldberg IH
    . 1989. Identification of 2-deoxyribonolactone at the site of neocarzinostatin-induced cytosine release in the sequence d(AGC). Biochemistry 28:1027–1032.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Lerma-Ortiz C,
    2. Jeffryes JG,
    3. Cooper AJL,
    4. Niehaus TD,
    5. Thamm AMK,
    6. Frelin O,
    7. Aunins T,
    8. Fiehn O,
    9. de Crecy-Lagard V,
    10. Henry CS,
    11. Hanson AD
    . 2016. “Nothing of chemistry disappears in biology”: the top 30 damage-prone endogenous metabolites. Biochem Soc Trans 44:961–971. doi:10.1042/BST20160073.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Lehmann M,
    2. Tshisuaka B,
    3. Fetzner S,
    4. Lingens F
    . 1995. Molecular cloning of the isoquinoline 1-oxidoreductase genes from Pseudomonas diminuta 7, structural analysis of iorA and iorB, and sequence comparisons with other molybdenum-containing hydroxylases. J Biol Chem 270:14420–14429.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Bastard K,
    2. Smith AA,
    3. Vergne-Vaxelaire C,
    4. Perret A,
    5. Zaparucha A,
    6. De Melo-Minardi R,
    7. Mariage A,
    8. Boutard M,
    9. Debard A,
    10. Lechaplais C,
    11. Pelle C,
    12. Pellouin V,
    13. Perchat N,
    14. Petit JL,
    15. Kreimeyer A,
    16. Medigue C,
    17. Weissenbach J,
    18. Artiguenave F,
    19. De Berardinis V,
    20. Vallenet D,
    21. Salanoubat M
    . 2014. Revealing the hidden functional diversity of an enzyme family. Nat Chem Biol 10:42–49. doi:10.1038/nchembio.1387.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Lemieux RU,
    2. Anderson L,
    3. Conner AH
    . 1971. The mutarotation of 2-deoxy-d-erythro-pentose (“2-deoxy-d-ribose”). Conformations, kinetics, and equilibria. Carbohydr Res 20:59–72.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Bartsch O,
    2. Hagemann M,
    3. Bauwe H
    . 2008. Only plant-type (GLYK) glycerate kinases produce d-glycerate 3-phosphate. FEBS Lett 582:3025–3028. doi:10.1016/j.febslet.2008.07.038.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Zelcbuch L,
    2. Razo-Mejia M,
    3. Herz E,
    4. Yahav S,
    5. Antonovsky N,
    6. Kroytoro H,
    7. Milo R,
    8. Bar-Even A
    . 2015. An in vivo metabolic approach for deciphering the product specificity of glycerate kinase proves that both E. coli’s glycerate kinases generate 2-phosphoglycerate. PLoS One 10:e0122957. doi:10.1371/journal.pone.0122957.
    OpenUrlCrossRef
  16. 16.↵
    1. Yu NY,
    2. Wagner JR,
    3. Laird MR,
    4. Melli G,
    5. Rey S,
    6. Lo R,
    7. Dao P,
    8. Sahinalp SC,
    9. Ester M,
    10. Foster LJ,
    11. Brinkman FS
    . 2010. PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics 26:1608–1615. doi:10.1093/bioinformatics/btq249.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Imam S,
    2. Noguera DR,
    3. Donohue TJ
    . 2015. CceR and AkgR regulate central carbon and energy metabolism in alphaproteobacteria. mBio 6:e02461-14. doi:10.1128/mBio.02461-14.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Novichkov PS,
    2. Kazakov AE,
    3. Ravcheev DA,
    4. Leyn SA,
    5. Kovaleva GY,
    6. Sutormin RA,
    7. Kazanov MD,
    8. Riehl W,
    9. Arkin AP,
    10. Dubchak I,
    11. Rodionov DA
    . 2013. RegPrecise 3.0—a resource for genome-scale exploration of transcriptional regulation in bacteria. BMC Genomics 14:745. doi:10.1186/1471-2164-14-745.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Dobbek H,
    2. Gremer L,
    3. Meyer O,
    4. Huber R
    . 1999. Crystal structure and mechanism of CO dehydrogenase, a molybdo iron-sulfur flavoprotein containing S-selanylcysteine. Proc Natl Acad Sci U S A 96:8884–8889.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Simon O,
    2. Klaiber I,
    3. Huber A,
    4. Pfannstiel J
    . 2014. Comprehensive proteome analysis of the response of Pseudomonas putida KT2440 to the flavor compound vanillin. J Proteomics 109:212–227. doi:10.1016/j.jprot.2014.07.006.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Graf N,
    2. Altenbuchner J
    . 2014. Genetic engineering of Pseudomonas putida KT2440 for rapid and high-yield production of vanillin from ferulic acid. Appl Microbiol Biotechnol 98:137–149. doi:10.1007/s00253-013-5303-1.
    OpenUrlCrossRef
  22. 22.↵
    1. Eichhorn MM,
    2. Cynkin MA
    . 1965. Microbial metabolism of 2-deoxyglucose; 2-deoxygluconic acid dehydrogenase. Biochemistry 4:159–165.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Watanabe S,
    2. Shimada N,
    3. Tajima K,
    4. Kodaki T,
    5. Makino K
    . 2006. Identification and characterization of l-arabonate dehydratase, l-2-keto-3-deoxyarabonate dehydratase, and l-arabinolactonase involved in an alternative pathway of l-arabinose metabolism. Novel evolutionary insight into sugar metabolism. J Biol Chem 281:33521–33536. doi:10.1074/jbc.M606727200.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Aro-Kärkkäinen N,
    2. Toivari M,
    3. Maaheimo H,
    4. Ylilauri M,
    5. Pentikäinen OT,
    6. Andberg M,
    7. Oja M,
    8. Penttilä M,
    9. Wiebe MG,
    10. Ruohonen L,
    11. Koivula A
    . 2014. l-Arabinose/d-galactose 1-dehydrogenase of Rhizobium leguminosarum bv. trifolii characterised and applied for bioconversion of l-arabinose to l-arabonate with Saccharomyces cerevisiae. Appl Microbiol Biotechnol 98:9653–9665. doi:10.1007/s00253-014-6039-2.
    OpenUrlCrossRef
  25. 25.↵
    1. Wattam AR,
    2. Davis JJ,
    3. Assaf R,
    4. Boisvert S,
    5. Brettin T,
    6. Bun C,
    7. Conrad N,
    8. Dietrich EM,
    9. Disz T,
    10. Gabbard JL,
    11. Gerdes S,
    12. Henry CS,
    13. Kenyon RW,
    14. Machi D,
    15. Mao C,
    16. Nordberg EK,
    17. Olsen GJ,
    18. Murphy-Olson DE,
    19. Olson R,
    20. Overbeek R,
    21. Parrello B,
    22. Pusch GD,
    23. Shukla M,
    24. Vonstein V,
    25. Warren A,
    26. Xia F,
    27. Yoo H,
    28. Stevens RL
    . 2017. Improvements to PATRIC, the all-bacterial Bioinformatics Database and Analysis Resource Center. Nucleic Acids Res 45:D535–D542. doi:10.1093/nar/gkw1017.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Landenmark HKE,
    2. Forgan DH,
    3. Cockell CS
    . 2015. An estimate of the total DNA in the biosphere. PLoS Biol 13:e1002168. doi:10.1371/journal.pbio.1002168.
    OpenUrlCrossRef
  27. 27.↵
    1. Rand JM,
    2. Pisithkul T,
    3. Clark RL,
    4. Thiede JM,
    5. Mehrer CR,
    6. Agnew DE,
    7. Campbell CE,
    8. Markley AL,
    9. Price MN,
    10. Ray J,
    11. Wetmore KM,
    12. Suh Y,
    13. Arkin AP,
    14. Deutschbauer AM,
    15. Amador-Noguez D,
    16. Pfleger BF
    . 2017. A metabolic pathway for catabolizing levulinic acid in bacteria. Nat Microbiol 2:1624–1634. doi:10.1038/s41564-017-0028-z.
    OpenUrlCrossRef
  28. 28.↵
    1. Cole BJ,
    2. Feltcher ME,
    3. Waters RJ,
    4. Wetmore KM,
    5. Mucyn TS,
    6. Ryan EM,
    7. Wang G,
    8. Ul-Hasan S,
    9. McDonald M,
    10. Yoshikuni Y,
    11. Malmstrom RR,
    12. Deutschbauer AM,
    13. Dangl JL,
    14. Visel A
    . 2017. Genome-wide identification of bacterial plant colonization genes. PLoS Biol 15:e2002860. doi:10.1371/journal.pbio.2002860.
    OpenUrlCrossRef
  29. 29.↵
    1. Sinha R,
    2. Stanley G,
    3. Gulati GS,
    4. Ezran C,
    5. Travaglini KJ,
    6. Wei E,
    7. Chan CKF,
    8. Nabhan AN,
    9. Su T,
    10. Morganti RM,
    11. Conley SD,
    12. Chaib H,
    13. Red-Horse K,
    14. Longaker MT,
    15. Snyder MP,
    16. Krasnow MA,
    17. Weissman IL
    . 2017. Index switching causes “spreading-of-signal” among multiplexed samples in Illumina HiSeq 4000 DNA sequencing. bioRxiv https://doi.org/10.1101/125724.
  30. 30.↵
    1. Price MN,
    2. Arkin AP
    . 2017. PaperBLAST: text mining papers for information about homologs. mSystems 2:e00039-17. doi:10.1128/mSystems.00039-17.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Finn RD,
    2. Bateman A,
    3. Clements J,
    4. Coggill P,
    5. Eberhardt RY,
    6. Eddy SR,
    7. Heger A,
    8. Hetherington K,
    9. Holm L,
    10. Mistry J,
    11. Sonnhammer EL,
    12. Tate J,
    13. Punta M
    . 2014. Pfam: the Protein Families Database. Nucleic Acids Res 42:D222–D230. doi:10.1093/nar/gkt1223.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Haft DH,
    2. Selengut JD,
    3. Richter RA,
    4. Harkins D,
    5. Basu MK,
    6. Beck E
    . 2012. Tigrfams and genome properties in 2013. Nucleic Acids Res 41:D387–D395. doi:10.1093/nar/gks1234.
    OpenUrlCrossRef
  33. 33.↵
    1. Kanehisa M,
    2. Sato Y,
    3. Kawashima M,
    4. Furumichi M,
    5. Tanabe M
    . 2016. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 44:D457–D462. doi:10.1093/nar/gkv1070.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Overbeek R,
    2. Olson R,
    3. Pusch GD,
    4. Olsen GJ,
    5. Davis JJ,
    6. Disz T,
    7. Edwards RA,
    8. Gerdes S,
    9. Parrello B,
    10. Shukla M,
    11. Vonstein V,
    12. Wattam AR,
    13. Xia F,
    14. Stevens R
    . 2014. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res 42:D206–D214. doi:10.1093/nar/gkt1226.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Marchler-Bauer A,
    2. Derbyshire MK,
    3. Gonzales NR,
    4. Lu S,
    5. Chitsaz F,
    6. Geer LY,
    7. Geer RC,
    8. He J,
    9. Gwadz M,
    10. Hurwitz DI,
    11. Lanczycki CJ,
    12. Lu F,
    13. Marchler GH,
    14. Song JS,
    15. Thanki N,
    16. Wang Z,
    17. Yamashita RA,
    18. Zhang D,
    19. Zheng C,
    20. Bryant SH
    . 2015. CDD: NCBI’s conserved domain database. Nucleic Acids Res 43:D222–D226. doi:10.1093/nar/gku1221.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Dehal PS,
    2. Joachimiak MP,
    3. Price MN,
    4. Bates JT,
    5. Baumohl JK,
    6. Chivian D,
    7. Friedland GD,
    8. Huang KH,
    9. Keller K,
    10. Novichkov PS,
    11. Dubchak IL,
    12. Alm EJ,
    13. Arkin AP
    . 2010. MicrobesOnline: an integrated portal for comparative and functional genomics. Nucleic Acids Res 38:D396–D400. doi:10.1093/nar/gkp919.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Bailey TL,
    2. Elkan C
    . 1994. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol 2:28–36.
    OpenUrlPubMed
  38. 38.↵
    1. Baran R,
    2. Bowen BP,
    3. Price MN,
    4. Arkin AP,
    5. Deutschbauer AM,
    6. Northen TR
    . 2013. Metabolic footprinting of mutant libraries to map metabolite utilization to genotype. ACS Chem Biol 8:189–199. doi:10.1021/cb300477w.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Gibson DG
    . 2011. Enzymatic assembly of overlapping DNA fragments. Methods Enzymol 498:349–361. doi:10.1016/B978-0-12-385120-8.00015-2.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    1. Plumb RS,
    2. Johnson KA,
    3. Rainville P,
    4. Smith BW,
    5. Wilson ID,
    6. Castro-Perez JM,
    7. Nicholson JK
    . 2006. UPLC/MS(E); a new approach for generating molecular fragment information for biomarker structure elucidation. Rapid Commun Mass Spectrom 20:1989–1994. doi:10.1002/rcm.2550.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Shliaha PV,
    2. Bond NJ,
    3. Gatto L,
    4. Lilley KS
    . 2013. Effects of traveling wave ion mobility separation on data independent acquisition in proteomics studies. J Proteome Res 12:2323–2339. doi:10.1021/pr300775k.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Distler U,
    2. Kuharev J,
    3. Navarro P,
    4. Levin Y,
    5. Schild H,
    6. Tenzer S
    . 2014. Drift time-specific collision energies enable deep-coverage data-independent acquisition proteomics. Nat Methods 11:167–170. doi:10.1038/nmeth.2767.
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.↵
    1. Neilson KA,
    2. Ali NA,
    3. Muralidharan S,
    4. Mirzaei M,
    5. Mariani M,
    6. Assadourian G,
    7. Lee A,
    8. van Sluyter SC,
    9. Haynes PA
    . 2011. Less label, more free: approaches in label-free quantitative mass spectrometry. Proteomics 11:535–553. doi:10.1002/pmic.201000553.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Nahnsen S,
    2. Bielow C,
    3. Reinert K,
    4. Kohlbacher O
    . 2013. Tools for label-free peptide quantification. Mol Cell Proteomics 12:549–556. doi:10.1074/mcp.R112.025163.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Oxidative Pathways of Deoxyribose and Deoxyribonate Catabolism
Morgan N. Price, Jayashree Ray, Anthony T. Iavarone, Hans K. Carlson, Elizabeth M. Ryan, Rex R. Malmstrom, Adam P. Arkin, Adam M. Deutschbauer
mSystems Feb 2019, 4 (1) e00297-18; DOI: 10.1128/mSystems.00297-18

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Oxidative Pathways of Deoxyribose and Deoxyribonate Catabolism
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
Share
Oxidative Pathways of Deoxyribose and Deoxyribonate Catabolism
Morgan N. Price, Jayashree Ray, Anthony T. Iavarone, Hans K. Carlson, Elizabeth M. Ryan, Rex R. Malmstrom, Adam P. Arkin, Adam M. Deutschbauer
mSystems Feb 2019, 4 (1) e00297-18; DOI: 10.1128/mSystems.00297-18
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • APPENDIX
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

deoxyribonate catabolism
deoxyribose catabolism
high-throughput genetics

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2019 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077