Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication

RiboGrove: a database of full-length prokaryotic 16S rRNA genes derived from completely assembled genomes

16S rRNA gene is frequently used for the identification of prokaryotic organisms and for phylogeny inference. Several specialized public databases exist that contain complete and partial sequences of 16S rRNA genes. In this paper, we present RiboGrove: the first publicly available database that comprises only full-length sequences of 16S rRNA genes originating from completely assembled prokaryotic genomes deposited in RefSeq. Despite being strongly biased towards frequently sequenced genomes, RiboGrove is a useful complement to existing 16S rRNA resources and allows for analyses that would not be possible using amplicon-derived gene sequences. For instance, the absence of partial gene sequences in RiboGrove allowed us to make a summary of prokaryotic organisms, which lack core anti-Shine-Dalgarno sequence in their 16S rRNA genes. In this study, we describe the collected sequence data and present the results of exploratory data analysis of 16S rRNA gene sequences.

Transcription-Replication Collisions and Chromosome Fragility

Accurate replication of the entire genome is critical for cell division and propagation. Certain regions in the genome, such as fragile sites (common fragile sites, rare fragile sites, early replicating fragile sites), rDNA and telomeres, are intrinsically difficult to replicate, especially in the presence of replication stress caused by, for example, oncogene activation during tumor development. Therefore, these regions are particularly prone to deletions and chromosome rearrangements during tumorigenesis, rendering chromosome fragility. Although, the mechanism underlying their “difficult-to-replicate” nature and genomic instability is still not fully understood, accumulating evidence suggests transcription might be a major source of endogenous replication stress (RS) leading to chromosome fragility. Here, we provide an updated overview of how transcription affects chromosome fragility. Furthermore, we will use the well characterized common fragile sites (CFSs) as a model to discuss pathways involved in offsetting transcription-induced RS at these loci with a focus on the recently discovered atypical DNA synthesis repair pathway Mitotic DNA Synthesis (MiDAS).

Recent Advancements in DNA Damage-Transcription Crosstalk and High-Resolution Mapping of DNA Breaks

Until recently, DNA damage arising from physiological DNA metabolism was considered a detrimental by-product for cells. However, an increasing amount of evidence has shown that DNA damage could have a positive role in transcription activation. In particular, DNA damage has been detected in transcriptional elements following different stimuli. These physiological DNA breaks are thought to be instrumental for the correct expression of genomic loci through different mechanisms. In this regard, although a plethora of methods are available to precisely map transcribed regions and transcription start sites, commonly used techniques for mapping DNA breaks lack sufficient resolution and sensitivity to draw a robust correlation between DNA damage generation and transcription. Recently, however, several methods have been developed to map DNA damage at single-nucleotide resolution, thus providing a new set of tools to correlate DNA damage and transcription. Here, we review how DNA damage can positively regulate transcription initiation, the current techniques for mapping DNA breaks at high resolution, and how these techniques can benefit future studies of DNA damage and transcription.

Conflict Resolution in the Genome: How Transcription and Replication Make It Work

The complex machineries involved in replication and transcription translocate along the same DNA template, often in opposing directions and at different rates. These processes routinely interfere with each other in prokaryotes, and mounting evidence now suggests that RNA polymerase complexes also encounter replication forks in higher eukaryotes. Indeed, cells rely on numerous mechanisms to avoid, tolerate, and resolve such transcription-replication conflicts, and the absence of these mechanisms can lead to catastrophic effects on genome stability and cell viability. In this article, we review the cellular responses to transcription-replication conflicts and highlight how these inevitable encounters shape the genome and impact diverse cellular processes.

Regulation of DNA Replication through Natural Impediments in the Eukaryotic Genome

All living organisms need to duplicate their genetic information while protecting it from unwanted mutations, which can lead to genetic disorders and cancer development. Inaccuracies during DNA replication are the major cause of genomic instability, as replication forks are prone to stalling and collapse, resulting in DNA damage. The presence of exogenous DNA damaging agents as well as endogenous difficult-to-replicate DNA regions containing DNA–protein complexes, repetitive DNA, secondary DNA structures, or transcribing RNA polymerases, increases the risk of genomic instability and thus threatens cell survival. Therefore, understanding the cellular mechanisms required to preserve the genetic information during S phase is of paramount importance. In this review, we will discuss our current understanding of how cells cope with these natural impediments in order to prevent DNA damage and genomic instability during DNA replication.

