- Peer Review Report
- 10.7554/elife.83845.sa1
Decision letter: Characterisation of an Escherichia coli line that completely lacks ribonucleotide reduction yields insights into the evolution of parasitism and endosymbiosis
- Feb 02, 2023
- Anders Hofer
Article Figures and data Abstract Editor's evaluation Introduction Results Discussion Materials and methods Data availability References Decision letter Author response Article and author information Metrics Abstract Life requires ribonucleotide reduction for de novo synthesis of deoxyribonucleotides. As ribonucleotide reduction has on occasion been lost in parasites and endosymbionts, which are instead dependent on their host for deoxyribonucleotide synthesis, it should in principle be possible to knock this process out if growth media are supplemented with deoxyribonucleosides. We report the creation of a strain of Escherichia coli where all three ribonucleotide reductase operons have been deleted following introduction of a broad spectrum deoxyribonucleoside kinase from Mycoplasma mycoides. Our strain shows slowed but substantial growth in the presence of deoxyribonucleosides. Under limiting deoxyribonucleoside levels, we observe a distinctive filamentous cell morphology, where cells grow but do not appear to divide regularly. Finally, we examined whether our lines can adapt to limited supplies of deoxyribonucleosides, as might occur in the switch from de novo synthesis to dependence on host production during the evolution of parasitism or endosymbiosis. Over the course of an evolution experiment, we observe a 25-fold reduction in the minimum concentration of exogenous deoxyribonucleosides necessary for growth. Genome analysis reveals that several replicate lines carry mutations in deoB and cdd. deoB codes for phosphopentomutase, a key part of the deoxyriboaldolase pathway, which has been hypothesised as an alternative to ribonucleotide reduction for deoxyribonucleotide synthesis. Rather than complementing the loss of ribonucleotide reduction, our experiments reveal that mutations appear that reduce or eliminate the capacity for this pathway to catabolise deoxyribonucleotides, thus preventing their loss via central metabolism. Mutational inactivation of both deoB and cdd is also observed in a number of obligate intracellular bacteria that have lost ribonucleotide reduction. We conclude that our experiments recapitulate key evolutionary steps in the adaptation to life without ribonucleotide reduction. Editor's evaluation Nearly all organisms require a ribonucleotide reductase (RNR) to convert ribonucleotides to their deoxyribonucleotide counterparts. In this important study, the reader learns how the model organism Escherichia coli can adapt to survive without any of its three RNRs. Compelling microbiology experiments to develop this model and analysis of compensatory mutations reveals patterns that are conserved in the few known pathogens that have also eliminated their dependence on an RNR. The manuscript will be of interest to microbiologists, biochemists, and those who work on the evolution of microbial metabolism. https://doi.org/10.7554/eLife.83845.sa0 Decision letter Reviews on Sciety eLife's review process Introduction All life on our planet requires ribonucleotide reduction for de novo synthesis of deoxyribonucleotides, from ribonucleotides (Nordlund and Reichard, 2006; Lundin et al., 2009). The emergence of ribonucleotide reduction likely drove the transition from RNA to DNA early in the evolution of life (Lundin et al., 2015; Reichard, 1993; Torrents et al., 2002; Poole et al., 2002) and its centrality is reflected in its near ubiquity in organismal genomes (Lundin et al., 2009). Structural analyses indicate that the catalytic cores of ribonucleotide reductases (RNRs) are homologous (Lundin et al., 2015; Sintchak et al., 2002), despite sequence similarities between them being extremely limited (Tauer and Benner, 1997). RNRs catalyse deoxyribonucleotide synthesis via a common mechanism wherein a cysteinyl-free radical is generated in the active site. While this similarity in catalysis underscores a common evolutionary origin, RNRs have diverse mechanisms for radical generation (Nordlund and Reichard, 2006). Three broad, evolutionarily distinct, families exist which are also divided on the basis of their radical generation mechanism (Nordlund and Reichard, 2006; Lundin et al., 2015; Lundin et al., 2010). This divergence likely coincides with adaptation to varied environments (Lundin et al., 2010; Hofer et al., 2012). Class I enzymes are strictly aerobic in that they require oxygen or superoxide for activation (Rose et al., 2018). This typically, though not exclusively (Blaesi et al., 2018; Srinivas et al., 2018), occurs through oxidation of a dimetal ion centre from which a stable cysteine radical is generated, either directly, or indirectly through formation of a tyrosyl radical in the small subunit. Class III enzymes generate a stable glycine radical via cleavage of S-adenosyl methionine, and are strictly anaerobic, while class II enzymes generate a 5'-deoxyadenosyl radical through cleavage of adenosylcobalamin and operate irrespective of oxygen presence or absence. To date, very few species have dispensed with ribonucleotide reduction, albeit indirectly, in that they instead rely on their hosts for deoxyribonucleotides (Lundin et al., 2009). Examples include a strain of Buchnera aphidicola str. Cc from the Cedar bark aphid, Cinara cedri (Pérez-Brocal et al., 2006). B. aphidicola are maternally inherited aphid endosymbionts that synthesise amino acids in short supply in the aphid diet, are obligately intracellular, and appear to be on the path to transitioning to an organelle (Andersson, 2000; Tamas et al., 2002). Two bacterial pathogens that live in very close association with their hosts—Ureaplasma urealyticum (Glass et al., 2000) and Borrelia burgdorferi (Fraser et al., 1997)—have also lost genes for ribonucleotide reduction, as have two eukaryotic parasites, Giardia lamblia and Entamoeba histolytica (Lundin et al., 2009). While ribonucleotide reduction is essential, we reasoned that it should be possible to dispense with genes for ribonucleotide reduction if deoxyribonucleotides or their precursors are available via growth media. We report the successful elimination of all three RNR operons from the bacterium Escherichia coli. The knockout strain was created by introducing a broad-spectrum deoxyribonucleoside kinase from Mycoplasma mycoides (Wang et al., 2001), which we predicted would enable use of media-supplied deoxyribonucleosides for deoxyribonucleotide synthesis in the absence of ribonucleotide reduction. As predicted, the resulting strain is dependent on deoxyribonucleosides in the growth medium but cannot grow if the media are supplemented with deoxyribose (dR) plus the four nucleobases (A, G, C, T). To understand the impact of limited deoxyribonucleoside availability, we subjected our knockout line to experimental evolution followed by genome sequencing. We were interested in establishing whether, under such conditions, our lines would utilise the reverse deoxyriboaldolase (DERA) pathway for deoxyribonucleotide synthesis. This pathway is reversible in vitro (Horinouchi et al., 2006a; Ogawa et al., 2003), and has long been considered a plausible alternative to ribonucleotide reduction (Racker, 1951; Racker, 1952) for dNTP synthesis, with some considering it to be a plausible ancestral route for the origin of DNA (Benner et al., 1989; Poole et al., 2014). Our evolution experiments reveal that this pathway is not coopted for deoxyribonucleotide synthesis following loss of ribonucleotide reduction, with cells grown under limiting deoxyribonucleoside levels exhibiting a filamentous morphology indicative of stress. Instead, it appears that the DERA pathway is a liability under conditions of limited deoxyribonucleotide availability; we observe loss-of-function mutations in deoB, which are predicted to prevent catabolism of 2-deoxy-D-ribose-1-phosphate. Available genome sequence data from obligate intracellular bacteria that lack ribonucleotide reduction indicate that deoB has in fact been lost on multiple occasions, suggesting that recycling of dR is disadvantageous under conditions where this sugar is in limited supply; deletion of deoB would enable this sugar to be rerouted for deoxyribonucleotide synthesis but would also preclude de novo synthesis via the reverse DERA reaction. Our results thus illuminate a key adaptive step taken by obligate intracellular and pathogenic species to mitigate the loss of ribonucleotide reduction. Results Creation of an E. coli line lacking ribonuclease reduction E. coli carries genes for three RNRs: aerobic class Ia (encoded by the nrdAB operon), class Ib (nrdHIEF), and anaerobic class III (nrdDG) (Figure 1). Under aerobic growth, ribonucleotide reduction is primarily performed by the iron-dependent class Ia enzyme (Monje-Casas et al., 2001), while class Ib has a manganese metal centre (Cotruvo and Stubbe, 2011) and is therefore able to support growth when iron is scarce (Martin and Imlay, 2011). Figure 1 with 2 supplements see all Download asset Open asset Deletion of ribonucleotide reductase (RNR) genes in E. coli. (A) Schematic of each of the three RNR operons in E. coli and the genomic regions following scarless deletion. Arrows indicate locations of PCR primers (Supplementary file 1B) used to confirm RNR operon presence/absence. (B) Results of PCR with primers external to all three nrd operons in ΔRNR compared with the wild-type progenitor strain (REL606). Lane 1, ladder; lanes 2 and 6, no template. Lanes 3–5, band sizes consistent with operon deletion (ΔnrdAB 1645 bp, ΔnrdDG 1620 bp and ΔnrdHIEF 1477 bp) in the ΔRNR line. Lanes 6–9, band sizes for all three operons in wild-type E. coli REL606 (nrdAB 5303 bp, nrdDG 4382 bp, and nrdHIEF 5191 bp). Figure 1—source data 1 This file contains unlabelled and labelled full gel pictures for Figure 1. https://cdn.elifesciences.org/articles/83845/elife-83845-fig1-data1-v2.zip Download elife-83845-fig1-data1-v2.zip We reasoned that, while ribonucleotide reduction is an essential process, it should be possible to compensate for its loss via media supplementation. To eliminate ribonucleotide reduction, we sequentially deleted the RNR operons using a scarless deletion protocol (Materials and methods). We first deleted nrdDG, as this operon is known not to be expressed under aerobic conditions. We next deleted the nrdHIEF operon from the resulting ΔnrdDG line. Finally, we attempted to delete the nrdAB operon the ΔnrdHIEF ΔnrdDG line, grown in the presence of deoxyribonucleosides (dNs). We were initially unsuccessful in deleting nrdAB. This is not unexpected as, while dNs are known to be taken up, some (deoxyguanine [dG], deoxycytidine [dC], deoxyadenine [dA]) are unavailable for DNA synthesis, even when genes in the DERA pathway are mutated, preventing their catabolism (Karlström, 1970). We reasoned that may be because of the absence of a suitable deoxynucleoside kinase activity for converting those dNs to deoxyribonucleotides. We therefore tried knocking out nrdAB in the presence of a heterologously expressed deoxyadenosine kinase (dAK) gene from M. mycoides. Mycoplasmas are known to possess deoxyadenosine kinases genes that permit utilisation of dNs for DNA synthesis (Wang et al., 2001). To establish if resulting lines lacked all three nrd operons, we first PCR screened for evidence of chromosomal nrd operon deletions. Our results indicate that all three operons were successfully deleted (Figure 1). As deletion does not exclude relocation of functional gene copies to another genomic location, we confirmed gene absence via PCR using primers internal to nrd genes and performed RT-PCR to confirm absence of gene expression (Figure 1—figure supplement 1). Finally, knockout status was confirmed with whole genome sequencing. This confirmed deletion of all genes for ribonucleotide reduction, with the nrdAB operon successfully deleted under heterologous expression of dAK from M. mycoides (Figure 1—figure supplement 2). Of four knockout lines, we selected one isolate (hereafter called ΔRNR) for all subsequent work. In addition to lacking all nrd genes, genome sequencing revealed 23 SNPs that likely appeared during the creation of ΔRNR from wild-type progenitor line (REL606) (Supplementary file 1C). ΔRNR is dependent on dN supplementation Our ΔRNR line was created in the presence of dN supplementation, but this does not mean it is dependent on dN supplementation for growth. We therefore sought to understand what supplementation, if any, our ΔRNR line requires. To determine the lowest concentration of dNs that permits growth, we generated a series of growth curves for differing dN concentrations. At high dN concentrations (1 mg/mL), ΔRNR grows favourably, though shows a clear lag compared to wild-type on equivalent media (Figure 2A). At lower dN concentrations, ΔRNR growth is impaired, but there is still discernible growth at 0.05 mg/mL (Figure 2C). When this is dropped to 0.01 mg/mL we observed only marginal growth of ΔRNR (Figure 2D). As expected, in the absence of dN supplementation no growth is observed, while wild-type lines are unaffected (Figure 2E). Finally, we tried growing ΔRNR on dR plus the four bases (A, G, C, T). No growth was observed (Figure 2F). These results indicate that deletion of the three nrd operons from ΔRNR has completely eliminated the capacity for de novo deoxyribonucleotide synthesis; no other genes appear able to compensate for this deficiency. Figure 2 Download asset Open asset ΔRNR shows growth across a range deoxyribonucleoside (dN) concentrations but cannot grow in the absence of dN supplementation. Growth was monitored for 25 hr in 1× MOPS media + 1% glucose with supplementation, as indicated. Wild-type (REL606, black), ΔRNR (orange). Curves show mean OD595, error bars show SEM, all experiments performed in triplicate. (A) 1 mg/mL, (B) 0.25 mg/mL, (C) 0.05 mg/mL, (D) 0.01 mg/mL, and (E) 0 mg/mL dNs. (F) 0.1 mg/mL deoxyribose (dR) + 0.25 mg/mL of adenine (A), cytosine (C), thymine (T), and guanine (G). We next sought to establish whether ΔRNR requires supplementation of all four dNs, or whether any are dispensable. Strains were once again grown in 1× MOPS + 1% glucose, with individual dNs (dA, dG, deoxythymine [dT], or dC) added at a concentration of 0.25 mg/mL. This revealed that ΔRNR is unable to grow in the presence of dA, dG, or dT (either alone [Figure 3A–C] or in combination [Figure 3E]), even after 44 hr of monitoring (data not shown). However, ΔRNR can grow on dC alone (Figure 3D). Together these data indicate dC is the sole dN required for ΔRNR growth in minimal media. Figure 3 Download asset Open asset ΔRNR requires only deoxycytidine (dC) for growth in minimal media. Growth was monitored for wild-type (REL606, black) and ΔRNR (orange) in 1× MOPS media + 1% glucose, with deoxyribonucleoside (dN) supplementation, as indicated. Growth was monitored for 25 hr. Curves show mean OD595, error bars show SEM, all experiments performed in triplicate. (A) Deoxyadenine (dA) (0.25 mg/mL), (B) deoxyguanine (dG) (0.25 mg/mL), (C) deoxythymine (dT) (0.25 mg/mL), (D) deoxycytidine (dC) (0.25 mg/mL), (E) dA, dG, dT, each at 0.25 mg/mL. ΔRNR exhibits a filamentous cell morphology when grown under limiting dNs During our growth assays we observed a 'clumping' phenotype when ΔRNR is grown in liquid media at concentrations of dNs that limit growth. This contrasts with the uniform cloudy appearance of wild-type E. coli (Figure 4C, top left panel). Examination of cells at ×100 magnification revealed that ΔRNR cells are elongated and filamentous at very low levels of dNs, whereas ΔRNR cells have a similar morphology to wild-type at higher dN concentrations (Figure 4A). Some cells reached lengths several times that of wild-type (Supplementary file 1D). Figure 4 Download asset Open asset ΔRNR exhibit a reversible elongated, filamentous cell morphology with multiple nucleoids when grown in limiting deoxyribonucleoside (dN) concentrations. (A) ×100 Brightfield photographs of REL606 and ∆RNR at 1 and 0.01 mg/mL. (B) FITC and DAPI staining at ×100 of REL606 and ∆RNR at 0.01 mg/mL. (C) ∆RNR phenotype is reversible if dN concentration is increased. Loss of clumpy phenotype occurs when ∆RNR cells are transferred from low (0.01 mg/mL) to high (1 mg/mL) dNs. Brightfield microscopy images reveal that this change is accompanied by a reduction in elongated cell morphology upon transfer to high [dN] (magnification ×1000; scale bar 20 μM). This phenotype appears most pronounced at low dN concentrations. We therefore sought to establish if this phenotype is the result of environmental conditions. We transferred a population of ΔRNR cells exhibiting the filamentous cell phenotype from low dN (0.01 mg/mL) media to a less restrictive environment ([dN]=1 mg/mL). Following transfer to higher [dN], the filamentous phenotype is heavily diminished (Figure 4C). This demonstrates that the phenotype is environmental and can be reversed if cells are moved to a high [dN] environment. One possibility is that, under low [dN], cells are growing but unable to complete cell division. If so, this might be reflected by the presence of multiple DNA-dense regions across the length of the cells. To visualise DNA within filamentous cells, we stained ΔRNR cells grown in low [dN] (0.01 mg/mL) with FITC and DAPI. This revealed the presence of multiple DAPI-stained regions across the length of the cells (Figure 4B), suggestive of the presence of multiple DNA nucleoids in ΔRNR cells grown at low [dN]. Evolution of ∆RNR lines under restricted dN supplementation The loss of ribonucleotide reduction in obligate intracellular lifestyles presumably resulted from relaxed selection on deoxyribonucleotide production when deoxyribonucleotides are available from the environment. However, for both parasitism and endosymbiosis, the host dNTP pool must be shared. We were therefore interested to assess if the ∆RNR line adapts to a reduction in dN availability, as might occur in the evolutionary switch from de novo synthesis to dependency on host production. In order to allow the strains to adapt to a lower concentration of dN supplementation in the growth media, we initiated evolution experiments at dN concentrations where ∆RNR lines grow, but not as well as wild-type (1 mg/mL and 0.25 mg/mL). Five independent lines of ∆RNR were serially passaged in each condition by growing the culture to stationary phase, then transferring 50 μL into a new six-well plate containing 5 mL of fresh MOPS media + 1% glucose and dNs (Figure 5). We also established three control lines of the wild-type progenitor strain (REL606) at each dN concentration. The initial experiment was run for a total of 30 transfers. Glycerol stocks were created every five transfers for each line, with contamination checks performed concurrently. Figure 5 Download asset Open asset Overview of the ΔRNR evolution experiment. Five lines of ΔRNR and three lines of wild-type progenitor (REL606) were established at one of two conditions (1 mg/mL or 0.25 mg/mL deoxyribonucleosides [dNs] in MOPS + 1% glucose), and serially passaged for 30 transfers. Genomic material from each line was then extracted from a single colony and sent for sequencing giving us clonal-level genome information. To further investigate adaptation to low concentrations of dNs, the 'fittest' ΔRNR line grown at 0.25 mg/mL dNs (ΔRNR_250_T30_1) was used to seed a subsequent experiment. Ten replicate lines of ΔRNR_250_T30 _1 were serially passaged for an additional 10 transfers in MOPS + 1% glucose and 0.01 mg/mL dNs. DNA from a single colony of each of these 10 lines (ΔRNR_10_T40_1–10) was then extracted and sent for whole genome sequencing. Elongate cell morphology in ∆RNR lines diminishes over the course of the evolution experiment Our initial observations of unevolved ∆RNR (henceforth ∆RNR_T0) revealed an elongated filamentous cellular phenotype at lower concentrations of dNs (Figure 4). We therefore monitored cell morphology during our evolution experiment. Every five transfers, we examined length of ∆RNR and REL606 cells from each condition (1 mg/mL or 0.25 mg/mL dNs; cells from well position and or were used for all and mean length At the of the experiment, ∆RNR cells are than at 0.25 mg/mL dNs (Figure We cell lengths of 20 at this lower concentration of a times the length of wild-type cells. the the experiment reached 30 transfers, cell length was to wild-type at both dN concentrations (Figure These in morphology cells to adapt to restricted dN Figure Download asset Open asset The ∆RNR cell phenotype diminishes over the course of an evolution experiment. (A) cell length during the course of a evolution experiment. length was for 20 cells for replicate line 1 of ∆RNR and REL606 in either 1 mg/mL or 0.25 mg/mL deoxyribonucleosides (dNs). were taken every five transfers. Over the course of the evolution experiment, ∆RNR cell length at both dN concentrations. (B) of wild-type and ∆RNR cells at transfer 1 and transfer 30 at each dN concentration (1 mg/mL or 0.25 mg/mL) (magnification ×1000; scale bar 2 μM). ∆RNR lines exhibit growth following evolution under restricted dN availability 30 transfers, we sought to determine if ∆RNR lines their capacity to grow under restricted (either 1 mg/mL or 0.25 mg/mL) dN We use the following for our evolution is replicate line of ∆RNR at 1 mg/mL dNs for 30 transfers. At 1 mg/mL, through all growth to (Figure compared to which no growth at 0.01 mg/mL dNs (Figure all lines growth at this low concentration (Figure similar of was also for the five ∆RNR replicate lines in 0.25 mg/mL (Figure supplement 1). that cell morphology and preclude use of for of or cell et al., in the error to REL606 can only be used to a of growth. All lines to grow in the absence of dN supplementation (Figure Figure supplement they on media supplementation. Figure with 1 supplement see all Download asset Open asset Growth of ∆RNR lines at 1 mg/mL deoxyribonucleosides for 30 transfers. Growth was monitored for wild-type and lines in 1 mg/mL dNs for 30 Growth experiments were performed in 1× MOPS media + 1% glucose with the addition of dNs at concentrations. Growth was monitored for 25 hr. Curves show mean OD595, error bars show SEM, all experiments performed in triplicate. (A) 1 mg/mL dNs, (B) 0.25 mg/mL dNs, (C) mg/mL dNs, (D) mg/mL dNs, (E) 0 mg/mL dNs. Evolution of ∆RNR at 0.25 mg/mL dNs results in Following 30 transfers, we the genomes of our lines in 0.25 mg/mL dNs and 1 mg/mL dNs and each of three wild-type control lines and from single The filamentous phenotype of our ∆RNR lines under restricted dN supplement of the number of which in of We therefore report total observed mutations plus for each experimental line (Figure file Our knockout lines and all mutations than wild-type at the dN concentration (Figure file Our lines the of mutations lines exhibit a small but in total SNPs to wild-type (Figure supplement 1 and file was no in total mutations between wild-type lines in either 1 mg/mL or 0.25 mg/mL dNs Figure with 4 supplements see all Download asset Open asset observed following experimental evolution of ∆RNR for 30 transfers. single plus in each line at evolution in 1 mg/mL deoxyribonucleosides RNR evolution in 0.25 mg/mL, RNR is well that RNRs deoxyribonucleotide pool sizes through et al., 2012). We were therefore interested to see whether loss of ribonucleotide reduction resulted in in E. coli. Our was that there would not be any as our lines are supplemented with concentrations of all four dNs. If there was we this to be most in our lines, as these the number of of SNPs in our lines are of one (Figure Figure supplement 1). ∆RNR lines of our genome data from lines at transfer 30 revealed that one which codes for was in of 10 lines from transfer 30 (Supplementary file and of and dC to and of the data revealed that of lines an deletion (Supplementary file Figure supplement while one line an that was predicted to result in a As the are we that this deletion may have in the ancestral but was not at the we initiated our evolution experiment. These genome data were from colony it was not possible to determine if this was the the of this was in the of our evolution experiment. To determine whether cdd is lost in response to the experimental conditions, we screened our ∆RNR to establish if any cdd. We PCR of stocks for two knockout lines This revealed that the that we used to establish our evolution experiments was a with both and cdd deletion (Figure supplement it appears that, after 30 transfers, all 10 lines were for the deleted this the deletion is very likely to have in the ancestral and that our lines were from a despite the knockout lines being established from individual following scarless deletion of nrdAB and methods). of revealed no evidence of a deletion (Figure supplement We therefore established a short evolution experiment with five independent lines, in 0.25 mg/mL dNs, the early of our evolution experiment with the two transfers, we observed the emergence of a cdd deletion in several replicate lines (Figure supplement In similar experiments not and using these lines, we see loss of the cdd under similar growth conditions (data not shown). Finally, that an in cdd but not the suggesting that it is deletion of the cdd that is selected under these conditions. ∆RNR lines in 0.01 mg/mL dNs the DERA pathway transfer one of our lines growth similar to wild-type when grown in 0.25 mg/mL dNs (Figure supplement 1). this line as we initiated a evolution experiment using 10 lines, grown in 0.01 mg/mL dNs (Figure 5). 10 transfers, we these lines All a of with SNPs (Figure supplement 4). As with the lines from our first 30 transfers, the of mutations were Genome analysis of lines revealed mutations (Supplementary file of which were in than one line. Of these were in two lines, in three lines, and 4 in four lines (Supplementary file One deoB, was in lines and is in the of our experiment in that it codes for This enzyme the transfer of a between the and of and is for the step in deoxyribonucleotide (Figure analysis revealed that several of the mutations to deoB were is possible that
Read more