Genome plasticity and ori-ter rebalancing in Salmonella typhi

Genome plasticity resulting from frequent rearrangement of the bacterial genome is a fascinating but poorly understood phenomenon. First reported in Salmonella typhi, it has been observed only in a small number of Salmonella serovars, although the over 2,500 known Salmonella serovars are all very closely related. To gain insights into this phenomenon and elucidate its roles in bacterial evolution, especially those involved in the formation of particular pathogens, we systematically analyzed the genomes of 127 wild-type S. typhi strains isolated from many places of the world and compared them with the two sequenced strains, Ty2 and CT18, attempting to find possible associations between genome rearrangement and other significant genomic features. Like other host-adapted Salmonella serovars, S. typhi contained large genome insertions, including the 134 kb Salmonella pathogenicity island, SPI7. Our analyses showed that SPI7 disrupted the physical balance of the bacterial genome between the replication origin (ori) and terminus (ter) when this DNA segment was inserted into the genome, and rearrangement in individual strains further changed the genome balance status, with a general tendency toward a better balanced genome structure. In a given S. typhi strain, genome diversification occurred and resulted in different structures among cells in the culture. Under a stressed condition, bacterial cells with better balanced genome structures were selected to greatly increase in proportion; in such cases, bacteria with better balanced genomes formed larger colonies and grew with shorter generation times. Our results support the hypothesis that genome plasticity as a result of frequent rearrangement provides the opportunity for the bacterial genome to adopt a better balanced structure and thus eventually stabilizes the genome during evolution.     

Surrogate origins of replication in the mitochondrial genomes of ori-zero petite mutants of yeast.

We have investigated the mitochondrial genome of eight ori-zero spontaneous petite mutants of Saccharomyces cerevisiae. The tandem repeat units of these genomes do not contain any of the seven canonical ori sequences of the wild-type genome. Instead, they contain one, or more, ori-S sequences. These 44-nucleotide long surrogate origins of replication are a subset of GC clusters characterized by a potential secondary fold with two sequences ATAG and GGAG , inserted in AT spacers, two AT base pairs just following them, a GC stem (broken in the middle, and, in most cases also near the base, by non-paired nucleotides), and a terminal loop. This structure is reminiscent of that of GC clusters A and B from canonical ori sequences and supports the view (Bernardi, 1982a ) that the GC clusters of the mitochondrial genome arose, by an expansion process, from the canonical ori sequences. Like the latter, ori-S sequences are present in both orientations, are located in intergenic regions, and can be used as excision sequences when tandemly oriented. Again as in the case of canonical ori sequences, the density of ori-S sequences on the repeat units of petite genomes are correlated with the replication efficiency of the latter, as assessed by the outcome of crosses with wild-type or petite tester strains.     

The Z curve database: a graphic representation of genome sequences

MOTIVATION: Genome projects for many prokaryotic and eukaryotic species have been completed and more new genome projects are being underway currently. The availability of a large number of genomic sequences for researchers creates a need to find graphic tools to study genomes in a perceivable form. The Z curve is one of such tools available for visualizing genomes. The Z curve is a unique three-dimensional curve representation for a given DNA sequence in the sense that each can be uniquely reconstructed given the other. The Z curve database for more than 1000 genomes have been established here. RESULTS: The database contains the Z curves for archaea, bacteria, eukaryota, organelles, phages, plasmids, viroids and viruses, whose genomic sequences are currently available. All the 3-dimensional Z curves and their three component curves are stored in the database. The applications of the Z curve database on comparative genomics, gene prediction, computation of G+C content with a windowless technique, prediction of replication origins and terminations of bacterial and archaeal genomes and study of local deviations from the Chargaff Parity Rule 2 etc. are presented in detail. The Z curve database reported here is a treasure trove in which biologists could find useful biological knowledge.     

Regions flanking ori sequences affect the replication efficiency of the mitochondrial genome of ori+ petite mutants from yeast

