The replication-related organization of bacterial genomes

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

SpoIIIE strips proteins off the DNA during chromosome translocation

The FtsK/SpoIIIE family of DNA transporters are responsible for translocating missegregated chromosomes after the completion of cell division. An extreme example of this post-cytokinetic DNA segregation occurs during spore formation in the bacterium Bacillus subtilis, where SpoIIIE pumps three-quarters of the chromosome (>3 megabases) into one of the two daughter cells. Here, we investigate the fate of the proteins associated with the translocated DNA. Taking advantage of several unique features of Bacillus sporulation, we demonstrate that RNA polymerase, transcription factors, and chromosome remodeling proteins are stripped off the DNA during translocation of the chromosome into the forespore compartment. Furthermore, we show that in vitro the soluble ATPase domain of SpoIIIE can displace RNA polymerase bound to DNA, suggesting that SpoIIIE alone is capable of this wire-stripping activity. Our data suggest that the bulk of the forespore chromosome is translocated naked into the forespore compartment. We propose that the translocation-stripping activity of SpoIIIE plays a key role in reprogramming developmental gene expression in the forespore.

Identification of replication origins in prokaryotic genomes

The availability of hundreds of complete bacterial genomes has created new challenges and simultaneously opportunities for bioinformatics. In the area of statistical analysis of genomic sequences, the studies of nucleotide compositional bias and gene bias between strands and replichores paved way to the development of tools for prediction of bacterial replication origins. Only a few (about 20) origin regions for eubacteria and archaea have been proven experimentally. One reason for that may be that this is now considered as an essentially bioinformatics problem, where predictions are sufficiently reliable not to run labor-intensive experiments, unless specifically needed. Here we describe the main existing approaches to the identification of replication origin (oriC) and termination (terC) loci in prokaryotic chromosomes and characterize a number of computational tools based on various skew types and other types of evidence. We also classify the eubacterial and archaeal chromosomes by predictability of their replication origins using skew plots. Finally, we discuss possible combined approaches to the identification of the oriC sites that may be used to improve the prediction tools, in particular, the analysis of DnaA binding sites using the comparative genomic methods.

Plastic architecture of bacterial genome revealed by comparative genomics of Photorhabdus variants

BACKGROUND

The phenotypic consequences of large genomic architecture modifications within a clonal bacterial population are rarely evaluated because of the difficulties associated with using molecular approaches in a mixed population. Bacterial variants frequently arise among Photorhabdus luminescens, a nematode-symbiotic and insect-pathogenic bacterium. We therefore studied genome plasticity within Photorhabdus variants.

RESULTS

We used a combination of macrorestriction and DNA microarray experiments to perform a comparative genomic study of different P. luminescens TT01 variants. Prolonged culturing of TT01 strain and a genomic variant, collected from the laboratory-maintained symbiotic nematode, generated bacterial lineages composed of primary and secondary phenotypic variants and colonial variants. The primary phenotypic variants exhibit several characteristics that are absent from the secondary forms. We identify substantial plasticity of the genome architecture of some variants, mediated mainly by deletions in the 'flexible' gene pool of the TT01 reference genome and also by genomic amplification. We show that the primary or secondary phenotypic variant status is independent from global genomic architecture and that the bacterial lineages are genomic lineages. We focused on two unusual genomic changes: a deletion at a new recombination hotspot composed of long approximate repeats; and a 275 kilobase single block duplication belonging to a new class of genomic duplications.

CONCLUSION

Our findings demonstrate that major genomic variations occur in Photorhabdus clonal populations. The phenotypic consequences of these genomic changes are cryptic. This study provides insight into the field of bacterial genome architecture and further elucidates the role played by clonal genomic variation in bacterial genome evolution.

Dynamics of genome rearrangement in bacterial populations

Genome structure variation has profound impacts on phenotype in organisms ranging from microbes to humans, yet little is known about how natural selection acts on genome arrangement. Pathogenic bacteria such as Yersinia pestis, which causes bubonic and pneumonic plague, often exhibit a high degree of genomic rearrangement. The recent availability of several Yersinia genomes offers an unprecedented opportunity to study the evolution of genome structure and arrangement. We introduce a set of statistical methods to study patterns of rearrangement in circular chromosomes and apply them to the Yersinia. We constructed a multiple alignment of eight Yersinia genomes using Mauve software to identify 78 conserved segments that are internally free from genome rearrangement. Based on the alignment, we applied Bayesian statistical methods to infer the phylogenetic inversion history of Yersinia. The sampling of genome arrangement reconstructions contains seven parsimonious tree topologies, each having different histories of 79 inversions. Topologies with a greater number of inversions also exist, but were sampled less frequently. The inversion phylogenies agree with results suggested by SNP patterns. We then analyzed reconstructed inversion histories to identify patterns of rearrangement. We confirm an over-representation of "symmetric inversions"-inversions with endpoints that are equally distant from the origin of chromosomal replication. Ancestral genome arrangements demonstrate moderate preference for replichore balance in Yersinia. We found that all inversions are shorter than expected under a neutral model, whereas inversions acting within a single replichore are much shorter than expected. We also found evidence for a canonical configuration of the origin and terminus of replication. Finally, breakpoint reuse analysis reveals that inversions with endpoints proximal to the origin of DNA replication are nearly three times more frequent. Our findings represent the first characterization of genome arrangement evolution in a bacterial population evolving outside laboratory conditions. Insight into the process of genomic rearrangement may further the understanding of pathogen population dynamics and selection on the architecture of circular bacterial chromosomes.

Highly plastic genome of Microcystis aeruginosa PCC 7806, a ubiquitous toxic freshwater cyanobacterium

Background

The colonial cyanobacterium Microcystis proliferates in a wide range of freshwater ecosystems and is exposed to changing environmental factors during its life cycle. Microcystis blooms are often toxic, potentially fatal to animals and humans, and may cause environmental problems. There has been little investigation of the genomics of these cyanobacteria.

Results

Deciphering the 5,172,804 bp sequence of Microcystis aeruginosa PCC 7806 has revealed the high plasticity of its genome: 11.7% DNA repeats containing more than 1,000 bases, 6.8% putative transposases and 21 putative restriction enzymes. Compared to the genomes of other cyanobacterial lineages, strain PCC 7806 contains a large number of atypical genes that may have been acquired by lateral transfers. Metabolic pathways, such as fermentation and a methionine salvage pathway, have been identified, as have genes for programmed cell death that may be related to the rapid disappearance of Microcystis blooms in nature. Analysis of the PCC 7806 genome also reveals striking novel biosynthetic features that might help to elucidate the ecological impact of secondary metabolites and lead to the discovery of novel metabolites for new biotechnological applications. M. aeruginosa and other large cyanobacterial genomes exhibit a rapid loss of synteny in contrast to other microbial genomes.

