The genomic pattern of tDNA operon expression in E. coli

The layered costs and benefits of translational redundancy

The rate and accuracy of translation hinges upon multiple components - including transfer RNA (tRNA) pools, tRNA modifying enzymes, and rRNA molecules - many of which are redundant in terms of gene copy number or function. It has been hypothesized that the redundancy evolves under selection, driven by its impacts on growth rate. However, we lack empirical measurements of the fitness costs and benefits of redundancy, and we have poor a understanding of how this redundancy is organized across components. We manipulated redundancy in multiple translation components of Escherichia coli by deleting 28 tRNA genes, 3 tRNA modifying systems, and 4 rRNA operons in various combinations. We find that redundancy in tRNA pools is beneficial when nutrients are plentiful and costly under nutrient limitation. This nutrient-dependent cost of redundant tRNA genes stems from upper limits to translation capacity and growth rate, and therefore varies as a function of the maximum growth rate attainable in a given nutrient niche. The loss of redundancy in rRNA genes and tRNA modifying enzymes had similar nutrient-dependent fitness consequences. Importantly, these effects are also contingent upon interactions across translation components, indicating a layered hierarchy from copy number of tRNA and rRNA genes to their expression and downstream processing. Overall, our results indicate both positive and negative selection on redundancy in translation components, depending on a species' evolutionary history with feasts and famines.

An optimal growth law for RNA composition and its partial implementation through ribosomal and tRNA gene locations in bacterial genomes

The distribution of cellular resources across bacterial proteins has been quantified through phenomenological growth laws. Here, we describe a complementary bacterial growth law for RNA composition, emerging from optimal cellular resource allocation into ribosomes and ternary complexes. The predicted decline of the tRNA/rRNA ratio with growth rate agrees quantitatively with experimental data. Its regulation appears to be implemented in part through chromosomal localization, as rRNA genes are typically closer to the origin of replication than tRNA genes and thus have increasingly higher gene dosage at faster growth. At the highest growth rates in E. coli, the tRNA/rRNA gene dosage ratio based on chromosomal positions is almost identical to the observed and theoretically optimal tRNA/rRNA expression ratio, indicating that the chromosomal arrangement has evolved to favor maximal transcription of both types of genes at this condition.

Unlike the proteome composition, RNA composition is often assumed to be independent of growth rate in bacteria, despite experimental evidence for a growth rate dependence in many microbes. In this work, we derived a growth-rate dependent optimal tRNA/rRNA concentration ratio by minimizing the combined costs of ribosome and ternary complex at the required protein production rate. The predicted optimal tRNA/rRNA expression ratio, which is a monotonically decreasing function of growth rate, agrees with experimental data for E. coli and other fast-growing microbes. This indicates the existing of an RNA composition growth law. Due to the presence of partially replicated chromosomes, gene dosage is higher for those genes whose DNA is replicated earlier, an effect that becomes stronger at higher growth rates. Because rRNA genes are located closer to origin of replication than tRNA genes in fast-growing species, the tRNA/rRNA gene dosage ratio scales with growth rate in the same direction as the optimal tRNA/rRNA expression ratio. Thus, it appears that the RNA growth law is–at least in part–implemented simply through the genomic positions of tRNA and rRNA genes. This finding indicates that growth rate-dependent optimal resource allocation can influence the genomic organization in bacteria.

Vibrionaceae core, shell and cloud genes are non-randomly distributed on Chr 1: An hypothesis that links the genomic location of genes with their intracellular placement

Background

The genome of Vibrionaceae bacteria, which consists of two circular chromosomes, is replicated in a highly ordered fashion. In fast-growing bacteria, multifork replication results in higher gene copy numbers and increased expression of genes located close to the origin of replication of Chr 1 (ori1). This is believed to be a growth optimization strategy to satisfy the high demand of essential growth factors during fast growth. The relationship between ori1-proximate growth-related genes and gene expression during fast growth has been investigated by many researchers. However, it remains unclear which other gene categories that are present close to ori1 and if expression of all ori1-proximate genes is increased during fast growth, or if expression is selectively elevated for certain gene categories.

Results