A few sequence polymorphisms among isolates of Maize bushy stunt phytoplasma associate with organ proliferation symptoms of infected maize plants

Background and Aims

Maize bushy stunt phytoplasma (MBSP) is a bacterial pathogen of maize ( Zea mays L.) across Latin America. MBSP belongs to the 16SrI-B sub-group within the genus ' Candidatus Phytoplasma'. MBSP and its insect vector Dalbulus maidis (Hemiptera: Cicadellidae) are restricted to maize; both are thought to have coevolved with maize during its domestication from a teosinte-like ancestor. MBSP-infected maize plants show a diversity of symptoms. and it is likely that MBSP is under strong selection for increased virulence and insect transmission on maize hybrids that are widely grown in Brazil. In this study it was investigated whether the differences in genome sequences of MBSP isolates from two maize-growing regions in South-east Brazil explain variations in symptom severity of the MBSP isolates on various maize genotypes.

Methods

MBSP isolates were collected from maize production fields in Guaíra and Piracicaba in South-east Brazil for infection assays. One representative isolate was chosen for de novo whole-genome assembly and for the alignment of sequence reads from the genomes of other phytoplasma isolates to detect polymorphisms. Statistical methods were applied to investigate the correlation between variations in disease symptoms of infected maize plants and MBSP sequence polymorphisms.

Key Results

MBSP isolates contributed consistently to organ proliferation symptoms and maize genotype to leaf necrosis, reddening and yellowing of infected maize plants. The symptom differences are associated with polymorphisms in a phase-variable lipoprotein, which is a candidate effector, and an ATP-dependent lipoprotein ABC export protein, whereas no polymorphisms were observed in other candidate effector genes. Lipoproteins and ABC export proteins activate host defence responses, regulate pathogen attachment to host cells and activate effector secretion systems in other pathogens.

Conclusions

Polymorphisms in two putative virulence genes among MBSP isolates from maize-growing regions in South-east Brazil are associated with variations in organ proliferation symptoms of MBSP-infected maize plants.

Coordinating Replication with Transcription

DNA topological transitions occur when replication forks encounter other DNA transactions such as transcription. Failure in resolving such conflicts leads to generation of aberrant replication and transcription intermediates that might have adverse effects on genome stability. Cells have evolved numerous surveillance mechanisms to avoid, tolerate, and resolve such replication-transcription conflicts. Defects or non-coordination in such cellular mechanisms might have catastrophic effect on cell viability. In this chapter, we review consequences of replication encounters with transcription and its associated events, topological challenges, and how these inevitable conflicts alter the genome structure and functions.

Function of the CRISPR-Cas System of the Human Pathogen Clostridium difficile

Clostridium difficile is the cause of most frequently occurring nosocomial diarrhea worldwide. As an enteropathogen, C. difficile must be exposed to multiple exogenous genetic elements in bacteriophage-rich gut communities. CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) systems allow bacteria to adapt to foreign genetic invaders. Our recent data revealed active expression and processing of CRISPR RNAs from multiple type I-B CRISPR arrays in C. difficile reference strain 630. Here, we demonstrate active expression of CRISPR arrays in strain R20291, an epidemic C. difficile strain. Through genome sequencing and host range analysis of several new C. difficile phages and plasmid conjugation experiments, we provide evidence of defensive function of the CRISPR-Cas system in both C. difficile strains. We further demonstrate that C. difficile Cas proteins are capable of interference in a heterologous host, Escherichia coli. These data set the stage for mechanistic and physiological analyses of CRISPR-Cas-mediated interactions of important global human pathogen with its genetic parasites.

Clostridium difficile is the major cause of nosocomial infections associated with antibiotic therapy worldwide. To survive in bacteriophage-rich gut communities, enteropathogens must develop efficient systems for defense against foreign DNA elements. CRISPR-Cas systems have recently taken center stage among various anti-invader bacterial defense systems. We provide experimental evidence for the function of the C. difficile CRISPR system against plasmid DNA and bacteriophages. These data demonstrate the original features of active C. difficile CRISPR system and bring important insights into the interactions of this major enteropathogen with foreign DNA invaders during its infection cycle.

The B. subtilis Accessory Helicase PcrA Facilitates DNA Replication through Transcription Units