Conclusion

Microcystis aeruginosa PCC 7806 appears to have adopted an evolutionary strategy relying on unusual genome plasticity to adapt to eutrophic freshwater ecosystems, a property shared by another strain of M. aeruginosa (NIES-843). Comparisons of the genomes of PCC 7806 and other cyanobacterial strains indicate that a similar strategy may have also been used by the marine strain Crocosphaera watsonii WH8501 to adapt to other ecological niches, such as oligotrophic open oceans.

Investigations of oligonucleotide usage variance within and between prokaryotes

Oligonucleotide usage in archaeal and bacterial genomes can be linked to a number of properties, including codon usage (trinucleotides), DNA base-stacking energy (dinucleotides), and DNA structural conformation (di- to tetranucleotides). We wanted to assess the statistical information potential of different DNA 'word-sizes' and explore how oligonucleotide frequencies differ in coding and non-coding regions. In addition, we used oligonucleotide frequencies to investigate DNA composition and how DNA sequence patterns change within and between prokaryotic organisms. Among the results found was that prokaryotic chromosomes can be described by hexanucleotide frequencies, suggesting that prokaryotic DNA is predominantly short range correlated, i.e., information in prokaryotic genomes is encoded in short oligonucleotides. Oligonucleotide usage varied more within AT-rich and host-associated genomes than in GC-rich and free-living genomes, and this variation was mainly located in non-coding regions. Bias (selectional pressure) in tetranucleotide usage correlated with GC content, and coding regions were more biased than non-coding regions. Non-coding regions were also found to be approximately 5.5% more AT-rich than coding regions, on average, in the 402 chromosomes examined. Pronounced DNA compositional differences were found both within and between AT-rich and GC-rich genomes. GC-rich genomes were more similar and biased in terms of tetranucleotide usage in non-coding regions than AT-rich genomes. The differences found between AT-rich and GC-rich genomes may possibly be attributed to lifestyle, since tetranucleotide usage within host-associated bacteria was, on average, more dissimilar and less biased than free-living archaea and bacteria.

The lifestyle of Corynebacterium urealyticum derived from its complete genome sequence established by pyrosequencing

Corynebacterium urealyticum is a lipid-requiring, urealytic bacterium of the human skin flora that has been recognized as causative agent of urinary tract infections. We report the analysis of the complete genome sequence of C. urealyticum DSM7109, which was initially recovered from a patient with alkaline-encrusted cystitis. The genome sequence was determined by a combination of pyrosequencing and Sanger technology. The chromosome of C. urealyticum DSM7109 has a size of 2,369,219bp and contains 2024 predicted coding sequences, of which 78% were considered as orthologous with genes in the Corynebacterium jeikeium K411 genome. Metabolic analysis of the lipid-requiring phenotype revealed the absence of a fatty acid synthase gene and the presence of a beta-oxidation pathway along with a large repertoire of auxillary genes for the degradation of exogenous fatty acids. A urease locus with the gene order ureABCEFGD may play a pivotal role in virulence of C. urealyticum by the alkalinization of human urine and the formation of struvite stones. Multidrug resistance of C. urealyticum DSM7109 is mediated by transposable elements, conferring resistances to macrolides, lincosamides, ketolides, aminoglycosides, chloramphenicol, and tetracycline. The complete genome sequence of C. urealyticum revealed a detailed picture of the lifestyle of this opportunistic human pathogen.

Improved determination of plasmid copy number using quantitative real-time PCR for monitoring fermentation processes

BACKGROUND

Recombinant protein production in Escherichia coli cells is a complex process, where among other parameters, plasmid copy number, structural and segregational stability of plasmid have an important impact on the success of productivity. It was recognised that a method for accurate and rapid quantification of plasmid copy number is necessary for optimization and better understanding of this process. Lately, qPCR is becoming the method of choice for this purpose. In the presented work, an improved qPCR method adopted for PCN determination in various fermentation processes was developed.

RESULTS

To avoid experimental errors arising from irreproducible DNA isolation, whole cells, treated by heating at 95 degrees C for 10 minutes prior to storage at -20 degrees C, were used as a template source. Relative quantification, taking into account different amplification efficiencies of amplicons for chromosome and plasmid, was used in the PCN calculation. The best reproducibility was achieved when the efficiency estimated for specific amplicon, obtained within one run, was averaged. It was demonstrated that the quantification range of 2 log units (100 to 10000 bacteria per well) enable quantification in each time point during fermentation. The method was applied to study PCN variation in fermentation at 25 degrees C and the correlation between PCN and protein accumulation was established.

CONCLUSION

Using whole cells as a template source and relative quantification considering different PCR amplification efficiencies are significant improvements of the qPCR method for PCN determination. Due to the approaches used, the method is suitable for PCN determination in fermentation processes using various media and conditions.

Reliability and applications of statistical methods based on oligonucleotide frequencies in bacterial and archaeal genomes

Background

The increasing number of sequenced prokaryotic genomes contains a wealth of genomic data that needs to be effectively analysed. A set of statistical tools exists for such analysis, but their strengths and weaknesses have not been fully explored. The statistical methods we are concerned with here are mainly used to examine similarities between archaeal and bacterial DNA from different genomes. These methods compare observed genomic frequencies of fixed-sized oligonucleotides with expected values, which can be determined by genomic nucleotide content, smaller oligonucleotide frequencies, or be based on specific statistical distributions. Advantages with these statistical methods include measurements of phylogenetic relationship with relatively small pieces of DNA sampled from almost anywhere within genomes, detection of foreign/conserved DNA, and homology searches. Our aim was to explore the reliability and best suited applications for some popular methods, which include relative oligonucleotide frequencies (ROF), di- to hexanucleotide zero'th order Markov methods (ZOM) and 2.order Markov chain Method (MCM). Tests were performed on distant homology searches with large DNA sequences, detection of foreign/conserved DNA, and plasmid-host similarity comparisons. Additionally, the reliability of the methods was tested by comparing both real and random genomic DNA.

Results

Our findings show that the optimal method is context dependent. ROFs were best suited for distant homology searches, whilst the hexanucleotide ZOM and MCM measures were more reliable measures in terms of phylogeny. The dinucleotide ZOM method produced high correlation values when used to compare real genomes to an artificially constructed random genome with similar %GC, and should therefore be used with care. The tetranucleotide ZOM measure was a good measure to detect horizontally transferred regions, and when used to compare the phylogenetic relationships between plasmids and hosts, significant correlation (R² = 0.4) was found with genomic GC content and intra-chromosomal homogeneity.