We calculated the pangenome of all complete genomes from the Vibrionaceae family and mapped the four pangene categories, core, softcore, shell and cloud, to their chromosomal positions. This revealed that core and softcore genes were found heavily biased towards ori1, while shell genes were overrepresented at the opposite part of Chr 1 (i.e., close to ter1). RNA-seq of Aliivibrio salmonicida and Vibrio natriegens showed global gene expression patterns that consistently correlated with chromosomal distance to ori1. Despite a biased gene distribution pattern, all pangene categories contributed to a skewed expression pattern at fast-growing conditions, whereas at slow-growing conditions, softcore, shell and cloud genes were responsible for elevated expression.

Conclusion

The pangene categories were non-randomly organized on Chr 1, with an overrepresentation of core and softcore genes around ori1, and overrepresentation of shell and cloud genes around ter1. Furthermore, we mapped our gene distribution data on to the intracellular positioning of chromatin described for V. cholerae, and found that core/softcore and shell/cloud genes appear enriched at two spatially separated intracellular regions. Based on these observations, we hypothesize that there is a link between the genomic location of genes and their cellular placement.

The rare lncRNA GOLLD is widespread and structurally conserved among Mycobacterium tRNA arrays

Noncoding RNA (ncRNA) genes produce transcripts involved in a wide range of functions, including catalytic and regulatory functions. Besides, some transcripts have highly complex structures that may impact their activities. Among the largest bacterial ncRNAs, there is the rare GOLLD RNA, which is associated with tRNA genes and supposed to be chromosome- and phage-encoded in specialized groups of bacteria, including those from Lactobacillales and Actinomycetales orders. The only GOLLD structure was inferred from a variety of sequences, including many marine metagenomes. To explore GOLLD RNA in bacterial genomes, we mined the GOLLD gene in thousands of Mycobacterium and virus genomes using Infernal software. We identified this gene in 350 mycobacteria, including megaplasmids, and 39 bacteriophages, mainly in the genomic context of tRNA arrays. Mycobacterium GOLLD genes presented a high diversity and were distributed in three phylogenetic groups: (i) Mycobacterium exclusive; (ii) Mycobacterium and mycobacteriophages; and (iii) mycobacteriophage exclusive. We also determined the GOLLD secondary structure of each group using R2 R software based on GOLLD alignments generated by Infernal software. All GOLLD groups displayed a 3' half conserved structure, including utter E-loops pseudoknots substructures, also shared by non-Mycobacterium GOLLD while the 5' half motif was different among the groups. Here, we showed that the lncRNA GOLLD is widespread in Mycobacterium within tRNA arrays and corroborated the previously predicted GOLLD secondary structure.

Significance and roles of synonymous codon usage in the evolutionary process of Proteus

Proteus spp. bacteria frequently serve as opportunistic pathogens that can infect many animals and show positive survival and existence in various natural environments. The evolutionary pattern of Proteus spp. is an unknown topic, which benefits understanding the different evolutionary dynamics for excellent bacterial adaptation to various environments. Here, the eight whole genomes of different Proteus species were analyzed for the interplay between nucleotide usage and synonymous codon usage. Although the orthologous average nucleotide identity and average nucleotide identity display the genetic diversity of these Proteus species at the genome level, the principal component analysis further shows that these species sustain the specific genetic niche at the aspect of synonymous codon usage patterns. Interestingly, although these Proteus species have A/T rich genes with underrepresented G (guanine) or C (cytosine) at the third codon positions and overrepresented A or T at these positions, some synonymous codons with A or T end are obviously suppressed in usage. The overall codon usage pattern reflected by the effective number of codons (ENC) has a significantly positive correlation with GC3 content (GC content at the third codon position), and ENC has a significantly negative correlation with the adaptation index for these species. These results suggest that the mutation pressure caused by nucleotide composition constraint serves as a dominant evolutionary dynamic driving evolutionary trend of Proteus spp., along with other selections related to natural selection, replication and fine-tune translation, and so on. Taken together, the analyses help to understand the evolutionary interplay between nucleotide and codon usage at the gene level of Proteus.

Exploring tRNA gene cluster in archaea

BACKGROUND

Shared traits between prokaryotes and eukaryotes are helpful in the understanding of the tree of life evolution. In bacteria and eukaryotes, it has been shown a particular organisation of tRNA genes as clusters, but this trait has not been explored in the archaea domain.

