Conservation of the links between gene transcription and chromosomal organization in the highly reduced genome of Buchnera aphidicola
© Viñuelas et al; licensee BioMed Central Ltd. 2007
Received: 07 February 2007
Accepted: 04 June 2007
Published: 04 June 2007
Genomic studies on bacteria have clearly shown the existence of chromosomal organization as regards, for example, to gene localization, order and orientation. Moreover, transcriptomic analyses have demonstrated that, in free-living bacteria, gene transcription levels and chromosomal organization are mutually influenced. We have explored the possible conservation of relationships between mRNA abundances and chromosomal organization in the highly reduced genome of Buchnera aphidicola, the primary endosymbiont of the aphids, and a close relative to Escherichia coli.
Using an oligonucleotide-based microarray, we normalized the transcriptomic data by genomic DNA signals in order to have access to inter-gene comparison data. Our analysis showed that mRNA abundances, gene organization (operon) and gene essentiality are correlated in Buchnera (i.e., the most expressed genes are essential genes organized in operons) whereas no link between mRNA abundances and gene strand bias was found. The effect of Buchnera genome evolution on gene expression levels has also been analysed in order to assess the constraints imposed by the obligate symbiosis with aphids, underlining the importance of some gene sets for the survival of the two partners. Finally, our results show the existence of spatial periodic transcriptional patterns in the genome of Buchnera.
Despite an important reduction in its genome size and an apparent decay of its capacity for regulating transcription, this work reveals a significant correlation between mRNA abundances and chromosomal organization of the aphid-symbiont Buchnera.
Past genomic studies have comprehensively described the organization of the bacterial chromosome, for example in terms of gene localization, order and orientation. The degree of organization has been shown to increase with genome size, overall GC composition and the presence of nucleoid-binding proteins . This organization of the chromosome can be described as being an adaptive and functional tool, essential for the survival of the bacterial cell. More precisely, several studies have identified strand asymmetries in the distribution of genes between the leading and the lagging strand of DNA (for a review see Rocha ). Indeed, bacterial genomes carry, on average, from 78% (for genomes containing the polymerase PolC) to 58% (for the other genomes) of their genes on the leading strand . This bias is even more important when the essentiality of genes is taken into account, and essential genes distribution bias reaches 76% and 94% in Escherichia coli and Bacillus subtilis respectively . The asymmetry of the distribution of the genes between the two strands of DNA is explained as being a means of minimizing interruptions in gene transcription due to collisions between the DNA and RNA polymerases. The most generally accepted hypothesis is that co-directional collisions on the leading strand have a weaker effect on RNA polymerase processivity than the head-on collisions occurring on the lagging strand. This "replicational selection" should guarantee that genes on the leading strand, particularly the essential ones, are efficiently transcribed [5–8].
Gene organization in the bacterial chromosome has very well known effects on transcription and a standard example is the organization of genes into operons, which allows for a sophisticated regulation of gene expression . It has also been shown that neighbouring genes in a bacterial chromosome tend to be co-expressed, even if they are not in the same operon , suggesting the existence of a "supra-operonic" organization . This is a perfect illustration of the mutual influence between gene expression levels and chromosomal organization in bacteria. Another example is the significant effect of DNA supercoiling on transcription. Recent studies have suggested that the dynamical structure of the nucleoid acts as a "transcription factor" in E. coli , and is probably responsible for the presence of spatial transcriptional patterns in free-living bacteria [13–15].
While these observations were made on free-living bacteria, very few data exist on symbiotic bacteria characterized by reduced genomes, shaped by their adaptation to the host metabolic requirements and by their high evolution rate. Buchnera aphidicola, the endosymbiont of the aphids, is one of the best studied examples of an intracellular bacterium with a reduced genome [16, 17]. The genomes of four Buchnera from different aphid species have been sequenced so far, with sizes ranging from 0.42 to 0.65 Mb [18–21]. In all the Buchnera genomes, most of the genes involved in the biosynthesis of essential amino acids (EAAs) that the insect cannot synthesize, or find in sufficient amounts in the phloem sap, were retained, whereas almost all genes regulating their expression were lost . In Buchnera from the pea aphid Acyrthosiphon pisum (BAp), among the 608 chromosomal genes, 56% are located on the leading strand DNA. This observation reveals a small gene-strand bias, although equivalent to that of its closely related free-living bacteria:E. coli (55%) .
Gene essentiality has been defined by experimental techniques in cultivable bacteria, such as E. coli for which a large repertoire of null mutants is available. For a small number of genes, even in the absence of mutants, essentiality is defined on the basis of their function. For example, genes encoding ribosomal proteins are usually described as essential, whereas genes involved in the flagellar apparatus are classified as non-essential in bacteria . Due to the specific lifestyle and symbiotic functions of Buchnera, the definition of essentials genes for this bacterium is particularly difficult. For this reason, we considered the minimal gene set for supporting bacterial life proposed by Gil et al.  as the most appropriate for our study. Indeed, the determination of this set is based on comparative analyses of five bacterial insect endosymbionts with small genomes (among them three Buchnera genomes), one of the smallest bacterial genome sequenced (Mycoplasma genitalium), and the essential gene lists available for B. subtilis and E. coli in the literature. Based on this minimal set of essential genes, BAp harbour 190 essential encoding-protein genes and 60% of them are located on the leading strand. This indicates a strand bias for essentiality, albeit lower than that observed in free-living bacteria. Recent transcriptomic analyses have shown that, despite a low response to different physiological conditions imposed on the aphid host, Buchnera retained the capability to express genes at different levels in basal conditions. This expression seems to be correlated to genome organization in the putative transcription units (pTU), even in the absence of many of the known transcriptional regulatory proteins .
This work focuses on studying the relationships between gene transcription and chromosome organization in the highly reduced genome of Buchnera, by coupling genomic information with transcriptomic data. For that, mRNA abundances were measured with a full-genome oligo-array dedicated to BAp. The basal transcriptome of Buchnera was then analysed regarding (i) the location of the genes on the leading/lagging strand, the gene essentiality and operon organization, (ii) the gene evolution rate (as measured by the GC rate), and (iii) the spatial location of genes along the chromosome. One crucial point in this study was the choice of the data normalization procedure. Since 2002, authors have described the possible advantages of genomic DNA (gDNA) standard for microarray data normalization in the investigation of microbial genomes [27–29]. This useful attractive standard is readily available, inexpensive, and invariant over time and from laboratory to laboratory. Moreover, the main advantage of the normalization by gDNA is that it takes into account the probe/target affinity in order to allow for inter-gene calibration and expression level comparisons. For these reasons, we have developed and validated in this study a "spot by spot" microarray data normalization method based on gDNA.
Effectiveness of using genomic DNA for microarray data normalization
Finally, we took advantage of the design of several probes (2 or 3) per gene in our microarray to validate the gDNA normalization procedure. Hence, we compared the within-gene fluorescence signal standard deviations of all the 3-probe genes, before and after normalization. The results presented in Figure 1 show, after normalization, an overall reduction of the fluorescent signals variability for about 79% of the genes. Close results were obtained for the 2-probe genes (65%, data not shown). Although low mRNA abundances seemed to be the major source of failure for gDNA normalization (data not shown) no intensity correction was included because the estimation of variability was not enough powerful with only 2 or 3 different probes per gene.
Analysis of the links between Buchnera normalized mRNA abundances and chromosome organization
DNA strand position, putative operons and gene essentiality
In this part of the work, we analysed the effect of gene DNA strand position, operon organization and gene essentiality on the variability of the basal mRNA abundances in Buchnera.
Analysis of the pTU factor alone revealed that normalized basal gene transcription (nBGT) variation was always smaller within pTU than between them, indicating that genes belonging to the same pTU tend to be transcribed at similar rates (ANOVA F-test, R2 = 0.640, P-value < 10-4). Analysis of essentiality alone revealed that the Buchnera genome contains 32% of essential genes versus only 6% for the E. coli genome .
Effects and interactions of the global ANOVA model fitted on the Buchnera log2 nBGT levels
Sum of Squares‡
< 0.0001 *
< 0.0001 *
Essentiality × Operon
Essentiality × Strand
Strand × Operon
Essentiality × Strand × Operon
Gene evolution rate
In the second part of this work we have analysed the relationship between nBGT levels and gene evolution rates in Buchnera. Gene evolution rates were estimated here by the GC content of genes instead of by non-synonymous substitution rates (Ka), since both parameters are highly correlated in Buchnera. Figure 3 shows that mRNA abundances and GC content are correlated (R2 = 0.151, P-value < 10-4): highly expressed genes are the most GC rich (the most highly conserved in Buchnera). More precisely, areas 1 and 2 in Figure 3 correspond to genes that are slowly evolving and that are weakly or highly expressed respectively, whereas areas 3 and 4 correspond to genes that are rapidly evolving and that are weakly or highly expressed respectively. In the second area we found almost all genes encoding ribosomal proteins (among them: rplL, rplP, rpmB, rpsJ, rpsK) but also genes encoding chaperone proteins (dnaK, dnaJ and mopA,mopB). An important proportion of the flagellar genes seem to evolve rapidly: 16 out of 26 in Buchnera (among them 9/12 fli genes) and almost all of them are highly expressed (11 out of 16). Moreover, 11 out of these 16 highly evolving flagellar genes are located on the leading strand (among them 8 fli genes). Interestingly, the only two Buchnera orphan genes (yba3 and yba4), are highly expressed and show a low GC content (area 4). The genes encoding the enzymes involved in the biosynthesis of EAAs generally evolve slowly and show high mRNA abundances. Among them, we have observed that the genes belonging to the two operons involved in the biosynthesis of isoleucine and valine are particularly highly expressed and well conserved in Buchnera (area 2). Moreover, ileS, valS, and six other genes encoding aminoacyl-tRNA synthetase for EAAs are weakly expressed, whereas genes encoding aminoacyl-tRNA synthetase for non-EAAs are moderately or highly expressed (Median test, P-value = 1.4 × 10-2). Finally, area 3 includes pseudogenes, such as apbE, cvpA, and yigL.
Spatial location of genes along the chromosome
Some recent papers [32–34] reported that spurious periodicities in transcriptomic data can be attributed to the spatial arrangements of the probes on an oligo-array. Even if by construction our gDNA normalization procedure is appropriate to take into account this bias, we checked the possible presence of periodic components in gDNA signals. The absence of significant periodicities for gDNA (black curve of Figure 5A) supports the idea that the observed spatial patterns for mRNA signals are proper to the Buchnera transcriptome.
Autocorrelation values of Buchnera nBGT levels for the different permutations of gene positions
Results of the spectral analysis testing the periodicity of Buchnera log2 nBGT
To our knowledge, we have presented here the first analysis showing i) a representation of whole-genome transcriptional data in Buchnera and ii) the links between mRNA abundances and chromosomal organization in a highly reduced bacterial genome, on the basis of experimental data. The large number of prokaryotic genome sequences available in databases has made it possible to study these links in many bacterial chromosomes. Using codon adaptation index (CAI) values computed from ribosomal proteins, Rocha and Danchin have shown that, in B. subtilis and in E. coli, the importance of the expressiveness in determining the localization of the genes on the leading strand is negligible, or even absent, when essentiality is taken into account . They have confirmed these results for other sequenced genomes, with the exception of some non-free-living bacteria among which Buchnera was one of the most important . However, the authors underlined the difficulty of obtaining correct CAI values in the genomes of intracellular bacteria, which generally do not show sufficient codon usage bias. Moreover, another explanation of the exceptions they found for Buchnera can be attributed to authors' assignment of gene essentiality, which was exclusively based on the homology with E. coli. Indeed, while Rocha et al.  found no essential genes strand distribution bias in Buchnera, we found a bias of 60% using the minimal gene set proposed by Gil et al. . More recently, Price et al.  have revisited the question of gene-strand bias in bacteria, using gene expression microarray data, and they have shown that, in B. subtilis and E. coli, the genes in operons on the leading strand DNA are more highly expressed than genes in operons on the lagging strand. This observation was true for both essential and non-essential genes.
In Buchnera, we found that highly expressed genes are generally the essential genes within pTU. Independently of the essentiality factor, the genes within pTU are more highly expressed than singletons. We also found that for genes within pTU the essential genes are more highly expressed than non essential ones whereas this effect was not observed for singletons. These results underline the conservation of a coherent relationship between mRNA abundances and gene essential functions in the reduced genome of Buchnera.
In this work we also analysed the relationship between mRNA abundances and genes GC ratio in Buchnera (taken as an indicator of gene evolution rate). By combining transcription analysis, evolution rates and comparative genomics, we were able to define new candidates for the essential gene set of Buchnera. In bacteria, genes encoding ribosomal, cell division and chaperone/protease proteins are considered as essential and they are also known to be constitutively highly expressed. As we expected, our data showed that almost all genes encoding ribosomal proteins, and genes encoding chaperonins, are relatively well conserved and also highly expressed in Buchnera. Interestingly, among the most highly conserved and expressed genes, we found genes involved in the biosynthesis of EAAs and, in particular, all the genes of the isoleucine and valine pathway. Also remarkable is that the genes encoding the aminoacyl-tRNA synthetases for 8 out of the 10 aphid EAAs are weakly expressed, whereas genes encoding aminoacyl-tRNA synthetases for non-EAAs are either moderately or highly expressed. A possible explanation for this observation is that Buchnera, which is a relatively slow-growing bacterium, does not necessitate high rates of protein production and constitutively synthesizes EAAs in order to furnish them to the aphid. Reducing the abundance of specific aminoacyl-tRNA synthetases might increase the concentration of free EAAs in Buchnera cells facilitating the transport of these amino acids to the aphid host cell. This observation is reminiscent of a similar putative adaptive response of Buchnera, which selectively underexpresses pheT under aromatic EEAs shortage .
In our study, we also found that the orphan genes yba3 and yba4 seem to evolve rapidly and are highly expressed, yba4 being one of the ten highest expressed genes in Buchnera. The conservation of these genes in BAp and in the Buchnera harboured by the aphid Schizaphis graminum (BSg), coupled with their high expression level, suggest that they could be of particular relevance in the symbiosis of Buchnera with the aphids S. graminum and A. pisum, as it has also been proposed by Shimomura et al.  for A. pisum.
Among the genes rapidly evolving (this work, Tamas et al.  and Reymond et al. ) and highly expressed in Buchnera, we found most of the flagellar genes. These data, taken together with recent experimental evidence that the Buchnera incomplete flagellar apparatus can function as a "protein transporter" , support the idea that flagellar genes are taking on new important functions in the symbiotic context. The flagellar gene set of BAp is composed of 26 genes of which fliEFGHIJKMNPQR and flhAB genes are located on the leading strand (except fliE), and flgNABCDEFGHIJK genes are located on the lagging strand DNA (except flgN and flgA). These flagellar genes are also conserved in BSg. However, Buchnera from a distant aphid lineage, Baizongia pistaciae (BBp), has lost five flagellar genes (flgA, flgD, flgE, flgK, flgN) but not any of the fli ones. Finally, Buchnera with the most dramatically reduced genome, from the aphid Cinara cedri (BCc), has lost the same genes as BBp but also four other flg genes (flgB,flgC, flgG, flgJ) and four fli genes (fliE, fliJ, fliK, fliM), hence preserving only a minimal type III virulence secretion system . The evolutionary selection of the majority of the fli genes also suggests the possible importance of the new putative transport function of these genes in Buchnera. It is to note, that the most conserved fli genes among the Buchnera lineages, probably involved in the new function, are located on the leading strand, whereas most of the gene losses occurred on the lagging strand. Previous comparative genomics analyses [40, 41] had tempted to dissect the evolutionary forces driving the genome organization in several Buchnera lineages (i.e., gene strand bias, gene and protein composition, gene expressivity, gene evolution rate and gene loss). Our results are partly consistent with these previous analyses as we found (1) highly conserved genes are highly expressed, (2) essential genes in pTU are highly expressed and probably preserved from mutations by purifying selection, and (3) positive selection may shape new "symbiotic" functions for some genes highly expressed and highly evolving. However, we reject the former idea, based on CAI analyses, that expressiveness is a factor driving gene strand bias in Buchnera.
Finally, an important result of this study was the discovery of spatial patterns of transcriptional activity in the chromosome of Buchnera: i.e., the transcription of the genes along the chromosome is determined according to spatial constraints. From autocorrelation and spectral analysis, four groups of spatial patterns can be defined: (i) autocorrelated short-range (between 2 and 8 genes), (ii) periodic short-range (up to 17 genes), (iii) periodic medium-range (between 23 and 61 genes) and (iv) periodic long-range (over 87 genes) structural components.
Autocorrelated short-range patterns, determined by the autocorrelation function of gene transcription levels, showed that genes spaced by a gene-to-gene distance of less than 8 have highly correlated expressions. As has been suggested for E. coli and B. subtilis, we propose for Buchnera that these correlations reflect the co-ordinated transcription of genes within operons [14, 15]. This observation reinforces the result mentioned above concerning the conservation of functional transcription units in Buchnera. However, by permuting gene position along the chromosome, we observed that the organization of genes into putative operons is not sufficient to fully explain the observed periodic spatial patterns of transcription. Indeed, if the presence of spatial periodic components in Buchnera gene expression was only due to the conservation of operon structures, the modification of the order and/or the number of the singletons located between the pTU should not affect the periodicities and the autocorrelation values. We have shown, however, that these modifications reduced the spatial patterns, alleviating the importance of a high-order arrangement of all the genes along the chromosome of the endosymbiont.
Jeong et al.  have classified the transcriptional periods that they have found in E. coli into three categories: short-range (up to 16 genes); medium-range (100–125 genes); and long-range (600–800 genes). The existence of short periods, up to 17 genes in Buchnera, allows us to corroborate the hypothesis, proposed in a previous study on B. subtilis and on E. coli, that this short range element could be a property of the structural nucleoid common to other bacteria, corresponding to large DNA spirals on the nucleoid surface . The medium and the long-range periods are shorter in Buchnera than those identified for free-living bacteria, which is probably due to the greatly reduced size of its genome. However, these two kinds of periods are not yet understood and do not correspond to the domains identified so far in the nucleoid . The effect of the second and third simulated gene permutations (Perm. 2 and Perm. 3 on Figure 6) on the medium and long-range periods respectively, could be explained by the importance of the spatial location of the operons along the chromosome and by the neighbourhood of singletons that form "supra-operonic" structures in Buchnera. Moreover, the decrease of the maximum size of the autocorrelated groups of genes (Table 2), for the different permutations of gene positions, corroborates the hypothesis of "supra-operonic" structures in Buchnera. But this speculation needs to be more studied and experimentally confirmed. Finally, the observation of transcriptional periodic patterns, coupled with the conservation in its genome of some Nucleoid Associated Proteins (NAP) such as H-NS, IHF and Fis, suggest that Buchnera has maintained a nucleoid structure responsible for the differences in gene transcription levels in basal conditions. These three NAPs were previously found to be differentially expressed in Buchnera  facing nutritional constraints. However, their role in transcriptional regulation remains presently speculative in the aphid endosymbiont.
In conclusion, the analysis of mRNA abundances with regard to chromosomal organization in Buchnera has shown, despite an important reduction in its genome size and an apparent decay of its transcriptional regulatory capacity, a conservation of the relationship between these two parameters in the aphid symbiont. Our work shows that the organization of the genes into operons and their essential functions influence transcription in Buchnera, whereas no significant leading/lagging strand bias was observed. This work also underlines the difficulties in defining gene essentiality in intracellular symbiotic bacteria and the importance of an additional specific gene set for BAp. Finally, we showed the existence of a transcriptional periodicity along the chromosome of Buchnera. From these results, comparative analyses performed on Buchnera taken from other species of aphids, but also on other endosymbiotic bacteria, such as Blochmannia spp. from the carpenter ant or Wigglesworthia spp. from the tsetse fly would establish whether the results presented here are a common trait of the insect endosymbionts.
A long-established parthenogenetic clone (LL01) of A. pisum was maintained at 21°C, with a 16 hour light photoperiod, on Vicia fabae. It was shown to be free of any of the five taxa of secondary endosymbionts identified to date. For these experiments, we used aphids which were reared from birth to day-7 on the plant (basal conditions of gene expression in Buchnera) and Buchnera were purified from their host as described by Charles and Ishikawa .
Nucleic acids purification and labelling
The protocols for RNA extraction and labelling have been previously described by Calevro et al. . Briefly, total RNA was purified using the Trizol method and any possible gDNA contaminants were removed using DNase RQ1 RNase-free (Promega, Madison, WI, USA). Total RNA was subsequently purified on an RNeasy® column (Qiagen, Hilden, Germany). Then, 15 μg of total RNA were indirectly labelled by incorporating aminoallyl-dUTP into reverse transcript cDNA (Amersham Biosciences, Piscataway, NJ, USA). Finally, Cy3 fluorescent dyes (Amersham) were coupled to the targets prior to purification on an Autoseq™ G-50 column (Amersham).
Genomic DNA was treated with RNase A (Sigma-Aldrich, St Quentin Fallavier, France) and classically purified with phenol/chloroform/isoamyl alcohol (25:24:1, v/v/v) (Q-Biogene, Irvine, CA, USA). After testing different concentrations of gDNA (from 1 to 6 μg), a concentration of 5 μg was chosen for the normalization of mRNA signals (data not shown). gDNA labelling was performed using the Nick Translation Kit (Amersham) and by directly incorporating dUTP-Cy3 (Amersham). Non-incorporated fluorescent dyes were eliminated by purification on an Autoseq™ G-50 column.
Slide preparation and hybridization
The Buchnera oligo-array is composed of 6144 spots, including two probes per Buchnera gene that were designed by ROSO software . A third additional probe was designed for specific genes chosen for their function in the symbiotic relationships (amino acid biosynthetic genes and putative transporters). The 35-mer oligonucleotide probes were printed, as quadruplets, by two different pins onto Quantifoil Micro Tools aldehyde slides (Interchim, Montluçon, France). Positive and negative controls were spread all over the slide.
The slides were hybridized using an automated Ventana Discovery station (Ventana Medical Systems, Illkirch, France) at 45°C for 8 hours. Following hybridization, the slides were washed with solutions of variable stringency, dried by centrifugation and scanned for fluorescence using a GeneTac LSIV scanner (Genomic Solutions, Huntingdon, UK). The signal intensity value for each spot (pixel medians) was recorded and qualitatively analysed with GenePix Pro 4.1 software (Axon Instruments, Foster City, CA, USA).
Four slides were hybridized: 2 with Buchnera gDNA and 2 with cDNA. Microarray data are available in the Array Express database, accession E-TABM-193.
Genomic DNA normalization
For the validation of the normalization procedure, we used two specific data sets of 43 and 371 genes with 3 and 2 different probes per gene respectively. These probes sets were characterized by high quality signals in the two types of hybridization (2 slides with gDNA and 2 slides with cDNA)
Buchnera genomic data
Sequence data for the Buchnera aphidicola genome, and the corresponding annotations, were retrieved from GenBank . The definition of essential genes was based on the minimal gene set for supporting bacterial life (see introduction part) described in Gil et al. . Putative operons in Buchnera have been manually identified by comparison with the E. coli genome. For that, each E. coli orthologous gene in Buchnera was searched for in EcoCyc database , and groups of Buchnera genes corresponding to transcription units in E. coli were collected. After removing the singletons, we determined 82 pTUs in Buchnera consisting of 2 to 13 genes. In this study, we estimated gene evolution rates from the GC content of genes instead of from the non-synonymous substitution rates (Ka) as the two parameters are highly correlated in the AT-rich genome of Buchnera, and because the GC content allows for estimations of the two orphan genes of Buchnera.
Spectral analysis and estimation of the periodicities
Note that the distance between two genes used in this article is the difference of their ranks on the chromosome (approximately equivalent to the number of kb).
In order to analyse the dependence of transcriptional activity as a function of chromosomal position, we also used the spectral density of the autocorrelation function, estimated by Fourier transform techniques. The periodogram obtained by this method illustrates the periodic components, and their corresponding intensities, of the transcriptomic signals along the chromosome of Buchnera. The periodic signals of the main periods in Buchnera are drawn by averaging and smoothing the signals of all the segments making up each period. Smoothing of the signals was performed by the "Fit Spline" function of the JMP software.
Statistical analyses of microarray data, mean comparisons (ANOVA F-test) and distribution comparisons (Likelihood ratio test) were performed using JMP software. For distribution comparisons, we previously split Buchnera gene expression levels into three bins: the 100 most highly expressed genes (labelled as "highly expressed"); the 100 least expressed genes (labelled as "weakly expressed"); and the others (labelled as "moderately expressed").
- BAp :
Buchnera from the aphid Acyrthosiphon pisum
- BBp :
Buchnera from the aphid Baizongia pistaciae
- BCc :
Buchnera from the aphid Cinara cedri
- BSg :
Buchnera from the aphid Schizaphis graminum
- EAAs :
essential amino acids
- gDNA :
- nBGT :
normalized basal gene transcription
- pTU :
putative transcription units
We thank G. Duport for the technical assistance, N. Gadot from the ANIPATH platform (Faculty of Medicine, RTH Laennec, Lyon) for helping with the hybridization protocol, and V. James for the critical revision of the English language. We also thank the Rhône-Alpes Genopole transcriptomic platform DTAMB for slide spotting. This work was supported by the "Programme Fédérateur INRA de Biologie Intégrative AgroBI 2006", and by the "action bioinformatique INSA de Lyon 2003–2006".
- Allen TE, Price ND, Joyce AR, Palsson BO: Long-range periodic patterns in microbial genomes indicate significant multi-scale chromosomal organization. PLoS Comput Biol. 2006, 2: e2-PubMed CentralPubMedView ArticleGoogle Scholar
- Rocha EP: The replication-related organization of bacterial genomes. Microbiology. 2004, 150: 1609-1627.PubMedView ArticleGoogle Scholar
- Rocha E: Is there a role for replication fork asymmetry in the distribution of genes in bacterial genomes?. Trends Microbiol. 2002, 10: 393-395.PubMedView ArticleGoogle Scholar
- Rocha EP, Danchin A: Gene essentiality determines chromosome organisation in bacteria. Nucleic Acids Res. 2003, 31: 6570-6577.PubMed CentralPubMedView ArticleGoogle Scholar
- Brewer BJ: When polymerases collide: replication and the transcriptional organization of the E. coli chromosome. Cell. 1988, 53: 679-686.PubMedView ArticleGoogle Scholar
- French S: Consequences of replication fork movement through transcription units in vivo. Science. 1992, 258: 1362-1365.PubMedView ArticleGoogle Scholar
- Omont N, Kepes F: Transcription/replication collisions cause bacterial transcription units to be longer on the leading strand of replication. Bioinformatics. 2004, 20: 2719-2725.PubMedView ArticleGoogle Scholar
- Mirkin EV, Mirkin SM: Mechanisms of transcription-replication collisions in bacteria. Mol Cell Biol. 2005, 25: 888-895.PubMed CentralPubMedView ArticleGoogle Scholar
- Lawrence JG: Gene organization: selection, selfishness, and serendipity. Annu Rev Microbiol. 2003, 57: 419-440.PubMedView ArticleGoogle Scholar
- Korbel JO, Jensen LJ, von Mering C, Bork P: Analysis of genomic context: prediction of functional associations from conserved bidirectionally transcribed gene pairs. Nat Biotechnol. 2004, 22: 911-917.PubMedView ArticleGoogle Scholar
- Audit B, Ouzounis CA: From genes to genomes: universal scale-invariant properties of microbial chromosome organisation. J Mol Biol. 2003, 332: 617-633.PubMedView ArticleGoogle Scholar
- Peter BJ, Arsuaga J, Breier AM, Khodursky AB, Brown PO, Cozzarelli NR: Genomic transcriptional response to loss of chromosomal supercoiling in Escherichia coli. Genome Biol. 2004, 5: R87-PubMed CentralPubMedView ArticleGoogle Scholar
- Kepes F: Periodic transcriptional organization of the E. coli genome. J Mol Biol. 2004, 340: 957-964.PubMedView ArticleGoogle Scholar
- Jeong KS, Ahn J, Khodursky AB: Spatial patterns of transcriptional activity in the chromosome of Escherichia coli. Genome Biol. 2004, 5: R86-PubMed CentralPubMedView ArticleGoogle Scholar
- Carpentier AS, Torresani B, Grossmann A, Henaut A: Decoding the nucleoid organisation of Bacillus subtilis and Escherichia coli through gene expression data. BMC Genomics. 2005, 6: 84-PubMed CentralPubMedView ArticleGoogle Scholar
- Charles H, Ishikawa H: Physical and genetic map of the genome of Buchnera, the primary endosymbiont of the pea aphid Acyrthosiphon pisum. J Mol Evol. 1999, 48: 142-150.PubMedView ArticleGoogle Scholar
- 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 USA. 2002, 99: 4454-4458.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Tamas I, Klasson L, Canback B, Naslund AK, Eriksson AS, Wernegreen JJ, Sandstrom JP, Moran NA, Andersson SG: 50 million years of genomic stasis in endosymbiotic bacteria. Science. 2002, 296: 2376-2379.PubMedView ArticleGoogle Scholar
- van Ham RC, Kamerbeek J, Palacios C, Rausell C, Abascal F, Bastolla U, Fernandez JM, Jimenez L, Postigo M, Silva FJ, Tamames J, Viguera E, Latorre A, Valencia A, Moran F, Moya A: Reductive genome evolution in Buchnera aphidicola. Proc Natl Acad Sci USA. 2003, 100: 581-586.PubMed CentralPubMedView ArticleGoogle Scholar
- Perez-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.PubMedView ArticleGoogle Scholar
- Douglas AE: The nutritional physiology of aphids. Adv Insect Physiol. 2003, 31: 73-140.View ArticleGoogle Scholar
- Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF, Gregor J, Davis NW, Kirkpatrick HA, Goeden MA, Rose DJ, Mau B, Shao Y: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277: 1453-1474.PubMedView ArticleGoogle Scholar
- Jordan IK, Rogozin IB, Wolf YI, Koonin EV: Essential genes are more evolutionarily conserved than are nonessential genes in bacteria. Genome Res. 2002, 12: 962-968.PubMed CentralPubMedView ArticleGoogle Scholar
- Gil R, Silva FJ, Pereto J, Moya A: Determination of the core of a minimal bacterial gene set. Microbiol Mol Biol Rev. 2004, 68: 518-537.PubMed CentralPubMedView ArticleGoogle Scholar
- Reymond N, Calevro F, Vinuelas J, Morin N, Rahbe Y, Febvay G, Laugier C, Douglas A, Fayard JM, Charles H: Different levels of transcriptional regulation due to trophic constraints in the reduced genome of Buchnera aphidicola APS. Appl Environ Microbiol. 2006, 72: 7760-7766.PubMed CentralPubMedView ArticleGoogle Scholar
- Talaat AM, Howard ST, Hale W, Lyons R, Garner H, Johnston SA: Genomic DNA standards for gene expression profiling in Mycobacterium tuberculosis. Nucleic Acids Res. 2002, 30: e104-PubMed CentralPubMedView ArticleGoogle Scholar
- Weil MR, Macatee T, Garner HR: Toward a universal standard: comparing two methods for standardizing spotted microarray data. Biotechniques. 2002, 32: 1310-1314.PubMedGoogle Scholar
- Kim H, Zhao B, Snesrud EC, Haas BJ, Town CD, Quackenbush J: Use of RNA and genomic DNA references for inferred comparisons in DNA microarray analyses. Biotechniques. 2002, 33: 924-930.PubMedGoogle Scholar
- Reymond N, Charles H, Duret L, Calevro F, Beslon G, Fayard JM: ROSO: optimizing oligonucleotide probes for microarrays. Bioinformatics. 2004, 20: 271-273.PubMedView ArticleGoogle Scholar
- The Profiling of Escherichia coli Chromosome database. [http://www.shigen.nig.ac.jp/ecoli/pec]
- Balazsi G, Kay KA, Barabasi AL, Oltvai ZN: Spurious spatial periodicity of co-expression in microarray data due to printing design. Nucleic Acids Res. 2003, 31: 4425-4433.PubMed CentralPubMedView ArticleGoogle Scholar
- Kluger Y, Yu H, Qian J, Gerstein M: Relationship between gene co-expression and probe localization on microarray slides. BMC Genomics. 2003, 4: 49-PubMed CentralPubMedView ArticleGoogle Scholar
- Lercher MJ, Hurst LD: Co-expressed yeast genes cluster over a long range but are not regularly spaced. J Mol Biol. 2006, 359: 825-831.PubMedView ArticleGoogle Scholar
- Moran NA, Mira A: The process of genome shrinkage in the obligate symbiont Buchnera aphidicola. Genome Biol. 2001, 2: RESEARCH0054-PubMed CentralPubMedView ArticleGoogle Scholar
- Rocha EP, Danchin A: Essentiality, not expressiveness, drives gene-strand bias in bacteria. Nat Genet. 2003, 34: 377-378.PubMedView ArticleGoogle Scholar
- Price MN, Alm EJ, Arkin AP: Interruptions in gene expression drive highly expressed operons to the leading strand of DNA replication. Nucleic Acids Res. 2005, 33: 3224-3234.PubMed CentralPubMedView ArticleGoogle Scholar
- Shimomura S, Shigenobu S, Morioka M, Ishikawa H: An experimental validation of orphan genes of Buchnera, a symbiont of aphids. Biochem Biophys Res Commun. 2002, 292: 263-267.PubMedView ArticleGoogle Scholar
- Maezawa K, Shigenobu S, Taniguchi H, Kubo T, Aizawa S, Morioka M: Hundreds of flagellar basal bodies cover the cell surface of the endosymbiotic bacterium Buchnera aphidicola sp. strain APS. J Bacteriol. 2006, 188: 6539-6543.PubMed CentralPubMedView ArticleGoogle Scholar
- Delmotte F, Rispe C, Schaber J, Silva FJ, Moya A: Tempo and mode of early gene loss in endosymbiotic bacteria from insects. BMC Evol Biol. 2006, 6: 56-PubMed CentralPubMedView ArticleGoogle Scholar
- Schaber J, Rispe C, Wernegreen J, Buness A, Delmotte F, Silva FJ, Moya A: Gene expression levels influence amino acid usage and evolutionary rates in endosymbiotic bacteria. Gene. 2005, 352: 109-117.PubMedView ArticleGoogle Scholar
- Calevro F, Charles H, Reymond N, Dugas V, Cloarec JP, Bernillon J, Rahbe Y, Febvay G, Fayard JM: Assessment of 35mer amino-modified oligonucleotide based microarray with bacterial samples. J Microbiol Methods. 2004, 57: 207-218.PubMedView ArticleGoogle Scholar
- Charles H, Calevro F, Vinuelas J, Fayard JM, Rahbe Y: Codon usage bias and tRNA over-expression in Buchnera aphidicola after aromatic amino acid nutritional stress on its host Acyrthosiphon pisum. Nucleic Acids Res. 2006, 34: 4583-4592.PubMed CentralPubMedView ArticleGoogle Scholar
- The GenBank database. [http://www.ncbi.nlm.nih.gov]
- The Encyclopedia of Escherichia coli K-12 Genes and Metabolism. [http://ecocyc.org]
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.