Characterisation of an Escherichia coli line that completely lacks ribonucleotide reduction yields insights into the evolution of parasitism and endosymbiosis
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.sa0Introduction
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).

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
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.

Δ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.

Δ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).

Δ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.

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 A1 [∆RNR] and A1 or B1 [REL606] were used for all measurements, and mean length calculated). At the beginning of the experiment, ∆RNR cells are substantially longer than wild-type, particularly at 0.25 mg/mL dNs (Figure 6). We recorded cell lengths of 20 μM at this lower concentration (average of 6 μM), a staggering 10–20 times the length of wild-type cells. By the time the experiment had reached 30 transfers, average cell length was comparable to wild-type at both dN concentrations (Figure 6A). These changes in gross morphology suggest cells had begun to adapt to restricted dN availability.

The ∆RNR elongate cell phenotype diminishes over the course of an evolution experiment.
(A) Mean cell length during the course of a 30-transfer evolution experiment. Cell length was determined for 20 cells for replicate line 1 of ∆RNR and REL606 evolved in either 1 mg/mL or 0.25 mg/mL deoxyribonucleosides (dNs). Cell measurements were taken every five transfers. Over the course of the evolution experiment, ∆RNR cell length gradually reduces at both dN concentrations. (B) Morphology of wild-type (top row) and ∆RNR (bottom row) 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 improved growth following evolution under restricted dN availability
After 30 transfers, we sought to determine if ∆RNR lines had improved their capacity to grow under restricted (either 1 mg/mL or 0.25 mg/mL) dN availability. We use the following nomenclature for our evolution lines: ∆RNR_[dN] _transfer#_line#, so ∆RNR1000_T30_L1 is replicate line #1 of ∆RNR evolved at 1 mg/mL (i.e. 1000 μg/mL) dNs for 30 transfers. At 1 mg/mL, ∆RNR1000_T30_L1 through L5 all showed improved growth relative to ∆RNR_T0 (Figure 7A). Moreover, compared to ∆RNR_T0, which exhibited almost no growth at 0.01 mg/mL dNs (Figure 2D), all evolved lines (∆RNR1000_T30_L1–5) showed improved growth at this low concentration (Figure 7D). A similar overall pattern of improvement was also seen for the five ∆RNR replicate lines evolved in 0.25 mg/mL (∆RNR250_T30_L1–5) (Figure 7—figure supplement 1). Note however that elongate cell morphology and clumping preclude use of OD measurements for accurate estimation of doubling time or cell counts (Stevenson et al., 2016) (reflected in the large standard error relative to REL606 controls), so can only be used to give a general indication of growth. All lines failed to grow in the absence of dN supplementation (Figure 7E, Figure 7—figure supplement 1), indicating they remain reliant on media supplementation.

Growth characterisation of ∆RNR lines evolved at 1 mg/mL deoxyribonucleosides (dNs) for 30 transfers.
Growth was monitored for wild-type (REL606), ancestor (∆RNR_T0), and evolved lines (∆RNR1000_T30_L1–5, evolved in 1 mg/mL dNs for 30 transfers). Growth experiments were performed in 1× MOPS media + 1% glucose with the addition of dNs at indicated 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) 0.055 mg/mL dNs, (D) 0.011 mg/mL dNs, (E) 0 mg/mL dNs.
Evolution of ∆RNR at 0.25 mg/mL dNs results in A:T → G:C mutational skew
Following 30 transfers, we sequenced the genomes of ∆RNR_T0, our lines evolved in 0.25 mg/mL dNs (∆RNR250_T30_L1–5) and 1 mg/mL dNs (∆RNR1000_T30_L1–5), and each set of three wild-type control lines (REL606_1000_T30_L1–3 and REL606_250_T30_L1–3) from single colonies. The filamentous phenotype of our ∆RNR lines under restricted dN supplement precluded accurate estimation of the number of generations per transfer, which in turn precluded reliable calculation of mutation rates. We therefore report total observed mutations (single nucleotide substitutions plus indels) for each experimental line (Figure 8, Supplementary file 1E). Our evolved knockout lines (∆RNR250_T30_L1–5 and ∆RNR1000_T30_L1–5) all accumulated substantially more mutations than wild-type controls evolved at the same dN concentration (Figure 8, Supplementary file 1E). Our ∆RNR250_T30 lines accumulated the greatest numbers of mutations (range: 39–118 SNPs). ∆RNR1000_T30 lines did exhibit a small but significant increase in total SNPs relative to wild-type (Figure 8—figure supplement 1 and Supplementary file 1E; p=0.01, unpaired t-test). There was no significant difference in total mutations between wild-type lines evolved in either 1 mg/mL or 0.25 mg/mL dNs (p=0.10; unpaired t-test).

Mutations observed following experimental evolution of ∆RNR for 30 transfers.
Total single nucleotide substitutions plus indels present in each line at T30. For evolution in 1 mg/mL deoxyribonucleosides (dNs), REL606=control, RNR = ∆RNR1000_T30_L1–5. For evolution in 0.25 mg/mL, REL 606=control, RNR = ∆RNR250_T30_L1–5.
It is well documented that RNRs keep deoxyribonucleotide pool sizes balanced through allosteric regulation (Hofer et al., 2012). We were therefore interested to see whether loss of ribonucleotide reduction resulted in mutational skew in E. coli. Our null expectation was that there would not be any detectable skew as our lines are supplemented with equal concentrations of all four dNs. If there was skew, we expected this to be most evident in our ∆RNR250_T30 lines, as these accumulated the greatest number of mutations. Indeed, 88% of SNPs in our ∆RNR250_T30 lines are of one type: A:T→G:C (Figure 8, Figure 8—figure supplement 1).
∆RNR lines lose cytidine deaminase
Analysis of our genome data from lines at transfer 30 revealed that one gene, cdd, which codes for cytidine deaminase, was mutated in 8 of 10 lines from transfer 30 (Supplementary file 1F and G). Cytidine deaminase catalyses deamination of cytidine and dC to uridine and deoxyuridine. Analysis of the data revealed that seven of eight lines carried an identical deletion (Supplementary file 1F, G; Figure 8—figure supplement 2), while one line possessed an SNP that was predicted to result in a truncation. As the deletions are identical, we suspected that this deletion may have occurred in the ancestral population, but was not fixed at the time we initiated our evolution experiment. These genome data were derived from colony isolates, meaning it was not possible to determine if this was the case. Thus, the significance of this mutation was unclear in the context of our evolution experiment. To determine whether cdd is lost in response to the experimental conditions, we screened our original ∆RNR isolates to establish if any carried intact cdd. We undertook PCR screening of glycerol stocks for two knockout lines (∆RNR31, ∆RNR34). This revealed that the glycerol stock that we used to establish our evolution experiments (∆RNR31) was a mixed population, with both intact and cdd deletion present (Figure 8—figure supplement 3). Moreover, it appears that, after 30 transfers, all 10 lines (∆RNR_T30_1000_L1–5, ∆RNR_T30_250_L1–5) were polymorphic for the deleted locus. Together, this indicates the deletion is very likely to have occurred in the ancestral population, and that our lines were derived from a genetically heterogeneous population, despite the knockout lines being established from individual colonies following scarless deletion of nrdAB (Methods and methods). Screening of ∆RNR34 glycerol stock however revealed no evidence of a deletion (Figure 8—figure supplement 3). We therefore established a short evolution experiment with five independent ∆RNR34 lines, in 0.25 mg/mL dNs, mimicking the early stages of our evolution experiment with the ∆RNR31 isolate. After two transfers, we observed the emergence of a cdd deletion in several replicate lines (Figure 8—figure supplement 3). In similar experiments not reported here and using these lines, we repeatedly see loss of the cdd locus under similar growth conditions (data not shown). Finally, that ∆RNR_T30_250_L5 carried an inactivating mutation in cdd but not the upstream locus (yohk), suggesting that it is deletion of the cdd locus that is selected under these conditions.
∆RNR lines evolved in 0.01 mg/mL dNs mutate the DERA salvage pathway
By transfer 30, one of our lines (∆RNR_T30_250_L1) exhibited growth similar to wild-type when grown in 0.25 mg/mL dNs (Figure 7—figure supplement 1). Using this line as progenitor, we initiated a second evolution experiment using 10 parallel lines, grown in 0.01 mg/mL dNs (Figure 5). After 10 transfers, we sequenced these lines (∆RNR_T40_10_L1.1–1.10). All exhibited a marked accumulation of mutations, with L1.7 accumulating 220 SNPs (Figure 8—figure supplement 4). As with the lines from our first 30 transfers, the majority of mutations were A:T→G:C substitutions.
Genome analysis of T40 lines revealed 648 mutations (Supplementary file 1H), of which 80 were found in more than one line. Of these 59 were present in two lines, 16 in three lines, and 4 in four lines (Supplementary file 1I).
One gene, deoB, was mutated in eight lines (Table 1), and is particularly noteworthy in the context of our experiment in that it codes for phosphopentomutase. This enzyme catalyses the transfer of a phosphate group between the C1 and C5 carbon atoms of deoxribose and is responsible for the committed step in deoxyribonucleotide salvage (Figure 9). Closer analysis revealed that several of the mutations to deoB were identical. It is possible that these mutations have occurred independently, but it may also be possible that these were the result of cross-contamination between lines during the evolution experiment. To assess this, we examined whether the affected lines carried a large number of shared mutations. Lines L1.1 and L1.8 had 47 and 157 unique mutations, with only the deoB mutation shared. For this pattern to be explained by cross-contamination, this mutation would have had to have been the very first mutation to appear, and cross-contamination would have had to have occurred before any other mutation occurred. With no evidence of deoB mutation in the progenitor line at transfer 30 (∆RNR250_T30_L1), we consider it unlikely that cross-contamination between these lines has occurred. In the second case, L1.3 (53 SNPs), L1.5 (40 SNPs), and L1.10 (123 SNPs) share a common deoB mutation along with 13 other mutations. On these data, it is harder to rule out cross-contamination. Conservatively, it thus appears that deoB has independently mutated five or six times under the conditions of our experiment.