Conclusion

The statistical methods examined are fast, easy to implement, and powerful for a number of different applications involving genomic sequence comparisons. However, none of the measures examined were superior in all tests, and therefore the choice of the statistical method should depend on the task at hand.

The genome of Pelotomaculum thermopropionicum reveals niche-associated evolution in anaerobic microbiota

The anaerobic biodegradation of organic matter is accomplished by sequential syntrophic catabolism by microbes in different niches. Pelotomaculum thermopropionicum is a representative syntrophic bacterium that catalyzes the intermediate bottleneck step in the anaerobic-biodegradation process, whereby volatile fatty acids (VFAs) and alcohols produced by upstream fermenting bacteria are converted to acetate, hydrogen, and carbon dioxide (substrates for downstream methanogenic archaea). To reveal genomic features that contribute to our understanding of the ecological niche and evolution of P. thermopropionicum, we sequenced its 3,025,375-bp genome and performed comparative analyses with genomes of other community members available in the databases. In the genome, 2920 coding sequences (CDSs) were identified. These CDSs showed a distinct distribution pattern in the functional categories of the Clusters of Orthologous Groups database, which is considered to reflect the niche of this organism. P. thermopropionicum has simple catabolic pathways, in which the propionate-oxidizing methylmalonyl-CoA pathway constitutes the backbone and is linked to several peripheral pathways. Genes for most of the important catabolic enzymes are physically linked to those for PAS-domain-containing regulators, suggesting that the catabolic pathways are regulated in response to environmental conditions and/or global cellular situations rather than specific substrates. Comparative analyses of codon usages revealed close evolutionary relationships between P. thermopropionicum and other niche members, while it was distant from phylogenetically related sugar-fermenting bacteria. These analyses suggest that P. thermopropionicum has evolved as a syntrophy specialist by interacting with niche-associated microbes.

Statistical evaluation of genetic footprinting data

As transposomics is extended to genome scale, appropriate statistical methods need to be developed to assign significance to gene essentiality. In this chapter, the author presents a set of steps that, together with genome-scale insertion data and the complete genome sequence of a prokaryote, can be used to classify the genes of the organism as either "essential" or "nonessential."

Selective chromosome amplification in Vibrio cholerae

Most bacteria have one chromosome but some have more than one, as is common in eukaryotes. How multiple chromosomes are maintained in bacteria remains largely obscure. Here we have examined the behaviour of the two Vibrio cholerae chromosomes as a function of growth rate. At slow growth rates, both chromosomes were maintained at copy numbers of one to two per cell. Increasing the growth rate by nutritional shift-up amplified the origin-proximal DNA of the larger chromosome (chrI) to four copies per cell, but not that of the smaller chrII. The latter was amplified when its specific initiator was supplied in excess or a specific negative regulator was deleted. The growth rate-insensitive behaviour of chrII, whose origin is similar to origins of members of a major class of plasmids, was shared by some but not all of several representative plasmids tested in V. cholerae. Also, unlike plasmid replication, chrII replication is known to be initiated at a specific stage of the cell cycle. Raising chrII copy number decreased growth rate, suggesting that this chromosome might serve as a repository for necessary but potentially deleterious genes.

Separating the effects of mutation and selection in producing DNA skew in bacterial chromosomes

Background

Many bacterial chromosomes display nucleotide asymmetry, or skew, between the leading and lagging strands of replication. Mutational differences between these strands result in an overall pattern of skew that is centered about the origin of replication. Such a pattern could also arise from selection coupled with a bias for genes coded on the leading strand. The relative contributions of selection and mutation in producing compositional skew are largely unknown.

Results

We describe a model to quantify the contribution of mutational differences between the leading and lagging strands in producing replication-induced skew. When the origin and terminus of replication are known, the model can be used to estimate the relative accumulation of G over C and of A over T on the leading strand due to replication effects in a chromosome with bidirectional replication arms. The model may also be implemented in a maximum likelihood framework to estimate the locations of origin and terminus. We find that our estimations for the origin and terminus agree very well with the location of genes that are thought to be associated with the replication origin. This indicates that our model provides an accurate, objective method of determining the replication arms and also provides support for the hypothesis that these genes represent an ancestral cluster of origin-associated genes.

Conclusion

The model has several advantages over other methods of analyzing genome skew. First, it quantifies the role of mutation in generating skew so that its effect on composition, for example codon bias, can be assessed. Second, it provides an objective method for locating origin and terminus, one that is based on chromosome-wide accumulation of leading vs lagging strand nucleotide differences. Finally, the model has the potential to be utilized in a maximum likelihood framework in order to analyze the effect of chromosome rearrangements on nucleotide composition.

From GC skews to wavelets: a gentle guide to the analysis of compositional asymmetries in genomic data

Compositional asymmetries are pervasive in DNA sequences. They are the result of the asymmetric interactions between DNA and cellular mechanisms such as replication and transcription. Here, we review many of the methods that have been proposed over the years to analyse compositional asymmetries in DNA sequences. Among these we list GC skews, oligonucleotide skews and wavelets, which among other uses have been extensively employed to delimitate origins and termini of replication in genomes. We also review the use of multivariate methods, such as factorial correspondence analysis, discriminant analysis and analysis of variance, which allow assigning compositional strand asymmetries to the different biological processes shaping sequence composition. Finally, we review methods that have been used to infer substitution matrices and allow understanding the mutational processes underlying strand asymmetry. We focus on replication asymmetries because they have been more thoroughly studied, but the methods may be adapted, and often are, to other problems. Although strand asymmetry has been studied more frequently through compositional skews of nucleotides or oligonucleotides, we recall that, depending on the goal of the analysis, other methods may be more appropriate to answer certain biological questions. We also refer to programs freely available to analyse strand asymmetry.

Quantitative analysis of mutation and selection pressures on base composition skews in bacterial chromosomes

Background

Most bacterial chromosomes exhibit asymmetry of base composition with respect to leading vs. lagging strands (GC and AT skews). These skews reflect mainly those in protein coding sequences, which are driven by asymmetric mutation pressures during replication and transcription (notably asymmetric cytosine deamination) plus subsequent selection for preferred structures, signals, amino acid or codons. The transcription-associated effects but not the replication-associated effects contribute to the overall skews through the uneven distribution of the coding sequences on the leading and lagging strands.

Results