The mitochondrial genomes of progenies from 26 crosses between 17 cytoplasmic, spontaneous, suppressive, ori+ petite mutants of Saccharomyces cerevisiae have been studied by electrophoresis of restriction fragments. Only parental genomes (or occasionally, genomes derived from them by secondary excisions) were found in the progenies of the almost 500 diploids investigated; no evidence for illegitimate, site-specific mitochondrial recombination was detected. One of the parental genomes was always found to be predominate over the other one, although to different extents in different crosses. This predominance appears to be due to a higher replication efficiency, which is correlated with a greater density of ori sequences on the mitochondrial genome (and with a shorter repeat unit size of the latter). Exceptions to the 'repeat-unit-size rule' were found, however, even when the parental mitochondrial genomes carried the same ori sequence. This indicates that noncoding, intergenic sequences outside ori sequences also play a role in modulating replication efficiency. Since in different petites such sequences differ in primary structure, size, and position relative to ori sequences, this modulation is likely to take place through an indirect effect on DNA and nucleoid structure.     

Conversion at large intergenic regions of mitochondrial DNA in Saccharomyces cerevisiae.

Saccharomyces cerevisiae mitochondrial DNA deletion mutants have been used to examine whether base-biased intergenic regions of the genome influence mitochondrial biogenesis. One strain (delta 5.0) lacks a 5-kilobase (kb) segment extending from the proline tRNA gene to the small rRNA gene that includes ori1, while a second strain (delta 3.7) is missing a 3.7-kb region between the genes for ATPase subunit 6 and glutamic acid tRNA that encompasses ori7 plus ori2. Growth of these strains on both fermentable and nonfermentable substrates does not differ from growth of the wild-type strain, indicating that the deletable regions of the genome do not play a direct role in the expression of mitochondrial genes. Examination of whether the 5- or 3.7-kb regions influence mitochondrial DNA transmission was undertaken by crossing strains and examining mitochondrial genotypes in zygotic colonies. In a cross between strain delta 5.0, harboring three active ori elements (ori2, ori3, and ori5), and strain delta 3.7, containing only two active ori elements (ori3 and ori5), there is a preferential recovery of the genome containing two active ori elements (37% of progeny) over that containing three active elements (20%). This unexpected result, suggesting that active ori elements do not influence transmission of respiratory-competent genomes, is interpreted to reflect a preferential conversion of the delta 5.0 genome to the wild type (41% of progeny). Supporting evidence for conversion over biased transmission is shown by preferential recovery of a nonparental genome in the progeny of a heterozygous cross in which both parental molecules can be identified by size polymorphisms.     

Two types of recombination hotspots in bacteriophage T4: one requires DNA damage and a replication origin and the other does not

Recombination hotspots have previously been discovered in bacteriophage T4 by two different approaches, marker rescue recombination from heavily damaged phage genomes and recombination during co-infection by two undamaged phage genomes. The phage replication origin ori(34) is located in a region that has a hotspot in both assays. To determine the relationship between the origin and the two kinds of hotspots, we generated phage carrying point mutations that should inactivate ori(34) but not affect the gene 34 reading frame (within which ori(34) is located). The mutations eliminated the function of the origin, as judged by both autonomous replication of plasmids during T4 infection and two-dimensional gel analysis of phage genomic replication intermediates. As expected from past studies, the ori(34) mutations also eliminated the hotspot for marker rescue recombination from UV-irradiated genomes. However, the origin mutations had no effect on the recombination hotspot that is observed with co-infecting undamaged phage genomes, demonstrating that some DNA sequence other than the origin is responsible for inflated recombination between undamaged genomes. The hotspots for marker rescue recombination may result from a replication fork restart process that acts upon origin-initiated replication forks that become blocked at nearby DNA damage. The two-dimensional gel analysis also revealed phage T4 replication intermediates not previously detected by this method, including origin theta forms.     

Temperature can reversibly modify the structure and the functional efficiency of ori sequences of the yeast mitochondrial genome

