Skip to main content


Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

A substrate ambiguous enzyme facilitates genome reduction in an intracellular symbiont



Genome evolution in intracellular microbial symbionts is characterized by gene loss, generating some of the smallest and most gene-poor genomes known. As a result of gene loss these genomes commonly contain metabolic pathways that are fragmented relative to their free-living relatives. The evolutionary retention of fragmented metabolic pathways in the gene-poor genomes of endosymbionts suggests that they are functional. However, it is not always clear how they maintain functionality. To date, the fragmented metabolic pathways of endosymbionts have been shown to maintain functionality through complementation by host genes, complementation by genes of another endosymbiont and complementation by genes in host genomes that have been horizontally acquired from a microbial source that is not the endosymbiont. Here, we demonstrate a fourth mechanism.


We investigate the evolutionary retention of a fragmented pathway for the essential nutrient pantothenate (vitamin B5) in the pea aphid, Acyrthosiphon pisum endosymbiosis with Buchnera aphidicola. Using quantitative analysis of gene expression we present evidence for complementation of the Buchnera pantothenate biosynthesis pathway by host genes. Further, using complementation assays in an Escherichia coli mutant we demonstrate functional replacement of a pantothenate biosynthesis enzyme, 2-dehydropantoate 2-reductase (E.C., by an endosymbiont gene, ilvC, encoding a substrate ambiguous enzyme.


Earlier studies have speculated that missing enzyme steps in fragmented endosymbiont metabolic pathways are completed by adaptable endosymbiont enzymes from other pathways. Here, we experimentally demonstrate completion of a fragmented endosymbiont vitamin biosynthesis pathway by recruitment of a substrate ambiguous enzyme from another pathway. In addition, this work extends host/symbiont metabolic collaboration in the aphid/Buchnera symbiosis from amino acid metabolism to include vitamin biosynthesis.


Extensive genome reduction in obligate endosymbiotic bacteria is a hallmark of intracellularity [1]-[3]. Recent genome sequencing of diverse bacterial endosymbionts of insects reveals ongoing symbiont genome reduction, resulting in small gene dense genomes that contain a subset of genes compared to their closest free-living relatives [4]-[7]. For example, the genome of the gammaproteobacterium Buchnera from the aphid, Cinara cedri, is only 416 kbp, approximately 11-fold smaller that its closest free-living relative Escherichia coli [4],[8]. Currently, the smallest known bacterial genomes, containing fewer than 200 protein-coding genes, belong to insect associated obligate endosymbionts [6],[7],[9],[10]. Eliminated genes include those that are either functionally redundant, or not essential for maintenance of the symbiosis [1]-[3]. Many of the genes evolutionarily retained in endosymbiont genomes are those required for biosynthesis of nutrients essential to the host. However, even endosymbiont essential nutrient biosynthesis pathways experience gene loss such that it is common to find fragmented, but evolutionary retained, pathways that appear incomplete and non-functional [11],[12].

Table 1 Pantothenate biosynthesis genes and enzyme activities

Recent genome sequencing of holosymbionts (hosts + their symbionts) reveals that fragmented biosynthesis pathways can be completed by three processes: (1) complementation by host genes with overlapping function [13]-[15]; (2) complementation by genes encoded by another symbiont [4]-[6]; and (3) complementation by genes horizontally transferred to the host nuclear genome from other bacterial genomes [11],[15]. All three of these processes can contribute to functional maintenance of metabolic pathways in a single endosymbiont and yet, even accounting for these three independent processes, some evolutionarily retained symbiont pathways remain fragmented making it unclear if they are functional [11],[12].

A fourth process for completion of fragmented biosynthesis pathways was proposed by Shigenobu et al. [12] when they hypothesized that missing enzyme steps in fragmented Buchnera amino acid biosynthesis pathways are mediated by other (unidentified) Buchnera enzymes with overlapping activity [12],[14]. In addition to well-defined primary roles within linear metabolic pathways, it is common for enzymes to catalyze side reactions [16]-[18]. Such multifunctional enzymes that are capable of catalyzing weak side reactions, exhibit either catalytic promiscuity, catalysis of a different reaction type on a single or range of substrates; or substrate ambiguity, catalysis of the same reaction type on a range of substrates [16]-[18]. Extraordinary levels of metabolic plasticity within organisms are facilitated by weak side reactions [19]-[24]. Importantly, weak side activities provide precursor catalytic activity for the evolution of new enzyme function [16]-[18],[21],[22]. It is not without precedent that multifunctional enzymes with overlapping activity ‘fill holes’ in fragmented essential metabolic pathways [25]-[27]. For example, the fragmented folate biosynthesis pathway from the obligate intracellular pathogen Chlamydia is functional because of the recruitment of promiscuous and adaptable enzymes from other biosynthetic pathways [25].

Here, working with the best-characterized holosymbiosis, that of the pea aphid Acyrthosiphon pisum with the gammaproteobacterium Buchnera, we reveal that the fragmented Buchnera pantothenate (vitamin B5) biosynthesis pathway can be completed using a Buchnera substrate ambiguous enzyme recruited from another pathway. We demonstrate using complementation assays in E. coli mutants that retention of a substrate ambiguous enzyme in the massively reduced genome of Buchnera facilitates the collaborative biosynthesis of pantothenate (vitamin B5) in the A. pisum/Buchnera holosymbiont.

Results and discussion

Buchnerahas a fragmented but evolutionary retained pantothenate biosynthesis pathway

Compared to E. coli, the obligate endosymbiont Buchnera from the pea aphid, A. pisum, has a fragmented pantothenate biosynthesis pathway that is missing panD, panE and ilvE (Figure 1). Gene identifiers and enzyme activities are provided in Table 1. Buchnera with fully sequenced genomes fall into two groups: (1) those from species such as Baizongia pistaciae, C. cedri and Uroleucon ambrosiae that have lost the whole pantothenate biosynthesis pathway [4],[28],[29]; and (2) those from species such as A. pisum, Myzus persicae and Schizaphis graminum that retain identically fragmented pathways characterized by loss of panD, panE and ilvE [12],[30],[31] (Figure 1). Evolutionary retention by multiple Buchnera lineages of identically fragmented pantothenate biosynthesis pathways that include two genes (panB and panC) that only function in pantothenate biosynthesis, led us to hypothesize that the pantothenate biosynthesis pathway in these species is functional.

Figure 1

Comparative analysis of the Buchnera pantothenate biosynthesis pathway and complementary A. pisum β-alanine biosynthesis pathways. Metabolic pathways were reconstructed using KEGG pantothenate biosynthesis pathway from E. coli K-12 (eco00770) and from Buchnera from host aphids that included: A. pisum (Ap, buc00770), B. pistaciae (Bp, bab00770), C. cedri (Cc, bcc00770), M. persicae (Mp, [30]), S. graminum (Sg, bas00770) and U. ambrosiae (Ua, bua00770) [34],[35]. Complementary β-alanine biosynthesis pathways of A. pisum (api00410) [34],[35] are shown. The presence/absence of biosynthesis genes in the Buchnera, E. coli and A. pisum genomes is indicated. Full details of A. pisum, Buchnera and E. coli gene identifiers and enzyme activities are provided in Table 1. KEGG, Kyoto Encyclopedia of Genes and Genomes.

To begin the investigation of our hypothesis that evolutionary retention of identically fragmented pantothenate biosynthesis pathways indicates retention of pantothenate biosynthesis capacity we characterized expression of ilvH, ilvC, ilvD, panB and panC in the specialized aphid bacteriocyte cells that house Buchnera. All five of these Buchnera genes were highly expressed in A. pisum bacteriocytes and consistently highly expressed across three genetically independent A. pisum lines (Figure 2A). Functionality of these fragmented pathways requires (1) host supply of β-alanine and (2) complementation of panE encoded 2-dehydropantoate 2-reductase (E.C.

Figure 2

Quantitative expression analyses of Buchnera and A. pisum pantothenoate biosynthesis genes. Expression profiles for (A) Buchnera pantothenate biosynthesis genes and (B) complementary A. pisum β-alanine biosynthesis genes were generated for three independent A. pisum lineages (5A, 9-2-1 and LSR1). Expression of Buchnera genes in bacteriocytes (Ba) was normalized to atpE (BU003) while expression of A. pisum genes in whole insects (WI), bacteriocyte (Ba) and gut (G) tissues was normalized to GAPDH (ACYPI009769). Relative abundance of each gene was compiled into a heat map showing high expression (yellow) to low expression (blue). The box shows gene expression fold-change in A. pisum LSR1 bacteriocyte (Ba) tissues relative to whole insect (WI) expression. Fold-change was calculated according to [36]; for each aphid lineage n = 3.

A. pisumsynthesizes β-alanine via the uracil degradation pathway in bacteriocytes

A. pisum and other insects synthesize β-alanine via two main pathways: (1) the uracil degradation pathway [32] and (2) via decarboxylation of L-aspartate [33] (Figure 1). Gene identifiers and enzyme activities are provided in Table 1. While synthesis of β-alanine is evolutionarily conserved in the major insect orders (Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway map 00410), we propose that β-alanine biosynthesis in A. pisum is now additionally required to supply Buchnera.

Gene expression analysis demonstrates that β-alanine biosynthesis via the A. pisum uracil degradation pathway is upregulated in bacteriocytes relative to whole insect expression (Figure 2B). Consistent with uracil degradation being the major β-alanine biosynthesis route in bacteriocytes, the alternate pathway that synthesizes β-alanine via decarboxylation of L-aspartate is down regulated in bacteriocytes (Figure 2B). Aspartate decarboxylase (ADC) is 4.3-fold down in LSR1 bacteriocytes, 25-fold down in 5A bacteriocytes and 100-fold down in 9-2-1 bacteriocytes relative to whole insect expression.

Buchnera ilvC complements panE to complete the pantothenate biosynthesis pathway in A. pisum

Given host supply of β-alanine, functional pantothenate biosynthesis in Buchnera with the fragmented pantothenate biosynthesis pathway requires complementation of lost panE encoded 2-dehydropantoate 2-reductase (E.C. Previous work in Salmonella typhimurium [40] and E. coli [41] has shown that ilvC, a gene encoding ketol-acid reductoisomerase (E.C. that mediates reactions in the biosynthesis of the essential amino acids isoleucine and valine and production of 2-oxoisovalerate from pyruvate in pantothenate biosynthesis, also has limited 2-dehydropantoate 2-reductase (E.C. activity. While ilvC can reduce 2-dehydropantoate to pantoate, the activity of ilvC in S. typhimurium to function in pantothenate biosynthesis is 20-fold lower than its activity as ketol-acid reductoisomerase [40] and so low in E. coli that the ilvC knockout mutant was indistinguishable from wild-type K12 in its reduction of 2-dehydropantoate to pantoate [41]. However, given that generation of E. coli pantoate auxotrophs requires generation of panE, ilvC double mutants, we hypothesized that Buchnera ilvC has sufficient 2-dehydropantoate 2-reductase activity and, thereby, critically maintains functional pantothenate biosynthesis in Buchnera that retain fragmented pantothenate biosysnthesis pathways.

To test the hypothesis that Buchnera ilvC has 2-dehydropantoate 2-reductase activity sufficient to complete pantoate synthesis we utilized a 2-dehydropantoate 2-reductase negative strain of E. coli (strain FE8: ilvC - and panE - [41]) and functionally complemented E. coli FE8 with Buchnera ilvC (Figure 3A). When cells were plated on selective plates (minimal media plates, supplemented with L-isoleucine and 3-methyl-2-oxobutanoic acid to complement deletion of ilvC but without pantothenate supplementation), rescue of E. coli FE8 was only achieved by complementation with E. coli panE and Buchnera ilvC (Figure 3A). In contrast, under selective conditions, very limited growth of FE8 cells transformed with E. coli ilvC was observed, suggesting that E. coli ilvC does not contain sufficient 2-dehydropantoate 2-reductase activity to rescue 2-dehydropantoate 2-reductase negative E. coli cells, a result consistent with the earlier work of Elischewski et al., [41]. As a negative control, E. coli FE8 was transformed with pUC19 empty expression vector. As expected, these pUC19 empty vector transformed cells did not grow without pantothenate supplementation (Figure 3A). When grown in non-selective conditions (minimal media plates with L-isoleucine, 3-methyl-2-oxobutanoic acid and pantothenate supplementation) FE8 cells transformed with E. coli panE, Buchnera ilvC, E. coli ilvC and pUC19 empty expression vector (negative control) established robust colonies (Figure 3A).

Figure 3

Functional complementation of E. coli pantothenate biosynthesis auxotrophs with Buchnera pantothenate biosynthesis genes. (A) E. coli FE8 cells (ilvC - and panE -) were transformed with expression plasmids containing E. coli panE, Buchnera ilvC, E. coli ilvC and negative control pUC19 empty vector (- control). For complementation assays 5 μl of cells (from OD 600 = 1.0) was plated on minimal media plates with pantothenate (+B5, non-selective) and without pantothenate (-B5, selective). Recovery of E. coli mutant growth was assessed after four days at 26°C. (B) E. coli AT1371 cells (panC -) were transformed with expression plasmids containing Buchnera panC, E. coli panC and negative control pUC19 empty vector (- control). Complementation assays were completed as described in (A).

Functional complementation assays were also completed to test the pantoate-β-alanine ligase (E.C. activity of Buchnera panC. Buchnera panC shares only 48% amino acid sequence identity with E. coli panC. Despite this high level of sequence divergence, Buchnera panC retains pantoate-β-alanine ligase activity (Figure 3B). Both Buchnera APS panC and E. coli panC (positive control) functionally rescue E. coli AT1371 (panC - strain) when grown on a minimal media without pantothenate supplementation. Negative control E. coli AT1371 transformed with pUC19 empty expression vector was unable to grow on minimal media without pantothenate supplementation. With pantothenate supplementation, growth of AT1371 transformed with Buchnera panC, E. coli panC and negative control cells is observed (Figure 3B).

Loss of Buchnera ilvEfacilitates holosymbiont integration

Loss of Buchnera ilvE (encoding branched-chain-amino-acid transaminase, E.C. is not critical to the integrity of the pantothenate biosynthesis pathway in Buchnera because Buchnera with fragmented pathways retain ilvH/I, ilvC and ilvD and are predicted to synthesize 2-oxoisovalerate from pyruvate (Figure 1). However, the loss of Buchnera ilvE has important coevolutionary implications for the biosynthesis of the branched chain amino acids (isoleucine, leucine and valine). Symbiont loss of ilvE was first reported in Buchnera APS from A. pisum [12] a decade before the genome of A. pisum was published [38]. While Shigenobu et al., [12] had speculated that Buchnera APS likely utilized the activity of another transaminase to complete branched chain amino acid biosynthesis, members of the Pea Aphid Genome Consortium offered an alternate explanation [42], that Buchnera loss of ilvE encoded branched-chain-amino-acid transaminase is complemented by the host genome encoded branched-chain-amino-acid transaminase (BCAT, ACYPI008371, [42]). Recently, the collaborative role of BCAT in the biosynthesis of leucine, one of the three branched chain amino acids, has been confirmed experimentally [43]. Identical patterns of holosymbiont coevolution involving symbiont loss of ilvE have likely also coevolved in the citrus mealybug, Planococcus citri [11], and the whitefly, Bemisia tabaci [31].

Host/symbiont metabolic collaboration and a substrate ambiguous enzyme facilitate maintenance of functional pantothenate biosynthesis pathways in some Buchneralineages

On the basis of holosymbiont metabolic pathway analysis, gene expression analysis and complementation tests in E. coli pantothenate biosynthesis mutants, we propose a model of collaborative pantothenate biosynthesis that critically depends on host supply of β-alanine, and the substrate ambiguous activity of ketol-acid reductoisomerase (E.C. encoded by ilvC as 2-dehydropantoate 2-reductase (E.C. (Figure 4).

Figure 4

Pantothenate (vitamin B5) biosynthesis in the A. pisum / Buchnera holosymbiont.

A. pisum synthesize β-alanine in bacteriocytes via the uracil degradation pathway (Figure 4). To complete pantothenate biosynthesis, β-alanine is transported from the bacteriocyte cytoplasm to Buchnera via a currently unidentified transporter (Figure 4). In Buchnera’s closest free-living relative E. coli, β-alanine is actively transported via amino acid transporter cycA [44]. However, due to extensive gene loss, Buchnera has lost almost all substrate specific transporters, including cycA [12],[45]. Our current working hypothesis, as previously proposed by Shigenobu et al., [12], is that Buchnera utilizes a flagellum-derived general substrate transport system for import of essential nutrients, including β-alanine.

Notably, we demonstrate that panE encoding 2-dehydropantoate 2-reductase has been lost from Buchnera’s genome and functionally replaced with ilvC, an enzyme that encodes both ketol-acid reductoisomerase (E.C. and 2-dehydropantoate 2-reductase (E.C. activity. The gene product of ilvC functions as a substrate ambiguous reductase in branched chain amino acid biosynthesis and pantothenate biosynthesis. Functional replacement of panE encoded ketopantoate reductase with a substrate ambiguous ilvC encoded ketopantoate reductase is not without precedent since loss of panE and functional replacement by ilvC has previously been demonstrated in the gram-positive soil bacterium Corynebacterium glutamicum [46]. In the context of holosymbiont evolution, we hypothesize that retention of multifunctional endosymbiont enzymes facilitates genome reduction while maintaining the functional integrity of essential biosynthetic pathways, a process that facilitates genome reduction but results in increased functional complexity of symbiont-retained genes. Consistent with our hypothesis, in genome-wide analyses, it has recently been demonstrated that bacterial species with small, reduced genomes have increased protein functional complexity, leading the authors to suggest that an increase in functional complexity of retained genes in small genomes compensates for gene losses [47].


There are currently four mechanisms that facilitate retention of metabolic functionality concomitant with gene loss in the primary endosymbionts of insects. All four mechanisms can operate within a single genome and all four contribute to explaining why fragmented and disconnected biosynthesis pathways, that on first observation appear incomplete and non-functional, are retained across diverse obligate endosymbiont taxa. As more holosymbiont genomes become available, metabolic reconstruction of obligate intracellular symbionts must consider a complex evolutionary patchwork that comprises host metabolic complementarity, symbiont replacement, horizontal transfer of complementary biosynthesis genes to the host genome and now also protein multifunctionality in endosymbionts, all of which contribute to the maintenance of endosymbiont metabolic function.


Aphid culture

Parthenogenetic females of the pea aphid A. pisum, lineage LSR1 [38], 5A and 9-2-1 [48] were maintained on pre-flowering Vicia faba cv. The Sutton. All aphids were kept at 20°C in a long day regime of 16 hours of light and 8 hours of dark.

Comparative analysis of vitamin B5 biosynthesis

Metabolic pathways for pantothenate (vitamin B5) biosynthesis were analyzed using the KEGG pathway database [34],[35]. These included the pantothenate biosynthesis pathways from E. coli K-12 MG1655 (eco00770) and B. aphidicola from host aphids, including: A. pisum (Buchnera Ap, buc00770 [12]), B. pistaciae (Buchnera Bp, bab00770 [28]), C. cedri (Buchnera Cc, bcc00770 [4]), M. persicae (Buchnera Mp [30]), S. graminum (Buchnera Sg, bas00770 [31]) and U. ambrosiae (Ua, bua00770 [29]). The complementary A. pisum β-alanine metabolism pathway was identified in A. pisum LSR1 (api00410 [38]).

Real-time quantitative RT-PCR

Real-time quantitative PCR (QPCR) was used to compare A. pisum and Buchnera pantothenate biosynthesis gene expression in whole adult female aphids and isolated bacteriocytes using 2-ΔΔCT methodology [36]. Methods are identical to those described in detail by [49]. Primers were designed using Primer3Plus software for real-time PCR [50], and specificity checked using Primer-BLAST against NCBI A. pisum RefSeq gene models (primer sequences are shown in Additional file 1: Table S1). For comparative gene expression analysis, all QPCR primers had amplification efficiencies >90% and <105%. A. pisum gene expression was normalized to glyceraldehyde-3-dehydrogenase (GAPDH, ACYPI009769) and Buchnera gene expression was normalized to atpE ATP synthase F0F1 subunit C (atpE, BU003). The normalized expression value for each gene was compiled into a heat map (z-score x 10; z = (each value - average)/standard deviation). Yellow: z >0, blue: z <0, black: z = 0.

ilvC, panC and panEexpression constructs

Full-length Buchnera ilvC (BU599) and panC (BU196) coding sequences were amplified from Buchnera LSR1 gDNA, and full-length E. coli ilvC (EG10495), panC (EG11746) and panE (EG13271) coding sequences were amplified from E. coli K-12 MG1655 gDNA; all coding sequences were amplified using Phusion high-fidelity DNA polymerase (Thermo Scientific, Waltham, MA, USA). Primers were designed to incorporate restriction enzyme sites to allow directional cloning of ilvC, panC and panE coding sequences into pUC19, in frame and downstream of the lacZ start codon; thus, generating pUC19 expression constructs containing ilvC, panC and panE coding sequences under the control of the lacZ promoter (primer sequences are shown in Additional file 1: Table S1). Amplification reactions were completed on a Mastercycler PCR system (Eppendorf, Hauppauge, NY, USA) with the following conditions: 98°C for 30 seconds, followed by 30 cycles at 98°C for 10 seconds, 55°C for 30 seconds and 72°C for 1 minute. Amplified coding sequences were digested with either HindIII and BamHI, or BamHI and EcoRI (see Additional file 1: Table S1) and cloned into the respective sites of pUC19. Expression constructs containing Buchnera and E. coli ilvC coding sequences were named pBUC-ilvC and pECO-ilvC, respectively. Expression constructs containing Buchnera and E. coli panC coding sequences were named pBUC-panC and pECO-panC, respectively. The expression construct containing E. coli panE coding sequence was named pECO-panE. All expression constructs were fully sequenced using Sanger sequencing in standard Big Dye terminator v3.1 reactions (Life Technologies, Green Island, NY, USA). Reaction products were analyzed on 3130xl genetic analyzer (Life Technologies), and sequence data assembled into contiguous sequences using Sequencher v4.9 (Gene Codes Corp., Ann Arbor, MI, USA).

Functional complementation of E. coliauxotrophs

Functional complement experiments were carried out using the E. coli ilvC/panE double mutant FE8 (ilvC::aacC1, panE::tetA; [41]), obtained as a kind gift from Jörn Kalinowski (Institute for Genome Research and Systems Biology, Bielefield University, Germany); and E. coli panC mutant AT1371 (F-, panC4, Δ(gpt-proA)62, lacY1, tsx-29, glnV44(AS), galK2(Oc), λ -, Rac-0, hisG4(Oc), rfbC1, xylA5, mtl-1, argE3(Oc), thiE1 [51],[52]), obtained from the Yale University E.coli Genetic Stock Center (CGSC, Yale University, New Haven, CT, USA.).

E. coli FE8 (ilvC - and panE -) cells were transformed with plasmids pBUC-ilvC, pECO-ilvC, pECO-panE and pUC19 empty vector (negative control) according to [53]. Transformants were grown overnight in Luria-Bertani (LB) media with 50 μg/ml ampicillin, 15 μg/ml gentamicin and 10 μg/ml tetracycline and washed twice in sterile distilled water and finally resuspended to an attenuance of OD 600 = 1.0. Functional complementation of E. coli FE8 pantothenate auxotrophs was carried out by plating aliquots (5 μl) of cell suspensions on M63 minimal media agar plates [54] containing 13.6 g/L KH2PO4 (pH 7.0), 2 g/L (NH4)2SO4, 0.4% (w/v) glycerol, 0.25 g/L MgSO4.7H2O, 50 mg/ml L-isoleucine, 50 mg/ml 3-methyl-2-oxobutanoic acid (Acros Organics, Waltham, MA, USA), 0.5 mg/L FeSO4.7H2O, 0.5 mg/L thiamine, 1.5% (w/v) agar, 50 μg/ml ampicillin, 15 μg/ml gentamicin and 10 μg/ml tetracycline. Aliquots of cell suspensions were also plated on M63 minimal media agar plates supplemented with 50 μg/ml D-pantothenic acid (non-selective) to assess cell viability. E. coli FE8 cell growth on selective plates (without vitamin B5) and non-selective plates (with vitamin B5) was assessed after four days of growth at 26°C.

E. coli AT1371 (panC -) cells were transformed with plasmids pBUC-panC, pECO-panC (positive control) and pUC19 empty vector (negative control) according to [53]. Transformants were grown overnight in LB media with 50 μg/ml ampicillin and washed twice in sterile distilled water and finally resuspended to an attenuance of OD600 = 1.0. Functional complementation of AT1371 pantothenate auxotrophs was carried out by plating aliquots (5 μl) of cell suspensions on M63 minimal media agar plates [54] containing 13.6 g/L KH2PO4 (pH 7.0), 2 g/L (NH4)2SO4, 0.4 % (w/v) glycerol, 0.25 g/L MgSO4.7H2O, 0.5 mg/LFeSO4.7H2O, 0.5 mg/L thiamine, 1.5% (w/v) agar and 100 μg/ml ampicillin. Aliquots of cell suspensions were also plated on M63 minimal media agar plates supplemented with 50 μg/ml D-pantothenic acid (non-selective) to assess cell viability. E. coli FE8 cell growth on selective plates (without pantothenate) and non-selective plates (with pantothenate) was assessed after four days of growth at 26°C.

Authors’ information

DRGP is a molecular biologist who uses functional genomic approaches to characterize the interaction between phloem-feeding insects and their obligate endosymbiotic bacteria. ACCW studies genome coevolution between plant sap-feeding insects and their obligate intracellular bacterial symbionts. Recent work in Wilson’s group focuses on functional characterization of amino acid transporters in sternorrhynchan insects.

Additional file


  1. 1.

    McCutcheon JP: The bacterial essence of tiny symbiont genomes. Curr Opin Microbiol. 2010, 13: 73-78. 10.1016/j.mib.2009.12.002.

  2. 2.

    Moran NA, Bennett GM: The tiniest tiny genomes. Annu Rev Microbiol. 2014, 68: 195-215. 10.1146/annurev-micro-091213-112901.

  3. 3.

    McCutcheon JP, Moran NA: Extreme genome reduction in symbiotic bacteria. Nat Rev Microbiol. 2012, 10: 13-26.

  4. 4.

    Pérez-Brocal V, Gil R, Ramos S, Lamelas A, Postigo M, Michelena JM, Silva FJ, Moya A, Latorre A: A small microbial genome: the end of a long symbiotic relationship?. Science. 2006, 314: 312-313. 10.1126/science.1130441.

  5. 5.

    Lamelas A, Gosalbes MJ, Manzano-Marín A, Peretó J, Moya A, Latorre A: Serratia symbiotica from the aphid Cinara cedri: a missing link from facultative to obligate insect endosymbiont. PLoS Genet. 2011, 7: 1-11. 10.1371/journal.pgen.1002357.

  6. 6.

    McCutcheon JP, von Dohlen CD: An interdependent metabolic patchwork in the nested symbiosis of mealybugs. Curr Biol. 2011, 21: 1366-1372. 10.1016/j.cub.2011.06.051.

  7. 7.

    Bennett GM, Moran NA: Small, smaller, smallest: the origins and evolution of ancient dual symbioses in a phloem-feeding insect. Genome Biol Evol. 2013, 5: 1675-1688. 10.1093/gbe/evt118.

  8. 8.

    Gil R, Sabater-Munoz B, Latorre A, Silva FJ, Moya A: Extreme genome reduction in Buchnera spp.: toward the minimal genome needed for symbiotic life. Proc Natl Acad Sci U S A. 2002, 99: 4454-4458. 10.1073/pnas.062067299.

  9. 9.

    McCutcheon JP, McDonald BR, Moran NA: Origin of an alternative genetic code in the extremely small and GC-rich genome of a bacterial symbiont. PLoS Genet. 2009, 5: e1000565-10.1371/journal.pgen.1000565.

  10. 10.

    Nakabachi A, Yamashita A, Toh H, Ishikawa H, Dunbar HE, Moran NA, Hattori M: The 160-kilobase genome of the bacterial endosymbiont Carsonella . Science. 2006, 314: 267-10.1126/science.1134196.

  11. 11.

    Husnik F, Nikoh N, Koga R, Ross L, Duncan RP, Fujie M, Tanaka M, Satoh N, Bachtrog D, Wilson ACC, von Dohlen CD, Fukatsu T, McCutcheon JP: Horizontal gene transfer from diverse bacteria to an insect genome enables a tripartite nested mealybug symbiosis. Cell. 2013, 153: 1567-1578. 10.1016/j.cell.2013.05.040.

  12. 12.

    Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H: Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature. 2000, 407: 81-86. 10.1038/35024074.

  13. 13.

    Hansen AK, Moran NA: Aphid genome expression reveals host-symbiont cooperation in the production of amino acids. Proc Natl Acad Sci U S A. 2011, 108: 2849-2854. 10.1073/pnas.1013465108.

  14. 14.

    Shigenobu S, Wilson ACC: Genomic revelations of a mutualism: the pea aphid and its obligate bacterial symbiont. Cell Mol Life Sci. 2011, 68: 1297-1309. 10.1007/s00018-011-0645-2.

  15. 15.

    Sloan DB, Nakabachi A, Richards S, Qu J, Murali SC, Gibbs RA, Moran NA: Parallel histories of horizontal gene transfer facilitated extreme reduction of endosymbiont genomes in sap-feeding insects. Mol Biol Evol. 2014, 31: 857-871. 10.1093/molbev/msu004.

  16. 16.

    Jensen RA: Enzyme recruitment in evolution of new function. Annu Rev Microbiol. 1976, 30: 409-425. 10.1146/annurev.mi.30.100176.002205.

  17. 17.

    Khersonsky O, Tawfik DS: Enzyme promiscuity: a mechanistic and evolutionary perspective. Annu Rev Biochem. 2010, 79: 471-505. 10.1146/annurev-biochem-030409-143718.

  18. 18.

    O’Brien PJ, Herschlag D: Catalytic promiscuity and the evolution of new enzymatic activities. Chem Biol. 1999, 6: R91-R105. 10.1016/S1074-5521(99)80033-7.

  19. 19.

    D'Ari R, Casadesus J: Underground metabolism. BioEssays. 1998, 20: 181-186. 10.1002/(SICI)1521-1878(199802)20:2<181::AID-BIES10>3.0.CO;2-0.

  20. 20.

    Carbonell P, Lecointre G, Faulon JL: Origins of specificity and promiscuity in metabolic networks. J Biol Chem. 2011, 286: 43994-44004. 10.1074/jbc.M111.274050.

  21. 21.

    Nam H, Lewis NE, Lerman JA, Lee DH, Chang RL, Kim D, Palsson BO: Network context and selection in the evolution to enzyme specificity. Science. 2012, 337: 1101-1104. 10.1126/science.1216861.

  22. 22.

    Notebaart RA, Szappanos B, Kintses B, Pál F, Györkei A, Bogos B, Lázár V, Spohn R, Csörgő B, Wagner A, Ruppin E, Pál C, Papp B: Network-level architecture and the evolutionary potential of underground metabolism. Proc Natl Acad Sci U S A. 2014, 111: 11762-11767. 10.1073/pnas.1406102111.

  23. 23.

    Yip SH, Matsumura I: Substrate ambiguous enzymes within the Escherichia coli proteome offer different evolutionary solutions to the same problem. Mol Biol Evol. 2013, 30: 2001-2012. 10.1093/molbev/mst105.

  24. 24.

    Greenspan RJ: The flexible genome. Nat Rev Genet. 2001, 2: 383-387. 10.1038/35072018.

  25. 25.

    Adams NE, Thiaville JJ, Proestos J, Juárez-Vázquez AL, McCoy AJ, Barona-Gómez F, Iwata-Reuyl D, de Crécy-Lagard V, Maurelli AT: Promiscuous and adaptable enzymes fill “holes” in the tetrahydrofolate pathway in Chlamydia species. MBio. 2014, 5: e01378-01314. 10.1128/mBio.01378-14.

  26. 26.

    Tian J, Bryk R, Shi S, Erdjument-Bromage H, Tempst P, Nathan C: Mycobacterium tuberculosis appears to lack α-ketoglutarate dehydrogenase and encodes pyruvate dehydrogenase in widely separated genes. Mol Microbiol 2005, 57:859–868.,

  27. 27.

    Wagner T, Bellinzoni M, Wehenkel A, O’Hare HM, Alzari PM: Functional plasticity and allosteric regulation of α-ketoglutarate decarboxylase in central mycobacterial metabolism. Chem Biol. 2011, 18: 1011-1020. 10.1016/j.chembiol.2011.06.004.

  28. 28.

    van Ham RC, Kamerbeek J, Palacios C, Rausell C, Abascal F, Bastolla U, Fernández JM, Jiménez L, Postigo M, Silva FJ, Tamames J, Viguera E, Latorre A, Valencia A, Morán F, Moya A: Reductive genome evolution in Buchnera aphidicola . Proc Natl Acad Sci U S A. 2003, 100: 581-586. 10.1073/pnas.0235981100.

  29. 29.

    Degnan PH, Ochman H, Moran NA: Sequence conservation and functional constraint on intergenic spacers in reduced genomes of the obligate symbiont Buchnera . PLoS Genet. 2011, 7: e1002252-10.1371/journal.pgen.1002252.

  30. 30.

    Jiang Z, Jones DH, Khuri S, Tsinoremas NF, Wyss T, Jander G, Wilson ACC: Comparative analysis of genome sequences from four strains of the Buchnera aphidicola Mp endosymbion of the green peach aphid. Myzus persicae. BMC Genomics. 2013, 14: 917-10.1186/1471-2164-14-917.

  31. 31.

    Tamas I, Klasson L, Canbäck B, Näslund AK, Eriksson AS, Wernegreen JJ, Sandström JP, Moran NA, Andersson SG: 50 million years of genomic stasis in endosymbiotic bacteria. Science. 2002, 296: 2376-2379. 10.1126/science.1071278.

  32. 32.

    Piskur J, Schnackerz KD, Andersen G, Bjornberg O: Comparative genomics reveals novel biochemical pathways. Trends Genet. 2007, 23: 369-372. 10.1016/j.tig.2007.05.007.

  33. 33.

    Richardson G, Ding H, Rocheleau T, Mayhew G, Reddy E, Han Q, Christensen BM, Li J: An examination of aspartate decarboxylase and glutamate decarboxylase activity in mosquitoes. Mol Biol Rep. 2010, 37: 3199-3205. 10.1007/s11033-009-9902-y.

  34. 34.

    Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.

  35. 35.

    Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40: D109-D114. 10.1093/nar/gkr988.

  36. 36.

    Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

  37. 37.

    Zhou J, Rudd KE: EcoGene 3.0. Nucleic Acids Res. 2013, 41: D613-D624. 10.1093/nar/gks1235.

  38. 38.

    International Aphid Genomics Consortium: Genome sequence of the pea aphid Acyrthosiphon pisum . PLoS Biol 2010, 8:e1000313.,

  39. 39.

    Schomburg I, Chang A, Placzek S, Söhngen C, Rother M, Lang M, Munaretto C, Ulas S, Stelzer M, Grote A, Scheer M, Schomburg D: BRENDA in 2013: integrated reactions, kinetic data, enzyme function data, improved disease classification: new options and contents in BRENDA. Nucleic Acids Res. 2013, 41: D764-D772. 10.1093/nar/gks1049.

  40. 40.

    Primerano DA, Burns RO: Role of acetohydroxy acid isomeroreductase in biosynthesis of pantothenic acid in Salmonella typhimurium . J Bacteriol. 1983, 153: 259-269.

  41. 41.

    Elischewski F, Puhler A, Kalinowski J: Pantothenate production in Escherichia coli K12 by enhanced expression of the panE gene encoding ketopantoate reductase. J Biotechnol. 1999, 75: 135-146. 10.1016/S0168-1656(99)00153-4.

  42. 42.

    Wilson ACC, Ashton PD, Calevro F, Charles H, Colella S, Febvay G, Jander G, Kushlan PF, Macdonald SJ, Schwartz JF, Thomas GH, Douglas AE: Genomic insight into the amino acid relations of the pea aphid, Acyrthosiphon pisum, with its symbiotic bacterium Buchnera aphidicola . Insect Mol Biol. 2010, 19: 249-258. 10.1111/j.1365-2583.2009.00942.x.

  43. 43.

    Russell CW, Bouvaine S, Newell PD, Douglas AE: Shared metabolic pathways in a coevolved insect-bacterial symbiosis. Appl Environ Microbiol. 2013, 79: 6117-6123. 10.1128/AEM.01543-13.

  44. 44.

    Schneider F, Kramer R, Burkovski A: Identification and characterization of the main β-alanine uptake system in Escherichia coli . Appl Microbiol Biotechnol. 2004, 65: 576-582. 10.1007/s00253-004-1636-0.

  45. 45.

    Charles H, Balmand S, Lamelas A, Cottret L, Pérez-Brocal V, Burdin B, Latorre A, Febvay G, Colella S, Calevro F, Rahbé Y: A genomic reappraisal of symbiotic function in the aphid/Buchnera symbiosis: reduced transporter sets and variable membrane organisations. PloS One. 2011, 6: e29096-10.1371/journal.pone.0029096.

  46. 46.

    Merkamm M, Chassagnole C, Lindley ND, Guyonvarch A: Ketopantoate reductase activity is only encoded by ilvC in Corynebacterium glutamicum . J Biotechnol. 2003, 104: 253-260. 10.1016/S0168-1656(03)00145-7.

  47. 47.

    Kelkar YD, Ochman H: Genome reduction promotes increase in protein functional complexity in bacteria. Genetics. 2013, 193: 303-307. 10.1534/genetics.112.145656.

  48. 48.

    Dunbar HE, Wilson ACC, Ferguson NR, Moran NA: Aphid thermal tolerance is governed by a point mutation in bacterial symbionts. PLoS Biol. 2007, 5: 1006-1015. 10.1371/journal.pbio.0050096.

  49. 49.

    Price DRG, Duncan RP, Shigenobu S, Wilson ACC: Genome expansion and differential expression of amino acid transporters at the aphid/Buchnera symbiotic interface. Mol Biol Evol. 2011, 28: 3113-3126. 10.1093/molbev/msr140.

  50. 50.

    Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM: Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007, 35: W71-W74. 10.1093/nar/gkm306.

  51. 51.

    Cronan JE: β-alanine synthesis in Escherichia coli . J Bacteriol. 1980, 141: 1291-1297.

  52. 52.

    Bukhari AI, Taylor AL: Genetic analysis of diaminopimelic acid- and lysine-requiring mutants of Escherichia coli . J Bacteriol. 1971, 105: 844-854.

  53. 53.

    Tauch A, Kirchner O, Wehmeier L, Kalinowski J, Puhler A: Corynebacterium glutamicum DNA is subjected to methylation-restriction in Escherichia coli . FEMS Microbiol Lett. 1994, 123: 343-347. 10.1111/j.1574-6968.1994.tb07246.x.

  54. 54.

    Elbing K, Brent R: Media preparation and bacteriological tools. Curr Protoc Mol Biol 2002, Chapter 1:Unit 1.1. doi:10.1002/0471142727.mb0101s59.

Download references


We are grateful to Philip Kushlan for useful discussion. This work was supported by National Science Foundation Award 1121847 (to ACCW and DRGP).

Author information

Correspondence to Alex CC Wilson.

Additional information

Competing interests

Both authors declare that they have no competing interests.

Authors’ contributions

DRGP and ACCW conceived of and designed the study and wrote the paper. DRGP performed all experimental work. Both authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Price, D.R., Wilson, A.C. A substrate ambiguous enzyme facilitates genome reduction in an intracellular symbiont. BMC Biol 12, 110 (2014).

Download citation


  • Bacteriocyte
  • Symbiosis
  • Endosymbiosis
  • Co-evolution
  • Gammaproteobacteria