Analysis of 185 representative bacterial chromosomes showed diverse and characteristic patterns of skews among different clades. The base composition skews in the coding sequences were used to derive quantitatively the effect of replication-driven mutation plus subsequent selection ('replication-associated pressure', RAP), and the effect of transcription-driven mutation plus subsequent selection at translation level ('transcription-associate pressure', TAP). While different clades exhibit distinct patterns of RAP and TAP, RAP is absent or nearly absent in some bacteria, but TAP is present in all. The selection pressure at the translation level is evident in all bacteria based on the analysis of the skews at the three codon positions. Contribution of asymmetric cytosine deamination was found to be weak to TAP in most phyla, and strong to RAP in all the Proteobacteria but weak in most of the Firmicutes. This possibly reflects the differences in their chromosomal replication machineries. A strong negative correlation between TAP and G+C content and between TAP and chromosomal size were also revealed.

Conclusion

The study reveals the diverse mutation and selection forces associated with replication and transcription in various groups of bacteria that shape the distinct patterns of base composition skews in the chromosomes during evolution. Some closely relative species with distinct base composition parameters are uncovered in this study, which also provides opportunities for comparative bioinformatic and genetic investigations to uncover the underlying principles for mutation and selection.

Replication-associated purine asymmetry may contribute to strand-biased gene distribution

Among prokaryotic genomes, the distribution of genes on the leading and lagging strands of the replication fork is known to be biased. Several hypotheses explaining this strand-biased gene distribution (SGD) have been proposed, but none have been tested or supported by sufficient data analyses. In this work we have analyzed 211 prokaryotic genomes in terms of compositional strand asymmetries and the presence or absence of polC and have found that SGD correlates not only with polC, but also with purine asymmetry (PAS). Furthermore, SGD, PAS, and polC are all features associated with a group of low-GC, gram-positive bacteria (Firmicutes). We conclude that PAS is a characteristic of organisms with a heterodimeric DNA polymerase III alpha-subunit constituted by polC and dnaE, which may play a direct role in the maintenance of SGD.

The unconventional Xer recombination machinery of Streptococci/Lactococci

Homologous recombination between circular sister chromosomes during DNA replication in bacteria can generate chromosome dimers that must be resolved into monomers prior to cell division. In Escherichia coli, dimer resolution is achieved by site-specific recombination, Xer recombination, involving two paralogous tyrosine recombinases, XerC and XerD, and a 28-bp recombination site (dif) located at the junction of the two replication arms. Xer recombination is tightly controlled by the septal protein FtsK. XerCD recombinases and FtsK are found on most sequenced eubacterial genomes, suggesting that the Xer recombination system as described in E. coli is highly conserved among prokaryotes. We show here that Streptococci and Lactococci carry an alternative Xer recombination machinery, organized in a single recombination module. This corresponds to an atypical 31-bp recombination site (dif(SL)) associated with a dedicated tyrosine recombinase (XerS). In contrast to the E. coli Xer system, only a single recombinase is required to recombine dif(SL), suggesting a different mechanism in the recombination process. Despite this important difference, XerS can only perform efficient recombination when dif(SL) sites are located on chromosome dimers. Moreover, the XerS/dif(SL) recombination requires the streptococcal protein FtsK(SL), probably without the need for direct protein-protein interaction, which we demonstrated to be located at the division septum of Lactococcus lactis. Acquisition of the XerS recombination module can be considered as a landmark of the separation of Streptococci/Lactococci from other firmicutes and support the view that Xer recombination is a conserved cellular function in bacteria, but that can be achieved by functional analogs.

New criteria for selecting the origin of DNA replication in Wolbachia and closely related bacteria

Background

The annotated genomes of two closely related strains of the intracellular bacterium Wolbachia pipientis have been reported without the identifications of the putative origin of replication (ori). Identifying the ori of these bacteria and related alpha-Proteobacteria as well as their patterns of sequence evolution will aid studies of cell replication and cell density, as well as the potential genetic manipulation of these widespread intracellular bacteria.

Results

Using features that have been previously experimentally verified in the alpha-Proteobacterium Caulobacter crescentus, the origin of DNA replication (ori) regions were identified in silico for Wolbachia strains and eleven other related bacteria belonging to Ehrlichia, Anaplasma, and Rickettsia genera. These features include DnaA-, CtrA- and IHF-binding sites as well as the flanking genes in C. crescentus. The Wolbachia ori boundary genes were found to be hemE and COG1253 protein (CBS domain protein). Comparisons of the putative ori region among related Wolbachia strains showed higher conservation of bases within binding sites.

Conclusion

The sequences of the ori regions described here are only similar among closely related bacteria while fundamental characteristics like presence of DnaA and IHF binding sites as well as the boundary genes are more widely conserved. The relative paucity of CtrA binding sites in the ori regions, as well as the absence of key enzymes associated with DNA replication in the respective genomes, suggest that several of these obligate intracellular bacteria may have altered replication mechanisms. Based on these analyses, criteria are set forth for identifying the ori region in genome sequencing projects.

Being pathogenic, plastic, and sexual while living with a nearly minimal bacterial genome

Mycoplasmas are commonly described as the simplest self-replicating organisms, whose evolution was mainly characterized by genome downsizing with a proposed evolutionary scenario similar to that of obligate intracellular bacteria such as insect endosymbionts. Thus far, analysis of mycoplasma genomes indicates a low level of horizontal gene transfer (HGT) implying that DNA acquisition is strongly limited in these minimal bacteria. In this study, the genome of the ruminant pathogen Mycoplasma agalactiae was sequenced. Comparative genomic data and phylogenetic tree reconstruction revealed that ∼18% of its small genome (877,438 bp) has undergone HGT with the phylogenetically distinct mycoides cluster, which is composed of significant ruminant pathogens. HGT involves genes often found as clusters, several of which encode lipoproteins that usually play an important role in mycoplasma–host interaction. A decayed form of a conjugative element also described in a member of the mycoides cluster was found in the M. agalactiae genome, suggesting that HGT may have occurred by mobilizing a related genetic element. The possibility of HGT events among other mycoplasmas was evaluated with the available sequenced genomes. Our data indicate marginal levels of HGT among Mycoplasma species except for those described above and, to a lesser extent, for those observed in between the two bird pathogens, M. gallisepticum and M. synoviae. This first description of large-scale HGT among mycoplasmas sharing the same ecological niche challenges the generally accepted evolutionary scenario in which gene loss is the main driving force of mycoplasma evolution. The latter clearly differs from that of other bacteria with small genomes, particularly obligate intracellular bacteria that are isolated within host cells. Consequently, mycoplasmas are not only able to subvert complex hosts but presumably have retained sexual competence, a trait that may prevent them from genome stasis and contribute to adaptation to new hosts.

