Translational interference at overlapping reading frames in prokaryotic messenger RNA

In overlapping reading frames of prokaryotic mRNA, the ribosome-binding site (RBS) of the downstream cistron is part of the coding sequence of the upstream message. We have examined whether the rate of translation in Escherichia coli can be sufficiently high to preclude the use of an RBS in initiation of protein synthesis when it is part of an actively decoded reading frame. The two sets of gene overlap present in the RNA phage MS2 are used as a model system. We find that translation of an upstream cistron can fully block initiation of protein synthesis at the overlapping RBS of the downstream cistron. Nonsense mutations in the upstream gene restore the translation of the downstream gene.     

First committed step of lipid A biosynthesis in Escherichia coli: sequence of the lpxA gene.

The min 4 region of the Escherichia coli genome contains genes (lpxA and lpxB) that encode proteins involved in lipid A biosynthesis. We have determined the sequence of 1,350 base pairs of DNA upstream of the lpxB gene. This fragment of DNA contains the complete coding sequence for the 28.0-kilodalton lpxA gene product and an upstream open reading frame capable of encoding a 17-kilodalton protein (ORF17). In addition there appears to be an additional open reading frame (ORF?) immediately upstream of ORF17. The initiation codon for lpxA is a GUG codon, and the start codon for ORF17 is apparently a UUG codon. The start and stop codons overlap between ORF? and ORF17, ORF17 and lpxA, and lpxA and lpxB. This overlap is suggestive of translational coupling and argues that the genes are cotranscribed. Crowell et al. (D.N. Crowell, W.S. Reznikoff, and C.R.H. Raetz, J. Bacteriol. 169:5727-5734, 1987) and Tomasiewicz and McHenry (H.G. Tomasiewicz and C.S. McHenry, J. Bacteriol. 169:5735-5744, 1987) have demonstrated that there are three similarly overlapping coding regions downstream of lpxB including dnaE, suggesting the existence of a complex operon of at least seven genes: 5'-ORF?-ORF17-lpxA-lpxB-ORF23-dnaE-ORF37-3 '.     

Complete nucleotide sequence of the Escherichia coli recB gene.

The complete nucleotide sequence of the Escherichia coli recB gene which encodes a subunit of the ATP-dependent DNase, Exonuclease V, has been determined. The proposed coding region for the RecB protein is 3543 nucleotides long and would encode a polypeptide of 1180 amino acids with a calculated molecular weight of 133,973. The start of the recB coding sequence overlaps the 3' end of the upstream ptr gene, and the recB termination codon overlaps the initiation codon of the downstream recD gene, suggesting that these genes may form an operon. No sequences which reasonably fit the consensus for an E. coli promoter could be identified upstream of the proposed recB translational start. The predicted RecB amino acid sequence contains regions of homology with ATPases, DNA binding proteins and DNA repair enzymes.     

Engineering gene overlaps to sustain genetic constructs in vivo

Evolution is often an obstacle to the engineering of stable biological systems due to the selection of mutations inactivating costly gene circuits. Gene overlaps induce important constraints on sequences and their evolution. We show that these constraints can be harnessed to increase the stability of costly genes by purging loss-of-function mutations. We combine computational and synthetic biology approaches to rationally design an overlapping reading frame expressing an essential gene within an existing gene to protect. Our algorithm succeeded in creating overlapping reading frames in 80% of E. coli genes. Experimentally, scoring mutations in both genes of such overlapping construct, we found that a significant fraction of mutations impacting the gene to protect have a deleterious effect on the essential gene. Such an overlap thus protects a costly gene from removal by natural selection by associating the benefit of this removal with a larger or even lethal cost. In our synthetic constructs, the overlap converts many of the possible mutants into evolutionary dead-ends, reducing the evolutionary potential of the system and thus increasing its stability over time.     

Purifying and directional selection in overlapping prokaryotic genes

In overlapping genes, the same DNA sequence codes for two proteins using different reading frames. Analysis of overlapping genes can help in understanding the mode of evolution of a coding region from noncoding DNA. We identified 71 pairs of convergent genes, with overlapping 3' ends longer than 15 nucleotides, that are conserved in at least two prokaryotic genomes. Among the overlap regions, we observed a statistically significant bias towards the 123:132 phase (i.e. the second codon base in one gene facing the degenerate third position in the second gene). This phase ensures the least mutual constraint on nonconservative amino acid replacements in both overlapping coding sequences. The excess of this phase is compatible with directional (positive) selection acting on the overlapping coding regions. This could be a general evolutionary mode for genes emerging from noncoding sequences, in which the protein sequence has not been subject to selection.     