In bacteria the concurrence of DNA replication and transcription leads to potentially deleterious encounters between the two machineries, which can occur in either the head-on (lagging strand genes) or co-directional (leading strand genes) orientations. These conflicts lead to replication fork stalling and can destabilize the genome. Both eukaryotic and prokaryotic cells possess resolution factors that reduce the severity of these encounters. Though Escherichia coli accessory helicases have been implicated in the mitigation of head-on conflicts, direct evidence of these proteins mitigating co-directional conflicts is lacking. Furthermore, the endogenous chromosomal regions where these helicases act, and the mechanism of recruitment, have not been identified. We show that the essential Bacillus subtilis accessory helicase PcrA aids replication progression through protein coding genes of both head-on and co-directional orientations, as well as rRNA and tRNA genes. ChIP-Seq experiments show that co-directional conflicts at highly transcribed rRNA, tRNA, and head-on protein coding genes are major targets of PcrA activity on the chromosome. Partial depletion of PcrA renders cells extremely sensitive to head-on conflicts, linking the essential function of PcrA to conflict resolution. Furthermore, ablating PcrA’s ATPase/helicase activity simultaneously increases its association with conflict regions, while incapacitating its ability to mitigate conflicts, and leads to cell death. In contrast, disruption of PcrA’s C-terminal RNA polymerase interaction domain does not impact its ability to mitigate conflicts between replication and transcription, its association with conflict regions, or cell survival. Altogether, this work establishes PcrA as an essential factor involved in mitigating transcription-replication conflicts and identifies chromosomal regions where it routinely acts. As both conflicts and accessory helicases are found in all domains of life, these results are broadly relevant.

In bacteria the concurrence of DNA replication and transcription leads to potentially deleterious encounters between the two machineries. These encounters can destabilize the genome and lead to mutations. Both eukaryotic and prokaryotic cells possess conflict resolution factors that reduce the detrimental effects of these collisions. In this study we show that without the essential Bacillus subtilis accessory DNA helicase, PcrA, the replication machinery slows down at certain regions of the chromosome in a transcription-dependent manner. PcrA is essential to life but incomplete depletion of PcrA only partially inhibits cell survival. We find that, under these conditions, partial survival defects are significantly exacerbated in the presence of a single severe conflict. In summary our work identifies a high degree of conservation for accessory helicase function in conflict resolution, directly establishes PcrA’s role in co-directional conflict resolution, and maps the natural chromosomal regions where such activities are routinely needed. Because both conflicts and accessory helicases are found in all domains of life, the results of this work are broadly relevant.

Engineered ribosomal RNA operon copy-number variants of E. coli reveal the evolutionary trade-offs shaping rRNA operon number

Ribosomal RNA (rrn) operons, characteristically present in several copies in bacterial genomes (7 in E. coli), play a central role in cellular physiology. We investigated the factors determining the optimal number of rrn operons in E. coli by constructing isogenic variants with 5–10 operons. We found that the total RNA and protein content, as well as the size of the cells reflected the number of rrn operons. While growth parameters showed only minor differences, competition experiments revealed a clear pattern: 7–8 copies were optimal under conditions of fluctuating, occasionally rich nutrient influx and lower numbers were favored in stable, nutrient-limited environments. We found that the advantages of quick adjustment to nutrient availability, rapid growth and economic regulation of ribosome number all contribute to the selection of the optimal rrn operon number. Our results suggest that the wt rrn operon number of E. coli reflects the natural, ‘feast and famine’ life-style of the bacterium, however, different copy numbers might be beneficial under different environmental conditions. Understanding the impact of the copy number of rrn operons on the fitness of the cell is an important step towards the creation of functional and robust genomes, the ultimate goal of synthetic biology.

The GC content of LSU rRNA evolves across topological and functional regions of the ribosome in all three domains of life