We have compared the suppressibility of three isonuclear spontaneous, cytoplasmic petite mutants of Saccharomyces cerevisiae, as measured at three temperatures, 23 degrees C, 28 degrees C and 33 degrees C. The three petites have mitochondrial genomes made up of repeat units which are about 400 bp in size, and carry an origin of replication, ori1. This ori sequence is intact in petite Z1, whereas it lacks GC cluster A in petite 26 and cluster A plus some contiguous nucleotides in petite 14. These deletions lead to the impossibility to form a stem-and-loop structure of the ori sequence, the 'A-B fold', which involves two GC clusters, A and B, and the nucleotides in between. Instead, a 'replacement fold', only involving AT base pairs, is feasible. In petites 14 and 26, suppressivity decreases when the temperature is raised from 28 degrees C to 33 degrees C, and increases when the temperature is lowered from 28 degrees C to 23 degrees C. In contrast, no changes are seen in petite Z1. These temperature effects correlate with the stability of the 'A-B fold' and the instability of the 'replacement folds'. Since suppressibility measures the replicative competitiveness of the petite genome relative to the wild-type genome, these results indicate that an environmental parameter, temperature, can reversibly affect the structure and the functional efficiency of ori sequences in vivo. The evolutionary implications of these findings are discussed.     

Causal analysis of CpG suppression in the Mycoplasma genome

Some bacterial genomes are known to have low CpG dinucleotide frequencies. While their causes are not clearly understood, the frequency of CpG is suppressed significantly in the genome of Mycoplasma genitalium, but not in that of Mycoplasma pneumoniae. We compared orthologous gene pairs of the two closely related species to analyze CpG substitution patterns between these two genomes. We also divided genome sequences into three regions: protein-coding, noncoding, and RNA-coding, and obtained the CpG frequencies for each region for each organism. It was found that the observed/expected ratio of CpG dinucleotides is low in both the protein-coding and noncoding regions; while that ratio is in the normal range in the RNA-coding region. Our results indicate that CpG suppression of the Mycoplasma genome is not caused by (1) biased usage amino acid; (2) biased usage of synonymous codon; or (3) methylation effects by the CpG methyltransferase in the genomes of their hosts. Instead, we consider it likely that a certain global pressure, such as genome-wide pressure for the advantages of DNA stability or replication, has the effect of decreasing CpG over the entire genome, which, in turn, resulted in the biased codon usage.     

Analysis of distribution indicates diverse functions of simple sequence repeats in Mycoplasma genomes

Simple sequence repeats (SSRs) composed of extensive tandem iterations of a single nucleotide or a short oligonucleotide are rare in most bacterial genomes, but they are common among Mycoplasma. Some of these repeats act as contingency loci in association with families of surface antigens. By contraction or expansion during replication, these SSRs increase genetic variance of the population and facilitate avoidance of the immune response of the host. Occurrence and distribution of SSRs are analyzed in complete genomes of 11 Mycoplasma and 3 related Mollicutes in order to gain insights into functional and evolutionary diversity of the SSRs in Mycoplasma. The results revealed an unexpected variety of SSRs with respect to their distribution and composition and suggest that it is unlikely that all SSRs function as contingency loci or recombination hot spots. Various types of SSRs are most abundant in Mycoplasma hyopneumoniae, whereas Mycoplasma penetrans, Mycoplasma mobile, and Mycoplasma synoviae do not contain unusually long SSRs. Mycoplasma hyopneumoniae and Mycoplasma pulmonis feature abundant short adenine and thymine runs periodically spaced at 11 and 12 bp, respectively, which likely affect the supercoiling propensities of the DNA molecule. Physiological roles of long adenine and thymine runs in M. hyopneumoniae appear independent of location upstream or downstream of genes, unlike contingency loci that are typically located in protein-coding regions or upstream regulatory regions. Comparisons among 3 M. hyopneumoniae strains suggest that the adenine and thymine runs are rarely involved in genome rearrangements. The results indicate that the SSRs in the Mycoplasma genomes play diverse roles, including modulating gene expression as contingency loci, facilitating genome rearrangements via recombination, affecting protein structure and possibly protein-protein interactions, and contributing to the organization of the DNA molecule in the cell.     

Single replication origin of the archaeon Methanosarcina mazei revealed by the Z curve method

