Previous Article | Next Article ![]()
Journal of Bacteriology, June 2005, p. 4229-4237, Vol. 187, No. 12
0021-9193/05/$08.00+0 doi:10.1128/JB.187.12.4229-4237.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona 85721
Received 21 January 2005/ Accepted 7 March 2005
|
|
|---|
|
|
|---|
-Proteobacteria provide strong evidence that B. aphidicola genomes are derived from larger genomes through loss of most genes (12, 19). Despite this massive gene loss, the three completed genome sequences for B. aphidicola of different aphid species show that each retains the enzyme-encoding genes for the biosynthesis of all or most essential amino acids. These 45 to 55 genes comprise approximately 10% of the genome. Retention of these genes fits with a major nutritional role of B. aphidicola, as hypothesized by early investigators (3) and as supported by experimental evidence (5-7, 26). Using relatively abundant nonessential amino acids and sugar derived from the aphid's diet of phloem sap, B. aphidicola synthesizes essential amino acids required by animals but rare in phloem (23). The mutualistic role of these pathways is highlighted by the observations that genes underlying synthesis of the nonessential amino acids have been eliminated from B. aphidicola genomes and that all amino acid biosynthetic pathways are typically missing from similarly small genomes of pathogenic bacteria, which obtain the needed nutrients from hosts. Since persistence of B. aphidicola depends on host reproduction, natural selection should favor the ability to regulate these genes appropriately in the context of the symbiosis.
Reduced bacterial genomes are particularly lacking in genes regulating transcription. For example, the E. coli genome encodes 233 experimentally verified regulatory proteins, including those targeting individual pathways and those acting as global regulatory proteins (27), but B. aphidicola contains only about five such genes, as inferred from homology. Most of the difference reflects gene loss in the B. aphidicola lineage. With regard to amino acid biosynthesis, B. aphidicola retains many of the structural genes encoding the enzymes but has lost almost all genes regulating their expression. For example, in E. coli and some other bacteria, transcription of argA, underlying the rate-limiting step in arginine biosynthesis, is repressed by argR, with repression enhanced in the presence of arginine (14). Two of the sequenced B. aphidicola genomes retain argA (and other genes for arginine biosynthesis) but have lost argR. Of the 16 E. coli genes known to encode regulators affecting transcription of amino acid biosynthetic genes (Table 1), only two are present in the sequenced genome of B. aphidicola of Schizaphis graminum [Buchnera(Sg)]. Only one of these, metR, has regulation of amino acid synthesis as a primary role; MetR regulates methionine production in E. coli (15, 32). The other, himA, encodes a protein that affects transcription of numerous E. coli genes, including those encoding carbamoyl phosphate synthetase, which underlies steps in both arginine and pyrimidine biosynthesis (9). In E. coli, most regulatory genes affecting amino acid biosynthesis act by repressing transcription in the presence of the amino acid end product. In the evolution of the Buchnera-aphid symbiosis, selection for overproduction of amino acids may have favored the inactivation or elimination of repressors. Indeed, metR, the only gene encoding a specific regulator of amino acid biosynthesis that is retained in Buchnera(Sg), is unusual in acting as a transcriptional activator. Evidence from E. coli and other bacteria shows that MetR enhances the transcription of metE especially when bound to the MetE substrate, homocysteine (32, 35).
|
View this table: [in a new window] |
TABLE 1. Regulatory genes for essential amino acid biosynthetic pathways in E. coli, with those retained in the genomes of Buchnera(Sg) and Buchnera(Ap) indicated
|
Although B. aphidicola has clearly lost most ancestral regulatory elements, inspection of the genome sequence cannot rule out the possibility that novel mechanisms for transcriptional regulation have evolved. For example, modifications of individual genes or binding sites may yield novel control mechanisms suited to the symbiotic lifestyle. That B. aphidicola may have novel mechanisms for controlling amino acid biosynthetic genes is strongly suggested by the relocation of genes for biosynthesis of leucine and tryptophan to plasmids (2, 11, 22).
In this paper, we present a first attempt to determine the capacity of B. aphidicola for controlling the expression of amino acid biosynthetic genes in response to changes in nutrition. Initially we employed a full-genome microarray based on the genome sequence of Buchnera(Sg) (20, 30, 36) to identify transcripts that respond to experimentally induced dietary changes. Quantitative reverse transcriptase PCR (qRT-PCR) was used to verify microarray findings and to examine consequences of interspecies differences in the presence of intact metR.
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Aphids, symbiotic bacteria, and corresponding food plants for organisms used in this study
|
Two additional sets of experiments were performed for use in qRT-PCR assays aimed at examining relative expression of metE. Treatments consisted of a control (tap water), tap water supplemented with glutamine (100 mM), and tap water supplemented with homocysteine (20 mM), with each solution placed in the reservoirs containing plant leaves, as described above. As in the treatments for the microarray experiments, plants were left in the treatment conditions for 24 h, and then aphids were added and allowed to feed for 24 to 30 h. In these treatments, aphids were mixed-age individuals transferred from a large colony. Each condition was replicated four times within each experimental set in separate chambers.
To compare the changes in metE transcript levels of Buchnera from A. pisum [Buchnera(Ap)] to those of Buchnera(Sg), we performed similar amino acid manipulations for A. pisum using cut fava bean seedlings in tap water and glutamine and homocysteine solutions (same concentrations as above). Aphids were mixed-age individuals from a large colony on fava bean. These manipulations were performed synchronously with the treatments for the qRT-PCR experiments for S. graminum. These samples were analyzed using qRT-PCR only.
When S. graminum was harvested from barley leaves, each sample was separated into two subsamples. One was used for microarray hybridizations, with RNA extraction and labeling protocols as described by Wilcox et al. (36). The other subsample was used to verify that the treatments did affect amino acid profiles in aphid body fluids. Free amino acids were extracted and individual concentrations determined using high-performance liquid chromatography, as described by Telang et al. (31). In addition, for some of the amino acid treatments and controls, we collected leaf exudates in order to assess treatment effects on amino acid profiles in phloem sap, as described by Telang et al. (31).
As for S. graminum, portions of the A. pisum experimental samples were also reserved for free amino acid determinations of aphid body fluids.
Microarrays and hybridization. The design and construction of the microarrays are described by Wilcox et al. (36) and only briefly summarized here. Individual genes were represented by PCR-amplified fragments obtained using primers designed using the Buchnera(Sg) genome sequence (30). The amplicons corresponded to 95% of coding genes, plus numerous spacers that were longer than 100 nucleotides, stamped onto glass slides, with each gene fragment present as eight copies per slide, as paired spots within four subgrids. All amino acid biosynthetic genes were represented. RNA extraction, labeling of cDNA with Cy3/Cy5, hybridization, and image scoring were carried out as described by Wilcox et al. (36). Total extracted RNA from the aphid sample was used, including host RNA and stable RNAs. Each slide was hybridized with transcriptomes resulting from two treatments, labeled with separate dyes. Each treatment was hybridized a total of two times, labeled once with Cy3 and once with Cy5 (arginine, glutamic acid, and leucine treatments), or three times, including labeling with both dyes (control, glutamine, tryptophan, and starvation treatments). Different hybridizations represented separate treatment sets; therefore, hybridizations are both biological and technical replicates. Slides were scanned using an Applied Precision ArrayWorx scanner, and spot intensity data were extracted using Digital Genome Version 2.19 (MolecularWare). Hybridization with extracted DNA was also carried out to determine differences in spot intensities arising from differential binding or fragment length of individual genes.
Microarray data analyses. Spot signal data were imported into SAS version 9 (25). Analyses followed the procedure described by Wilcox et al. (36) for assessing signals for individual genes. Briefly, after log2 transformation of the raw signal data for individual spots, we applied a mixed-model ANOVA to remove systematic effects attributable to individual slides, to dye, and to interactions. Average spot intensities from the DNA hybridization were used to remove gene-specific effects such as those reflecting individual fragment length, DNA concentration in the probe spots, or sequence features. The signal was normalized for each slide, and we obtained the residuals for each spot. We used the microarrays to explore the overall patterns and magnitudes of response, and we did not select a particular P value as a threshold for significance. However, we considered those with P < 0.0001 to be candidates for further examination. Our focus in these experiments was on the 55 genes underlying amino acid biosynthesis.
qRT-PCR.
To verify expression changes in metE, the gene appearing to be most responsive to any treatments based on the microarray hybridizations, we performed qRT-PCR. These tests were performed for glutamine and homocysteine treatments relative to the control and were conducted for both Buchnera(Sg) and Buchnera(Ap). Primers were designed to amplify short regions of these genes, using the determined genome sequences for these two B. aphidicola isolates (28, 30). Oligonucleotide primers (5'
3') for metE were AAAAGCGATGCTAACTGGACCA and CCAACCCGATTTGTTTAGCAA for Buchnera(Sg) and TAAAGGGATGTTAACAGGACCG and CTAAAGCAATTTGTGTTGCAA for Buchnera(Ap), with both reactions yielding 95-bp products. In addition to the target gene, we included two control genes for each comparison. The first, rpsL, encodes a ribosomal protein; the second, atpE, encodes a subunit of ATP synthase. Both control genes showed near-constant relative transcript abundance under different treatments based on the microarray results, and neither is directly related to amino acid synthesis. Primers (5'
3') for rpsL wereGGCCACGGTTAATCAATTGGT and AACGCCTCTTTTTTGCGGAC[Buchnera(Sg)] and CTACACCTAAAAAGCCTAATTCAGC and CCTGCACAATCTAAGGAGCC [Buchnera(Ap)], giving products of 97 and 199, bp respectively. Primers (5'
3') for rpsL were TAGAAGGAGCTGCAAGGC and AGCCTAATCCAACTGCAATC [Buchnera(Sg)] and AAGTTTTTAGAGGGAGCCGC and CAGCAATCATTGGAATCGC [Buchnera(Ap)], yielding products of 108 and 103 bp, respectively. qRT-PCR was carried out with a Roche Lightcycler and Roche reagents. We used the same "touchdown" reaction protocol for all reactions, as follows: 95°C for 10 min and then 40 cycles of 95°C for 5 s, variable annealing temperature for 15 s, and 72°C for 5 s. The annealing temperature was 68°C for the first three cycles and was then lowered by 1°C each cycle until it reached 55°C, with remaining cycles at 55°C. An internal standard curve was generated for each of the reactions, using serial dilutions of the cloned template. Each experimental sample consisted of several pooled aphids from a treatment replicate. The same RNA extraction was used for control and test genes, and template copy number was estimated from the standard curve for each gene. Therefore, variation in quantity of sample or in efficiency of an extraction would affect control and test genes equally.
For metE, the ratio of its copy number to the copy number for each of the two control genes was calculated to yield a measure of relative transcript level for a particular treatment. Within each experimental set, all treatments were run synchronously, and analyses were performed separately for each set. Data were categorized according to treatment (control, glutamine, homocysteine). The ratios of transcript copy numbers (metE/control gene) were normalized by a factor that set the lowest value at 1. For parametric tests (one-way ANOVA and t test), results were log transformed to meet criteria of normality and equality of variances; nonparametric tests were also performed (Mann-Whitney for set 1, Kruskal-Wallis for sets 2 and 3). Statistical tests were carried out using the InStat (GraphPad) software package.
|
|
|---|
![]() View larger version (25K): [in a new window] |
FIG. 1. Effects of homocysteine (Hcn) and glutamine (Gln) treatments on relative concentrations of amino acids in phloem sap of barley (A), body fluids of S. granimum (B), and body fluids of A. pisum (C). Data are presented as proportions of the total amino acid pool.
|
The results in Fig. 1 depict relative abundances. We presume that adding an amino acid did not result in absolute decreases in concentrations of any amino acids (i.e., decreased proportions of some amino acids in treatment samples must reflect elevation of other amino acids). Overall, the evidence consistently showed that adding amino acids to the phloem resulted in large effects on amino acid content of the aphid diet and on free amino acid levels in aphid body fluids. Similar changes in body fluid composition occurred in both S. graminum and A. pisum (Fig. 1B and C).
Microarray results. Overviews of gene expression changes in response to dietary manipulations are shown in Fig. 2, which contrasts the control treatments with the various experimental treatments, in the form of "volcano" plots of fold change and statistical significance, and Fig. 3, which shows the average fold changes for genes underlying different enzymatic steps in amino acid biosynthesis. Only one gene, metE, shows a strong response to any of the treatments; transcripts of metE changed >fourfold in response to two nonessential amino acids (glutamic acid and glutamine), >fourfold in response to the essential amino acid arginine, and about twofold in response to tryptophan (Fig. 2A to D and 3). In contrast, metE showed no response to leucine or to starvation (Fig. 2E and F). After metE, the most dramatic change in an amino acid biosynthetic gene involved argA, which showed a slight but statistically robust increase in relative transcript abundance in response to glutamine and arginine and a decrease in response to starvation but no response to leucine or glutamic acid (Fig. 2 and 3). In all cases, changes in argA transcript levels were less than twofold. The only other amino acid biosynthetic gene that showed a possibly significant response was trpA, which showed a minor decrease in tryptophan-treated samples (Fig. 2D). The starvation treatment, compared to controls, showed a decrease in transcripts of argA and some increase in the plasmid-encoded genes (trpEG and leuABCD). Other than the change in metE transcript abundance, all of the responses were small in magnitude (<twofold), and most did not approach statistical significance.
![]() View larger version (38K): [in a new window] |
FIG. 2. Volcano plots showing the statistical significance (y axis) and the fold change (x axis) for differences between controls and treatments involving the experimental addition of single amino acids to the host plants of S. graminum. Statistical significance is based on a mixed-model ANOVA and reflects both biological and technical replication. Significance levels less than 0.0001 were rounded to that value.
|
![]() View larger version (80K): [in a new window] |
FIG. 3. Average effects of treatments on transcript levels of amino acid biosynthetic genes of Buchnera(Sg), based on microarray results.
|
![]() View larger version (20K): [in a new window] |
FIG. 4. Transcript abundance of metE following control, glutamine, and homocysteine treatments. Values are based on estimates from qRT-PCR. The three sets for Buchnera(Sg) were performed at separate times, with the third set of bars based on the glutamine and control treatments performed for the microarrays. Each value represents four replicates, with replicates derived from aphids on different plants in different growth chambers.
|
Buchnera(Ap), which diverged from Buchnera(Sg) about 60 million years ago (18), lacks metR and retains metE. Thus, if MetR is the basis for the major changes in metE transcript abundance in Buchnera(Sg), Buchnera(Ap) is predicted to lack the metE response. Treatments with glutamine and homocysteine had similar effects on amino acid profiles of aphid body fluids in both species (Fig. 1). But, as predicted if MetR is involved in metE repression, these treatments failed to induce a metE response in Buchnera(Ap) (Fig. 4, P > 0.8 for both Kruskal-Wallis and ANOVA on log-transformed data). Thus, MetR appears to be the major regulator underlying the large responses of metE to changes in amino acid pools in Buchnera(Sg), and loss of metR in the lineage leading to Buchnera(Ap) corresponds to a lack of regulation of metE expression.
|
|
|---|
These results for transcript abundances under different treatments would reflect any regulatory responses affecting either transcription rates or transcript longevity. Based on observations of other bacteria, transcriptional control is expected to be central in regulating these pathways, but our study does not exclude the possibility of regulatory changes involving the translation, persistence, or activity of enzymes.
Our results indicate that strong shifts in transcriptional patterns occur solely when the ancestral regulator is retained: metE is the only gene showing a strong response in response to the numerous dietary manipulations that we applied. Furthermore, this response occurs only when metR, the known regulator of metE, is present in the genome. Buchnera(Ap), lacking metR, shows no such response. Interestingly, metR is one of the few regulatory genes that acts as an activator of transcription rather than as a repressor. Therefore, the immediate effect of its loss would be to lower rather than increase transcription of its target gene.
Another complicating factor in understanding the nutrient exchange underlying this symbiosis is that not all amino acids enter the bacteriome and the symbiont cells. Compartmentalization of the symbionts, with barriers to the free exchange of materials with the host, means that excess amounts of a particular essential amino acid in hemolymph or in host cells may not be detectable or useable by the symbionts. Previous nutritional studies using labeled N sources provided evidence that many amino acids are not taken up by B. aphidicola (6, 7, 13). Possibly, only nonessential amino acids are imported; two of these, glutamic acid and aspartic acid, serve as substrates for the pathways producing the essential amino acids. Arginine is readily broken down into glutamic acid plus glutamine in animal tissues, probably explaining why arginine and glutamine treatments had similar effects on transcript profiles (Figs. 2 and 3). In our treatments, tryptophan was also metabolized, at least partly in the plant, to give an increase in glutamine in hemolymph (data not shown), consistent with the elevated expression of metE and mildly elevated argA, as observed for glutamine and glutamic acid treatments. In contrast, the addition of leucine, which appears not to be catabolized in plants or aphids, did not result in elevated glutamine (Fig. 1) and had no effect on metE transcript levels (Fig. 2 and 3). It is possible that manipulation of glutamine, glutamic acid, arginine, and tryptophan all increased homocysteine levels within B. aphidicola cells. In plant and insect samples, the levels of homocysteine were always low relative to other amino acids; we did not assay absolute concentrations.
This study, as well as a previous study on transcriptome responses to heat stress (36), point to a lack of novel mechanisms for regulating transcript levels in B. aphidicola. Buchnera(Sg) does mount a heat shock response, but compared to free-living bacteria, it is relatively feeble in terms of numbers of loci involved and magnitude of changes (36). For variability in amino acid pools and in temperature, the only strong responses are those expected on the basis of the few ancestral regulatory mechanisms that are retained. Furthermore, these appear to function in a fashion similar to that in ancestors. For example, the consensus sequence of the global heat shock promoter is exactly the same in Buchnera(Sg) and in E. coli, and only the loci retaining this promoter show a response to heat stress (36). Likewise, homocysteine appears to be the main molecule affecting MetR activation of metE promoters, as in E. coli and other enterics. Most "heat shock" genes and most genes for amino acid biosynthesis show little or no transcriptional responses to these variables, in contrast to the case for free-living, related bacteria.
The difference between Buchnera(Sg) and Buchnera(Ap) in the presence and action of metR indicates that loss of regulatory genes is ongoing during Buchnera evolution. It is possible that differences in content of sulfur-containing organic compounds between grasses (the host plants of S. graminum) and herbaceous legumes (used by A. pisum) affect the need for regulating methionine production and thus underlie the differential retention of metR by these two species. Based on the genome sequences, the two species also differ in ability to utilize sulfate; the genes underlying reduction to sulfide are intact in Buchnera(Ap) but are inactivated as pseudogenes in Buchnera(Sg), possibly due to the presence of more sulfur-containing organic compounds in grasses (30).
Distantly related bacteria often display completely different regulatory mechanisms for the expression of amino acid biosynthetic pathways, indicating that new regulatory systems often can replace old ones while the structural genes encoding the biosynthetic enzymes are continuously maintained (e.g., see references 37 and 38). In light of this lability in regulatory mechanisms and of the great age of the Buchnera symbiosis (>100 million years [18]), the lack of novel transcriptional responses distinctive to the symbiotic lifestyle is somewhat surprising. Among the possible reasons that B. aphidicola lacks novel regulatory systems is the remarkable stasis of its genome contents and arrangement (30), which may have limited its ability to acquire novel regulatory elements. Another explanation may lie in its small genetic population size, which may limit the effectiveness of selection for more fine-tuned transcriptional regulation (8, 16, 21).
Funding was from an NSF postdoctoral fellowship to J.L.W. and from grant NSF-0313737 to N.A.M.
|
|
|---|
-Proteobacteria. PLOS Biol. 1:101-108.[CrossRef]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»