Mycoplasmas are cell wall–lacking prokaryotes that evolved from ancestors common to Gram-positive bacteria by way of massive losses of genetic material. With their minimal genome, mycoplasmas are considered to be the simplest free-living organisms, yet several species are successful pathogens of man and animal. In this study, we challenged the commonly accepted view in which mycoplasma evolution is driven only by genome down-sizing. Indeed, we showed that a significant amount of genes underwent horizontal transfer among different mycoplasma species that share the same ruminant hosts. In these species, the occurrence of a genetic element that can promote DNA transfer via cell-to-cell contact suggests that some mycoplasmas may have retained or acquired sexual competence. Transferred genes were found to encode proteins that are likely to be associated with mycoplasma–host interactions. Sharing genetic resources via horizontal gene transfer may provide mycoplasmas with a means for adapting to new niches or to new hosts and for avoiding irreversible genome erosion.

Conservation of the links between gene transcription and chromosomal organization in the highly reduced genome of Buchnera aphidicola

BACKGROUND

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.

RESULTS

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.

CONCLUSION

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.

Solving ambiguities in contig assembly of Idiomarina loihiensis L2TR chromosome by in silico analyses

Nucleotide composition analyses of bacterial genomes such as cumulative GC skew highlight the atypical, strongly asymmetric architecture of the recently published chromosome of Idiomarina loihiensis L2TR, suggesting that an inversion of a 600-kb chromosomal segment occurred. The presence of 3.4-kb inverted repeated sequences at the borders of the putative rearrangement supports this hypothesis. Reverting in silico this segment restores (1) a symmetric chromosome architecture; (2) the co-orientation of transcription of all rRNA operons with DNA replication; and (3) a better conservation of gene order between this chromosome and other gamma-proteobacterial ones. Finally, long-range PCRs encompassing the ends of the 600-kb segment reveal the existence of the reverted configuration but not of the published one. This demonstrates how cumulative nucleotide-skew analyses can validate genome assemblies.

Genome-wide coorientation of replication and transcription reduces adverse effects on replication in Bacillus subtilis

In many bacteria, there is a strong bias for genes to be encoded on the leading strand of DNA, resulting in coorientation of replication and transcription. In Bacillus subtilis, transcription of the majority of genes (75%) is cooriented with replication. By using genome-wide profiling of replication with DNA microarrays, we found that this coorientation bias reduces adverse effects of transcription on replication. We found that in wild-type cells, transcription did not appear to affect the rate of replication elongation. However, in mutants with reversed transcription bias for an extended region of the chromosome, replication elongation was slower. This reduced replication rate depended on transcription and was limited to the region in which the directions of replication and transcription are opposed. These results support the hypothesis that the strong bias to coorient transcription and replication is due to selective pressure for processive, efficient, and accurate replication.

Mutational bias suggests that replication termination occurs near the dif site, not at Ter sites

In bacteria, Ter sites bound to Tus/Rtp proteins halt replication forks moving only in one direction, providing a convenient mechanism to terminate them once the chromosome had been replicated. Considering the importance of replication termination and its position as a checkpoint in cell division, the accumulated knowledge on these systems has not dispelled fundamental questions regarding its role in cell biology: why are there so many copies of Ter, why are they distributed over such a large portion of the chromosome, why is the tus gene not conserved among bacteria, and why do tus mutants lack measurable phenotypes? Here we examine bacterial genomes using bioinformatics techniques to identify the region(s) where DNA polymerase III-mediated replication has historically been terminated. We find that in both Escherichia coli and Bacillus subtilis, changes in mutational bias patterns indicate that replication termination most likely occurs at or near the dif site. More importantly, there is no evidence from mutational bias signatures that replication forks originating at oriC have terminated at Ter sites. We propose that Ter sites participate in halting replication forks originating from DNA repair events, and not those originating at the chromosomal origin of replication.

Validating the significance of genomic properties of Chi sites from the distribution of all octamers in Escherichia coli

Chi sites (5'-GCTGGTGG-3') are homologous recombinational hotspot octamer sequences, which attenuate the exonuclease activity of RecBCD in Escherichia coli. They are overrepresented in the genome (1008 occurrences), preferentially located within coding regions (98%), oriented in the direction of replication (75%), and occur most commonly on the mRNA-synonymous sense strand of the double helix (79%). Previous statistical studies of the genome sequence suggested that these genomic properties of Chi sites appear to be related to their role in recombinational repair and therefore to replication and transcription. In this study, we employ three mathematical models to predict the properties of Chi sites from single nucleotide and multi-nucleotide compositions, and validate them statistically using the distribution of all octamer sequences in the entire genome, or exclusively within ORFs. The model based on the overall distribution of all octamers provided better predictions than the single nucleotide composition model, and the ORF and sense strand preference of Chi sites were shown to be within the standard deviation of all octamers. In contrast, the orientation bias of the Chi sites in the direction of replication was significant, although the bias was not as pronounced as with the single nucleotide composition model, suggesting a selective pressure related to the role of RecBCD in replication.

Noise-reduction filtering for accurate detection of replication termini in bacterial genomes

Bacterial chromosomes are highly polarized in their nucleotide composition through mutational selection related to replication. Using compositional skews such as the GC skew, replication origin and terminus can be predicted in silico by observing the shift points. However, the genome sequence is affected by myriad functional requirements and selection on numerous subgenomic features, and elimination of this "noise" should lead to better predictions. Here, we present a noise-reduction approach that uses low-pass filtering through Fast Fourier transform coupled with cumulative skew graphs. It increases the prediction accuracy of the replication termini compared with previously documented methods based on genomic base composition.

Similar compositional biases are caused by very different mutational effects

Compositional replication strand bias, commonly referred to as GC skew, is present in many genomes of prokaryotes, eukaryotes, and viruses. Although cytosine deamination in ssDNA (resulting in C-->T changes on the leading strand) is often invoked as its major cause, the precise contributions of this and other substitution types are currently unknown. It is also unclear if the underlying mutational asymmetries are the same among taxa, are stable over time, or how closely the observed biases are to mutational equilibrium. We analyzed nearly neutral sites of seven taxa each with between three and six complete bacterial genomes, and inferred the substitution spectra of fourfold degenerate positions in nonhighly expressed genes. Using a bootstrap procedure, we extracted compositional biases associated with replication and identified the significant asymmetries. Although all taxa showed an overrepresentation of G relative to C on the leading strand (and imbalances between A and T), widely variable substitution asymmetries are noted. Surprisingly, all substitution types show significant asymmetry in at least one taxon, but none were universally biased in all taxa. Notably, in the two most biased genomes, A-->G, rather than C-->T, shapes the compositional bias. Given the variability in these biases, we propose that the process is multifactorial. Finally, we also find that most genomes are not at compositional equilibrium, and suggest that mutational-based heterotachy is deeply imprinted in the history of biological macromolecules. This shows that similar compositional biases associated with the same essential well-conserved process, replication, do not reflect similar mutational processes in different genomes, and that caution is required in inferring the roles of specific mutational biases on the basis of contemporary patterns of sequence composition.