Polarity Effects in the Hisg Gene of Salmonella Require a Site within the Coding Sequence

A single site in the middle of the coding sequence of the hisG gene of Salmonella is required for most of the polar effect of mutations in this gene. Nonsense and insertion mutations mapping upstream of this point in the hisG gene all have strong polar effects on expression of downstream genes in the operon; mutations mapping promotor distal to this site have little or no polar effect. Two previously known hisG mutations, mapping in the region of the polarity site, abolish the polarity effect of insertion mutations mapping upstream of this region. New polarity site mutations have been selected which have lost the polar effect of upstream nonsense mutations. All mutations abolishing the function of the site are small deletions; three are identical, 28-bp deletions which have arisen independently. A fourth mutation is a deletion of 16 base pairs internal to the larger deletion. Several point mutations within this 16-bp region have no effect on the function of the polarity site. We believe that a small number of polarity sites of this type are responsible for polarity in all genes. The site in the hisG gene is more easily detected than most because it appears to be the only such site in the hisG gene and because it maps in the center of the coding sequence.     

Same-strand overlapping genes in bacteria: compositional determinants of phase bias

Background

Same-strand overlapping genes may occur in frameshifts of one (phase 1) or two nucleotides (phase 2). In previous studies of bacterial genomes, long phase-1 overlaps were found to be more numerous than long phase-2 overlaps. This bias was explained by either genomic location or an unspecified selection advantage. Models that focused on the ability of the two genes to evolve independently did not predict this phase bias. Here, we propose that a purely compositional model explains the phase bias in a more parsimonious manner. Same-strand overlapping genes may arise through either a mutation at the termination codon of the upstream gene or a mutation at the initiation codon of the downstream gene. We hypothesized that given these two scenarios, the frequencies of initiation and termination codons in the two phases may determine the number for overlapping genes.

Results

We examined the frequencies of initiation- and termination-codons in the two phases, and found that termination codons do not significantly differ between the two phases, whereas initiation codons are more abundant in phase 1. We found that the primary factors explaining the phase inequality are the frequencies of amino acids whose codons may combine to form start codons in the two phases. We show that the frequencies of start codons in each of the two phases, and, hence, the potential for the creation of overlapping genes, are determined by a universal amino-acid frequency and species-specific codon usage, leading to a correlation between long phase-1 overlaps and genomic GC content.

Conclusion

Our model explains the phase bias in same-strand overlapping genes by compositional factors without invoking selection. Therefore, it can be used as a null model of neutral evolution to test selection hypotheses concerning the evolution of overlapping genes.

Reviewers

This article was reviewed by Bill Martin, Itai Yanai, and Mikhail Gelfand.     

Stability and evolution of overlapping genes

Abstract.— When the same sequence of nucleotides codes for regions of more than one functional polypeptide, this sequence contains overlapping genes. Overlap is most common in rapidly evolving genomes with high mutation rates such as viruses, bacteria, and mitochondria. Overlap is thought to be important as: (1) a means of compressing a maximum amount of information into short sequences of structural genes; and (2) as a mechanism for regulating gene expression through translational coupling of functionally related polypeptides. The stability of overlapping codes is examined in relation to the information cost of overlap and the mutation rate of the genome. The degree of overlap in a given population will tend to become monomorphic. Evolution toward partial overlap of genes is shown to depend on a convex cost function of overlap. Overlap does not evolve when expression of overlapping genes is mutually exclusive and produced by rare mutations to the wild-type genome. Assuming overlap increases coupling between functionally related genes, the conditions favoring overlap are explored in relation to the kinetics of gene activation and decay. Coupling is most effective for genes in which the gene overlapping at its 5'end (leading gene) decays rapidly, while the gene overlapping at the 3'end (induced gene) decays slowly. If gene expression can feedback on itself (autocatalysis), then high rates of activation favor overlap.     

The role of population size, pleiotropy and fitness effects of mutations in the evolution of overlapping gene functions