Large-subunit rRNA is the ribozyme that catalyzes protein synthesis by translation, and many of its features vary along a deep-to-superficial gradient. By measuring the G+C proportions in this rRNA in all three domains of life (60 bacteria, 379 eukaryote, and 23 archaean sequences), we tested whether the proportion of GC nucleotides varies along this in-out gradient. The rRNA regions used were several zones identified by Bokov and Steinberg (2009) as being arranged from deep to superficial within the LSU. To the Bokov-Steinberg zones, we added the most superficial zone of all, the divergent domains (expansion segments), which are greatly enlarged in eukaryotes. Regression lines constructed from the hundreds of species of organisms revealed the expected in-out gradient, showing that species with high %GC (or high %AT) in their rRNA distribute more of these abundant nucleotides into the peripheral zones. This could be explained by the evolutionary rates of replacement of all nucleotides (A, C, G, T), because these latter rates are fastest at the periphery and slowest near the conserved core. As an overall explanation, we propose that when extrinsic factors (whole-genome nucleotide composition, or environmental temperature) demand the percentage of GC in the rRNA of a species be high or low, then the deep-lying zones are buffered against GC variation because they are the slowest to evolve. The deep, conserved zones are also the most involved in translation, hinting that stabilizing selection there prevents a high GC variability that would diminish LSU rRNA's core functions. We found only a few domain-specific trends in rRNA-GC distribution, which relate to many Archaea living at high temperatures or to the highly complex genes and adaptations of Eukaryota. Use of rRNA sequences in molecular phylogenetic studies, for reconstructing the relationships of organisms across the tree of life, requires accurate models of how rRNA evolves. The demonstration that GC distributes in regular patterns across rRNA regions can improve these tree-reconstruction models in the future and should yield phylogenies of greater accuracy.

Dealing with transcriptional outbursts during S phase to protect genomic integrity

Transcription during S phase needs to be spatially and temporally regulated to prevent collisions between the transcription and replication machineries. Cells have evolved a number of mechanisms to make both processes compatible under normal growth conditions. When conflict management fails, the head-on encounter between RNA and DNA polymerases results in genomic instability unless conflict resolution mechanisms are activated. Nevertheless, there are specific situations in which cells need to dramatically change their transcriptional landscape to adapt to environmental challenges. Signal transduction pathways, such as stress-activated protein kinases (SAPKs), serve to regulate gene expression in response to environmental insults. Prototypical members of SAPKs are the yeast Hog1 and mammalian p38. In response to stress, p38/Hog1 SAPKs control transcription and also regulate cell cycle progression. When yeast cells are stressed during S phase, Hog1 promotes gene induction and, remarkably, also delays replication by directly affecting early origin firing and fork progression. Therefore, by delaying replication, Hog1 plays a key role in preventing conflicts between RNA and DNA polymerases. In this review, we focus on the genomic determinants and mechanisms that make compatible transcription with replication during S phase to prevent genomic instability, especially in response to environmental changes.

Molecular traffic jams on DNA

All aspects of DNA metabolism-including transcription, replication, and repair-involve motor enzymes that move along genomic DNA. These processes must all take place on chromosomes that are occupied by a large number of other proteins. However, very little is known regarding how nucleic acid motor proteins move along the crowded DNA substrates that are likely to exist in physiological settings. This review summarizes recent progress in understanding how DNA-binding motor proteins respond to the presence of other proteins that lie in their paths. We highlight recent single-molecule biophysical experiments aimed at addressing this question, with an emphasis placed on analyzing the single-molecule, ensemble biochemical, and in vivo data from a mechanistic perspective.

RNA polymerase and the ribosome: the close relationship

In bacteria transcription and translation are linked in time and space. When coupled to RNA polymerase (RNAP), the translating ribosome ensures transcriptional processivity by preventing RNAP backtracking. Recent advances in the field have characterized important linker proteins that bridge the gap between transcription and translation: In particular, the NusE(S10):NusG complex and the NusG homolog, RfaH. The direct link between the moving ribosome and RNAP provides a basis for maintaining genomic integrity while enabling efficient transcription and timely translation of various genes within the bacterial cell.

Interference between DNA replication and transcription as a cause of genomic instability

Replication and transcription are key aspects of DNA metabolism that take place on the same template and potentially interfere with each other. Conflicts between these two activities include head-on or co-directional collisions between DNA and RNA polymerases, which can lead to the formation of DNA breaks and chromosome rearrangements. To avoid these deleterious consequences and prevent genomic instability, cells have evolved multiple mechanisms preventing replication forks from colliding with the transcription machinery. Yet, recent reports indicate that interference between replication and transcription is not limited to physical interactions between polymerases and that other cotranscriptional processes can interfere with DNA replication. These include DNA-RNA hybrids that assemble behind elongating RNA polymerases, impede fork progression and promote homologous recombination. Here, we discuss recent evidence indicating that R-loops represent a major source of genomic instability in all organisms, from bacteria to human, and are potentially implicated in cancer development.

CarD: a new RNA polymerase modulator in mycobacteria