Analysis of the terminus region of the Caulobacter crescentus chromosome and identification of the dif site

The terminus region of the Caulobacter crescentus chromosome and the dif chromosome dimer resolution site were characterized. The Caulobacter genome contains skewed sequences that abruptly switch strands at dif and may have roles in chromosome maintenance and segregation. Absence of dif or the XerCD recombinase results in a chromosome segregation defect. The Caulobacter terminus region is unusual, since it contains many essential or highly expressed genes.

Evolution of the Terminal Regions of the Streptomyces Linear Chromosome

Comparative analysis of the Streptomyces chromosome sequences, between Streptomyces coelicolor, Streptomyces avermitilis, and Streptomyces ambofaciens ATCC23877 (whose partial sequence is released in this study), revealed a highly compartmentalized genetic organization of their genome. Indeed, despite the presence of specific genomic islands, the central part of the chromosome appears highly syntenic. In contrast, the chromosome of each species exhibits large species-specific terminal regions (from 753 to 1,393 kb), even when considering closely related species (S. ambofaciens and S. coelicolor). Interestingly, the size of the central conserved region between species decreases as the phylogenetic distance between them increases, whereas the specific terminal fraction reciprocally increases in size. Between highly syntenic central regions and species-specific chromosomal parts, there is a notable degeneration of synteny due to frequent insertions/deletions. This reveals a massive and constant genomic flux (from lateral gene transfer and DNA rearrangements) affecting the terminal contingency regions. We speculate that a gradient of recombination rate (i.e., insertion/deletion events) toward the extremities is the force driving the exclusion of essential genes from the terminal regions (i.e., chromosome compartmentalization) and generating a fast gene turnover for strong adaptation capabilities.

An integrative method for accurate comparative genome mapping

We present MAGIC, an integrative and accurate method for comparative genome mapping. Our method consists of two phases: preprocessing for identifying "maximal similar segments," and mapping for clustering and classifying these segments. MAGIC's main novelty lies in its biologically intuitive clustering approach, which aims towards both calculating reorder-free segments and identifying orthologous segments. In the process, MAGIC efficiently handles ambiguities resulting from duplications that occurred before the speciation of the considered organisms from their most recent common ancestor. We demonstrate both MAGIC's robustness and scalability: the former is asserted with respect to its initial input and with respect to its parameters' values. The latter is asserted by applying MAGIC to distantly related organisms and to large genomes. We compare MAGIC to other comparative mapping methods and provide detailed analysis of the differences between them. Our improvements allow a comprehensive study of the diversity of genetic repertoires resulting from large-scale mutations, such as indels and duplications, including explicitly transposable and phagic elements. The strength of our method is demonstrated by detailed statistics computed for each type of these large-scale mutations. MAGIC enabled us to conduct a comprehensive analysis of the different forces shaping prokaryotic genomes from different clades, and to quantify the importance of novel gene content introduced by horizontal gene transfer relative to gene duplication in bacterial genome evolution. We use these results to investigate the breakpoint distribution in several prokaryotic genomes.

Staphylococcus aureus helicase but not Escherichia coli helicase stimulates S. aureus primase activity and maintains initiation specificity

Bacterial primases are essential for DNA replication due to their role in polymerizing the formation of short RNA primers repeatedly on the lagging-strand template and at least once on the leading-strand template. The ability of recombinant Staphylococcus aureus DnaG primase to utilize different single-stranded DNA templates was tested using oligonucleotides of the sequence 5'-CAGA (CA)5 XYZ (CA)3-3', where XYZ represented the variable trinucleotide. These experiments demonstrated that S. aureus primase synthesized RNA primers predominately on templates containing 5'-d(CTA)-3' or TTA and to a much lesser degree on GTA-containing templates, in contrast to results seen with the Escherichia coli DnaG primase recognition sequence 5'-d(CTG)-3'. Primer synthesis was initiated complementarily to the middle nucleotide of the recognition sequence, while the third nucleotide, an adenosine, was required to support primer synthesis but was not copied into the RNA primer. The replicative helicases from both S. aureus and E. coli were tested for their ability to stimulate either S. aureus or E. coli primase. Results showed that each bacterial helicase could only stimulate the cognate bacterial primase. In addition, S. aureus helicase stimulated the production of full-length primers, whereas E. coli helicase increased the synthesis of only short RNA polymers. These studies identified important differences between E. coli and S. aureus related to DNA replication and suggest that each bacterial primase and helicase may have adapted unique properties optimized for replication.

Tracking of controlled Escherichia coli replication fork stalling and restart at repressor-bound DNA in vivo

We report an efficient, controllable, site-specific replication roadblock that blocks cell proliferation, but which can be rapidly and efficiently reversed, leading to recovery of viability. Escherichia coli replication forks of both polarities stalled in vivo within the first 500 bp of a 10 kb repressor-bound array of operator DNA-binding sites. Controlled release of repressor binding led to rapid restart of the blocked replication fork without the participation of homologous recombination. Cytological tracking of fork stalling and restart showed that the replisome-associated SSB protein remains associated with the blocked fork for extended periods and that duplication of the fluorescent foci associated with the blocked operator array occurs immediately after restart, thereby demonstrating a lack of sister cohesion in the region of the array. Roadblocks positioned near oriC or the dif site did not prevent replication and segregation of the rest of the chromosome.

A novel FtsZ-like protein is involved in replication of the anthrax toxin-encoding pXO1 plasmid in Bacillus anthracis

Plasmid pXO1 encodes the tripartite anthrax toxin, which is the major virulence factor of Bacillus anthracis. In spite of the important role of pXO1 in anthrax pathogenesis, very little is known about its replication and maintenance in B. anthracis. We cloned a 5-kb region of the pXO1 plasmid into an Escherichia coli vector and showed that this plasmid can replicate when introduced into B. anthracis. Mutational analysis showed that open reading frame 45 (repX) of pXO1 was required for the replication of the miniplasmid in B. anthracis. Interestingly, repX showed limited homology to bacterial FtsZ proteins that are involved in cell division. A mutation in the predicted GTP binding domain of RepX abolished its replication activity. Genes almost identical to repX are contained on several megaplasmids in members of the Bacillus cereus group, including a B. cereus strain that causes an anthrax-like disease. Our results identify a novel group of FtsZ-related initiator proteins that are required for the replication of virulence plasmids in B. anthracis and possibly in related organisms. Such replication proteins may provide novel drug targets for the elimination of plasmids encoding the anthrax toxin and other virulence factors.

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.