OBJECTIVE

Explore the occurrence of tRNA gene clusters in archaea.

METHODS

In-silico analyses of complete and draft archaeal genomes based on tRNA gene isotype and synteny, tRNA gene cluster content and mobilome elements.

FINDINGS

We demonstrated the prevalence of tRNA gene clusters in archaea. tRNA gene clusters, composed of archaeal-type tRNAs, were identified in two Archaea class, Halobacteria and Methanobacteria from Euryarchaeota supergroup. Genomic analyses also revealed evidence of the association between tRNA gene clusters to mobile genetic elements and intra-domain horizontal gene transfer.

MAIN CONCLUSIONS

tRNA gene cluster occurs in the three domains of life, suggesting a role of this type of tRNA gene organisation in the biology of the living organisms.

Beyond the Limits: tRNA Array Units in Mycobacterium Genomes

tRNA array unit, a genomic region presenting an intriguing high tRNA gene number and density, was supposed to occur only in few bacteria phyla, particularly Firmicutes. Here, we identified and characterized an abundance and diversity of tRNA array units in Mycobacterium associated genomes. These genomes comprised chromosome, bacteriophages and plasmids from mycobacteria. Firstly, we had identified 32 tRNA genes organized in an array unit within a mycobacteria plasmid genome and therefore, we hypothesized the presence of such structures in Mycobacterium genus. However, at the time, bioinformatics tools only predict tRNA genes, not characterizing their arrangement as arrays. In order to test our hypothesis, we developed and applied an in-house Perl script that identified tRNA genes organization as an array unit. This survey included a total of 7,670 complete and drafts genomes of Mycobacterium genus, 4312 mycobacteriophage genomes and 40 mycobacteria plasmids. We showed that tRNA array units are abundant in genomes associated to the Mycobacterium genus, mainly in Mycobacterium abscessus complex species, being spread in chromosome, prophage, and plasmid genomes. Moreover, other non-coding RNA species (tmRNA and structured RNA) were also identified in these regions. Our results revealed that tRNA array units are not restrict, as previously assumed, to few bacteria phyla and genomes being present in one of the most diverse bacteria genus. We also provide a bioinformatics tool that allows further exploration of this issue in huge genomic databases. The presence of tRNA array units in plasmids and bacteriophages, associated with horizontal gene transfer, and in a bacteria genus that explores diverse niches, are indicatives that tRNA array units have impact in the bacteria biology.

Transcriptome analysis of Haloquadratum walsbyi: vanity is but the surface

Background

Haloquadratum walsbyi dominates saturated thalassic lakes worldwide where they can constitute up to 80-90% of the total prokaryotic community. Despite the abundance of the enigmatic square-flattened cells, only 7 isolates are currently known with 2 genomes fully sequenced and annotated due to difficulties to grow them under laboratory conditions. We have performed a transcriptomic analysis of one of these isolates, the Spanish strain HBSQ001 in order to investigate gene transcription under light and dark conditions.

Results

Despite a potential advantage for light as additional source of energy, no significant differences were found between light and dark expressed genes. Constitutive high gene expression was observed in genes encoding surface glycoproteins, light mediated proton pumping by bacteriorhodopsin, several nutrient uptake systems, buoyancy and storage of excess carbon. Two low expressed regions of the genome were characterized by a lower codon adaptation index, low GC content and high incidence of hypothetical genes.

Conclusions

Under the extant cultivation conditions, the square hyperhalophile devoted most of its transcriptome towards processes maintaining cell integrity and exploiting solar energy. Surface glycoproteins are essential for maintaining the large surface to volume ratio that facilitates light and organic nutrient harvesting whereas constitutive expression of bacteriorhodopsin warrants an immediate source of energy when light becomes available.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3892-2) contains supplementary material, which is available to authorized users.

Dealing with Gene-Dosage Imbalance during S Phase