Impact of mutations on deoxyribonucleotide synthesis.
(A) Deoxyribonucleotides in free-living bacteria are primarily synthesised intracellularly via ribonucleotide reduction (RNR), with some uptake of deoxyribonucleosides (dNs) from the environment also possible in principle. Deoxyribonucleotides are required for genome synthesis, but may also be metabolised via the DERA pathway, the end product (G3P) of which feeds into glycolysis (grey arrows). (B) In intracellular bacteria lacking RNR (crossed-out light green arrow), deoxyribonucleotides must instead be derived from the host as dNs (dark green arrow). Loss of RNR leads to loss of phosphopentomutase (PPMase) activity (crossed-out light blue arrow), which prevents loss of dNs via glycolysis. Abbreviations: RNR—ribonucleotide reductase; dNs—deoxyribonucleosides; DR1P—2-deoxyribose 1-phosphate; DR5P—2-deoxyribose 5-phosphate; G3P—D-glyceraldehyde 3-phosphate; DHAP—dihydroxyacetone phosphate; Glu—glucose; FDP—fructose 1,6-diphosphate; Pyr—pyruvate, dNKases—deoxynucleotide kinases; DERA—deoxyriboaldolase.
Genome sequencing of ∆RNR10_T40 lines identifies multiple mutations to deoB.
Mutation | Location in ORF | Effect | Lines impacted |
---|---|---|---|
Substitution | 1057 (A→C) | AA substitution (T353P) | L1 L8 |
Deletion | 583 (G(1)→G(0)) | ORF truncation (frameshift) | L2 |
Deletion | 926 (G(4)→G(3)) | ORF truncation (frameshift) | L3 L5 L10 |
Substitution | 932 (A→C) | AA substitution (H311P) | L4 |
Substitution | 1147 (A→C) | AA substitution (T383P) | L6 |
To next sought to assess the impact of the mutations in deoB, on phosphopentomutase function. We first generated an alignment of 100 homologous phosphopentomutases, including from E. coli (REL606) and Bacillus cereus (Supplementary file 3). The latter was included as protein structure and active site residues have been characterised (Panosian et al., 2011). Using the alignment as a guide, we next assessed functional impact of observed deoB mutations by mapping these onto the structure of B. cereus phosphopentomutase (PDB ID: C3M8). Of the five unique mutations observed in our experiment (Table 1), three directly impact conserved active site residues, either through truncation or mutation of active site residues (Table 2). The remaining two mutations are not directly associated with the active site but result in amino acid substitutions not observed in any sequences in our alignment. We therefore conclude that these mutations to the deoB gene are all likely to impact phosphopentomutase function, with some expected to completely inactivate the enzyme.
Predicted impact of observed deoB mutations on phosphopentomutase function.
B. cereus PPMase* | E. coli PPMase† | L1, L8 | L2 | L4 | L6 | L3, L5, L10 |
---|---|---|---|---|---|---|
Active site residues | ||||||
D13 | D10 | ✓ | ✓ | ✓ | ✓ | ✓ |
T85 | T98 | ✓ | ✓ | ✓ | ✓ | ✓ |
D156 | D173 | ✓ | ✓ | ✓ | ✓ | ✓ |
D286 | D306 | ✓ | - | ✓ | ✓ | ✓ |
H291 | H311 | ✓ | - | H311P | ✓ | - |
A327 | D347 | ✓ | - | ✓ | ✓ | - |
H328 | H348 | ✓ | - | ✓ | ✓ | - |
H339 | H359 | ✓ | - | ✓ | ✓ | - |
Non-active site residues | ||||||
I333 | T353 | T353P | - | ✓ | ✓ | - |
T365 | T383P | ✓ | - | ✓ | T383P | - |
-
*
PDB ID: 3M8W.
-
†
Homologous residues derived from multiple sequence alignment (Supplementary file 3).
Loss of ribonucleotide reduction among natural obligate intracellular bacteria is frequently associated with loss of deoB and cdd coding potential
Loss of ribonucleotide reduction genes has previously been noted in a small number species (Lundin et al., 2009), including the bacteria B. burgdorferi, U. urealyticum, and B. aphidicola str. Cc. We were therefore interested in examining if loss of deoB and cdd genes observed in our experiments extends to these species and the genera to which they belong. Sequences used for blast searches are given in Supplementary file 1J. For the three specific species above, blastp searches (Supplementary file 1K) confirmed absence of RNRs and phosphopentomutase (encoded by deoB). A hit to cytidine deaminase (encoded by cdd) was identified for B. burgdorferi. To examine if our results extend to other members of each genus, we repeated our searches for all members of each genus for which there are genome assemblies in the NCBI genome database (https://www.ncbi.nlm.nih.gov/genome/). For Ureaplasma (taxid:2129), of 58 genome assemblies, only one, Candidatus Ureaplasma intestinipullorum (a MAG from a chicken gut microbiome study; Gilroy et al., 2021) returned hits to RNR query sequences (Supplementary file 1L). No other Ureaplasma assemblies revealed evidence of RNR, PPMase, CDD, or YohK coding potential. Thus, with the possible exception of Candidatus Ureaplasma intestinipullorum, no members of the genus Ureaplasma carry genes for RNR, PPM, or CDD. These results align well with our experimental evolution study, where loss of ribonucleotide reduction resulted in subsequent loss of deoB and cdd genes.
The situation in Borrelia (taxid:138) is more complex. We screened 77 available genome assemblies, recovering 73 hits to NrdE and 67 to NrdF, indicating that the majority of Borrelia genome species code for class Ib RNR. At the time we performed our analysis, we noticed that Borrelia burgdorferi has been renamed Boreliella burgdorferi, following a proposal from Adeolu and Gupta, 2014, to split Borrelia into two genera, Borrelia and Borreliella. While this reclassification is the subject of ongoing debate (Barbour et al., 2017; Margos et al., 2017; Gupta and Bergström, 2019; Margos et al., 2018; Margos et al., 2020), we expanded our search to all members of Borreliella (taxid: 64895, 17 species, 607 genome assemblies), and use the taxonomic distinction here so as to enable replication of our analyses. Only one member of the Borelliella returned any hits to RNR: Candidatus Borreliella tachyglossi appears to have a class Ib RNR (Supplementary file 1M). Given that the placement of Candidatus Borreliella tachyglossi is deep but appears unresolved, this may help to pinpoint the timing of RNR Ib loss from the clade/genus Borelliella.
For CDD, our search against Borrelia (taxid:138) returned 20 hits, all annotated as cytidine deaminases but these are approximately half the length of the query sequence. This is not unexpected as the E. coli protein is known to carry two cytidine deaminase domains and is a homodimer, whereas the shorter (homotetrameric) cytidine deaminase proteins from Borrelia are common across bacteria. For the same search against Borreliella (taxid:64895), 33 hits were returned, all annotated as cytidine deaminases. For all genome assemblies from both candidate genera, no significant hits to either PPM (deoB) or YohK were returned. These data indicate that no members of either Borrelia or Borreliella carry deoB, with only the Borreliella apparently lacking ribonucleotide reduction.
Lundin et al., 2009, noted that only one specific strain of B. aphidicola (from the cedar aphid C. cedri) lacked RNR genes. At the time, this was the smallest Buchnera genome, and remains one of the smallest genomes reported for a bacterium (Pérez-Brocal et al., 2006). A search against B. aphidicola Str. Cc (taxid:372461) revealed no coding potential for RNR, PPMase, CDD, and YokH (Supplementary file 1K). We next examined all members of B. aphidicola (taxid:9). This revealed widespread presence of RNR Ia (encoded by nrdA, nrdB) and PPMase (deoB), but no evidence of either CDD or YohK coding potential. Of 74 assemblies, we found 22 strains that lack both ribonucleotide reduction and phosphopentomutase coding, two strains that possess a deoB gene but lack genes for ribonucleotide reduction, and two with the opposite pattern (Supplementary file 1N).
Discussion
While all life is dependent on ribonucleotide reduction for the synthesis of deoxyribonucleotides, a handful of species are known to lack RNR genes (Lundin et al., 2009) and presumably depend on their host for a source of dNTPs. In this work, we report the successful deletion of all RNR operons from E. coli (Figure 1). The resulting line (∆RNR) was sequenced to confirm genomic deletion of RNR genes. Our ∆RNR knockout line grows on MOPS minimal media + 1% glucose supplemented with dNs (Figure 2), but does not grow when provided with their constituent components (dR and the four bases A, G, C, T). Growth experiments reveal that, of the four dNs, only dC is essential for growth (Figure 3). Under limiting levels of dNs, our lines exhibit a filamentous phenotype that can be reversed by increasing dN supplementation (Figure 4). Filamentous phenotypes have been documented under a range of stress conditions, including in the presence of β-lactam antibiotics (Stevenson et al., 2016; Miller et al., 2004; Yao et al., 2012). As cell length changes during growth in our experiment, it is not possible to use OD measurements for quantification of either doubling time or cell number (Stevenson et al., 2016). Our experiments should thus be taken as a qualitative indication of growth. Interestingly, we found that this phenotype reduced over the course of the experiment, suggesting that our lines were adapting to the scarcity of dNs. It will be interesting to determine whether this change is associated specifically with disruption of deoB, the effect of which is predicted to make dNs more readily available for DNA synthesis.
In order to understand how obligate intracellular bacteria adapted to the loss of ribonucleotide reduction, we established an evolution experiment where we serially passaged replicate lines in either 1 mg/mL or 0.25 mg/mL dNs. The resulting lines were sequenced after 30 transfers and one line evolved in 0.25 mg/mL dNs (ΔRNR_250_T30_1) was selected for a second round of evolution at 0.01 mg/mL dNs through a further 10 transfers (Figure 5).
Genome sequencing revealed independent mutations to two genes across our replicate lines. One gene, cdd, carried mutations in 8/10 lines by transfer 30. The cdd locus codes for cytidine deaminase, which converts dC to deoxyuridine, which in turn will be available for dTTP production. Our data reveal two independent mutations (a segmental deletion in seven lines, and an inactivating SNP in the eighth). The segmental deletion eliminated part of the cdd ORF and part of ORF of the upstream gene, yohK, which, together with yohJ, has been shown to code a 3-hydroxypropionate transporter (Nguyen-Vo et al., 2020), and was positionally identical across all seven lines. This deletion likely occurred prior to the establishment of our replicate lines, but appeared not to have gone to fixation. An independent ∆RNR line (∆RNR34) that lacked this deletion rapidly lost cdd when passaged in the experimental growth media (Figure 8—figure supplement 3). This, together with the independent truncation mutation in ∆RNR_T30_250_L5, which impacts cdd but not yohK, suggests there is selection against cdd but not necessarily yohK. Given the fact that our ∆RNR lines can grow with dC as sole dN supplement but do not grow if dC is omitted from the dN mix, it may be that eliminating cytidine deamination prevents loss of this key dN through deamination. We have not assessed growth on dU, but one possibility is that deoxyuridine is less amenable as a substrate for production of all four deoxyribonucleotides. One intriguing observation is that, despite the rapid loss of cdd, this loss of function appears not to have occurred across all replicates by transfer 30. This may be a consequence of the length of the evolution experiment. Comparing the genomes of species known to lack genes for ribonucleotide reduction (Lundin et al., 2009) indicates that absence of cdd frequently coincides with this state. Our genome analyses show that all confirmed members of the genus Ureaplasma lack nrd genes and cdd, suggesting that loss of cdd function is selectively advantageous. Establishing whether loss of one gene drives the loss of the other is not straightforward however. An analysis of strains of B. aphidicola revealed that 48/74 strains retain ribonucleotide reduction, but 0/74 possess cdd (Supplementary file 1N), indicating that the latter gene can be lost in lines capable of synthesising their own deoxyribonucleotides. The picture is also not clear in Borrelia, with cdd found in species variously possessing or lacking ribonucleotide reduction (Supplementary file 1M).
Deletion of yohK and cdd results in a small, but intact ORF (Figure 8—figure supplement 2). While this seems unlikely to be functional, the creation of an ORF, particularly where the upstream gene codes for a transmembrane protein, may be consistent with the ‘Car Trunk’ hypothesis for the avoidance of mutations that would create a toxic protein (Omer et al., 2017). The fact that deletion of these genes is reproducible under defined experimental conditions (Figure 8—figure supplement 3) suggests that this region might be suitable for testing this hypothesis.
All lines grown in 0.25 mg/mL dNs showed a marked mutational skew (A:T→G:C) that differs from the G:C→A:T mutational skew expected for E. coli (Lee et al., 2012), and more generally for bacteria (Hershberg and Petrov, 2010). Knockout of cdd does not appear linked to this skew; 4/5 lines grown at 1 mg/mL dNs do not show an obvious A:T→G:C skew (Figure 8). For lines grown at 0.25 mg/mL dNs, two (L2, L3) carry intact cdd genes and both genomes show evidence of A:T→G:C mutational skew. Such skew has been previously shown to be associated with defective mismatch repair (Lee et al., 2012). By transfer 30, we observe only two lines with SNPs in mismatch repair genes, with only one mutation being nonsynonymous (Supplementary file 1). Thus, the skew we observe cannot be explained by widespread mutation to the mismatch repair pathway. The skew may instead be driven in part by the growth dependency on dC (Figure 3D), with the other dNs unable to compensate for lack of dC (Figure 3E) but appearing to provide minor improvements to growth (Figure 2B). It is not currently clear why this should be so. The two major nucleoside transporters in E. coli, NupC and NupG, are able to transport dNs (Patching et al., 2005), so import does not appear to be a barrier.
We also observed loss of function of deoB in our experiments (Table 1). This observation is interesting for several reasons. First, it indicates that loss of ribonucleotide reduction does not lead to compensatory deoxyribonucleotide synthesis via the DERA pathway, which has been speculated to be an ancestral route for dNTP synthesis, predating the advent of ribonucleotide reduction (Benner et al., 1989; Poole et al., 2014). While the enzymes from this pathway can operate biosynthetically (Horinouchi et al., 2006a; Ogawa et al., 2003; Horinouchi et al., 2006b), loss of phosphopentomutase function via mutation of deoB suggests that this pathway is not readily accessible for synthesis of deoxyribonucleotides in vivo, even though there would have been a strong selective advantage to such an alternative synthetic route during our experiment. We note however that both acetaldehyde toxicity (to both DERA and cells more generally) and low intracellular availability may preclude a simple switch to using the DERA pathway for dNTP production (Ogawa et al., 2003; Poole et al., 2014).
Indeed, under the conditions of our experiment, where dN supplementation is reduced by 25-fold, we observe strong selection to dispense with deoB. This suggests that not only is the reverse DERA pathway inaccessible for biosynthesis, but salvage is selected against when dNTPs are scarce. In the second phase of our experiment, dN supplementation is dropped from 0.25 mg/mL to 0.01 mg/mL, and it is under these conditions of severe dN restriction that we observe parallel loss-of function mutations at the deoB locus (Table 1). We conclude that, when deoxyribonucleotides are in short supply, deoB loss of function ensures that dNTPs are not lost to glycolysis via salvage (Figure 9). The biochemistry of the DERA pathway supports mutation of deoB as the most plausible target, since this step is responsible for interconversion of deoxyribose-5-phosphate and deoxyribose-1-phosphate. The latter is the substrate for phosphorylase-mediated addition of nucleobases, generating dNs, whereas the former is the substrate for DERA-mediated salvage. Thus, deoB loss-of-function mutations serve to prevent dR loss to central metabolism while still permitting phosphorylase-mediated dN resynthesis (Figure 9). Thus, while conversion of DERA to a synthetic pathway might be a better long-term solution to loss of ribonucleotide reduction, loss of function of this pathway is advantageous in the short term, precluding its cooption to deoxyribonucleotide synthesis. This interpretation is consistent with a recent experiment where mutation of dihydrofolate reductase led to deoB loss of function to prevent loss of deoxyribose-5-phosphate to glycolysis (Rodrigues and Shakhnovich, 2019).
Analysis of species that lack nrd genes reveals frequent loss of deoB, suggesting that loss of ribonucleotide reduction might drive loss of deoB, as observed in our experiments. For members of Ureaplasma, the pattern is clear (Supplementary file 1L): all members of this genus completely lack genes for ribonucleotide reduction and show no evidence of coding for either deoB or cdd, the two most mutated genes in our experiments. For B. aphidicola cdd appears to have been completely lost from all sequenced strains. Close to half of sequenced Buchnera strains show co-occurrence of ribonucleotide reduction and phosphopentomutase genes (48/74 strains), with 30% (22/74 strains) lacking genes for both functions (Supplementary file 1N), suggesting that loss of one function frequently results in loss of the other. However, two strains appear to carry deoB while lacking nrd genes, and two further strains show the opposite pattern. If confirmed, this would indicate loss of one function may not always result in loss of the other.
We observe a very clear pattern in members of the Borreliaceae. The NCBI Taxonomy Browser (Schoch et al., 2020) lists two genera with genome data, Borrelia and Borreliella, following a recent proposal (Adeolu and Gupta, 2014; Barbour et al., 2017; Gupta and Bergström, 2019) to split the genus Borrelia into two (though this proposal has been contested; Margos et al., 2017; Margos et al., 2018; Margos et al., 2020). Our analyses indicate that Borrelia (sensu Adeolu and Gupta, 2014) code class Ib RNRs but lack deoB, while Borrelliela lack both deoB and RNR genes. This pattern could be taken to indicate that members of the Borrelia lost deoB but retained ribonucleotide reduction. However, phylogenetic analyses (Zhong et al., 2006) and location on linear plasmids (Miller et al., 2013) indicate that the class Ib RNR genes (nrdEF) in Borrelia has been acquired via horizontal gene transfer, presumably following initial loss of ribonucleotide reduction in the Borreliaceae.
It is noteworthy that class Ib RNR, carried by members of the Borrelia, is manganese-dependent (Cotruvo and Stubbe, 2010), whereas class Ia is an iron-dependent enzyme (Torrents, 2014). Interestingly, the most well-studied member of the Borreliaceae, Borreliella burgdorferi (sensu Adeolu and Gupta) does not require iron (Posey and Gherardini, 2000). It is tempting to speculate that initial loss of ribonucleotide reduction was driven by selection for iron independence, with subsequent reacquisition of ribonucleotide reduction favouring the manganese-dependent enzyme.
More generally, the switch to an intracellular existence may result in relaxed selection on ribonucleotide reduction if host-derived deoxyribonucleotides are sufficient to replace endogenous production. This appears to be an infrequent outcome as most intracellular bacteria retain ribonucleotide reduction. This rare event might be driven by additional pressures, as evolution of iron independence in B. burgdorferi suggests. Our experiments indicate that, in the absence of endogenous ribonucleotide reduction, disruption of the DERA pathway via mutations to deoB may serve to minimise loss of these building blocks when dNTP availability is limited. Given the necessity of genome replication, we conclude that life without ribonucleotide reduction drives the loss of dNTP catabolism via DERA, enabling building block reuse for DNA synthesis.
It is interesting to note that loss of ribonucleotide reduction appears extremely rare in the evolution of endosymbiosis and parasitism, despite the seeming ubiquity of host-derived deoxyribonucleotides. Loss of production in favour of acquisition from the host ought perhaps to be a more common outcome. Acquisition from the host would appear straightforward as many bacteria are known to possess mechanisms for uptake of dNs, deoxyribonucleotides, and even DNA (Huang et al., 2022). With respect to dN uptake, E. coli for instance codes for two broad-spectrum nucleoside transporters, NupC and NupG, which have been reported to take up dNs (Patching et al., 2005). These broad specificities are noteworthy in the context of our experimental results showing ΔRNR is viable when the media is supplemented with dC alone (Figure 3). It may be that imported dNs are not readily utilised in DNA synthesis as they must first be converted to deoxyribonucleotides for utilisation in genome synthesis. In our experiments, it was only possible to eliminate ribonucleotide reduction after first introducing a broad-spectrum deoxyadenosine kinase gene from M. mycoides, presumably because it improves the opportunity for imported dNs to be coopted for deoxyribonucleotide synthesis instead of being catabolised. Indeed, a recent study shows that dNs, deoxyribonucleotides, and DNA can all be used as a nitrogen source by E. coli (Huang et al., 2022). Those authors noted elevated expression of genes involved in catabolism, including deoD, which codes nucleoside phosphorylase, part of the DERA pathway (Figure 9). Thus, under conditions where DNA and its monomeric building blocks are important sources of nitrogen and phosphorus, loss of ribonucleotide reduction would presumably have a deleterious impact on growth. This trade-off may explain the small number of instances of loss of ribonucleotide reduction.
It remains to be seen whether there are conditions wherein the reverse DERA pathway (Poole et al., 2014) can operate synthetically in vivo as an alternative to ribonucleotide reduction, but one possibility is that this may be favoured under conditions where other nitrogen and phosphorous sources are abundant.
Materials and methods
Strains and growth conditions
Request a detailed protocolAll E. coli strains used in this study are listed in Supplementary file 1A. E. coli B strain REL606 (Lenski et al., 1991) was used as the wild-type strain for all experiments. We elected to perform experiments using REL606 as it is well suited for experimental evolution, having been subjected to a long-running evolution experiment (Lenski et al., 1991). Two key features that make this strain convenient for such work are genomically coded resistance against streptomycin and a selectively neutral mutation in the araA gene, making it easy to detect via colony colour when grown on tetrazolium arabinose. This makes it suitable for tracing during competition assays.
Strains were grown in LB (1% tryptone, 1% NaCl, and 0.5% yeast extract ± 2% agar) with the addition of 1 mg/mL each of the four dNs (dG, dA, dC, and dT (AK Scientific)) for mutant strains, unless stated otherwise. For evolution experiments, strains were grown in 1× MOPS media (Neidhardt et al., 1974) supplemented with dNs, as specified. All experiments were conducted in the presence of 100 μg/mL streptomycin.
Creation of an RNR mutant
Request a detailed protocolE. coli has three RNR operons: nrdAB codes for class Ia RNR, nrdDG codes for class III RNR, and nrdHIEF codes for the class Ib enzyme. Deletion of RNR operons from the BL-21 E. coli line REL606 (Supplementary file 1A) was carried out in series using a scarless genome engineering method (Fehér, 2008): the nrdDG operon (class III RNR) was deleted first, followed by nrdHIEF (class IB) then nrdAB (class IA). For each operon, 500 bp upstream and downstream from the start and stop codon was PCR amplified. A second round of overlap PCR combined fragments together and the three resulting product was cloned into pST76a. The resulting construct was transformed into wild-type at 30°C, then with genome integration induced by increasing temperature to 42°C. Helper plasmid (pSTKST) was transformed into the resulting strain, and deletion was induced with chlortetracycline. Successful deletion was confirmed using PCR across the operon (primers used in this study are listed in Supplementary file 1B).
To knock out the nrdAB operon, we sought to identify conditions where de novo deoxyribonucleotide synthesis could be replaced by supplementation. We determined that dN supplementation could replace de novo synthesis provided that we first introduced dN kinase (dNKase) activity to make dNs available for dNTP synthesis. We therefore transformed our lines with a pBAD33 construct containing a dNKase gene from M. mycoides (Wang et al., 2001) that had been previously shown to phosphorylate dA, dG, and dC. As E. coli codes for thymidine kinase activity (tdk), we reasoned that expression of the M. mycoides dNKase (Mm-dAK—characterised as a dN adenosine kinase as it has highest affinity for this dN) should render the nrdAB operon non-essential and thus amenable to knock out. We cloned the Mm-dAK gene (codon-optimised for E. coli expression [GenScript; Supplementary file 2]) into pBAD33 under control of an IPTG-inducible tac promoter with rrnB terminator into our double knockout line (ΔnrdDGΔnrdHIEF). We used scarless deletion with lines grown in LB plus 0.5% wt/vol dA, dC, dG. We omitted dT as dTTP is synthesised via deamination of dCTP and dC to dUTP and dU, respectively.
The resulting strain, ΔRNR (Supplementary file 1A), was confirmed to lack all three operons coding for ribonucleotide reduction (ΔnrdAB, ΔnrdHIEF, ΔnrdDG) using PCR, RT-PCR (Figure 1—figure supplement 1), and whole genome sequencing.
Total RNA extraction and RT-PCR
Request a detailed protocolTo confirm absence of RNR expression, strains were grown overnight in LB and diluted to 1:100 in 10 mL fresh media after reaching stationary phase. The cultures were then grown for ~3 hr and total RNA was isolated from mid-log phase cultures using TRIzol Max Bacterial Isolation kit (Thermo Fisher, Catalogue# 16096020). RNA was quantified using a Qubit 4.0 fluorometer. Purified RNA was diluted to 300 ng/mL and treated with TURBO DNaseI (Ambion, Catalogue# AM2238). This DNA-free RNA was then subjected to RT-PCR using the SuperScript III One-Step RT-PCR system with Platinum Taq DNA Polymerase kit (Invitrogen, Catalogue# 12574018) with primers specific to the gene of interest (Supplementary file 1B). Mm-dAK expression was also confirmed using RT-PCR. We observed Mm-dAK expression in our ΔRNR line with and without IPTG induction (Figure 1—figure supplement 2). We therefore elected to perform evolution experiments without IPTG.
Growth assays
Request a detailed protocolBacterial strains were retrieved from –80°C stocks and grown overnight in LB (1 mg/mL dNs were added to media for growth of ΔRNR lines). Cultures were washed twice in 1× PBS and 10 μL equivalent was added to fresh 24-well plates containing 1 mL of 1× MOPS + 1% glucose. Differing concentrations of dNs were added from a 20 mg/mL stock solution to each individual well. Growth was monitored for 24–48 hr, taking measurements (OD595) every 6 min using a FLUOstar Omega Microplate Reader (BMG Labtech). All experiments were performed in triplicate.
Microscopy
Request a detailed protocolOvernight cultures were grown at 37°C in MOPS media (with dNs added as required). Cultures were twice spun down and washed in 1× PBS. Ten µl of culture was aliquoted onto a microscope slide. Brightfield images were taken using a LEICA ICC50 W microscope (Leica, v.3.2.0) and imported to Photoshop (v.22.4.2) for cell length measurement. Cell length measurements were determined for each strain by averaging from 20 observations (length of the first five cells from top left to right were counted from each of four images).
To visualise DNA, we used a modified version of a FITC/DRAQ-5 double-staining protocol (Silva et al., 2010) instead using DAPI in place of DRAQ-5. Coverslips were coated in poly-D-lysine and placed at 37°C overnight. Coverslips were washed twice with water, dried, and stored at 4°C until ready for use. Overnight cultures were washed in 1× PBS, resuspended, and placed on a poly-D-lysine-coated coverslip at 37°C for 1 hr. Cells were washed with PBS, and fixed using 4% paraformaldehyde at room temperature for 10 min. Coverslips were washed with 1× PBS, followed by 1% Triton X-100 for 5 min. Following a further wash with 1× PBS, coverslips were then incubated in 1× PBS containing 6 µg/mL FITC for 30 min at 37°C. The coverslips were washed again with 1× PBS, then suspended in 1× PBS containing 5 µg/mL DAPI and placed in the dark for 10 min. Coverslips were lastly washed twice more in 1× PBS, then were mounted on microscope slides. Strains were visualised on a Nikon Ni-E Fluorescence microscope using fluorescent filter cubes for DAPI and FITC, and a ×100 oil objective lens. Images were overlaid using Nikon NIS-Elements software.
Experimental evolution
Request a detailed protocolAll evolution experiments were performed at 37°C. Original bacterial strains were retrieved from –80°C stocks and grown overnight in LB (1 mg/mL dNs were added to media for growth of ΔRNR lines). Cultures were then washed twice in 1× PBS and 50 μL equivalent was added to fresh six-well plates containing 5 mL 1× MOPS media (1% glucose) supplemented with the dN concentration required for each condition (1 mg/mL, 0.25 mg/mL, 0.01 mg/mL). Five lineages of ΔRNR and 3 for WT (E. coli REL606; Lenski et al., 1991; Jeong et al., 2009) were passaged for each condition until transfer 30. Line 1 of ΔRNR grown at 0.25 mg/mL (ΔRNR_T30_250_1) was serially passaged as 10 independent lines for an additional 10 transfers in MOPS + 1% glucose and 0.01 mg/mL dNs (ΔRNR_T40_10_1–10). Plates were left to reach stationary with agitation at 120 rpm. Approximately 107–108 of wash cells were transferred to a fresh well. A glycerol stock was created for each line every five transfers. PCR contamination checks were performed every five transfers.
Sequencing and genome assembly
Request a detailed protocolAll strains and lineages required for sequencing were streaked for single colonies on LB media, with a single colony being used to inoculate an overnight culture to give clonal-level sequencing data. Genomic material was isolated using 20 μg genomic tips (QIAGEN) as per the manufacturer’s instructions and DNA quantified using a Qubit 4.0 fluorometer (Thermo Fisher). Libraries were generated using the NEXTFLEX rapid XP DNA-seq kit (PerkinElmer, Catalogue# NOVA-5149-21) and NEXTFLEX UDI Barcodes (PerkinElmer, Catalogue# NOVA-514150). Strains were sequenced to approximately 50-fold coverage at Auckland Genomics Facility using an Illumina MiSeq Platform with 2×150 bp paired-end reads. Reads were trimmed using the BBDuk (Bushnell, 2014) 1.0 plug-in for Geneious Prime 2021.1.1 (https://www.geneious.com). Mapping to the E. coli REL606 genome reference was performed using Bowtie2 (v.2.3.2) (Langmead and Salzberg, 2012) plug-in (v.7.2.1) and visualised using default settings in Geneious.
Sequence analyses
Request a detailed protocolFor creation of the phosphopentomutase alignment (Supplementary file 3), we first ran BLASTP with default settings (E-value threshold = 0.1, BLOSUM62, Filtering: none, Gapped: yes, Hits: 1000) against the Uniprot reference genomes plus Swiss-Prot database (https://www.uniprot.org/blast/) using the DeoB protein sequence from E. coli (REL606) as query. To create our alignment dataset, we filtered results for Swiss-Prot sequences and removed duplicates. We then added the E. coli REL606 query sequence and the phosphopentomutase from B. cereus, the crystal structure of which has been determined (Panosian et al., 2011). We generated our alignment with Muscle (Edgar, 2004) using the Geneious plugin (v.3.8.425) and default settings.
Three species of bacteria were noted by Lundin et al., 2009, as completely lacking ribonucleotide reduction: B. burgdorferi, U. urealyticum, and B. aphidicola str. Cc. To screen for presence/absence of ribonucleotide reduction, phosphopentomutase, and cytidine deaminases in these species, we ran blastp searches against the nr_protein database restricted to these taxa using query sequences listed in Supplementary file 1J. As the deletions of cdd involve loss of yohK in several instances, we also screened for presence/absence of the latter.
Searches were also performed genus/strain level in NCBI genome (https://www.ncbi.nlm.nih.gov/genome) for Ureaplasma (taxid:2129; 58 genome assemblies), B. aphidicola (taxid:9; 74 genome assemblies), and Borrelia (split into two genera), Borrelia (taxid:138, 77 genome assemblies) and Borreliella (taxid:64895, 17 species, 607 genome assemblies), following a proposal by Adeolu and Gupta, 2014. Under this proposed reclassification, Borellia burgdorferi has been renamed Borreliella burgdorferi. We note that this reclassification is the subject of ongoing debate (Barbour et al., 2017; Margos et al., 2017; Gupta and Bergström, 2019; Margos et al., 2018; Margos et al., 2020). As NCBI Taxonomy lists both genera and, we use the taxonomic distinction here so as to aid replication of our analyses.
Where data indicated gene absence, we confirmed this by examining the genome sequences for assembly issues, tblastn searches, and manual inspection of predicted ORFs from the genome assembly files (https://www.ncbi.nlm.nih.gov/genome/browse/#!/overview/). For cases where gene absence was clearly due to the assembly, these were excluded from further analysis. Reciprocal blast searches (blastp, nr_protein) were performed where necessary to confirm annotation of subject sequences was correct.
Data availability
All genome data have been deposited to the Single Read Archive (https://www.ncbi.nlm.nih.gov/sra) under accessions SRX17677859-SRX17677886. The following biosamples are associated with this project: SAMN30957267-SAMN30957273. The data have been deposited with links to BioProject accession number PRJNA882995 in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/).
-
NCBI BioProjectID PRJNA882995. Experimental evolution of an Escherichia coli line lacking ribonucleotide reduction.
References
-
A phylogenomic and molecular marker based proposal for the division of the genus Borrelia into two genera: The emended genus Borrelia containing only the members of the relapsing fever Borrelia, and the genus borreliella gen nov containing the members of the Lyme disease Borrelia (Borrelia burgdorferi sensu lato complex)Antonie van Leeuwenhoek 105:1049–1072.https://doi.org/10.1007/s10482-014-0164-x
-
Evolutionary genomics: Is buchnera a bacterium or an organelle?Current Biology 10:R866–R868.https://doi.org/10.1016/s0960-9822(00)00816-2
-
Division of the genus borrelia into two genera (corresponding to lyme disease and relapsing fever groups) reflects their genetic and phenotypic distinctiveness and will lead to a better understanding of these two groups of microbesInternational Journal of Systematic and Evolutionary Microbiology 67:2058–2067.https://doi.org/10.1099/ijsem.0.001815
-
Distinction between borrelia and borreliella is more robustly supported by molecular and phenotypic characteristics than all other neighbouring prokaryotic genera: response to margos’ et al. “the genus borrelia reloaded” (plos one 13(12): e0208432)PLOS ONE 14:e0221397.https://doi.org/10.1371/journal.pone.0221397
-
Dna building blocks: Keeping control of manufactureCritical Reviews in Biochemistry and Molecular Biology 47:50–63.https://doi.org/10.3109/10409238.2011.630372
-
Biochemical retrosynthesis of 2’-deoxyribonucleosides from glucose, acetaldehyde, and a nucleobaseApplied Microbiology and Biotechnology 71:615–621.https://doi.org/10.1007/s00253-005-0205-5
-
Escherichia coli can eat DNA as an excellent nitrogen source to grow quicklyFrontiers in Microbiology 13:894849.https://doi.org/10.3389/fmicb.2022.894849
-
Genome sequences of Escherichia coli B strains REL606 and BL21(DE3)Journal of Molecular Biology 394:644–652.https://doi.org/10.1016/j.jmb.2009.09.052
-
Inability of Escherichia coli B to incorporate added deoxycytidine, deoxyandenosine, and deoxyguanosine into DNAEuropean Journal of Biochemistry 17:68–71.https://doi.org/10.1111/j.1432-1033.1970.tb01135.x
-
Fast gapped-read alignment with bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Long-Term experimental evolution in Escherichia coli. I. adaptation and divergence during 2,000 generationsThe American Naturalist 138:1315–1341.https://doi.org/10.1086/285289
-
There is inadequate evidence to support the division of the genus BorreliaInternational Journal of Systematic and Evolutionary Microbiology 67:1081–1084.https://doi.org/10.1099/ijsem.0.001717
-
Comment on: gupta, 2019, distinction between borrelia and borreliella is more robustly supported by molecular and phenotypic characteristics than all other neighbouring prokaryotic genera: response to margos’ et al. “the genus borrelia reloaded” (plos one 13(12): e0208432). plos one 14(8):e0221397Ticks and Tick-Borne Diseases 11:101320.https://doi.org/10.1016/j.ttbdis.2019.101320
-
Large linear plasmids of Borrelia species that cause relapsing feverJournal of Bacteriology 195:3629–3639.https://doi.org/10.1128/JB.00347-13
-
Expression analysis of THE nrdhief operon from Escherichia coli conditions that trigger the transcript level in vivoThe Journal of Biological Chemistry 276:18031–18037.https://doi.org/10.1074/jbc.M011728200
-
Culture medium for enterobacteriaJournal of Bacteriology 119:736–747.https://doi.org/10.1128/jb.119.3.736-747.1974
-
Ribonucleotide reductasesAnnual Review of Biochemistry 75:681–706.https://doi.org/10.1146/annurev.biochem.75.103004.142443
-
Microbial production of 2-deoxyribose 5-phosphate from acetaldehyde and triosephosphate for the synthesis of 2’-deoxyribonucleosidesBioscience, Biotechnology, and Biochemistry 67:933–936.https://doi.org/10.1271/bbb.67.933
-
Does sequence conservation provide evidence for biological function?Trends in Microbiology 25:11–18.https://doi.org/10.1016/j.tim.2016.09.010
-
Bacillus cereus phosphopentomutase is an alkaline phosphatase family member that exhibits an altered entry point into the catalytic cycleThe Journal of Biological Chemistry 286:8043–8054.https://doi.org/10.1074/jbc.M110.201350
-
The evolution of the ribonucleotide reductases: much ado about oxygenJournal of Molecular Evolution 55:180–196.https://doi.org/10.1007/s00239-002-2315-3
-
The case for an early biological origin of DNAJournal of Molecular Evolution 79:204–212.https://doi.org/10.1007/s00239-014-9656-6
-
Enzymatic synthesis and breakdown of desoxyribose phosphateJournal of Biological Chemistry 196:347–365.https://doi.org/10.1016/S0021-9258(18)55739-2
-
The use of DRAQ5 to monitor intracellular DNA in Escherichia coli by flow cytometryJournal of Fluorescence 20:907–914.https://doi.org/10.1007/s10895-010-0636-y
-
General calibration of microbial growth in microplate readersScientific Reports 6:38828.https://doi.org/10.1038/srep38828
-
Ribonucleotide reductases: divergent evolution of an ancient enzymeJournal of Molecular Evolution 55:138–152.https://doi.org/10.1007/s00239-002-2311-7
-
Ribonucleotide reductases: essential enzymes for bacterial lifeFrontiers in Cellular and Infection Microbiology 4:52.https://doi.org/10.3389/fcimb.2014.00052
-
Function and evolution of plasmid-borne genes for pyrimidine biosynthesis in Borrelia sppJournal of Bacteriology 188:909–918.https://doi.org/10.1128/JB.188.3.909-918.2006
Decision letter
-
Amie K BoalReviewing Editor; Pennsylvania State University, United States
-
Christian R LandrySenior Editor; Université Laval, Canada
-
Anders HoferReviewer; Umeå University, Sweden
Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.
Decision letter after peer review:
Thank you for submitting your article "Characterisation of an Escherichia coli line that completely lacks ribonucleotide reduction yields insights into the evolution of obligate intracellularity." for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Christian Landry as the Senior Editor. The following individual involved in the review of your submission has agreed to reveal their identity: Anders Hofer (Reviewer #2).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1) Revise the introductory text about cofactor usage by class I RNRs recommended by reviewer #1. Also, modify the discussion of pathogens as detailed by reviewer #2 and comment on the selection of the starting strain as suggested by reviewer #3.
2) Make the changes to Figure 9 that are described by reviewer #2.
3) Comment further on how observation of cell elongation would impact OD measurements (reviewer #1 comment) and address the concerns of reviewer #3 about cell morphology changes and analysis of evolutionary trajectory.
4) Modify the abstract and introduction to include a statement about the need to introduce a nucleotide kinase gene.
5) Describe the control experiment to identify selective pressures induced by the antibiotic.
6) Clarify whether population-level sequencing experiments were performed as detailed by reviewer #3.
7) Revise Figure 1 to include the entire gel and include more effective annotations.
Reviewer #1 (Recommendations for the authors):
Line 72-73 – revise terminology referring to "an oxygen-containing metal ion centre" to describe Cys radical generation by class I RNRs. The Cys oxidant cofactor employed by class I enzymes does not contain oxygen. Class I RNR cofactors are assembled from either oxygen or an oxygen-derived oxidant (superoxide). But the important part about the oxygen dependence during assembly is not the incorporation of oxygen per se – it has not been established whether that happens in every case anyway – but the fact that oxygen provides oxidizing equivalents that can convert the precursor form of the cofactor (Mn(II), Fe(II), etc) to an oxidant with sufficient potency to oxidize a Cys side chain. Also, class Ie RNRs do not contain metal in their activated form – so it is also not accurate to say that all class I RNRs use a metal center to accomplish Cys oxidation.
Line 284 – there are a number of more recent studies that support the idea that class Ib RNR is used under conditions of iron limitation in E. coli.
Line 643 – missing word.
Reviewer #2 (Recommendations for the authors):
The broad questions are already described in the Public Review section (major findings and possible future experiments). Generally, I think the findings are very interesting but I have a few things that I would like to be corrected.
I find it a bit misleading that it is mentioned in the manuscript that the few organisms that lack RNR are obligately intracellular and Borrelia burgdorferi is mentioned as one of them. B. burgdorferi is mainly extracellular in the mammalian host, although it has been suggested that it can have intracellular stages. Another extracellular pathogen that lacks RNR is Giardia (several species). This is a protozoan parasite (eukaryote) that is extracellular. I think this organism should be mentioned as well in the manuscript. I am also surprised that it is stated that Mycoplasma is known to be dependent on salvage synthesis for dNTP production (row 300). They actually encode a novel subclass of class I RNRs that was published a few years ago (DOPA radical-dependent RNRs). Remove the sentence or explain this part further. The statement of lacking RNR is valid for Ureaplasma but not Mycoplasma.
I also think some adjustments need to be made in Figure 9. First of all, I think it would be good to use "ribonucleotides" and "deoxyribonucleotides" instead of NTPs and dNTPs in the top part of the figure showing the RNR reaction. This will both solve the problem that many RNRs use diphosphates as substrates and also that dNKs phosphorylate deoxyribonucleosides to the monophosphate level and not the triphosphate level. Another problem with the figure is that PPase normally means pyrophosphatase and because of that I think it is better to write "nucleoside phosphorylase" or alternatively just "phosphorylase" if the space is limited. A third problem is that it looks like the cells are taking up dNTPs, which is usually not the case, at least not in prokaryotes that do not have endocytosis. The uptake arrow from the outside should instead lead to dNs.
It should be stated in the abstract and introduction that they needed to introduce an exogenous deoxynucleoside kinase gene to the E. coli cells in order to get the system to work (E. coli can normally only phosphorylate thymidine). I would also like them to discuss more what is known about the substrate specificity of the kinase and what transporters E. coli has for the uptake of the deoxynucleosides.
Reviewer #3 (Recommendations for the authors):
Specific questions/comments for the authors:
– A general question, why did the authors choose REL606 as opposed to other strains? Is this due to evolutionary experiment convenience? If so, please briefly discuss this.
– Please clarify if a clonal and/or population-level sequence was performed, and how.
– What is the control experiment (i.e. a strain evolved in parallel) to detect any selective pressures that could be introduced by the strep in the evolution media?
– What is the significance of the cell morphology analyses for the evolutionary trajectory hypotheses presented in the study? An earlier justification or a better-justified discussion would be helpful here. To my knowledge, REL606 does not particularly exhibit a morphologically distinct response to many environments the bug has been tested/evolved before, or is this not the case? If so, then why is this a key observation attributable to the hypothesis tested here?
– Is there a correlation between the number of mutations vs. the morphological modifications of the evolved cell?
– Please clarify what is meant by this sentence (line 414): Note however that elongate cell morphology and clumping preclude use of OD measurements for accurate estimation of doubling time or cell counts (reflected in the large standard error relative to REL606 controls), so can only be used to give a general indication of growth.
– I'm intrigued by the cross-contamination speculation (and lack thereof) for the T30.L1 line. Did the authors perform population-level sequencing? Shouldn't that reveal whether early contamination took place?
– As a general comment, the second half of the Results section is rather difficult to follow. Some arguments are emphasized in this part (especially the loss of deoB and cdd genes section) and others can move to Discussion; there is a lot of back/forthing going on in this section, making it harder to grasp the sequence of events. Similarly, some paragraphs emphasize results are in the Discussion section (see, for instance, section starting line 653).
– Please include what blue/red lines correspond to on Figures 2 and 3 directly, in addition to the legend.
– Figure 1 should include the entire gel, not a trimmed version. Include an arrow pointing to exactly what we are looking at here.
– Overall summary figure (indicating the "Sequence" would be a good place to specify what you mean by sequence: and how (clonal, genome, how many replicates, etc.)).
– Specify what you mean by length (uM) on the figure in addition to the caption. The colors here do not indicate the specific strain on the Figure – is red the WT? Blue evolved WT? Unclear.
https://doi.org/10.7554/eLife.83845.sa1Author response
Essential revisions:
1) Revise the introductory text about cofactor usage by class I RNRs recommended by reviewer #1. Also, modify the discussion of pathogens as detailed by reviewer #2 and comment on the selection of the starting strain as suggested by reviewer #3.
We have now reworded the introductory text for accuracy and added references to recent papers reporting on the advances that reviewer #1 notes (ll 72-76). We also added references regarding the recent studies showing a likely role for Ib in replication during iron starvation (ll 295-296). Thanks for alerting us to these.
Thanks to Reviewer #2 for noting some inaccuracies in our descriptions of the key lineages that do not have RNRs. We have corrected these errors. We had originally omitted the eukaryote cases as our work is focused on bacteria, but agree it is helpful to mention these, so we have now added these in for completeness (ll 81-89).
We have also corrected the statement regarding the presence of ribonucleotide reductase in Mycoplasma (l 312). We now cite the discovery of DOPA radical-dependent RNRs in the introductory text. (ll 73-74)
Regarding the choice of starting strain, we did indeed choose this because of some of the convenient features of this strain for experimental evolution, including genomically coded antibiotic resistance and presence of a neutral SNP in a gene involved in arabinose utilisation, which makes this lineage suitable for competition experiments. We made use of the former in the project, but not the latter.
We have now made a statement regarding choice of strain under the methods subsection, ‘Strains and growth conditions’. (ll 125-131)
2) Make the changes to Figure 9 that are described by reviewer #2.
Figure 9 has been updated for clarity, following advice from Reviewer #2. We have replaced NTPs/dNTPs with ribonucleotides and deoxyribonucleotides, respectively. Thanks for this suggestion. While we and others have used the abbreviation PPase elsewhere, we do appreciate the point, and have altered this to remove any ambiguity. PPase is now replaced with “nucleoside phosphorylase". The arrow from outside the cell has been corrected to lead to the deoxyribonucleotides label inside the cell, and the cell boundary has been redrawn for clarity. We have also updated the figure legend (ll 834-845).
3) Comment further on how observation of cell elongation would impact OD measurements (reviewer #1 comment) and address the concerns of reviewer #3 about cell morphology changes and analysis of evolutionary trajectory.
We mentioned this in the abstract, but we have now added brief mention of this in the final paragraph of the introduction, where we summarise key findings (ll 109-111).
For establishing if there is either growth or no growth under various conditions, as we have done, a qualitative assessment such as the one presented in Figure 3 is sufficient. The issue of whether OD is impacted by cell elongation has been documented previously (Stevenson et al. https://www.nature.com/articles/srep38828), and is only a problem if trying to quantify parameters such as doubling time or when trying to estimate cell counts. We do not do either of these, as calculation of both requires an assumption of normal cell morphology in E. coli. We have added a note to clarify this in the first paragraph of the Discussion section, as per the suggestion from Reviewer #1. (ll 606-613)
4) Modify the abstract and introduction to include a statement about the need to introduce a nucleotide kinase gene.
We have now added statements to this effect in the abstract (ll 38-39) and in the final paragraph of the introduction (ll 94-98).
5) Describe the control experiment to identify selective pressures induced by the antibiotic.
Both the control and the experimental lines are resistant to streptomycin. The basis of the resistance is a genomic mutation in the rpsL gene which occurred following selection for resistance in the creation of REL606 (see Studier et al. 2009; https://doi.org/10.1016/j.jmb.2009.09.021). As streptomycin is included in the growth media of both the control and experimental lines, it should not affect interpretation of results.
If additional mutations do appear as a result of antibiotic resistance or some other form of media adaptation, they would be expected to appear in both control and experimental lines as both are evolved in the same antibiotic-containing growth media. Our genome sequencing indicated only one gene is mutated in both control and experimental lines. This gene, kdbD, encodes a putative histidine kinase. Across the control and experimental lines, we observe a range of mutations, including synonymous and nonsynonymous substitutions, plus one apparent frameshift in kdbD. Mutation of this gene cannot be specifically linked to any selection pressures associated with antibiotic resistance, and we note that this gene is not well studied. However, we plan to revisit this matter in a follow-up paper where we plan to present detailed population-level genomic data (see response to point 6, below).
Given that antibiotic presence is a constant across all lines, and given that we do not see evidence of specific adaptation to this antibiotic (as the lines are already preadapted to its presence) we do not believe it is necessary to discuss this, so we have not made any changes to the manuscript.
6) Clarify whether population-level sequencing experiments were performed as detailed by reviewer #3.
As noted in the Methods section (under Sequencing and genome assembly), genome data are derived from single colonies (ll 235-237). We have since performed population-level sequencing which we will present as part of a subsequent study, where we have extended the evolution experiment through to 100 transfers, and examined population dynamics.
7) Revise Figure 1 to include the entire gel and include more effective annotations.
Figure now edited and figure legend updated.
Reviewer #1 (Recommendations for the authors):
Line 643 – missing word.
Corrected, thanks.
Reviewer #3 (Recommendations for the authors):
– Please clarify if a clonal and/or population-level sequence was performed, and how.
Clonal – see above.
– What is the control experiment (i.e. a strain evolved in parallel) to detect any selective pressures that could be introduced by the strep in the evolution media?
Note that a control experiment is presented, but not for Streptomycin, for reasons noted above.
– What is the significance of the cell morphology analyses for the evolutionary trajectory hypotheses presented in the study? An earlier justification or a better-justified discussion would be helpful here. To my knowledge, REL606 does not particularly exhibit a morphologically distinct response to many environments the bug has been tested/evolved before, or is this not the case? If so, then why is this a key observation attributable to the hypothesis tested here?
See above. The cell morphology is reversible, so it is not significant in evolutionary terms. As now explained in the revised text, this type of morphology has been seen in a range of conditions, and is associated with stress. We show it is specifically a response to the environment (i.e. low [dN]) and is alleviated when cells are grown in higher [dN]. Perhaps what is interesting is that this environmental phenotype appears to be less severe as the experiment proceeds. This might be due to improved capacity to utilise dNs, perhaps via reduction of loss to catabolism following mutation of deoB. We now make this suggestion in the discussion (l l613-617).
– Is there a correlation between the number of mutations vs. the morphological modifications of the evolved cell?
No. We show that the morphology is environmental, not genetic – our cell lines lose this phenotype when grown in higher [dNs], though note above comment about mutations that might influence the availability of dNs for DNA synthesis.
– Please clarify what is meant by this sentence (line 414): Note however that elongate cell morphology and clumping preclude use of OD measurements for accurate estimation of doubling time or cell counts (reflected in the large standard error relative to REL606 controls), so can only be used to give a general indication of growth.
This is now discussed, as noted above. We cite a key paper (Stevenson et al.) which explains why these types of cellular morphology/behaviour preclude use of OD measurements for estimating doubling time/cell counts.
– I'm intrigued by the cross-contamination speculation (and lack thereof) for the T30.L1 line. Did the authors perform population-level sequencing? Shouldn't that reveal whether early contamination took place?
We have not performed population sequencing for the present work. However, we are currently undertaking this for a follow-up study which takes our evolution experiment through to 100 transfers. As the reviewer notes, with population data, we should indeed be able to confirm the status of the SNP in deoB in lines 1 and 8, and we look forward to presenting that work in our next article.
– As a general comment, the second half of the Results section is rather difficult to follow. Some arguments are emphasized in this part (especially the loss of deoB and cdd genes section) and others can move to Discussion; there is a lot of back/forthing going on in this section, making it harder to grasp the sequence of events. Similarly, some paragraphs emphasize results are in the Discussion section (see, for instance, section starting line 653).
We presume that, by the ‘second half’ of the results, the reviewer is referring to the bioinformatics-based analyses presented. We presented the results of screens for several genes for three genera, and elected to present results by genus, rather than by gene. We believe this is the more logical way to follow for understanding the results in biological context, but it does inevitably result in some back and forth, since we are discussing the same gene sets in several genera.
Regarding the results from line 653, we felt it useful to reemphasise the key results, but we have now added some further discussion around uptake of deoxyribonucleosides, which we hope puts this in clearer context. This is at ll 680-682 in the revised manuscript.
– Please include what blue/red lines correspond to on Figures 2 and 3 directly, in addition to the legend.
Corrected.
– Figure 1 should include the entire gel, not a trimmed version. Include an arrow pointing to exactly what we are looking at here.
Corrected.
– Overall summary figure (indicating the "Sequence" would be a good place to specify what you mean by sequence: and how (clonal, genome, how many replicates, etc.)).
Corrected. The legend to Figure 5 now reads:
“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 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”.
We have also indicated this in methods and in the main text, and we have also altered the text in the figure itself to make it clear we sequenced from single colonies.
– Specify what you mean by length (uM) on the figure in addition to the caption. The colors here do not indicate the specific strain on the Figure – is red the WT? Blue evolved WT? Unclear.
Corrected. Figure 6 has been updated to include the axis label to read “cell length (uM)”. The missing colour indicators in the key have been updated.
https://doi.org/10.7554/eLife.83845.sa2Article and author information
Author details
Funding
Royal Society Te Apārangi (17-UOA-257)
- Jun Ogawa
- Anthony M Poole
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Auckland Genomics for technical support with genome sequencing, and the Faculty of Science Imaging Centre (University of Auckland) for image analysis.
This work was supported by a Royal Society of New Zealand Marsden Fund grant to AMP and JO (17-UOA-257).
Senior Editor
- Christian R Landry, Université Laval, Canada
Reviewing Editor
- Amie K Boal, Pennsylvania State University, United States
Reviewer
- Anders Hofer, Umeå University, Sweden
Publication history
- Preprint posted: September 30, 2022 (view preprint)
- Received: September 30, 2022
- Accepted: April 5, 2023
- Accepted Manuscript published: April 6, 2023 (version 1)
- Version of Record published: April 21, 2023 (version 2)
Copyright
© 2023, Arras et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 508
- Page views
-
- 84
- Downloads
-
- 0
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Genetics and Genomics
Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated ‘elusive genes.’ These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.
-
- Evolutionary Biology
- Plant Biology
While the domestication process has been investigated in many crops, the detailed route of cultivation range expansion and factors governing this process received relatively little attention. Here using mungbean (Vigna radiata var. radiata) as a test case, we investigated the genomes of more than one thousand accessions to illustrate climatic adaptation’s role in dictating the unique routes of cultivation range expansion. Despite the geographical proximity between South and Central Asia, genetic evidence suggests mungbean cultivation first spread from South Asia to Southeast, East, and finally reached Central Asia. Combining evidence from demographic inference, climatic niche modeling, plant morphology, and records from ancient Chinese sources, we showed that the specific route was shaped by the unique combinations of climatic constraints and farmer practices across Asia, which imposed divergent selection favoring higher yield in the south but short-season and more drought-tolerant accessions in the north. Our results suggest that mungbean did not radiate from the domestication center as expected purely under human activity, but instead the spread of mungbean cultivation is highly constrained by climatic adaptation, echoing the idea that human commensals are more difficult to spread through the south-north axis of continents.