Selection for Chromosome Architecture in Bacteria

Bacterial chromosomes are immense polymers whose faithful replication and segregation are crucial to cell survival. The ability of proteins such as FtsK to move unidirectionally toward the replication terminus, and direct DNA translocation into the appropriate daughter cell during cell division, requires that bacterial genomes maintain an architecture for the orderly replication and segregation of chromosomes. We suggest that proteins that locate the replication terminus exploit strand-biased sequences that are overrepresented on one DNA strand, and that selection increases with decreased distance to the replication terminus. We report a generalized method for detecting these architecture imparting sequences (AIMS) and have identified AIMS in nearly all bacterial genomes. Their increased abundance on leading strands and decreased abundance on lagging strands toward replication termini are not the result of changes in mutational bias; rather, they reflect a gradient of long-term positive selection for AIMS. The maintenance of the pattern of AIMS across the genomes of related bacteria independent of their positions within individual genes suggests a well-conserved role in genome biology. The stable gradient of AIMS abundance from replication origin to terminus suggests that the replicore acts as a target of selection, where selection for chromosome architecture results in the maintenance of gene order and in the lack of high-frequency DNA inversion within replicores.

Strong asymmetric mutation bias in endosymbiont genomes coincide with loss of genes for replication restart pathways

A large majority of bacterial genomes show strand asymmetry, such that G and T preferentially accumulate on the leading strand. The mechanisms are unknown, but cytosine deaminations are thought to play an important role. Here, we have examined DNA strand asymmetry in three strains of the aphid endosymbiont Buchnera aphidicola. These are phylogenetically related, have similar genomic GC contents, and conserved gene order structures, yet B. aphidicola (Bp) shows a fourfold higher replication-induced strand bias than B. aphidicola (Sg) and (Ap). We rule out an increase in the overall substitution frequency as the major cause of the stronger strand bias in B. aphidicola (Bp). Instead, the results suggest that the higher GC skew in this species is caused by a different spectrum of mutations, including a relatively higher frequency of C to T mutations on the leading strand and/or of G to A mutations on the lagging strand. A comparative analysis of 20 gamma-proteobacterial genomes revealed that endosymbiont genomes lacking recA and other genes involved in replication restart processes, such as priA, which codes for primosomal helicase PriA, displayed the strongest strand bias. We hypothesize that cytosine deaminations accumulate during single-strand exposure at arrested replication forks and that inefficient restart mechanisms may lead to high DNA strand asymmetry in bacterial genomes.

Origin of replication in circular prokaryotic chromosomes

To predict origins of replication in prokaryotic chromosomes, we analyse the leading and lagging strands of 200 chromosomes for differences in oligomer composition and show that these correlate strongly with taxonomic grouping, lifestyle and molecular details of the replication process. While all bacteria have a preference for Gs over Cs on the leading strand, we discover that the direction of the A/T skew is determined by the polymerase-alpha subunit that replicates the leading strand. The strength of the strand bias varies greatly between both phyla and environments and appears to correlate with growth rate. Finally we observe much greater diversity of skew among archaea than among bacteria. We have developed a program that accurately locates the origins of replication by measuring the differences between leading and lagging strand of all oligonucleotides up to 8 bp in length. The program and results for all publicly available genomes are available from http://www.cbs.dtu.dk/services/GenomeAtlas/suppl/origin.

Genome evolution in bacteria: order beneath chaos

Bacterial genomes have been viewed as collections of genes, with each gene and genome evolving more-or-less independently through the acquisition of mutational changes. This historical view has been overturned by the finding that genomes of even closely-related taxa differ widely in gene content. Yet, genomes are more than ever-shuffling collections of genes. Some genes within a genome are more transient than others, conferring a layer of phenotypic lability over a core of genotypic stability; this core decreases in size as the taxa included become increasingly diverse. In addition, some lineages no longer experience high rates of gene turnover, and gene content alters primarily through slow rates of gene loss. More importantly, the cell and molecular biology of the bacterial cell imposes constraints on chromosome composition, maintaining a stable architecture in the face of gene turnover. As a result, genomes reflect the sum of processes that introduce variability, which is then arbitrated by processes that maintain stability.

Inference and analysis of the relative stability of bacterial chromosomes

The stability of genomes is highly variable, both in terms of gene content and gene order. Here I calibrate the loss of gene order conservation (GOC) through time by fitting a simple probabilistic model on pairwise comparisons involving 126 bacterial genomes. The model computes the probability of separation of pairs of contiguous genes per unit of time and fits the data better than previous ones while allowing a mechanistic interpretation for the loss of GOC with time. Although the information on operons is not used in the model, I observe, as expected, that most highly conserved pairs of genes are indeed within operons. However, even the other pairs are much more conserved than expected given the observed experimental rearrangement rates. After 500 Myr, about 50% of the originally contiguous orthologues remain so in the average genome. Hence, the large majority of rearrangements must be deleterious and random genome rearrangements are unlikely to provide for positively selected structural changes. I then use the deviations from the model to define an intrinsic measure of genome stability that allowed the comparison of distantly related genomes and the inference of ancestral states. This shows that clades differ in genome stability, with cyanobacteria being the least stable and gamma-proteobacteria the most stable. Without correction for phylogeny, free-living bacteria are the least stable group of genomes, followed by pathogens, and then endomutualists. However, after correction for phylogenetic inertia (or the removal of cyanobacteria from the analysis), there is no significant association between genome stability and lifestyle or genome size. Hence, although this method has allowed uncovering some of mechanisms leading to rearrangements, we still ignore the forces that differentially shape selection upon genome stability in different species.

Combining two genomes in one cell: stable cloning of the Synechocystis PCC6803 genome in the Bacillus subtilis 168 genome