Mycobacteria CarD is an essential RNAP binding protein that regulates many transcripts including rRNA. This article will review our present state of knowledge regarding CarD and compare the known functions of CarD with other RNAP binding proteins in E. coli, emphasizing how this information can guide future investigations.

Replication-transcription conflicts in bacteria

DNA replication and transcription use the same template and occur concurrently in bacteria. The lack of temporal and spatial separation of these two processes leads to their conflict, and failure to deal with this conflict can result in genome alterations and reduced fitness. In recent years major advances have been made in understanding how cells avoid conflicts between replication and transcription and how such conflicts are resolved when they do occur. In this Review, we summarize these findings, which shed light on the significance of the problem and on how bacterial cells deal with unwanted encounters between the replication and transcription machineries.

Transcription as a source of genome instability

Alterations in genome sequence and structure contribute to somatic disease, affect the fitness of subsequent generations and drive evolutionary processes. The crucial roles of highly accurate replication and efficient repair in maintaining overall genome integrity are well-known, but the more localized stability costs that are associated with transcribing DNA into RNA molecules are less appreciated. Here we review the diverse ways in which the essential process of transcription alters the underlying DNA template and thereby modifies the genetic landscape.

What happens when replication and transcription complexes collide?

The arrest of replication forks due to collisions with transcription complexes leads to genomic instability and cell death. Mechanisms that promote the progression of replication forks past transcription complexes are therefore essential for propagation and preservation of the genome. Recent studies of E. coli directly investigate the consequences of collisions of the replisome with RNAP polymerase (RNAP) in vitro and provide novel mechanisms by which these encounters may be resolved. Additionally, recent in vivo and in vitro studies support the longstanding hypothesis that auxiliary DNA helicases promote replication through roadblocks such as transcription complexes. Here we review past and recent advances that formulate our current understanding of how the bacterial replisome deals with transcription complexes along the path of chromosome duplication.

The chloroplast genome of the green alga Schizomeris leibleinii (Chlorophyceae) provides evidence for bidirectional DNA replication from a single origin in the chaetophorales

In the Chlorophyceae, the chloroplast genome is extraordinarily fluid in architecture and displays unique features relative to other groups of green algae. For the Chaetophorales, 1 of the 5 major lineages of the Chlorophyceae, it has been shown that the distinctive architecture of the 223,902-bp genome of Stigeoclonium helveticum is consistent with bidirectional DNA replication from a single origin. Here, we report the 182,759-bp chloroplast genome sequence of Schizomeris leibleinii, a member of the earliest diverging lineage of the Chaetophorales. Like its Stigeoclonium homolog, the Schizomeris genome lacks a large inverted repeat encoding the rRNA operon and displays a striking bias in coding regions that is associated with a bias in base composition along each strand. Our results support the notion that these two chaetophoralean genomes replicate bidirectionally from a putative origin located in the vicinity of the small subunit ribosomal RNA gene. Their shared structural characteristics were most probably inherited from the common ancestor of all chaetophoralean algae. Short dispersed repeats account for most of the 41-kb size variation between the Schizomeris and Stigeoclonium genomes, and there is no indication that homologous recombination between these repeated elements led to the observed gene rearrangements. A comparison of the extent of variation sustained by the Stigeoclonium and Schizomeris chloroplast DNAs (cpDNAs) with that observed for the cpDNAs of the chlamydomonadalean Chlamydomonas and Volvox suggests that gene rearrangements as well as changes in the abundance of intergenic and intron sequences occurred at a slower pace in the Chaetophorales than in the Chlamydomonadales.

Co-orientation of replication and transcription preserves genome integrity

In many bacteria, there is a genome-wide bias towards co-orientation of replication and transcription, with essential and/or highly-expressed genes further enriched co-directionally. We previously found that reversing this bias in the bacterium Bacillus subtilis slows replication elongation, and we proposed that this effect contributes to the evolutionary pressure selecting the transcription-replication co-orientation bias. This selection might have been based purely on selection for speedy replication; alternatively, the slowed replication might actually represent an average of individual replication-disruption events, each of which is counter-selected independently because genome integrity is selected. To differentiate these possibilities and define the precise forces driving this aspect of genome organization, we generated new strains with inversions either over ∼1/4 of the chromosome or at ribosomal RNA (rRNA) operons. Applying mathematical analysis to genomic microarray snapshots, we found that replication rates vary dramatically within the inverted genome. Replication is moderately impeded throughout the inverted region, which results in a small but significant competitive disadvantage in minimal medium. Importantly, replication is strongly obstructed at inverted rRNA loci in rich medium. This obstruction results in disruption of DNA replication, activation of DNA damage responses, loss of genome integrity, and cell death. Our results strongly suggest that preservation of genome integrity drives the evolution of co-orientation of replication and transcription, a conserved feature of genome organization.