The genomic sequence of the archaeon Methanosarcina mazei has been analyzed by the Z curve method. The Z curve is a three-dimensional curve that uniquely represents the given DNA sequence. The three-dimensional Z curve and its x and y components for the genome of M. mazei show a sharp peak and relatively broad peak, respectively. The cdc6 gene is located exactly at the position of the sharp peak. Based on the known behavior of the Z curves for the archaea whose replication origins have been identified, we hypothesize that the replication origin and termination sites correspond to the positions of the sharp peak and broad peak, respectively. We have located an intergenic region that is between the cdc6 gene (MM1314) and the gene for an adjacent protein (MM1315), which shows strong characteristics of the known replication origins. This region is highly rich in AT and contains multiple copies of consecutive repeats. Our results strongly suggest that the single replication origin of M. mazei is situated at the intergenic region between the cdc6 gene and the gene for the adjacent protein, from 1,564,657 to 1,566,241 bp of the genome.     

A region of extreme instability in the mitochondrial genome of yeast.

About half of the spontaneous petite mutants produced by wild-type Saccharomyces cerevisiae strain B (as well as by several other strains) have the same defective mitochondrial genome. Its repeat unit is a segment, 2200 base pairs (bp) long, which derives from an excision between the origins of replication ori 2 and ori 7 of the wild-type genome, and contains a hybrid ori 2-ori 7 sequence. The spontaneous petites carrying this defective ori.h genome are supersuppressive , i.e., they very rapidly compete out the wild-type genome in crosses. The main reasons for the exceptional frequency of ori.h petites are an extremely high excision frequency, due to the extended homology between the two tandemly oriented ori sequences 265 bp long and the short distance separating them. Such an excision frequency is very strongly increased in petite genomes encompassing the ori 2-ori 7 region, because of their higher concentration in these ori sequences.     

Transmission of the yeast mitochondrial genome to progeny: The impact of intergenic sequences

Summary In a previous publication it was shown that the output of yeast mitochondrial loci lacking nearby intergenic sequences (encompassing ori/rep elements) was reduced in crosses to strains with wild-type mtDNAs. In the present work, mitochondrial genomes carrying the intergenic deletions were marked at unlinked, loci by introducing specific antibiotic resistance mutations against erythromycin, oligomycin and paromomycin. These marked genomes were used to follow the output of unlinked regions of the genome from crosses between the intergenic deletion mutants and wild-type strains. Transmission of genetically unlinked markers in coding regions was substantially reduced when an intergenic deletion was present on the same genome. In general the transmission of the antibiotic markers was the same as or slightly higher than the corresponding intergenic marker. These results indicate that the presence of an intergenic deletion in the regions studied impairs the transmission to progeny of a mitochondrial genome as a whole. More specifically, the results suggest that ori/rep sequences, present in the regions that have been deleted, confer a competitive advantage over genomes lacking a full complement of such sequences. These results support the hypothesis that intergenic sequences, and specifically ori/rep elements, have a biological role in the mitochondrial genome. However, because of the exclusive presence of ori/rep sequences in the genus Saccharomyces, it may be that these sequences evolved in (or invaded) the mitochondrial genome relatively late in the evolution of the yeasts. Therefore, in a more general sense, variations in the amount and structure of intergenic sequences in various yeasts may reflect processes that have been of selective advantage in the metabolism of individual mitochondrial DNA in a particular environment and that have not drastically interrupted the respiratory phenotype.     

The orir to ori+ mutation in spontaneous yeast petites is accompanied by a drastic change in mitochondrial genome replication

The orir petite mutants of Saccharomyces cerevisiae show a very low level of suppressivity (5-12%; suppressivity is the percentage of diploid petites issued from a cross of the parental haploid petite with a wild-type cell), indicating a poor replication efficiency of their mitochondrial genome. The latter is made up of repeat units containing two inverted ori sequences and arranged as tandem pairs in inverted orientation relative to their nearest neighbors. After subcloning orir petites or crossing with wild-type cells a large number of ori+ petites are found in the progeny. In contrast to the orir petites, from which they are derived, these ori+ petites are characterized by high suppressivity levels (approx. 90%) and contain mitochondrial genomes made up of tandem repeat units containing single ori sequences. The structural changes underlying the orir to ori+ mutation are therefore accompanied by a dramatic increase in suppressivity, indicating that the elimination of inverted ori sequences causes a drastic change from very poor to very good replicative efficiency in the mitochondrial genome. Finally, crosses of ori0 petites with wild-type cells were also studied; the results obtained have clarified the reasons for the high frequency of petites having genomes similar to those of orir petites after mutagenesis with ethidium bromide.     

