Base pairing between Escherichia coli RNase P RNA and its substrate.

Base pairing between the substrate and the ribozyme has previously been shown to be essential for catalytic activity of most ribozymes, but not for RNase P RNA. By using compensatory mutations we have demonstrated the importance of Watson-Crick complementarity between two well-conserved residues in Escherichia coli RNase P RNA (M1 RNA), G292 and G293, and two residues in the substrate, +74C and +75C (the first and second C residues in CCA). We suggest that these nucleotides base pair (G292/+75C and G293/+74C) in the ribozyme-substrate complex and as a consequence the amino acid acceptor stem of the precursor is partly unfolded. Thus, a function of M1 RNA is to anchor the substrate through this base pairing, thereby exposing the cleavage site such that cleavage is accomplished at the correct position. Our data also suggest possible base pairing between U294 in M1 RNA and the discriminator base at position +73 of the precursor. Our findings are also discussed in terms of evolution.     

A comparison of the binding to polynucleotides of complementary and noncomplementary oligonucleotides.

The binding of oligonucleotides to synthetic polynucleotides has been studied as a control for investigations of the binding of oligonucleotides to natural RNA molecules. Only combinations that involved A-U, G-C, and G-G pairs were found to be significantly stable under the experimental conditions used here. The stability of the oligomerpolymer pairing increased with the length of the region paired and with its G + C content. Further, some different sequence isomers of the same G + C content exhibited quite different binding constants. This variability is consistent with certain sequence differences in the double-strand stacking interactions stabilizing the oligomer-polymer association. Oligomer binding was also shown to depend upon the identity of the polymer residues neighboring the binding site, indicating the effect upon oligomer binding of small changes in the single-strand conformation of the binding site. These observations validate the criteria that allow one to decide if an observed association constant of an oligomer to an RNA molecule reflects a complete complementarity between the two or not. This improves the basis for using oligonucleotide binding constants to RNA of known sequences to map secondary structure.     

DNA mismatch repair and synonymous codon evolution in mammals

It has been suggested that the differences in synonymous codon use between mammalian genes within a genome are due to differences in the efficiency of DNA mismatch repair. This hypothesis was tested by developing a model of mismatch repair, which was used to predict the expected relationship between the rate of substitution and G+C content at silent sites. It was found that the silent-substitution rate should decline with increasing G+C content over most of the G+C-content range, if it is assumed that mismatch repair is G+C biased, an assumption which is supported by data. This prediction was then tested on a set of 58 primate and artiodactyl genes. There was no evidence of a direct decline in substitution rate with increasing G+C content, for either twofold- or fourfold-degenerate sites. It was therefore concluded that variation in the efficiency of mismatch repair is not responsible for the differences in synonymous codon use between mammalian genes. In support of this conclusion, analysis of the model also showed that the parameter range over which mismatch repair can explain the differences in synonymous codon use between genes is very small.      

A novel method to calculate the G+C content of genomic DNA sequences.

The base composition of a DNA fragment or genome is usually measured by the proportion of A+T or G+C in the sequence. The G+C content along genomic sequences is usually calculated using an overlapping or non-overlapping sliding window method. The result and accuracy of such an approach depends on the size of the window and the moving distance adopted. In this paper, a novel windowless technique to calculate the G+C content of genomic sequences is proposed. By this method, the G+C content can be calculated at different "resolution". In an extreme case, the G+C content may be computed at a specific point, rather than in a window of finite size. This is particularly useful to analyze the fine variation of base composition along genomic sequences. As the first example, the variation of G+C content along each of 16 yeast chromosomes is analyzed. The G+C-rich regions with length larger than 5 kb sequences are detected and listed in details. It is found that each chromosome consists of several G+C-rich and G+C-poor regions alternatively, i.e., a mosaic structure. Another example is to analyze the G+C content for each of the two chromosomes of the Vibrio cholerae genome. Based on the variations of the G+C content in each chromosome, it is shown that some fragments in the Vibrio cholerae genome may have been transferred from other species. Especially, the position and size of the large integron island on the smaller chromosome was precisely predicted. This method would be a useful tool for analyzing genomic sequences.     