An important feature of genome organization is that transcription and replication are selectively co-oriented. This feature helps to avoid conflicts between head-on replication and transcription. The precise consequences of the conflict and how it affects genome organization remain to be understood. We previously found that reversing the transcription bias slows replication in the Bacillus subtilis genome. Here we engineered new inversions to avoid changes in other aspects of genome organization. We found that the reversed transcription bias is sufficient to decrease replication speed, and it results in lowered fitness of the inversion strains and a competitive disadvantage relative to wild-type cells in minimal medium. Further, by analyzing genomic copy-number snapshots to obtain replication speed as a function of genome position, we found that inversion of the strongly-transcribed rRNA genes obstructs replication during growth in rich medium. This confers a strong growth disadvantage to cells in rich medium, turns on DNA damage responses, and leads to cell death in a subpopulation of cells, while the surviving cells are more sensitive to genotoxic agents. Our results strongly support the hypothesis that evolution has favored co-orientation of transcription with replication, mainly to avoid these effects.

Cloning yeast actin cDNA leads to an investigative approach for the molecular biology laboratory

The emergence of molecular tools in multiple disciplines has elevated the importance of undergraduate laboratory courses that train students in molecular biology techniques. Although it would also be desirable to provide students with opportunities to apply these techniques in an investigative manner, this is generally not possible in the classroom because of the preparation, expense, and logistics involved in independent student projects. The authors have designed a 10-week lab series that mimics the research environment by tying separate fundamental lab techniques to a common goal: to build a plasmid with yeast actin cDNA cloned in a particular orientation. In the process of completing this goal, a problem arises in that students are unable to obtain the target plasmid and instead only recover the gene cloned in the opposite orientation. To address this problem, students identify four plausible hypotheses and work in teams to address them by designing and executing experiments. This project reinforces the utility and flexibility of techniques covered earlier in the class and serves to develop their skills in experimental design and analysis. As the project is focused on one problem, the diversity of experimental approaches is limited and may be prepared in advance with little additional expense in reagents or technical support.

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.

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.

Distinctive architecture of the chloroplast genome in the chlorophycean green alga Stigeoclonium helveticum

The chloroplast genome has experienced many architectural changes during the evolution of chlorophyte green algae, with the class Chlorophyceae displaying the lowest degree of ancestral traits. We have previously shown that the completely sequenced chloroplast DNAs (cpDNAs) of Chamydomonas reinhardtii (Chlamydomonadales) and Scenedesmus obliquus (Sphaeropleales) are highly scrambled in gene order relative to one another. Here, we report the complete cpDNA sequence of Stigeoclonium helveticum (Chaetophorales), a member of a third chlorophycean lineage. This genome, which encodes 97 genes and contains 21 introns (including four putatively trans-spliced group II introns inserted at novel sites), is remarkably rich in derived features and extremely rearranged relative to its chlorophycean counterparts. At 223,902 bp, Stigeoclonium cpDNA is the largest chloroplast genome sequenced thus far, and in contrast to those of Chlamydomonas and Scenedesmus, features no large inverted repeat. Interestingly, the pattern of gene distribution between the DNA strands and the bias in base composition along each strand suggest that the Stigeoclonium genome replicates bidirectionally from a single origin. Unlike most known trans-spliced group II introns, those of Stigeoclonium exhibit breaks in domains I and II. By placing our comparative genome analyses in a phylogenetic framework, we inferred an evolutionary scenario of the mutational events that led to changes in genome architecture in the Chlorophyceae.

Living with genome instability: the adaptation of phytoplasmas to diverse environments of their insect and plant hosts