Sequence organization of the mitochondrial genome of yeast--a review

We have compiled the available primary structural data for the mitochondrial genome of Saccharomyces cerevisiae and have estimated the size of the remaining gaps, which represent 12-13% of the genome. The lengths of sequenced regions and of gaps lead to a new assessment of genome sizes; these range (in round figures) from 85 000 bp for the long genomes, to 78 000 bp for the short genomes, to 74 000 bp for the supershort genome of Saccharomyces carlsbergensis. These values are 8-11% higher than those previously estimated from restriction fragments. Interstrain differences concern not only facultative intervening sequences (introns) and mini-inserts, but also insertions/deletions in intergenic sequences. The primary structure appears to be extremely conserved in genes and ori sequences, and highly conserved in intergenic sequences. Since coding sequences represent at most 33-35% of the genome, at least two thirds of the genome are formed by noncoding and yet highly conserved sequences. The G + C level of genes or exon is 25%, and that of intronic open reading frames (ORFs) 22%; increasingly lower values are shown by intronic closed reading frames (CRFs), 20%, ori sequences, 19%, intergenic ORFs, 17.5% and intergenic sequences, 15%.     

The GC clusters of the mitochondrial genome of yeast and their evolutionary origin

We have studied the primary and secondary structures, the location and the orientation of the 196 GC clusters present in the 90% of the mitochondrial genome of Saccharomyces cerevisiae which have already been sequenced. The vast majority of GC clusters is located in intergenic sequences (including ori sequences, intergenic open reading frames and the gene varl which probably arose from an intergenic spacer) and in intronic closed reading frames (CRF's); most of them can be folded into stem-and-loop systems; both orientations are equally frequent. The primary structures of GC clusters permit to group them into eight families, seven of which are related to the family formed by clusters A, B and C of the ori sequences. On the basis of the present work, we propose that the latter derive from a primitive ori sequence (probably made of only a monomeric cluster C and its flanking sequences r* and r) through (i) a series of duplication inversions generating clusters A and B; and (ii) an expansion process producing the AT stretches of ori sequences. Most GC clusters apparently originated from primary clusters also derived from the primitive ori sequence in the course of its evolution towards the present ori sequences. Finally, we propose that the function of GC clusters is predominantly, or entirely, associated with the structure and organization of the mitochondrial genome of yeast and, indirectly, with the regulation of its expression.     

Analysis of tandem repeats found in 44 prokaryotic genomes

Genomes contain various types of repetitive sequences. They may be used as probes for seeking genome rearrangements because they are rather free from the natural selection if they are located in the intergenic regions. In this study, we searched for tandem repeats (TRs) in 44 prokaryotic genomes by the color-coding method and sought the signs of genome rearrangements by detailed analysis of the detected TRs. We found 13,542 tandem repeats from 44 prokaryotic genomes in total ranging from several tens to one thousand per genome. The results of statistical analysis show that TRs tend to exist on high base composition bias regions in some genomes. Moreover, we recognized the characteristic distribution patterns of equivalent TR-pairs in 12 genomes, which are expected to indicate the occurrence of whole-genome duplication (WGD) on the genomes. It is demonstrated that TRs could indeed be used for seeking genome rearrangements. Although it has not been made clear at this time whether or not WGD had occurred in prokaryotic genomes, the results of the analyses of equivalent TR-pairs in this study are thought to be evidences of WGD in these genomes.     

Remarkable sequence signatures in archaeal genomes

Complete archaeal genomes were probed for the presence of long (> or = 25 bp) oligonucleotide repeats (words). We detected the presence of many words distributed in tandem with narrow ranges of periodicity (i.e., spacer length between repeats). Similar words were not identified in genomes of non-archaeal species, namely Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Mycoplasma genitalium and Mycoplasma pneumoniae. BLAST similarity searches against the GenBank nucleotide sequence database revealed that these words were archaeal species-specific, indicating that they are of a signature character. Sequence analysis and genome viewing tools showed these repeats to be restricted to non-coding regions. Thus, archaea appear to possess a non-coding genomic signature that is absent in bacterial species. The identification of a species-specific genomic signature would be of great value to archaeal genome mapping, evolutionary studies and analyses of genome complexity.