Sheltered from deleterious mutations, genes with overlapping or partially redundant functions may be important sources of novel gene functions. While most partially redundant genes originated in gene duplications, it is much less clear why genes with overlapping functions have been retained, in some cases for hundreds of millions of years. A case in point is the many partially redundant genes in vertebrates, the result of ancient gene duplications in primitive chordates. Their persistence and ubiquity become surprising when it is considered that duplicate and original genes often diversify very rapidly, especially if the action of natural selection is involved. Are overlapping gene functions perhaps maintained because of their protective role against otherwise deleterious mutations? There are two principal objections against this hypothesis, which are the main subject of this article. First, because overlapping gene functions are maintained in populations by a slow process of "second order" selection, population sizes need to be very high for this process to be effective. It is shown that even in small populations, pleiotropic mutations that affect more than one of a gene''s functions simultaneously can slow the mutational decay of functional overlap after a gene duplication by orders of magnitude. Furthermore, brief and transient increases in population size may be sufficient to maintain functional overlap. The second objection regards the fact that most naturally occurring mutations may have much weaker fitness effects than the rather drastic "knock-out" mutations that lead to detection of partially redundant functions. Given weak fitness effects of most mutations, is selection for the buffering effect of functional overlap strong enough to compensate for the diversifying force exerted by mutations? It is shown that the extent of functional overlap maintained in a population is not only independent of the mutation rate, but also independent of the average fitness effects of mutation. These results are discussed with respect to experimental evidence on redundant genes in organismal development.     

Detection of Functional Overlapping Genes: Simulation and Case Studies

As far as protein-coding genes are concerned, there is a non-zero probability that at least one of the five possible overlapping sequences of any gene will contain an open-reading frame (ORF) of a length that may be suitable for coding a functional protein. It is, however, very difficult to determine whether or not such an ORF is functional. Recently, we proposed a method that predicts functionality of an overlapping ORF if it can be shown that it has been subject to purifying selection during its evolution. Here, we use simulation to test this method under several conditions and compare it with the method of Firth and Brown. We found that under most conditions, our method detects functional overlapping genes with higher sensitivity than Firth and Brown’s method, while maintaining high specificity. Further, we tested the hypothesis that the two aminoacyl tRNA synthetase classes have originated from a pair of overlapping genes. A central piece of evidence ostensibly supporting this hypothesis is the assertion that an overlapping ORF of a heat-shock protein-70 gene, which exhibits some similarity to class 2 aminoacyl tRNA synthetases, is functional. We found signature of purifying selection only in highly divergent sequences, suggesting that the method yields false-positives in high sequence divergence and that the overlapping ORF is not a functional gene. Finally, we examined three cases of overlap in the human genome. We find varying signatures of purifying selection acting on these overlaps, raising the possibility that two of the overlapping genes may not be functional.     

Birth and death of gene overlaps in vertebrates

Background  

Between five and fourteen per cent of genes in the vertebrate genomes do overlap sharing some intronic and/or exonic sequence. It was observed that majority of these overlaps are not conserved among vertebrate lineages. Although several mechanisms have been proposed to explain gene overlap origination the evolutionary basis of these phenomenon are still not well understood. Here, we present results of the comparative analysis of several vertebrate genomes. The purpose of this study was to examine overlapping genes in the context of their evolution and mechanisms leading to their origin.     

Influences on gene expression in vivo by a Shine-Dalgarno sequence

The Shine-Dalgarno (SD+: 5'-AAGGAGG-3') sequence anchors the mRNA by base pairing to the 16S rRNA in the small ribosomal subunit during translation initiation. We have here compared how an SD+ sequence influences gene expression, if located upstream or downstream of an initiation codon. The positive effect of an upstream SD+ is confirmed. A downstream SD+ gives decreased gene expression. This effect is also valid for appropriately modified natural Escherichia coli genes. If an SD+ is placed between two potential initiation codons, initiation takes place predominantly at the second start site. The first start site is activated if the distance between this site and the downstream SD+ is enlarged and/or if the second start site is weakened. Upstream initiation is eliminated if a stable stem-loop structure is placed between this SD+ and the upstream start site. The results suggest that the two start sites compete for ribosomes that bind to an SD+ located between them. A minor positive contribution to upstream initiation resulting from 3' to 5' ribosomal diffusion along the mRNA is suggested. Analysis of the E. coli K12 genome suggests that the SD+ or SD-like sequences are systematically avoided in the early coding region suggesting an evolutionary significance.