Cloning the whole 3.5-megabase (Mb) genome of the photosynthetic bacterium Synechocystis PCC6803 into the 4.2-Mb genome of the mesophilic bacterium Bacillus subtilis 168 resulted in a 7.7-Mb composite genome. We succeeded in such unprecedented large-size cloning by progressively assembling and editing contiguous DNA regions that cover the entire Synechocystis genome. The strain containing the two sets of genome grew only in the B. subtilis culture medium where all of the cloning procedures were carried out. The high structural stability of the cloned Synechocystis genome was closely associated with the symmetry of the bacterial genome structure of the DNA replication origin (oriC) and its termination (terC) and the exclusivity of Synechocystis ribosomal RNA operon genes (rrnA and rrnB). Given the significant diversity in genome structure observed upon horizontal DNA transfer in nature, our stable laboratory-generated composite genome raised fundamental questions concerning two complete genomes in one cell. Our megasize DNA cloning method, designated megacloning, may be generally applicable to other genomes or genome loci of free-living organisms.

The genomic pattern of tDNA operon expression in E. coli

In fast-growing microorganisms, a tRNA concentration profile enriched in major isoacceptors selects for the biased usage of cognate codons. This optimizes translational rate for the least mass invested in the translational apparatus. Such translational streamlining is thought to be growth-regulated, but its genetic basis is poorly understood. First, we found in reanalysis of the E. coli tRNA profile that the degree to which it is translationally streamlined is nearly invariant with growth rate. Then, using least squares multiple regression, we partitioned tRNA isoacceptor pools to predicted tDNA operons from the E. coli K12 genome. Co-expression of tDNAs in operons explains the tRNA profile significantly better than tDNA gene dosage alone. Also, operon expression increases significantly with proximity to the origin of replication, oriC, at all growth rates. Genome location explains about 15% of expression variation in a form, at a given growth rate, that is consistent with replication-dependent gene concentration effects. Yet the change in the tRNA profile with growth rate is less than would be expected from such effects. We estimated per-copy expression rates for all tDNA operons that were consistent with independent estimates for rDNA operons. We also found that tDNA operon location, and the location dependence of expression, were significantly different in the leading and lagging strands. The operonic organization and genomic location of tDNA operons are significant factors influencing their expression. Nonrandom patterns of location and strandedness shown by tDNA operons in E. coli suggest that their genomic architecture may be under selection to satisfy physiological demand for tRNA expression at high growth rates.

The concentrations of tRNAs are co-adapted to codon usage frequencies in the transcriptomes of E. coli and other diverse organisms. But how are tRNA concentrations determined? Here, the researchers analyzed the E. coli tRNA concentration profile in its genomic context, using clustering and regression methods to partition tRNA concentration data to tDNA operons that were defined semi-automatically. They found that co-expression in operons explains the tRNA profile much better than tDNA gene dosage alone. Furthermore, they could significantly explain the total expression from tDNA operons by their distance from the genomic origin of replication. Per-copy transcription initiation rates from tDNA operons were also estimated. Although there is some evidence for replication-dependent effects on tDNA operon expression, this cannot explain how constant the tRNA profile is with growth rate. As a consequence, tDNA promoters are predicted to compensate for the location of their operons. Finally, the researchers found pronounced asymmetries between the leading and lagging genomic strands in the locations of tDNA operons, and on the effect of location on their expression. These nonrandom patterns suggest that the genomic location and strandedness of tDNA operons may be under some selection in E. coli to satisfy physiological demand for tRNAs at high growth rates.

Interruptions in gene expression drive highly expressed operons to the leading strand of DNA replication

In bacteria, most genes are on the leading strand of replication, a phenomenon attributed to collisions between the DNA and RNA polymerases. In Escherichia coli, these collisions slow the movement of the replication fork through actively transcribed genes only if they are coded on the lagging strand. For genes on both strands, however, these collisions sever nascent transcripts and interrupt gene expression. Based on these observations, we propose a new theory to explain strand bias: genes whose expression is important for fitness are selected to the leading strand because this reduces the duration of these interruptions. Our theory predicts that multi-gene operons, which are subject to longer interruptions, should be more strongly selected to the leading strand than singleton transcripts. We show that this is true even after controlling for the tendency for essential genes, which are strongly biased to the leading strand, to occur in operons. Our theory also predicts that other factors that are associated with strand bias should have stronger effects for genes that are in operons. We find that expression level and phylogenetic ubiquity are correlated with strand bias for both essential and non-essential genes, but only for genes in operons.

Real-time RT-PCR normalisation; strategies and considerations

Real-time RT-PCR has become a common technique, no longer limited to specialist core facilities. It is in many cases the only method for measuring mRNA levels of vivo low copy number targets of interest for which alternative assays either do not exist or lack the required sensitivity. Benefits of this procedure over conventional methods for measuring RNA include its sensitivity, large dynamic range, the potential for high throughout as well as accurate quantification. To achieve this, however, appropriate normalisation strategies are required to control for experimental error introduced during the multistage process required to extract and process the RNA. There are many strategies that can be chosen; these include normalisation to sample size, total RNA and the popular practice of measuring an internal reference or housekeeping gene. However, these methods are frequently applied without appropriate validation. In this review we discuss the relative merits of different normalisation strategies and suggest a method of validation that will enable the measurement of biologically meaningful results.

A bunch of fun-guys: the whole-genome view of yeast evolution

The resolution of the complete sequences of several hemiascomycete genomes provides new insights into the ways that yeast genomes change in size and in gene contents. These genomes provide evidence of whole-genome duplication occurring before the divergence of Saccharomyces cerevisiae and Candida glabrata, followed by massive gene loss that restored diploidy. The pattern of genome evolution in yeast differs from that in bacteria apparently as a result of stronger selective constraints on bacterial chromosomes.

Codon usage bias from tRNA's point of view: redundancy, specialization, and efficient decoding for translation optimization

The selection-mutation-drift theory of codon usage plays a major role in the theory of molecular evolution by explaining the co-evolution of codon usage bias and tRNA content in the framework of translation optimization. Because most studies have focused only on codon usage, we analyzed the tRNA gene pool of 102 bacterial species. We show that as minimal generation times get shorter, the genomes contain more tRNA genes, but fewer anticodon species. Surprisingly, despite the wide G+C variation of bacterial genomes these anticodons are the same in most genomes. This suggests an optimization of the translation machinery to use a small subset of optimal codons and anticodons in fast-growing bacteria and in highly expressed genes. As a result, the overrepresented codons in highly expressed genes tend to be the same in very different genomes to match the same most-frequent anticodons. This is particularly important in fast-growing bacteria, which have higher codon usage bias in these genes. Three models were tested to understand the choice of codons recognized by the same anticodons, all providing significant fit, but under different classes of genes and genomes. Thus, co-evolution of tRNA gene composition and codon usage bias in genomes seen from tRNA's point of view agrees with the selection-mutation-drift theory. However, it suggests a much more universal trend in the evolution of anticodon and codon choice than previously thought. It also provides new evidence that a selective force for the optimization of the translation machinery is the maximization of growth.