DNA replication perturbs the dosage balance between genes that replicate early during S phase and those that replicate late. If propagated to influence protein content, this dosage imbalance could influence cellular functions. In bacteria, mechanisms have evolved to use this imbalance to tune certain processes with the rate of cell growth. By contrast, eukaryotes buffer this dosage imbalance to ensure gene expression homeostasis also during S phase. Here, we outline classical and more recent studies describing how different organisms deal with this replication-dependent dosage imbalance, and describe recent results linking the eukaryotic buffering mechanism to replication-dependent histone acetylation. Finally, we discuss the possible implications of this buffering mechanism and speculate why it is specific to eukaryote cells.

A Comprehensive tRNA Genomic Survey Unravels the Evolutionary History of tRNA Arrays in Prokaryotes

Considering the importance of tRNAs in the translation machinery, scant attention has been paid to tRNA array units defined as genomic regions containing at least 20 tRNA genes with a minimal tRNA gene density of two tRNA genes per kilobase. Our analysis of Acidithiobacillus ferrivorans CF27 and Acidithiobacillus ferrooxidans ATCC 23270(T) genomes showed that both display a tRNA array unit with syntenic conservation which mainly contributed to the tRNA gene redundancy in these two organisms. Our investigations into the occurrence and distribution of tRNA array units revealed that 1) this tRNA organization is limited to few phyla and mainly found in Gram-positive bacteria; and 2) the presence of tRNA arrays favors the redundancy of tRNA genes, in particular those encoding the core tRNA isoacceptors. Finally, comparative array organization revealed that tRNA arrays were acquired through horizontal gene transfer (from Firmicutes or unknown donor), before being subjected to tRNA rearrangements, deletions, and duplications. In Bacilli, the most parsimonious evolutionary history involved two common ancestors and the acquisition of their arrays arose late in evolution, in the genera branches. Functional roles of the array units in organism lifestyle, selective genetic advantage and translation efficiency, as well as the evolutionary advantages of organisms harboring them were proposed. Our study offers new insight into the structural organization and evolution of tRNA arrays in prokaryotic organisms.

Auxiliary tRNAs: large-scale analysis of tRNA genes reveals patterns of tRNA repertoire dynamics

Decoding of all codons can be achieved by a subset of tRNAs. In bacteria, certain tRNA species are mandatory, while others are auxiliary and are variably used. It is currently unknown how this variability has evolved and whether it provides an adaptive advantage. Here we shed light on the subset of auxiliary tRNAs, using genomic data from 319 bacteria. By reconstructing the evolution of tRNAs we show that the auxiliary tRNAs are highly dynamic, being frequently gained and lost along the phylogenetic tree, with a clear dominance of loss events for most auxiliary tRNA species. We reveal distinct co-gain and co-loss patterns for subsets of the auxiliary tRNAs, suggesting that they are subjected to the same selection forces. Controlling for phylogenetic dependencies, we find that the usage of these tRNA species is positively correlated with GC content and may derive directly from nucleotide bias or from preference of Watson-Crick codon-anticodon interactions. Our results highlight the highly dynamic nature of these tRNAs and their complicated balance with codon usage.

Ancestral genome organization: an alignment approach

We present a comparative genomics approach for inferring ancestral genome organization and evolutionary scenarios, based on present-day genomes represented as ordered gene sequences with duplicates. We develop our methodology for a model of evolution restricted to duplication and loss, and then show how to extend it to other content-modifying operations, and to inversions. From a combinatorial point of view, the main consequence of ignoring rearrangements is the possibility of formulating the problem as an alignment problem. On the other hand, duplications and losses are asymmetric operations that are applicable to one of the two aligned sequences. Consequently, an ancestral genome can directly be inferred from a duplication-loss scenario attached to a given alignment. Although alignments are a priori simpler to handle than rearrangements, we show that a direct approach based on dynamic programming leads, at best, to an efficient heuristic. We present an exact pseudo-boolean linear programming algorithm to search for the optimal alignment along with an optimal scenario of duplications and losses. Although exponential in the worst case, we show low running times on real datasets as well as synthetic data. We apply our algorithm (*) in a phylogenetic context to the evolution of stable RNA (tRNA and rRNA) gene content and organization in Bacillus genomes. Our results lead to various biological insights, such as rates of ribosomal RNA proliferation among lineages, their role in altering tRNA gene content, and evidence of tRNA class conversion.

Back to the origin: reconsidering replication, transcription, epigenetics, and cell cycle control