Phytoplasmas ("Candidatus Phytoplasma," class Mollicutes) cause disease in hundreds of economically important plants and are obligately transmitted by sap-feeding insects of the order Hemiptera, mainly leafhoppers and psyllids. The 706,569-bp chromosome and four plasmids of aster yellows phytoplasma strain witches' broom (AY-WB) were sequenced and compared to the onion yellows phytoplasma strain M (OY-M) genome. The phytoplasmas have small repeat-rich genomes. This comparative analysis revealed that the repeated DNAs are organized into large clusters of potential mobile units (PMUs), which contain tra5 insertion sequences (ISs) and genes for specialized sigma factors and membrane proteins. So far, these PMUs appear to be unique to phytoplasmas. Compared to mycoplasmas, phytoplasmas lack several recombination and DNA modification functions, and therefore, phytoplasmas may use different mechanisms of recombination, likely involving PMUs, for the creation of variability, allowing phytoplasmas to adjust to the diverse environments of plants and insects. The irregular GC skews and the presence of ISs and large repeated sequences in the AY-WB and OY-M genomes are indicative of high genomic plasticity. Nevertheless, segments of approximately 250 kb located between the lplA and glnQ genes are syntenic between the two phytoplasmas and contain the majority of the metabolic genes and no ISs. AY-WB appears to be further along in the reductive evolution process than OY-M. The AY-WB genome is approximately 154 kb smaller than the OY-M genome, primarily as a result of fewer multicopy sequences, including PMUs. Furthermore, AY-WB lacks genes that are truncated and are part of incomplete pathways in OY-M.

The complete plastid genome sequence of the parasitic green alga Helicosporidium sp. is highly reduced and structured

BACKGROUND

Loss of photosynthesis has occurred independently in several plant and algal lineages, and represents a major metabolic shift with potential consequences for the content and structure of plastid genomes. To investigate such changes, we sequenced the complete plastid genome of the parasitic, non-photosynthetic green alga, Helicosporidium.

RESULTS

The Helicosporidium plastid genome is among the smallest known (37.5 kb), and like other plastids from non-photosynthetic organisms it lacks all genes for proteins that function in photosynthesis. Its reduced size results from more than just loss of genes, however; it has little non-coding DNA, with only one intron and tiny intergenic spaces, and no inverted repeat (no duplicated genes at all). It encodes precisely the minimal complement of tRNAs needed to translate the universal genetic code, and has eliminated all redundant isoacceptors. The Helicosporidium plastid genome is also highly structured, with each half of the circular genome containing nearly all genes on one strand. Helicosporidium is known to be related to trebouxiophyte green algae, but the genome is structured and compacted in a manner more reminiscent of the non-photosynthetic plastids of apicomplexan parasites.

CONCLUSION

Helicosporidium contributes significantly to our understanding of the evolution of plastid DNA because it illustrates the highly ordered reduction that occurred following the loss of a major metabolic function. The convergence of plastid genome structure in Helicosporidium and the Apicomplexa raises the interesting possibility that there are common forces that shape plastid genomes, subsequent to the loss of photosynthesis in an organism.

Natural selection is not required to explain universal compositional patterns in rRNA secondary structure categories

We have encountered an unexpected property of rRNA secondary structures that may generalize to all RNAs. Analysis of 8892 ribosomal RNA sequences and structures from a wide range of species revealed unexpected universal compositional trends. First, different categories of rRNA secondary structure (stems, loops, bulges, and junctions) have distinct, characteristic base compositions. Second, the observed patterns of variation are similar among sequences from large and small rRNA subunits and all domains of life, despite extensive evolutionary divergence. Surprisingly, these differences do not seem to be related to selection for different compositions in different structural categories, but rather relate to the overall composition of the molecule: Randomized RNAs with no evolutionary history show the same structure-dependent compositional biases as rRNAs. These compositional trends may improve the accuracy of RNA secondary structure prediction, because they allow us to compare predicted structures against known compositional preferences. They also suggest caution in interpreting differences in the rate of change of the GC content in different parts of the molecule as evidence of differential selection.

A study on the correlation of nucleotide skews and the positioning of the origin of replication: different modes of replication in bacterial species

Deviations from Chargaff's 2nd parity rule, according to which A approximately T and G approximately C in single stranded DNA, have been associated with replication as well as with transcription in prokaryotes. Based on observations regarding mainly the transcription-replication co-linearity in a large number of prokaryotic species, we formulate the hypothesis that the replication procedure may follow different modes between genomes throughout which the skews clearly follow different patterns. We draw the conclusion that multiple functional sites of origin of replication may exist in the genomes of most archaea and in some exceptional cases of eubacteria, while in the majority of eubacteria, replication occurs through a single fixed origin.

Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705

BACKGROUND

Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes.