Differences in (G+C) content between species: a commentary on Forsdyke's "chromosomal viewpoint" of speciation

Forsdyke (1999) has recently argued that differences in (G+C)%, or G+C content, may trigger new species formation. He further argues that the genic model has shortcomings that can be overcome by his "chromosomal" (hereafter, "G+C") model. We disagree on several counts. First, we do not accept that the genic model has the shortcomings suggested by Forsdyke. There is an abundance of empirical support for the contribution of individual genes, as well as of mapped chromosomal regions, to post-zygotic reproductive isolation (and Haldane's rule). Further, we argue that the G+C model suffers from the same theoretical difficulties as other speciation models based on underdominance. We also question the evidence Forsdyke uses to support his model. Finally, we describe analyses of G+C content in a well-studied model system of speciation (the Drosophila melanogaster species complex), the results of which are incompatible with the G+C model. Thus, while Forsdyke's G+C model cannot be explicitly ruled out, it is not directly supported by empirical data. In contrast, the genic model is well supported by empirical data, holds up on theoretical grounds, and does not require any assistance from the G+C model.     

Nucleotide sequence and presumed secondary structure of the 28S rRNA of pea aphid implication for diversification of insect rRNA

Determination of the entire nucleotide sequence of the aphid 28S ribosomal RNA gene (28S rDNA) revealed that it is 4,147 by in length with a G + C content of 60.3%. Based on the nucleotide sequence, we constructed a presumed secondary-structure model of the aphid 28S rRNA which indicated that the aphid 28S rRNA is characterized by the length and high G + C content of its variable regions. The G + C content of the aphid's variable regions was much higher than that of the entire sequence of the 28S rRNA, which formed a striking contrast to those ofDrosophila with the G + C content much lower than the entire 28S molecule. In this respect, the aphid 28S rRNA somewhat resembled those of vertebrates. This is the third report of a complete large-subunit rRNA sequence from an arthropod, and the first 28S rRNA sequence for a nondipterous insect. Correspondence to: H. Ishikawa     

Different conformational forms of Escherichia coli and rat liver 5S rRNA revealed by Pb(II)-induced hydrolysis.

Different stable forms of Escherichia coli and rat liver 5S rRNA have been probed by Pb(II)-induced hydrolysis. In the native A forms of 5S rRNA, Pb2+ reveal single-stranded RNA stretches and regions of increased conformational flexibility or distorted by the presence of bulged nucleotides. Hydrolysis of urea/EDTA-treated E. coli 5S rRNA (B form) shows the presence of two strong helical domains; helix A retained from the A form and a helix composed of RNA regions G33-C42 and G79-C88. Other RNA regions resistant to hydrolysis may be involved in alternative base pairing, causing conformational heterogeneity of that form. Pb(II)-induced hydrolysis distinguishes two different forms of rat liver 5S rRNA; the native A form and the form obtained by renaturation of 5S rRNA in the presence of EDTA. Pb(II)-hydrolysis data suggest that both forms are highly structured. In the latter form, the orientation of the bulged C66 is changed with respect to helix B. At the same time, a new helical segment is possibly formed, composed of nucleotides from helix C and loop c on one side and from helix E and loop d' on the other.     

Discovering isochores by least-squares optimal segmentation

The isochore structure of a genome is observable by variation in the G+C (guanine and cytosine) content within and between the chromosomes. Describing the isochore structure of vertebrate genomes is a challenging task, and many computational methods have been developed and applied to it. Here we apply a well-known least-squares optimal segmentation algorithm to isochore discovery. The algorithm finds the best division of the sequence into k pieces, such that the segments are internally as homogeneous as possible. We show how this simple segmentation method can be applied to isochore discovery using as input the G+C content of sliding windows on the sequence. To evaluate the performance of this segmentation technique on isochore detection, we present results from segmenting previously studied isochore regions of the human genome. Detailed results on the MHC locus, on parts of chromosomes 21 and 22, and on a 100 Mb region from chromosome 1 are similar to previously suggested isochore structures. We also give results on segmenting all 22 autosomal human chromosomes. An advantage of this technique is that oversegmentation of G+C rich regions can generally be avoided. This is because the technique concentrates on greater global, instead of smaller local, differences in the sequence composition. The effect is further emphasized by a log-transformation of the data that lowers the high variance that is observed in G+C rich regions. We conclude that the least-squares optimal segmentation method is computationally efficient and yields results close to previous biologically motivated isochore structures.     

Evidence of selection on silent site base composition in mammals: potential implications for the evolution of isochores and junk DNA.

It has been suggested that mutation bias is the major determinant of base composition bias at synonymous, intron, and flanking DNA sites in mammals. Here I test this hypothesis using population genetic data from the major histocompatibility genes of several mammalian species. The results of two tests are inconsistent with the mutation hypothesis in coding, noncoding, CpG-island, and non-CpG-island DNA, but are consistent with selection or biased gene conversion. It is argued that biased gene conversion is unlikely to affect silent site base composition in mammals. The results therefore suggest that selection is acting upon silent site G + C content. This may have broad implications, since silent site base composition reflects large-scale variation in G + C content along mammalian chromosomes. The results therefore suggest that selection may be acting upon the base composition of isochores and large sections of junk DNA.     

An unmethylated “3 SE” RNA in hamster mitochondria: A 5 S RNA-equivalent?

A novel low molecular weight (“3 S_E”) RNA associated with hamster cell mitochondria has been partially characterized. It was present at approx. 1:1 molar ratio with structural mitochondrial ribosomal RNA; it was unmethylated; and it resembled other mitochondrial RNA fractions in having a low content of G + C. These findings support the idea that 3 S_E RNA is a mitochondrial equivalent of 5 S ribosomal RNA.     

The genetic properties of the primary endosymbionts of mealybugs differ from those of other endosymbionts of plant sap-sucking insects

Mealybugs (Hemiptera, Coccoidea, Pseudococcidae), like aphids and psyllids, are plant sap-sucking insects that have an obligate association with prokaryotic endosymbionts that are acquired through vertical, maternal transmission. We sequenced two fragments of the genome of Tremblaya princeps, the endosymbiont of mealybugs, which is a member of the beta subdivision of the Proteobacteria. Each of the fragments (35 and 30 kb) contains a copy of 16S-23S-5S rRNA genes. A total of 37 open reading frames were detected, which corresponded to putative rRNA proteins, chaperones, and enzymes of branched-chain amino acid biosynthesis, DNA replication, protein translation, and RNA synthesis. The genome of T. princeps has a number of properties that distinguish it from the genomes of Buchnera aphidicola and Carsonella ruddii, the endosymbionts of aphids and psyllids, respectively. Among these properties are a high G+C content (57.1 mol%), the same G+C content in intergenic spaces and structural genes, and similar G+C contents of the genes encoding highly and poorly conserved proteins. The high G+C content has a substantial effect on protein composition; about one-third of the residues consist of four amino acids with high-G+C-content codons. Sequence analysis of DNA fragments containing the rRNA operon and adjacent regions from endosymbionts of several mealybug species suggested that there was a single duplication of the rRNA operon and the adjacent genes in an ancestor of the present T. princeps. Subsequently, in one mealybug lineage rpS15, one of the duplicated genes, was retained, while in another lineage it decayed. These results extend the diversity of the types of endosymbiotic associations found in plant sap-sucking insects.     

Mutation bias is the driving force of codon usage in the Gallus gallus genome

Synonymous codons are used with different frequencies both among species and among genes within the same genome and are controlled by neutral processes (such as mutation and drift) as well as by selection. Up to now, a systematic examination of the codon usage for the chicken genome has not been performed. Here, we carried out a whole genome analysis of the chicken genome by the use of the relative synonymous codon usage (RSCU) method and identified 11 putative optimal codons, all of them ending with uracil (U), which is significantly departing from the pattern observed in other eukaryotes. Optimal codons in the chicken genome are most likely the ones corresponding to highly expressed transfer RNA (tRNAs) or tRNA gene copy numbers in the cell. Codon bias, measured as the frequency of optimal codons (Fop), is negatively correlated with the G + C content, recombination rate, but positively correlated with gene expression, protein length, gene length and intron length. The positive correlation between codon bias and protein, gene and intron length is quite different from other multi-cellular organism, as this trend has been only found in unicellular organisms. Our data displayed that regional G + C content explains a large proportion of the variance of codon bias in chicken. Stepwise selection model analyses indicate that G + C content of coding sequence is the most important factor for codon bias. It appears that variation in the G + C content of CDSs accounts for over 60% of the variation of codon bias. This study suggests that both mutation bias and selection contribute to codon bias. However, mutation bias is the driving force of the codon usage in the Gallus gallus genome. Our data also provide evidence that the negative correlation between codon bias and recombination rates in G. gallus is determined mostly by recombination-dependent mutational patterns.     

An isochore map of the human genome based on the Z curve method

The distribution of the G+C content in the human genome has been studied by using a windowless technique derived from the Z curve method. The most important findings presented in this paper are twofold. First, abrupt variations of the G+C content along human chromosome sequences are the main variation patterns of G+C content. It is found that at some sites, the G+C content undergoes abrupt changes from a G+C-rich region to a G+C-poor region alternatively and vice versa. Second, it is shown that long domains with relatively homogeneous G+C content along each chromosome do exist. These domains are thought to be isochores, which usually have sharp boundaries. Consequently, 56 isochores longer than 3 Mb have been identified in chromosomes 1-22, X and Y. Boundaries, size and G+C content of each isochore identified are listed in detail. As an example to demonstrate the power of the method, the boundary between the Classes III and II isochores of the MHC sequence has been determined and found to be at 2,477,936, which is in good agreement with the experimental evidence. A homogeneity index is introduced to measure the homogeneity of G+C content in isochores. We emphasize that the homogeneity of G+C content is relative. The isochores in which the G+C content keeps absolutely constant do not exist. Isochore structures appear to be a basic organization of the human genome. Due to the relevance to many important biological functions, the clarification of isochore structures will provide much insight into the understanding of the human genome.     

Characterization of the 7S RNA and its gene from halobacteria.

The 7S RNA is an abundant nonribosomal RNA in H. halobium and other halobacteria. A specific 7S RNA gene probe shows high homology to genomic DNA of all halobacteria tested but not to those of several other archaebacteria, eubacteria and eukaryotes. All halobacterial genomes seem to carry a single copy of the 7S RNA gene. The coding region of the 7S RNA gene is highly G+C rich whereas the 5'- and 3'-noncoding regions possess a rather low G+C content. An extended double stranded structure for the 7S RNA is deduced from its nucleotide sequence. The 7S RNA of H. halobium (304 nucleotides) resembles in size and structure the 7S-L RNA from mammalian cells and shares with it a sequence homology of about 50% when arranged in a colinear fashion. The similarities in sequence are found particularly at the 3'- and 5'-termini. No similarity was detected between the 7S RNA from H. halobium and the nonribosomal 6S RNA from Escherichia coli.     

The covariation between TpA deficiency, CpG deficiency, and G+C content of human isochores is due to a mathematical artifact

CpG and TpA dinucleotides are underrepresented in the human genome. The CpG deficiency is due to the high mutation rate from C to T in methylated CpG's. The TpA suppression was thought to reflect a counterselection against TpA's destabilizing effect in RNA. Unexpectedly, the TpA and CpG deficiencies vary according to the G+C contents of sequences. It has been proposed that the variation in CpG suppression was correlated with a particular chromatin organization in G+C-rich isochores. Here, we present an improved model of dinucleotide evolution accounting for the overlap between successive dinucleotides. We show that an increased mutation rate from CpG to TpG or CpA induces both an apparent TpA deficiency and a correlation between CpG and TpA deficiencies and G+C content. Moreover, this model shows that the ratio of observed over expected CpG frequency underestimates the real CpG deficiency in G+C-rich sequences. The predictions of our model fit well with observed frequencies in human genomic data. This study suggests that previously published selectionist interpretations of patterns of dinucleotide frequencies should be taken with caution. Moreover, we propose new criteria to identify unmethylated CpG islands taking into account this bias in the measure of CpG depletion.     

A conserved hairpin structure in Alfamovirus and Bromovirus subgenomic promoters is required for efficient RNA synthesis in vitro

The coat protein gene in RNA 3 of alfalfa mosaic virus (AMV; genus Alfamovirus, family Bromoviridae) is translated from the subgenomic RNA 4. Analysis of the subgenomic promoter (sgp) in minus-strand RNA 3 showed that a sequence of 37 nt upstream of the RNA 4 start site (nt +1) was sufficient for full sgp activity in an in vitro assay with the purified viral RNA-dependent RNA-polymerase (RdRp). The sequence of nt -6 to -29 could be folded into a potential hairpin structure with a loop represented by nt -16, -17, and -18, and a bulge involving nt -23. By introducing mutations that disrupted base pairing and compensatory mutations that restored base pairing, it was shown that base pairing in the top half of the putative stem (between the loop and bulge) was essential for sgp activity, whereas base pairing in the bottom half of the stem was less stringently required. Deletion of the bulged residue A-23 or mutation of this residue into a C strongly reduced sgp activity, but mutation of A-23 into U or G had little effect on sgp activity. Mutation of loop residues A-16 and A-17 affected sgp activity, whereas mutation of U-18 did not. Using RNA templates corresponding to the sgp of brome mosaic virus (BMV; genus Bromovirus, family Bromoviridae) and purified BMV RdRp, evidence was obtained indicating that also in BMV RNA a triloop hairpin structure is required for sgp activity.     

Relationships Between Genomic G+C Content,RNA Secondary Structures,and Optimal Growth Temperature in Prokaryotes

G:C pairs are more stable than A:T pairs because they have an additional hydrogen bond. This has led to many studies on the correlation between the guanine+cytosine (G+C) content of nucleic acids and temperature over the last 20 years. We collected the optimal growth temperatures (T_opt) and the G+C contents of genomic DNA; 23S, 16S, and 5S ribosomal RNAs; and transfer RNAs for 764 prokaryotic species. No correlation was found between genomic G+C content and T_opt, but there were striking correlations between the G+C content of ribosomal and transfer RNA stems and T_opt. Two explanations have been proposed—neutral evolution and selection pressure—for the approximate equalities of G and C (respectively, A and T) contents within each strand of DNA molecules. Our results do not support the notion that selection pressure induces complementary oligonucleotides in close proximity and therefore numerous secondary structures in prokaryotic DNA, as the genomic G+C content does not behave in the same way as that of folded RNA with respect to optimal growth temperature. Received: 25 September 1996 / Accepted: 21 January 1997     

Gene structure of Enterococcus hirae (Streptococcus faecalis) F1F0-ATPase, which functions as a regulator of cytoplasmic pH.

Enterococcus hirae (formerly Streptococcus faecalis) ATCC 9790 has an F1F0-ATPase which functions as a regulator of the cytoplasmic pH but does not synthesize ATP. We isolated four clones which contained genes for c, b, delta, and alpha subunits of this enzyme but not for other subunit genes. It was revealed that two specific regions (upstream of the c-subunit gene and downstream of the gamma-subunit gene) were lost at a specific site in the clones we isolated, suggesting that these regions were unstable in Escherichia coli. The deleted regions were amplified by polymerase chain reaction, and the nucleotide sequences of these regions were determined. The results showed that eight genes for a, c, b, delta, alpha, gamma, beta, and epsilon subunits were present in this order. Northern (RNA) blot analysis showed that these eight genes were transcribed to one mRNA. The i gene was not found in the upper region of the a-subunit gene. Instead of the i gene, this operon contained a long untranslated region (240 bp) whose G + C content was only 30%. There was no typical promoter sequence such as was proposed for E. coli, suggesting that the promoter structure of this species is different from that of E. coli. Deduced amino acid sequences suggested that E. hirae H(+)-ATPase is a typical F1F0-type ATPase but that its gene structure is not identical to that of other bacterial F1F0-ATPases.