In bacteria, replication is a carefully orchestrated event that unfolds the same way for each bacterium and each cell division. The process of DNA replication in bacteria optimizes cell growth and coordinates high levels of simultaneous replication and transcription. In metazoans, the organization of replication is more enigmatic. The lack of a specific sequence that defines origins of replication has, until recently, severely limited our ability to define the organizing principles of DNA replication. This question is of particular importance as emerging data suggest that replication stress is an important contributor to inherited genetic damage and the genomic instability in tumors. We consider here the replication program in several different organisms including recent genome-wide analyses of replication origins in humans. We review recent studies on the role of cytosine methylation in replication origins, the role of transcriptional looping and gene gating in DNA replication, and the role of chromatin's 3-dimensional structure in DNA replication. We use these new findings to consider several questions surrounding DNA replication in metazoans: How are origins selected? What is the relationship between replication and transcription? How do checkpoints inhibit origin firing? Why are there early and late firing origins? We then discuss whether oncogenes promote cancer through a role in DNA replication and whether errors in DNA replication are important contributors to the genomic alterations and gene fusion events observed in cancer. We conclude with some important areas for future experimentation.

On ribosome load, codon bias and protein abundance

Different codons encoding the same amino acid are not used equally in protein-coding sequences. In bacteria, there is a bias towards codons with high translation rates. This bias is most pronounced in highly expressed proteins, but a recent study of synthetic GFP-coding sequences did not find a correlation between codon usage and GFP expression, suggesting that such correlation in natural sequences is not a simple property of translational mechanisms. Here, we investigate the effect of evolutionary forces on codon usage. The relation between codon bias and protein abundance is quantitatively analyzed based on the hypothesis that codon bias evolved to ensure the efficient usage of ribosomes, a precious commodity for fast growing cells. An explicit fitness landscape is formulated based on bacterial growth laws to relate protein abundance and ribosomal load. The model leads to a quantitative relation between codon bias and protein abundance, which accounts for a substantial part of the observed bias for E. coli. Moreover, by providing an evolutionary link, the ribosome load model resolves the apparent conflict between the observed relation of protein abundance and codon bias in natural sequences and the lack of such dependence in a synthetic gfp library. Finally, we show that the relation between codon usage and protein abundance can be used to predict protein abundance from genomic sequence data alone without adjustable parameters.

Search for poly(A) polymerase targets in E. coli reveals its implication in surveillance of Glu tRNA processing and degradation of stable RNAs

Polyadenylation is a universal post-transcriptional modification involved in degradation and quality control of bacterial RNAs. In Escherichia coli, it is admitted that any accessible RNA 3' end can be tagged by a poly(A) tail for decay. However, we do not have yet an overall view of the population of polyadenylated molecules. The sampling of polyadenylated RNAs presented here demonstrates that rRNA fragments and tRNA precursors originating from the internal spacer regions of the rrn operons, in particular, rrnB are abundant poly(A) polymerase targets. Focused analysis showed that Glu tRNA precursors originating from the rrnB and rrnG transcripts exhibit long 3' trailers that are primarily removed by PNPase and to a lesser extent by RNase II and poly(A) polymerase. Moreover, 3' trimming by exoribonucleases precedes 5' end maturation by RNase P. Interestingly, characterization of RNA fragments that accumulate in a PNPase deficient strain showed that Glu tRNA precursors still harbouring the 5' leader can be degraded by a 3' to 5' quality control pathway involving poly(A) polymerase. This demonstrates that the surveillance of tRNA maturation described for a defective tRNA also applies to a wild-type tRNA.

The influence of anticodon-codon interactions and modified bases on codon usage bias in bacteria