RESULTS

Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions.

CONCLUSION

Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

Mitochondrial-encoded membrane protein transcripts are pyrimidine-rich while soluble protein transcripts and ribosomal RNA are purine-rich

BACKGROUND

Eukaryotic organisms contain mitochondria, organelles capable of producing large amounts of ATP by oxidative phosphorylation. Each cell contains many mitochondria with many copies of mitochondrial DNA in each organelle. The mitochondrial DNA encodes a small but functionally critical portion of the oxidative phosphorylation machinery, a few other species-specific proteins, and the rRNA and tRNA used for the translation of these transcripts. Because the microenvironment of the mitochondrion is unique, mitochondrial genes may be subject to different selectional pressures than those affecting nuclear genes.

RESULTS

From an analysis of the mitochondrial genomes of a wide range of eukaryotic species we show that there are three simple rules for the pyrimidine and purine abundances in mitochondrial DNA transcripts. Mitochondrial membrane protein transcripts are pyrimidine rich, rRNA transcripts are purine-rich and the soluble protein transcripts are purine-rich. The transitions between pyrimidine and purine-rich regions of the genomes are rapid and are easily visible on a pyrimidine-purine walk graph. These rules are followed, with few exceptions, independent of which strand encodes the gene. Despite the robustness of these rules across a diverse set of species, the magnitude of the differences between the pyrimidine and purine content is fairly small. Typically, the mitochondrial membrane protein transcripts have a pyrimidine richness of 56%, the rRNA transcripts are 55% purine, and the soluble protein transcripts are only 53% purine.

CONCLUSION

The pyrimidine richness of mitochondrial-encoded membrane protein transcripts is partly driven by U nucleotides in the second codon position in all species, which yields hydrophobic amino acids. The purine-richness of soluble protein transcripts is mainly driven by A nucleotides in the first codon position. The purine-richness of rRNA is also due to an abundance of A nucleotides. Possible mechanisms as to how these trends are maintained in mtDNA genomes of such diverse ancestry, size and variability of A-T richness are discussed.

A genomic island present along the bacterial chromosome of the Parachlamydiaceae UWE25, an obligate amoebal endosymbiont, encodes a potentially functional F-like conjugative DNA transfer system

BACKGROUND

The genome of Protochlamydia amoebophila UWE25, a Parachlamydia-related endosymbiont of free-living amoebae, was recently published, providing the opportunity to search for genomic islands (GIs).

RESULTS

On the residual cumulative G+C content curve, a G+C-rich 19-kb region was observed. This sequence is part of a 100-kb chromosome region, containing 100 highly co-oriented ORFs, flanked by two 17-bp direct repeats. Two identical gly-tRNA genes in tandem are present at the proximal end of this genetic element. Several mobility genes encoding transposases and bacteriophage-related proteins are located within this chromosome region. Thus, this region largely fulfills the criteria of GIs. The G+C content analysis shows that several modules compose this GI. Surprisingly, one of them encodes all genes essential for F-like conjugative DNA transfer (traF, traG, traH, traN, traU, traW, and trbC), involved in sex pilus retraction and mating pair stabilization, strongly suggesting that, similarly to the other F-like operons, the parachlamydial tra unit is devoted to DNA transfer. A close relatedness of this tra unit to F-like tra operons involved in conjugative transfer is confirmed by phylogenetic analyses performed on concatenated genes and gene order conservation. These analyses and that of gly-tRNA distribution in 140 GIs suggest a proteobacterial origin of the parachlamydial tra unit.

CONCLUSIONS

A GI of the UWE25 chromosome encodes a potentially functional F-like DNA conjugative system. This is the first hint of a putative conjugative system in chlamydiae. Conjugation most probably occurs within free-living amoebae, that may contain hundreds of Parachlamydia bacteria tightly packed in vacuoles. Such a conjugative system might be involved in DNA transfer between internalized bacteria. Since this system is absent from the sequenced genomes of Chlamydiaceae, we hypothesize that it was acquired after the divergence between Parachlamydiaceae and Chlamydiaceae, when the Parachlamydia-related symbiont was an intracellular bacteria. It suggests that this heterologous DNA was acquired from a phylogenetically-distant bacteria sharing an amoebal vacuole. Since Parachlamydiaceae are emerging agents of pneumonia, this GI might be involved in pathogenicity. In future, conjugative systems might be developed as genetic tools for Chlamydiales.