One of the most abundant sources of organic carbon in the ocean is glycolate, the secretion of which by marine phytoplankton results in an estimated annual flux of one petagram of glycolate in marine environments1. Although it is generally accepted that glycolate is oxidized to glyoxylate by marine bacteria2,3,4, the further fate of this C2 metabolite is not well understood. Here we show that ubiquitous marine Proteobacteria are able to assimilate glyoxylate via the β-hydroxyaspartate cycle (BHAC) that was originally proposed 56 years ago5. We elucidate the biochemistry of the BHAC and describe the structure of its key enzymes, including a previously unknown primary imine reductase. Overall, the BHAC enables the direct production of oxaloacetate from glyoxylate through only four enzymatic steps, representing—to our knowledge—the most efficient glyoxylate assimilation route described to date. Analysis of marine metagenomes shows that the BHAC is globally distributed and on average 20-fold more abundant than the glycerate pathway, the only other known pathway for net glyoxylate assimilation. In a field study of a phytoplankton bloom, we show that glycolate is present in high nanomolar concentrations and taken up by prokaryotes at rates that allow a full turnover of the glycolate pool within one week. During the bloom, genes that encode BHAC key enzymes are present in up to 1.5% of the bacterial community and actively transcribed, supporting the role of the BHAC in glycolate assimilation and suggesting a previously undescribed trophic interaction between autotrophic phytoplankton and heterotrophic bacterioplankton.
This is a preview of subscription content
Subscription info for Chinese customers
We have a dedicated website for our Chinese customers. Please go to naturechina.com to subscribe to this journal.
Rent or Buy article
Get time limited or full article access on ReadCube.
All prices are NET prices.
The coordinates and structure factors of the crystal structures generated from this research are available at the PDB under accession numbers 6QKB and 6RQA. Mass spectrometry proteomics data are available via ProteomeXchange with the identifier PXD013274. MAGs are available under accession PRJEB28156 at the European Nucleotide Archive (ENA). All other relevant data are available in the Article and the Supplementary Information. Source Data for Figs. 2, 3 and Extended Data Fig. 1, 4, 5, 7–9 are provided with the paper.
Wright, R. T. & Shah, N. M. Trophic role of glycolic acid in coastal seawater. II. Seasonal changes in concentration and heterotrophic use in Ipswich Bay, Massachusetts, USA. Mar. Biol. 43, 257–263 (1977).
Wright, R. T. & Shah, N. M. Trophic role of glycolic acid in coastal seawater. I. Heterotrophic metabolism in seawater and bacterial cultures. Mar. Biol. 33, 175–183 (1975).
Fogg, G. E. The ecological significance of extracellular products of phytoplankton photosynthesis. Bot. Mar. 26, 3–14 (1983).
Lau, W. W., Keil, R. G. & Armbrust, E. V. Succession and diel transcriptional response of the glycolate-utilizing component of the bacterial community during a spring phytoplankton bloom. Appl. Environ. Microbiol. 73, 2440–2450 (2007).
Kornberg, H. L. & Morris, J. G. β-Hydroxyaspartate pathway: a new route for biosyntheses from glyoxylate. Nature 197, 456–457 (1963).
Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science 281, 237–240 (1998).
Duarte, C. M. & Cebrian, J. The fate of marine autotrophic production. Limnol. Oceanogr. 41, 1758–1766 (1996).
Hellebust, J. A. Excretion of some organic compounds by marine phytoplankton. Limnol. Oceanogr. 10, 192–206 (1965).
Tolbert, N. E. & Zill, L. P. Excretion of glycolic acid by algae during photosynthesis. J. Biol. Chem. 222, 895–906 (1956).
Leboulanger, C., Descolasgros, C. & Jupin, H. HPLC determination of glycolic acid in seawater. An estimation of phytoplankton photorespiration in the Gulf of Lions, western Mediterranean Sea. J. Plankton Res. 16, 897–903 (1994).
Leboulanger, C., Oriol, L., Jupin, H. & Descolas-Gros, C. Diel variability of glycolate in the eastern tropical Atlantic Ocean. Deep Sea Res. Part I Oceanogr. Res. Pap. 44, 2131–2139 (1997).
Casey, J. R., Ferrón, S. & Karl, D. M. Light-enhanced microbial organic carbon yield. Front. Microbiol. 8, 2157 (2017).
Lau, W. W. & Armbrust, E. V. Detection of glycolate oxidase gene glcD diversity among cultured and environmental marine bacteria. Environ. Microbiol. 8, 1688–1702 (2006).
Carini, P., Steindler, L., Beszteri, S. & Giovannoni, S. J. Nutrient requirements for growth of the extreme oligotroph ‘Candidatus Pelagibacter ubique’ HTCC1062 on a defined medium. ISME J. 7, 592–602 (2013).
Eiler, A. et al. Tuning fresh: radiation through rewiring of central metabolism in streamlined bacteria. ISME J. 10, 1902–1914 (2016).
Tripp, H. J. et al. Unique glycine-activated riboswitch linked to glycine–serine auxotrophy in SAR11. Environ. Microbiol. 11, 230–238 (2009).
Krakow, G. & Barkulis, S. S. Conversion of glyoxylate to hydroxypyruvate by extracts of Escherichia coli. Biochim. Biophys. Acta 21, 593–594 (1956).
Hansen, R. W. & Hayashi, J. A. Glycolate metabolism in Escherichia coli. J. Bacteriol. 83, 679–687 (1962).
Kornberg, H. L. & Morris, J. G. The utilization of glycollate by Micrococcus denitrificans: the β-hydroxyaspartate pathway. Biochem. J. 95, 577–586 (1965).
Gibbs, R. G. & Morris, J. G. Assay and properties of β-hydroxyaspartate aldolase from Micrococcus denitrificans. Biochim. Biophys. Acta 85, 501–503 (1964).
Gibbs, R. G. & Morris, J. G. Purification and properties of erythro-β-hydroxyaspartate dehydratase from Micrococcus denitrificans. Biochem. J. 97, 547–554 (1965).
Liu, J. Q., Dairi, T., Itoh, N., Kataoka, M. & Shimizu, S. A novel enzyme, d-3-hydroxyaspartate aldolase from Paracoccus denitrificans IFO 13301: purification, characterization, and gene cloning. Appl. Microbiol. Biotechnol. 62, 53–60 (2003).
Mortarino, M. et al. l-Aspartate oxidase from Escherichia coli. I. Characterization of coenzyme binding and product inhibition. Eur. J. Biochem. 239, 418–426 (1996).
Lenz, M., Borlinghaus, N., Weinmann, L. & Nestl, B. M. Recent advances in imine reductase-catalyzed reactions. World J. Microbiol. Biotechnol. 33, 199 (2017).
Hochreiter, M. C. & Schellenberg, K. A. α-Iminoglutarate formation by beef liver l-glutamate dehydrogenase. Detection by borohydride or dithionite reduction to glutamate. J. Am. Chem. Soc. 91, 6530–6531 (1969).
Hallen, A., Cooper, A. J., Smith, J. R., Jamie, J. F. & Karuso, P. Ketimine reductase/CRYM catalyzes reductive alkylamination of α-keto acids, confirming its function as an imine reductase. Amino Acids 47, 2457–2461 (2015).
Buchan, A., LeCleir, G. R., Gulvik, C. A. & González, J. M. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698 (2014).
Teeling, H. et al. Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms. eLife 5, e11888 (2016).
Teeling, H. et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 336, 608–611 (2012).
Shah, N. M. & Wright, R. T. Occurrence of glycolic acid in coastal sea-water. Mar. Biol. 24, 121–124 (1974).
Fogg, G. E., Burton, N. F. & Coughlan, S. J. The occurrence of glycollic acid in Antarctic waters. Br. Antarct. Surv. Bull. 41 & 42, 193–195 (1975).
Hasan-Al, R. H., Coughlan, S. J., Pant, A. & Fogg, G. E. Seasonal variations in phytoplankton and glycollate concentrations in the Menai Straits, Anglesey. J. Mar. Biol. Assoc. U.K. 55, 557–565 (1975).
Edenborn, H. M. & Litchfield, C. D. Glycolate turnover in the water column of the New York Bight apex. Mar. Biol. 95, 459–467 (1987).
Leboulanger, C., Serve, L., Comellas, L. & Jupin, H. Determination of glycolic acid released from marine phytoplankton by post-derivatization gas chromatography mass spectrometry. Phytochem. Anal. 9, 5–9 (1998).
Lord, J. M., Codd, G. A. & Merrett, M. J. The effect of light quality on glycolate formation and excretion in algae. Plant Physiol. 46, 855–856 (1970).
Smith, W. O. Extracellular release of glycolic acid by a marine diatom. J. Phycol. 10, 30–33 (1974).
Leboulanger, C., Martin-Jezequel, V., Descolas-Gros, C., Sciandra, A. & Jupin, H. J. Photorespiration in continuous culture of Dunaliella tertiolecta (Chlorophyta): relationships between serine, glycine, and extracellular glycolate. J. Phycol. 34, 651–654 (1998).
Schnitzler Parker, M., Armbrust, E. V., Piovia-Scott, J. & Keil, R. G. Induction of photorespiration by light in the centric diatom Thalassiosira weissflogii (Bacillariophyceae): molecular characterization and physiological consequences. J. Phycol. 40, 557–567 (2004).
Bertilsson, S., Berglund, O., Pullin, M. J. & Chisholm, S. W. Release of dissolved organic matter by Prochlorococcus. Vie Milieu 55, 225–231 (2005).
Simó, R. & Pedrós-Alió, C. Short-term variability in the open ocean cycle of dimethylsulfide. Glob. Biogeochem. Cycles 13, 1173–1181 (1999).
Bertani, G. Studies on lysogenesis. I. The mode of phage liberation by lysogenic Escherichia coli. J. Bacteriol. 62, 293–300 (1951).
Beijerinck, M. W. & Minkman, D. C. J. Bildung und Verbrauch von Stickoxydul durch Bakterien. Zentralbl. Bakteriol. Naturwiss. 25, 30–63 (1910).
Hahnke, S. M., Moosmann, P., Erb, T. J. & Strous, M. An improved medium for the anaerobic growth of Paracoccus denitrificans Pd1222. Front. Microbiol. 5, 18 (2014).
Laemmli, U. K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 227, 680–685 (1970).
Kitagawa, M. et al. Complete set of ORF clones of Escherichia coli ASKA library (a complete set of E. coli K-12 ORF archive): unique resources for biological research. DNA Res. 12, 291–299 (2005).
Lane, C. F. Sodium cyanoborohydride — a highly selective reducing agent for organic functional groups. Synthesis 1975, 135–146 (1975).
Kiefer, P., Schmitt, U. & Vorholt, J. A. eMZed: an open source framework in Python for rapid and interactive development of LC/MS data analysis workflows. Bioinformatics 29, 963–964 (2013).
Bradford, M. M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254 (1976).
Ledermann, R., Strebel, S., Kampik, C. & Fischer, H. M. Versatile vectors for efficient mutagenesis of Bradyrhizobium diazoefficiens and other Alphaproteobacteria. Appl. Environ. Microbiol. 82, 2791–2799 (2016).
Thoma, S. & Schobert, M. An improved Escherichia coli donor strain for diparental mating. FEMS Microbiol. Lett. 294, 127–132 (2009).
Moggridge, S., Sorensen, P. H., Morin, G. B. & Hughes, C. S. Extending the compatibility of the SP3 paramagnetic bead processing approach for proteomics. J. Proteome Res. 17, 1730–1740 (2018).
Glatter, T. et al. Large-scale quantitative assessment of different in-solution protein digestion protocols reveals superior cleavage efficiency of tandem Lys-C/trypsin proteolysis over trypsin digestion. J. Proteome Res. 11, 5145–5156 (2012).
Engilberge, S. et al. Crystallophore: a versatile lanthanide complex for protein crystallography combining nucleating effects, phasing properties, and luminescence. Chem. Sci. 8, 5909–5917 (2017).
Kabsch, W. XDS. Acta Crystallogr. D 66, 125–132 (2010).
Winn, M. D. et al. Overview of the CCP4 suite and current developments. Acta Crystallogr. D 67, 235–242 (2011).
Uhl, M. K. et al. The crystal structure of d-threonine aldolase from Alcaligenes xylosoxidans provides insight into a metal ion assisted PLP-dependent mechanism. PLoS ONE 10, e0124056 (2015).
Gallagher, D. T. et al. Structure of alanine dehydrogenase from Archaeoglobus: active site analysis and relation to bacterial cyclodeaminases and mammalian mu crystallin. J. Mol. Biol. 342, 119–130 (2004).
Waterhouse, A. et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 46, W296–W303 (2018).
Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. D 66, 213–221 (2010).
Emsley, P. & Cowtan, K. Coot: model-building tools for molecular graphics. Acta Crystallogr. D 60, 2126–2132 (2004).
Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60 (2015).
Francis, T. B., Krüger, K., Fuchs, B. M., Teeling, H. & Amann, R. I. Candidatus Prosiliicoccus vernus, a spring phytoplankton bloom associated member of the Flavobacteriaceae. Syst. Appl. Microbiol. 42, 41–53 (2019).
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055 (2015).
Parks, D. H. et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat. Biotechnol. 36, 996–1004 (2018).
Letunic, I. & Bork, P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245 (2016).
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004).
Silvestro, D. & Michalak, I. raxmlGUI: a graphical front-end for RAxML. Org. Divers. Evol. 12, 335–337 (2012).
Le, S. Q. & Gascuel, O. An improved general amino acid replacement matrix. Mol. Biol. Evol. 25, 1307–1320 (2008).
Gerlt, J. A. et al. Enzyme Function Initiative-Enzyme Similarity Tool (EFI-EST): A web tool for generating protein sequence similarity networks. Biochim. Biophys. Acta 1854, 1019–1037 (2015).
Villar, E. et al. The Ocean Gene Atlas: exploring the biogeography of plankton genes online. Nucleic Acids Res. 46, W289–W295 (2018).
Kaplun, A. et al. Glyoxylate carboligase lacks the canonical active site glutamate of thiamine-dependent enzymes. Nat. Chem. Biol. 4, 113–118 (2008).
Han, J., Gagnon, S., Eckle, T. & Borchers, C. H. Metabolomic analysis of key central carbon metabolism carboxylic acids as their 3-nitrophenylhydrazones by UPLC/ESI-MS. Electrophoresis 34, 2891–2900 (2013).
Rose, T. M. et al. Consensus-degenerate hybrid oligonucleotide primers for amplification of distantly related sequences. Nucleic Acids Res. 26, 1628–1635 (1998).
Rose, T. M., Henikoff, J. G. & Henikoff, S. CODEHOP (consensus-degenerate hybrid oligonucleotide primer) PCR primer design. Nucleic Acids Res. 31, 3763–3766 (2003).
Boyce, R., Chilana, P. & Rose, T. M. iCODEHOP: a new interactive program for designing consensus-degenerate hybrid oligonucleotide primers from multiply aligned protein sequences. Nucleic Acids Res. 37, W222–W228 (2009).
We thank L. Franzmeyer and K.-P. Rücknagel as well as the crew of the Aade and the staff at the Biological Station at Helgoland (BAH) for sample collection and processing; K. H. Wiltshire for providing Chl a data collected at BAH; S. Vidal-Melgosa, A. Bolte and J.-H. Hehemann for sharing samples; T. Ferdelman for help with radioactive tracer work; B. Vögeli, T. Schwander, G. Stoffel and S. Burgener for helpful discussions. We acknowledge the support from the staff scientists at the European Synchrotron Radiation Facility Grenoble, France (ESRF, beamlines ID29 & ID30B) as well as at the Deutsches Elektronen-Synchrotron Hamburg, Germany (DESY, beamline P13). Metagenome sequences were obtained within the COGITO (Coastal Microbe Genomic and Taxonomic Observatory) project granted to H. Teeling (MPI Bremen) as a community sequencing project by the Department of Energy’s Joint Genome Institute in Walnut Creek, CA, USA (Proposal ID 998; https://doi.org/10.25585/1488076). This study was funded by the Max-Planck-Society (R.I.A. and T.J.E.), FET-Open Grant 686330 (Future Agriculture) and the German Research Foundation (SFB987 ‘Microbial diversity in environmental signal response’ and FOR 2406 ‘Proteogenomics of marine polysaccharide utilization’).
The Max-Planck-Gesellschaft zur Förderung der Wissenschaften is the patent applicant for the following three patents. All patent applications are pending. L.S.v.B. and T.J.E. have filed European patent no. EP 19190404.4 for the production of plants with altered photorespiration due to implementation of the BHAC. L.S.v.B., J.Z. and T.J.E. have filed European patent no. EP 18167406.0 for the production of photoautotrophic organisms with altered photorespiration due to implementation of the BHAC. L.S.v.B. and T.J.E. have filed European patent no. 18211454.6 for the enantioselective preparation of primary amine compounds using the enzyme BhcD or its homologues.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Previously reported glycolate concentrations in environmental samples and cultures of photosynthetic organisms.
a, Bar diagram of glycolate concentrations as previously reported in environmental samples. For details on samples, replicates, and analytics see b and the literature cited therein. b, Table of glycolate concentrations as previously reported in environmental samples (E1, E2 and so on) and cultures of photosynthetic organisms (C1, C2 and so on). When reported in the reference1,2,4,9,10,11,30,31,32,33,34,35,36,37,38,39, the mean value ± error is given.
Extended Data Fig. 2 Crystal structure and phylogenetic analysis of the β-hydroxyaspartate aldolase BhcC.
a, Cartoon representation of the β-hydroxyaspartate aldolase homodimer (PDB 6QKB) with superimposed protein surface (left, side view; right, top view). b, Active site of β-hydroxyaspartate aldolase with covalently bound PLP (light cyan). Active site residues highlighted in pink (A160, A195 and S313) are completely conserved only among β-hydroxyaspartate aldolases, but differ in d-threonine aldolases. c, Active site of d-threonine aldolase (PDB 4V15). The corresponding conserved residues among d-threonine aldolases (Q155, S190 and C303) are highlighted as in b. d, Maximum likelihood phylogenetic tree of the type III PLP-dependent protein superfamily. Sequences of the β-hydroxyaspartate aldolase BhcC and its homologues form a distinct clade (blue) within the d-threonine aldolase branch of this superfamily. Bootstrap values of at least 50 are given on the respective nodes.
Extended Data Fig. 3 Crystal structure and phylogenetic analysis of the iminosuccinate reductase BhcD.
a, Cartoon representation of the iminosuccinate reductase homodimer (PDB 6RQA) with superimposed protein surface (left, side view; right, top view). b, Active site of BhcD with bound NAD+ (light cyan). Residues highlighted in pink (V39, R41, G52, K54 and H83) may contribute to substrate binding and are conserved among iminosuccinate reductases, but differ in l-alanine dehydrogenases. c, Active site of l-alanine dehydrogenase (PDB 1OMO). The corresponding conserved residues among l-alanine dehydrogenases (K41, Y43, R52, M54 and V81) are highlighted as in b. d, Maximum likelihood phylogenetic tree of the ornithine cyclodeaminase/µ-crystalline protein superfamily. Sequences of the iminosuccinate reductase BhcD and its homologues form a distinct clade (red) within this superfamily. Bootstrap values of at least 50 are given on the respective nodes. e, Sequence similarity network of 1,614 sequences from the ornithine cyclodeaminase/µ-crystalline protein superfamily. Connected sequences with more than 80% identity are clustered into nodes. The number in each node gives the number of sequences contained within. Nodes with more than 50% identity are connected by edges. Similar to the phylogenetic analysis shown in d, sequences of the iminosuccinate reductase BhcD and its homologues form a distinct clade (red) within this superfamily.
a, Michaelis–Menten kinetics for aspartate–glyoxylate aminotransferase (BhcA). b, Michaelis–Menten kinetics for β-hydroxyaspartate dehydratase (BhcB). c, Michaelis–Menten kinetics for β-hydroxyaspartate aldolase (BhcC). d, Michaelis–Menten kinetics for iminosuccinate reductase (BhcD). a–d, Data are shown from n = 3 independent experiments at different substrate concentrations. The data are summarized in Table 1.
a, Growth rate of wild-type P. denitrificans DSM 413 on the BHAC substrates glycolate and glyoxylate. The middle line and box are the median and interquartile range of n = 6 independent experiments and the whiskers indicate the maximum range of the dataset. b, c, Representative growth curves of wild-type P. denitrificans DSM 413 (grey) and bhc deletion strains (coloured) grown in the presence of 60 mM glycolate (b) or 60 mM glyoxylate (c). Deletion of any single gene in the bhc gene cluster is sufficient to completely abolish growth in the presence of glycolate and glyoxylate. These experiments were repeated three times independently with similar results. d–f, Growth rates (μ) of wild-type P. denitrificans DSM 413 (grey) and BHAC deletion strains (coloured) grown in the presence of 60 mM acetate (d), 30 mM succinate (e) or 20 mM glucose (f). Deletion of any single gene in the bhc gene cluster, or of the whole bhc gene cluster, still permits growth on acetate, succinate or glucose with comparable growth rates as for the wild type. Data are the mean ± s.d. of n = 3 independently grown cultures. g, Analysis of the proteome of glycolate-grown compared to succinate-grown P. denitrificans DSM 413. All proteins that were quantified by at least three unique peptides are shown. The 15 proteins that showed the strongest increase in abundance are marked in the volcano plot. The four enzymes of the BHAC are marked in red, the three subunits of glycolate oxidase in orange, the proteins of a putative operon for lactate utilization in white and the proteins directly downstream of the bhc gene cluster in light red. h, The abundance of these proteins, given as the percentage of the intensity-based absolute quantification (iBAQ) value. Data are the mean ± s.d. of n = 4 independently grown cultures. i, Specific activities of BHAC enzymes in cell-free extracts of glycolate-grown P. denitrificans DSM 413, as measured spectrophotometrically. Note that the activity of BhcD is plotted on the right y axis and consists of the actual iminosuccinate reductase activity (iminosuccinate to l-aspartate) as well as endogenous malate dehydrogenase activity (oxaloacetate to l-malate). j, Ratio of malate to aspartate determined by LC–MS during the enzyme assay for BhcD activity. The ratio remains approximately constant at 12:1, indicating that only approximately 8% of the activity (around 1.3 U mg−1) shown in i can be ascribed to iminosuccinate reductase. i, j, Data are the mean ± s.d. of n = 3 independently grown cultures; each data point represents the mean of n = 3 technical replicates. k, DNA-binding properties of BhcR. Left, a fluorescently labelled DNA fragment carrying the putative promoter region of the bhc gene cluster (Pbhc) was incubated with increasing amounts of purified BhcR protein and subsequently separated by electrophoresis to visualize DNA bound to BhcR and free DNA; a DNA fragment derived from the coding region of bhcA was used as a negative control. BhcR specifically forms a complex with the DNA fragment containing the putative promoter region of the bhc gene cluster. Right, the Pbhc–BhcR complex was incubated with increasing concentrations of glyoxylate and subsequently separated by electrophoresis to assess the effect of glyoxylate on complex formation; the bhcA DNA fragment together with BhcR was used as a negative control. Increasing concentrations of glyoxylate decrease the binding of BhcR to the Pbhc DNA fragment. For gel source data, see Supplementary Fig. 1.
a, Genome-based maximum likelihood phylogenetic tree of bacterial strains with the bhc gene cluster. The bhc gene cluster is found in Gammaproteobacteria (green), and in the alphaproteobacterial orders Rhizobiales (blue) and Rhodobacterales (red), as well as in one member each of Sphingomonadales and Kiloniellales. The phylogenetic tree is based on an alignment of 120 bacterial marker genes from 264 publicly available bacterial genomes and 5 MAGs and was calculated using GTDB-Tk64 (https://github.com/Ecogenomics/GtdbTk). If several strains from the same genus cluster together, nodes are collapsed at the genus level, and the size of the resulting circle corresponds to the respective number of strains. Loktanella*: collapsed node contains the MAGs 20110516_Bin_8_1 and 20110523_Bin_9_1; Planktotalea**: collapsed node contains the MAG 20110523_Bin_97_1; Litoricola***: collapsed node contains the MAG 20110526_Bin_19_1. b, Maximum likelihood phylogenetic tree of concatenated BHAC enzyme sequences. Colour code is the same as in a. Phylogenetic groups that were mostly isolated from terrestrial or freshwater habitats are marked with a black dot. Comparison with a reveals that the sequences of the BHAC enzymes are not phylogenetically representative, as, for example, alpha- and gammaproteobacterial sequences form a common branch and sequences from terrestrial or freshwater Rhizobiales and Rhodobacterales form another common branch. This suggests that the bhc gene cluster might have been subject to horizontal gene transfer between distantly related strains in shared habitats. The environmental bhc gene cluster sequence that could not be binned successfully is marked in bold and clusters together with isolated representatives of Pseudoruegeria, Litoreibacter and Pseudooceanicola. The phylogenetic tree is based on the concatenated alignments of the 4 enzymes (BhcA–BhcD) from 264 publicly available bacterial genomes and from 6 metagenome contigs. It was calculated using raxmlGUI67. Bootstrap values of at least 50 are given on the respective nodes; calculated branch lengths of the tree are ignored for the sake of better visualization. If several strains from the same genus cluster together, nodes are collapsed at the genus level, and the size of the resulting circle corresponds to the respective number of strains. If strains from the same genus cluster in more than one node, the respective branches are labelled as Genus_1, Genus_2, and so on, in a clockwise manner. Loktanella_2*: collapsed node contains the MAGs 20110516_Bin_8_1 and 20110523_Bin_9_1; Planktotalea**: collapsed node contains the MAG 20110523_Bin_97_1; Litoricola***: collapsed node contains the MAG 20110526_Bin_19_1. a, b, Taxonomy is based on GTDB (release 03-RS86; http://gtdb.ecogenomic.org/). All strains contained in the phylogenetic trees are listed in Supplementary Data 1.
a, Metagenomes collected during the Tara Oceans expedition were searched for the presence of BhcC as representative enzyme of the BHAC. Dots on the map mark sampling locations of metagenomes containing BhcC sequences; the colour of the dot corresponds to BhcC abundance in samples from surface water (0.22–3-µm size fraction), as shown in the legend. The map was made with Ocean Data View 5.1.5 (Schlitzer, R., Ocean Data View, odv.awi.de, 2018). b, Phylogenetic distribution of 104 unique BhcC sequences found in Tara Oceans metagenomes. c, Phylogenetic distribution of 32 unique Gcl (as representative enzyme of the glycerate pathway) sequences found in Tara Oceans metagenomes. Whereas BhcC is mainly found in Alphaproteobacteria, Gcl is largely restricted to Gammaproteobacteria. b, c, Taxonomy is based on GTDB64 (release 03-RS86; http://gtdb.ecogenomic.org/). d, Ratio of the abundances (in percentage of total reads) of BhcC to Gcl in Tara Oceans metagenomes. BhcC:Gcl ratios from n = 210 samples are plotted together (left) and clustered by sampling depth (SRF, upper layer zone (n = 101); DCM, deep chlorophyll maximum layer (n = 68); MES, mesopelagic zone (n = 41)). Samples from the 0.22–3-µm size fraction are denoted by a black dot; samples from the <0.22-µm size fraction are denoted by a blue dot. The median is shown in orange as centre value and error bars represent interquartile ranges. The median BhcC:Gcl ratio of all samples is 18.7. The highest BhcC:Gcl ratio is found in surface water samples (median = 41.8), with the ratio generally being higher in the 0.22–3-µm size fraction than in the <0.22-µm size fraction. Sequence IDs, abundances and BhcC:Gcl ratios are given in Supplementary Data 2.
a, The location of Helgoland Island approximately 40 km offshore the northern German coastline in the North Sea is marked with a red dot. The map was made with Ocean Data View 5.1.5 (R. Schlitzer, Ocean Data View, odv.awi.de, 2018). b, The long-term ecological research site ‘Kabeltonne’ (red dot: 54° 11.3′ N, 7° 54.0′ E) is located between Helgoland Island (left) and the small island Düne (right). Satellite image from WorldWind Explorer (B. Schubert, worldwind.earth/explorer, 2016–2018); the image was adapted to indicate the sampling site. c–e, Abundance of the bhc gene cluster (in RPKM) was calculated in 38 metagenomes from samples collected during the algal spring blooms of 2010 to 2012 in the North Sea close to Helgoland28. Six different sequences were investigated, five of which could be assigned to metagenome bins (Extended Data Fig. 6 and Supplementary Data 3), whereas the remaining, most abundant sequence (black) could not be binned successfully.
Degenerate primers for bhcC were used for qPCR with different amounts of genomic DNA from P. denitrificans DSM 413, Rhodobacter sphaeroides 2.4.1, and E. coli K-12 MG1655 as template. While the bhcC gene from P. denitrificans DSM 413 is amplified, genomic DNA from organisms that lack the bhc gene cluster does not result in reliable amplification. Data are mean ± s.d.; n = 3 independent experiments.
This file contains a figure showing an SDS-PAGE gel of all purified proteins of the bhc gene cluster, as well as uncropped gel data scans of this figure and Extended Data Figure 5k.
This file contains Supplementary Table 1 (an overview of known natural net glyoxylate assimilation pathways), Supplementary Table 2 (a list of the strains and plasmids used in this study) and Supplementary Table 3 (a list of the primers used in this study).
This dataset contains a list of bacterial isolates with the bhc gene cluster, as well as their taxonomic classification and the Gene IDs for their bhcC and glcD genes.
This dataset contains the Gene IDs, taxonomic classifications, abundances, and ratios of BhcC and Gcl sequences found in Tara Oceans metagenomes.
This dataset contains information on the quality of Helgoland metagenome bins containing the bhc gene cluster.
This dataset contains information on the yields of DNA and RNA extractions from North Sea water samples.
About this article
Cite this article
Schada von Borzyskowski, L., Severi, F., Krüger, K. et al. Marine Proteobacteria metabolize glycolate via the β-hydroxyaspartate cycle. Nature 575, 500–504 (2019). https://doi.org/10.1038/s41586-019-1748-4
Environmental Microbiome (2021)
Nature Metabolism (2021)
Nature Catalysis (2021)
Nature Reviews Earth & Environment (2021)