Most transfer RNAs (tRNAs) can translate more than one synonymous codon, and most codons can be translated by more than one isoacceptor tRNA. The rates of translation of synonymous codons are dependent on the concentrations of the tRNAs and on the rates of pairing of each anticodon-codon combination. Translational selection causes a significant bias in codon frequencies in highly expressed genes in most bacteria. By comparing codon frequencies in high and low-expression genes, we determine which codons are preferred for each amino acid in a large sample of bacterial genomes. We relate this to the number of copies of each tRNA gene in each genome. In two-codon families, preferred codons have Watson-Crick pairs (GC and AU) between the third codon base and the wobble base of the anticodon rather than GU pairs. This suggests that these combinations are more rapidly recognized by the ribosome. In contrast, in four-codon families, preferred codons do not correspond to Watson-Crick rules. In some cases, a wobble-U tRNA can pair with all four codons. In these cases, A and U codons are preferred over G and C. This indicates that the nonstandard UU combination appears to be translated surprisingly well. Differences in modified bases at the wobble position of the anticodon appear to be responsible for the differences in behavior of tRNAs in two- and four-codon families. We discuss the way changes in the bases in the anticodon influence both the speed and the accuracy of translation. The number of tRNA gene copies and the strength of translational selection correlate with the growth rate of the organism, as we would expect if the primary cause of translational selection in bacteria is the requirement to optimize the speed of protein production.

The systemic imprint of growth and its uses in ecological (meta)genomics

Microbial minimal generation times range from a few minutes to several weeks. They are evolutionarily determined by variables such as environment stability, nutrient availability, and community diversity. Selection for fast growth adaptively imprints genomes, resulting in gene amplification, adapted chromosomal organization, and biased codon usage. We found that these growth-related traits in 214 species of bacteria and archaea are highly correlated, suggesting they all result from growth optimization. While modeling their association with maximal growth rates in view of synthetic biology applications, we observed that codon usage biases are better correlates of growth rates than any other trait, including rRNA copy number. Systematic deviations to our model reveal two distinct evolutionary processes. First, genome organization shows more evolutionary inertia than growth rates. This results in over-representation of growth-related traits in fast degrading genomes. Second, selection for these traits depends on optimal growth temperature: for similar generation times purifying selection is stronger in psychrophiles, intermediate in mesophiles, and lower in thermophiles. Using this information, we created a predictor of maximal growth rate adapted to small genome fragments. We applied it to three metagenomic environmental samples to show that a transiently rich environment, as the human gut, selects for fast-growers, that a toxic environment, as the acid mine biofilm, selects for low growth rates, whereas a diverse environment, like the soil, shows all ranges of growth rates. We also demonstrate that microbial colonizers of babies gut grow faster than stabilized human adults gut communities. In conclusion, we show that one can predict maximal growth rates from sequence data alone, and we propose that such information can be used to facilitate the manipulation of generation times. Our predictor allows inferring growth rates in the vast majority of uncultivable prokaryotes and paves the way to the understanding of community dynamics from metagenomic data.

Microbial minimal generation times vary from a few minutes to several weeks. The reasons for this disparity have been thought to lie on different life-history strategies: fast-growing microbes grow extremely fast in rich media, but are less capable of dealing with stress and/or poor nutrient conditions. Prokaryotes have evolved a set of genomic traits to grow fast, including biased codon usage and transient or permanent gene multiplication for dosage effects. Here, we studied the relative role of these traits and show they can be used to predict minimal generation times from the genomic data of the vast majority of microbes that cannot be cultivated. We show that this inference can also be made with incomplete genomes and thus be applied to metagenomic data to test hypotheses about the biomass productivity of biotopes and the evolution of microbiota in the human gut after birth. Our results also allow a better understanding of the co-evolution between growth rates and genomic traits and how they can be manipulated in synthetic biology. Growth rates have been a key variable in microbial physiology studies in the last century, and we show how intimately they are linked with genome organization and prokaryotic ecology.

The relationship between DNA replication and human genome organization

Assessment of the impact of DNA replication on genome architecture in Eukaryotes has long been hampered by the scarcity of experimental data. Recent work, relying on computational predictions of origins of replication, suggested that replication might be a major determinant of gene organization in human (Huvet et al. 2007. Human gene organization driven by the coordination of replication and transcription. Genome Res. 17:1278-1285). Here, we address this question by analyzing the first large-scale data set of experimentally determined origins of replication in human: 283 origins identified in HeLa cells, in 1% of the genome covered by ENCODE regions (Cadoret et al. 2008. Genome-wide studies highlight indirect links between human replication origins and gene regulation. Proc Natl Acad Sci USA. 105:15837-15842). We show that origins of replication are not randomly distributed as they display significant overlap with promoter regions and CpG islands. The hypothesis of a selective pressure to avoid frontal collisions between replication and transcription polymerases is not supported by experimental data as we find no evidence for gene orientation bias in the proximity of origins of replication. The lack of a significant orientation bias remains manifest even when considering only genes expressed at a high rate, or in a wide number of tissues, and is not affected by the regional replication timing. Gene expression breadth does not appear to be correlated with the distance from the origins of replication. We conclude that the impact of DNA replication on human genome organization is considerably weaker than previously proposed.

tDNA locus polymorphism and ecto-chromosomal DNA insertion hot-spots are related to the phylogenetic group of Escherichia coli strains

tRNA-encoding genes (tDNA) are known hot-spots for the integration of ecto-chromosomal DNA (ECDNA) including genomic islands. However, only a few loci are currently known to be targeted by such insertions in Escherichia coli. A PCR-based screening of tDNA integrity was therefore performed on a collection of E. coli strains in order to identify tDNA loci that are most frequently intact and those that are preferred ECDNA insertion sites. It was shown that only a subset of tDNAs were hot-spots for ECDNA insertions, and that the majority of loci were never targeted by such insertions. Polycistronic tDNAs, highly transcribed tDNAs or tDNAs encoding tRNAs recognizing frequently used codons were generally not targeted by ECDNA insertions. Most interestingly, strains of different ECOR groups showed different patterns of tDNA loci polymorphism. More subtle differences were also observed between strains of different pathotypes.

Codon usage bias and tRNA over-expression in Buchnera aphidicola after aromatic amino acid nutritional stress on its host Acyrthosiphon pisum

Codon usage bias and relative abundances of tRNA isoacceptors were analysed in the obligate intracellular symbiotic bacterium, Buchnera aphidicola from the aphid Acyrthosiphon pisum, using a dedicated 35mer oligonucleotide microarray. Buchnera is archetypal of organisms living with minimal metabolic requirements and presents a reduced genome with high-evolutionary rate. Codonusage in Buchnera has been overcome by the high mutational bias towards AT bases. However, several lines of evidence for codon usage selection are given here. A significant correlation was found between tRNA relative abundances and codon composition of Buchnera genes. A significant codon usage bias was found for the choice of rare codons in Buchnera: C-ending codons are preferred in highly expressed genes, whereas G-ending codons are avoided. This bias is not explained by GC skew in the bacteria and might correspond to a selection for perfect matching between codon-anticodon pairs for some essential amino acids in Buchnera proteins. Nutritional stress applied to the aphid host induced a significant overexpression of most of the tRNA isoacceptors in bacteria. Although, molecular regulation of the tRNA operons in Buchnera was not investigated, a correlation between relative expression levels and organization in transcription unit was found in the genome of Buchnera.

Replication-associated gene dosage effects shape the genomes of fast-growing bacteria but only for transcription and translation genes

The bidirectional replication of bacterial genomes leads to transient gene dosage effects. Here, we show that such effects shape the chromosome organisation of fast-growing bacteria and that they correlate strongly with maximal growth rate. Surprisingly the predicted maximal number of replication rounds shows little if any phylogenetic inertia, suggesting that it is a very labile trait. Yet, a combination of theoretical and statistical analyses predicts that dozens of replication forks may be simultaneously present in the cells of certain species. This suggests a strikingly efficient management of the replication apparatus, of replication fork arrests and of chromosome segregation in such cells. Gene dosage effects strongly constrain the position of genes involved in translation and transcription, but not other highly expressed genes. The relative proximity of the former genes to the origin of replication follows the regulatory dependencies observed under exponential growth, as the bias is stronger for RNA polymerase, then rDNA, then ribosomal proteins and tDNA. Within tDNAs we find that only the positions of the previously proposed 'ubiquitous' tRNA, which translate the most frequent codons in highly expressed genes, show strong signs of selection for gene dosage effects. Finally, we provide evidence for selection acting upon genome organisation to take advantage of gene dosage effects by identifying a positive correlation between genome stability and the number of simultaneous replication rounds. We also show that gene dosage effects can explain the over-representation of highly expressed genes in the largest replichore of genomes containing more than one chromosome. Together, these results demonstrate that replication-associated gene dosage is an important determinant of chromosome organisation and dynamics, especially among fast-growing bacteria.