Transcriptional slippage occurs during elongation at runs of adenine or thymine in Escherichia coli

Avoidance of long mononucleotide repeats in codon pair usage

Protein is an essential component for life, and its synthesis is mediated by codons in any organisms on earth. While some codons encode the same amino acid, their usage is often highly biased. There are many factors that can cause the bias, but a potential effect of mononucleotide repeats, which are known to be highly mutable, on codon usage and codon pair preference is largely unknown. In this study we performed a genomic survey on the relationship between mononucleotide repeats and codon pair bias in 53 bacteria, 68 archaea, and 13 eukaryotes. By distinguishing the codon pair bias from the codon usage bias, four general patterns were revealed: strong avoidance of five or six mononucleotide repeats in codon pairs; lower observed/expected (o/e) ratio for codon pairs with C or G repeats (C/G pairs) than that with A or T repeats (A/T pairs); a negative correlation between genomic GC contents and the o/e ratios, particularly for C/G pairs; and avoidance of C/G pairs in highly conserved genes. These results support natural selection against long mononucleotide repeats, which could induce frameshift mutations in coding sequences. The fact that these patterns are found in all kingdoms of life suggests that this is a general phenomenon in living organisms. Thus, long mononucleotide repeats may play an important role in base composition and genetic stability of a gene and gene functions.

What is a microsatellite: a computational and experimental definition based upon repeat mutational behavior at A/T and GT/AC repeats

Microsatellites are abundant in eukaryotic genomes and have high rates of strand slippage-induced repeat number alterations. They are popular genetic markers, and their mutations are associated with numerous neurological diseases. However, the minimal number of repeats required to constitute a microsatellite has been debated, and a definition of a microsatellite that considers its mutational behavior has been lacking. To define a microsatellite, we investigated slippage dynamics for a range of repeat sizes, utilizing two approaches. Computationally, we assessed length polymorphism at repeat loci in ten ENCODE regions resequenced in four human populations, assuming that the occurrence of polymorphism reflects strand slippage rates. Experimentally, we determined the in vitro DNA polymerase-mediated strand slippage error rates as a function of repeat number. In both approaches, we compared strand slippage rates at tandem repeats with the background slippage rates. We observed two distinct modes of mutational behavior. At small repeat numbers, slippage rates were low and indistinguishable from background measurements. A marked transition in mutability was observed as the repeat array lengthened, such that slippage rates at large repeat numbers were significantly higher than the background rates. For both mononucleotide and dinucleotide microsatellites studied, the transition length corresponded to a similar number of nucleotides (approximately 10). Thus, microsatellite threshold is determined not by the presence/absence of strand slippage at repeats but by an abrupt alteration in slippage rates relative to background. These findings have implications for understanding microsatellite mutagenesis, standardization of genome-wide microsatellite analyses, and predicting polymorphism levels of individual microsatellite loci.

Transcription elongation complex stability: the topological lock

Transcription machinery from a variety of organisms shows striking mechanistic similarity. Both multi- and single subunit RNA polymerases have evolved an 8-10-base pair RNA-DNA hybrid as a part of a stably transcribing elongation complex. Through characterization of halted complexes that can readily carry out homopolymeric slippage synthesis, this study reveals that T7 RNA polymerase elongation complexes containing only a 4-base pair hybrid can nevertheless be more stable than those with the normal 8-base pair hybrid. We propose that a key feature of this stability is the topological threading of RNA through the complex and/or around the DNA template strand. The data are consistent with forward translocation as a mechanism to allow unthreading of the topological lock, as can occur during programmed termination of transcription.

In the non-insect-transmissible line of onion yellows phytoplasma (OY-NIM), the plasmid-encoded transmembrane protein ORF3 lacks the major promoter region

‘Candidatus Phytoplasma asteris’, onion yellows strain (OY), a mildly pathogenic line (OY-M), is a phytopathogenic bacterium transmitted by Macrosteles striifrons leafhoppers. OY-M contains two types of plasmids (EcOYM and pOYM), each of which possesses a gene encoding the putative transmembrane protein, ORF3. A non-insect-transmissible line of this phytoplasma (OY-NIM) has the corresponding plasmids (EcOYNIM and pOYNIM), but pOYNIM lacks orf3. Here we show that in OY-M, orf3 is transcribed from two putative promoters and that on EcOYNIM, one of the promoter sequences is mutated and the other deleted. We also show by immunohistochemical analysis that ORF3 is not expressed in OY-NIM-infected plants. Moreover, ORF3 protein seems to be preferentially expressed in OY-M-infected insects rather than in plants. We speculate that ORF3 may play a role in the interactions of OY with its insect host.

RNA polymerase elongation factors

The elongation phase of transcription by RNA polymerase is highly regulated and modulated. Both general and operon-specific elongation factors determine the local rate and extent of transcription to coordinate the appearance of transcript with its use as a messenger or functional ribonucleoprotein or regulatory element, as well as to provide operon-specific gene regulation.

Preferred and avoided codon pairs in three domains of life

BACKGROUND

Alternative synonymous codons are not used with equal frequencies. In addition, the contexts of codons - neighboring nucleotides and neighboring codons - can have certain patterns. The codon context can influence both translational accuracy and elongation rates. However, it is not known how strong or conserved the codon context preferences in different organisms are. We analyzed 138 organisms (bacteria, archaea and eukaryotes) to find conserved patterns of codon pairs.

RESULTS

After removing the effects of single codon usage and dipeptide biases we discovered a set of neighboring codons for which avoidances or preferences were conserved in all three domains of life. Such biased codon pairs could be divided into subtypes on the basis of the nucleotide patterns that influence the bias. The most frequently avoided type of codon pair was nnUAnn. We discovered that 95.7% of avoided nnUAnn type patterns contain out-frame UAA or UAG triplets on the sense and/or antisense strand. On average, nnUAnn codon pairs are more frequently avoided in ORFeomes than in genomes. Thus we assume that translational selection plays a major role in the avoidance of these codon pairs. Among the preferred codon pairs, nnGCnn was the major type.

CONCLUSION

Translational selection shapes codon pair usage in protein coding sequences by rules that are common to all three domains of life. The most frequently avoided codon pairs contain the patterns nnUAnn, nnGGnn, nnGnnC, nnCGCn, GUCCnn, CUCCnn, nnCnnA or UUCGnn. The most frequently preferred codon pairs contain the patterns nnGCnn, nnCAnn or nnUnCn.

Endosymbiont gene functions impaired and rescued by polymerase infidelity at poly(A) tracts

Among host-dependent bacteria that have evolved by extreme reductive genome evolution, long-term bacterial endosymbionts of insects have the smallest (160-790 kb) and most A + T-rich (>70%) bacterial genomes known to date. These genomes are riddled with poly(A) tracts, and 5-50% of genes contain tracts of 10 As or more. Here, we demonstrate transcriptional slippage at poly(A) tracts within genes of Buchnera aphidicola associated with aphids and Blochmannia pennsylvanicus associated with ants. Several tracts contain single frameshift deletions; these apparent pseudogenes showed patterns of constraint consistent with purifying selection on the encoded proteins. Transcriptional slippage yielded a heterogeneous population of transcripts with variable numbers of As in the tract. Across several frameshifted genes, including B. aphidicola cell wall biosynthesis genes and a B. pennsylvanicus histidine biosynthesis gene, 12-50% of transcripts contained corrected reading frames that could potentially yield full-length proteins. In situ immunostaining confirmed the production of the cell wall biosynthetic enzyme UDP-N-acetylmuramyl pentapeptide synthase encoded by the frameshifted murF gene. Simulation studies indicated an overrepresentation of poly(A) tracts in endosymbiont genomes relative to other A + T-rich bacterial genomes. Polymerase infidelity at poly(A) tracts rescues the functionality of genes with frameshift mutations and, conversely, reduces the efficiency of expression for in-frame genes carrying poly(A) regions. These features of homopolymeric tracts could be exploited to manipulate gene expression in small synthetic genomes.

Transcriptional slippage prompts recoding in alternate reading frames in the hepatitis C virus (HCV) core sequence from strain HCV-1

Since the first report of frameshifting in HCV-1, its sequence has been the paradigm for examining the mechanism that directs alternative translation of the hepatitis C virus (HCV) genome. The region encoding the core protein from this strain contains a cluster of 10 adenines at codons 8-11, which is thought to direct programmed ribosomal frameshifting (PRF), but formal evidence for this process has not been established unequivocally. To identify the mechanisms of frameshifting, this study used a bicistronic dual luciferase reporter system in a coupled transcription/translation in vitro assay. This approach revealed +1 as well as -1 frameshifting, whereas point mutations, selectively introduced between codons 8 and 11, demonstrated that PRF did not readily account for frameshifting in strain HCV-1. Sequence analysis of cDNAs derived from RNA transcribed by T7 RNA polymerase in the dual luciferase reporter system, as well as in both a subgenomic replicon and an infectious clone derived from strain JFH1, identified additions and deletions of adenines between codons 8 and 11 due to transcriptional slippage (TS). Moreover, RNA isolated from cells infected with virus generated by JFH1 containing the A-rich tract also contained heterogeneity in the adenine sequence, strongly suggesting TS by the NS5B viral polymerase. These findings have important implications for insight into frameshifting events in HCV-1 and demonstrate for the first time the involvement of transcriptional slippage in this recoding event.

Regulation of pyrimidine biosynthetic gene expression in bacteria: repression without repressors

SUMMARY

DNA-binding repressor proteins that govern transcription initiation in response to end products generally regulate bacterial biosynthetic genes, but this is rarely true for the pyrimidine biosynthetic (pyr) genes. Instead, bacterial pyr gene regulation generally involves mechanisms that rely only on regulatory sequences embedded in the leader region of the operon, which cause premature transcription termination or translation inhibition in response to nucleotide signals. Studies with Escherichia coli and Bacillus subtilis pyr genes reveal a variety of regulatory mechanisms. Transcription attenuation via UTP-sensitive coupled transcription and translation regulates expression of the pyrBI and pyrE operons in enteric bacteria, whereas nucleotide effects on binding of the PyrR protein to pyr mRNA attenuation sites control pyr operon expression in most gram-positive bacteria. Nucleotide-sensitive reiterative transcription underlies regulation of other pyr genes. With the E. coli pyrBI, carAB, codBA, and upp-uraA operons, UTP-sensitive reiterative transcription within the initially transcribed region (ITR) leads to nonproductive transcription initiation. CTP-sensitive reiterative transcription in the pyrG ITRs of gram-positive bacteria, which involves the addition of G residues, results in the formation of an antiterminator RNA hairpin and suppression of transcription attenuation. Some mechanisms involve regulation of translation rather than transcription. Expression of the pyrC and pyrD operons of enteric bacteria is controlled by nucleotide-sensitive transcription start switching that produces transcripts with different potentials for translation. In Mycobacterium smegmatis and other bacteria, PyrR modulates translation of pyr genes by binding to their ribosome binding site. Evidence supporting these conclusions, generalizations for other bacteria, and prospects for future research are presented.

Experimental evidence for negative selection in the evolution of a Yersinia pestis pseudogene

Yersinia pestis, the agent of bubonic plague, evolved from the enteric pathogen Yersinia pseudotuberculosis within the past 20,000 years. Because ancestor and descendant both exist, it is possible to infer steps in molecular evolution by direct experimental approaches. The Y. pestis life cycle includes establishment of a biofilm within its vector, the flea. Although Y. pseudotuberculosis makes biofilms in other environments, it fails to do so in the insect. We show that rcsA, a negative regulator of biofilms that is functional in Y. pseudotuberculosis, is a pseudogene in Y. pestis. Replacement of the pseudogene with the functional Y. pseudotuberculosis rcsA allele strongly represses biofilm formation and essentially abolishes flea biofilms. The conversion of rcsA to a pseudogene during Y. pestis evolution, therefore, was a case of negative selection rather than neutral genetic drift.

DNA polymerase kappa produces interrupted mutations and displays polar pausing within mononucleotide microsatellite sequences

Microsatellites are ubiquitously present in eukaryotic genomes and are implicated as positive factors in evolution. At the nucleotide level, microsatellites undergo slippage events that alter allele length and base changes that interrupt the repetitive tract. We examined DNA polymerase errors within a [T]₁₁ microsatellite using an in vitro assay that preferentially detects mutations other than unit changes. We observed that human DNA polymerase kappa (Pol κ) inserts dGMP and dCMP within the [T]₁₁ mononucleotide repeat, producing an interrupted 12-bp allele. Polymerase β produced such interruptions at a lower frequency. These data demonstrate that DNA polymerases are capable of directly producing base interruptions within microsatellites. At the molecular level, expanded microsatellites have been implicated in DNA replication fork stalling. Using an in vitro primer extension assay, we observed sequence-specific synthesis termination by DNA polymerases within mononucleotides. Quantitatively, intense, polar pausing was observed for both pol κ and polymerase α-primase within a [T]₁₁ allele. A mechanism is proposed in which pausing results from DNA bending within the duplex stem of the nascent DNA. Our data support the concept of a microsatellite life-cycle, and are consistent with the models in which DNA sequence or secondary structures contributes to non-uniform rates of replication fork progression.

Transcriptional slippage controls production of type III secretion apparatus components in Shigella flexneri

During transcription, series of approximately 9 As or Ts can direct RNA polymerase to incorporate into the mRNA nucleotides not encoded by the DNA, changing the reading frame downstream from the slippage site. We detected series of 9 or 10 As in spa13, spa33 and mxiA encoding type III secretion apparatus components. Analysis of cDNAs indicated that transcriptional slippage occurs in spa13, mxiA and spa33. Changes in the reading frame were confirmed by using plasmids carrying slippage sites in the 5' part of lacZ. Slippage is required for production of Spa13 from two overlapping reading frames and should lead to production of truncated MxiA and Spa33 proteins. Complementation of spa13 and mxiA mutants with plasmids carrying altered sites indicated that slippage in spa13 is required for assembly of the secretion apparatus and that slippage sites in spa13 and mxiA have not been selected to encode Lys residues or to produce two proteins endowed with different activities. The presence of slippage sites decreases production of Spa13 by 70%, of MxiA and Spa33 by 15% and of Spa32 (encoded downstream from spa13) by 50%. These results suggest that transcriptional slippage controls protein production by reducing the proportion of mRNA translated into functional proteins.

The rggC locus, with a frameshift mutation, is involved in oxidative stress response by Streptococcus thermophilus

In Streptococcus thermophilus, the locus rggC contains a frameshift mutation and thus consists of two open reading frames (ORFs), rggC (1) and rggC (2), which encode proteins exhibiting similarity with the Rgg transcriptional regulator family. In this work, mutants showing a partial deletion of rggC (1) and rggC (2 )were constructed and their response to menadione, a superoxide-generating compound, was analysed. These mutants exhibited different behaviour to this oxidative stress compared with the wild-type strain. Analysis of this locus among 21 strains of S. thermophilus showed a polythymine tract length variability and a strain-dependant adenine residue could be found upstream of this repeat. This interstrain polymorphism supports evidence for the hypothesis that the rggC locus is phase variable.

Transcriptional slippage in mxiE controls transcription and translation of the downstream mxiD gene, which encodes a component of the Shigella flexneri type III secretion apparatus

The Shigella flexneri transcription activator MxiE is produced by transcriptional slippage from two overlapping open reading frames. By using plasmids encoding a mxiD-lacZ translational fusion, we showed that transcriptional slippage in mxiE influences both transcription and translation of the downstream mxiD gene encoding an essential component of the type III secretion apparatus.

DNA sequences shaped by selection for stability

The sequence of a stretch of nucleotides affects its propensity for errors during replication and expression. Are proteins encoded by stable or unstable nucleotide sequences? If selection for variability is prevalent, one could expect an excess of unstable sequences. Alternatively, if selection against targets for errors were substantial, an excess of stable sequences would be expected. We screened the genome sequences of different organisms for an important determinant of stability, the presence of mononucleotide repeats. We find that codons are used to encode proteins in a way that avoids the emergence of mononucleotide repeats, and we can attribute this bias to selection rather than a neutral process. This indicates that selection for stability, rather than for the generation of variation, substantially influences how information is encoded in the genome.

Mutations are a double-edged sword. Most mutations are deleterious to an organism's fitness. On the other hand, without mutation, evolutionary change cannot occur. The rate of mutation is partially controlled by the organism, and one determinant of the mutation rate is the DNA sequence itself. Some DNA sequences are prone to mutations and errors during gene expression, whereas other sequences are more stable. Do organisms typically use stable or unstable DNA sequences in their genes? Both possibilities might seem plausible, and both have been postulated. To answer this question, the authors studied whether organisms' DNA sequences are more or less stable than expected by chance. Analyzing the genomes of a bacterium, a yeast, and a nematode, they find a overwhelming prevalence of stable DNA sequences, suggesting that selection for genetic stability is more important than selection for the generation of variation.

Vibrio cholerae ToxT independently activates the divergently transcribed aldA and tagA genes

The Vibrio cholerae ToxT regulon includes the genes encoding cholera toxin (CT) and the toxin-coregulated pilus (TCP), which are the major virulence factors required for causing cholera disease and colonizing the upper small intestine of the host, respectively. The genes encoding CT, ctxAB, and the genes encoding the components of the TCP, tcpA to tcpJ, are organized within operons, upstream of which are DNA binding sites for the transcriptional activator ToxT. ToxT is a member of the large AraC/XylS family of transcriptional regulators and also activates transcription of five other genes whose roles in V. cholerae pathogenesis, if any, are poorly understood. acfA and acfD are divergently transcribed genes required for efficient colonization of the intestine. Transcriptional activation of acfA and acfD requires a pair of central ToxT binding sites in an inverted-repeat configuration for ToxT-directed transcription of both genes. tcpI has an unknown role in pathogenesis. aldA and tagA are divergently transcribed genes that also have unknown roles in pathogenesis. In this study, we map the aldA and tagA promoters and identify the ToxT binding sites upstream of each gene. Our results suggest that two ToxT binding sites in an inverted-repeat configuration are required for ToxT-directed transcription of tagA and that a single ToxT binding site is required for ToxT-directed transcription of aldA. Furthermore, to direct transcription of tagA and aldA, ToxT uses independent binding regions upstream of each gene, in contrast to what we previously found for the divergently transcribed acfA and acfD genes, which share ToxT binding sites between the two genes.

Genome sequence of Blochmannia pennsylvanicus indicates parallel evolutionary trends among bacterial mutualists of insects

The distinct lifestyle of obligately intracellular bacteria can alter fundamental forces that drive and constrain genome change. In this study, sequencing the 792-kb genome of Blochmannia pennsylvanicus, an obligate endosymbiont of Camponotus pennsylvanicus, enabled us to trace evolutionary changes that occurred in the context of a bacterial-ant association. Comparison to the genome of Blochmannia floridanus reveals differential loss of genes involved in cofactor biosynthesis, the composition and structure of the cell wall and membrane, gene regulation, and DNA replication. However, the two Blochmannia species show complete conservation in the order and strand orientation of shared genes. This finding of extreme stasis in genome architecture, also reported previously for the aphid endosymbiont Buchnera, suggests that genome stability characterizes long-term bacterial mutualists of insects and constrains their evolutionary potential. Genome-wide analyses of protein divergences reveal 10- to 50-fold faster amino acid substitution rates in Blochmannia compared to related bacteria. Despite these varying features of genome evolution, a striking correlation in the relative divergences of proteins indicates parallel functional constraints on gene functions across ecologically distinct bacterial groups. Furthermore, the increased rates of amino acid substitution and gene loss in Blochmannia have occurred in a lineage-specific fashion, which may reflect life history differences of their ant hosts.

Frameshifting by transcriptional slippage is involved in production of MxiE, the transcription activator regulated by the activity of the type III secretion apparatus in Shigella flexneri

Bacteria of Shigella spp. are responsible for shigellosis in humans. They use a type III secretion (TTS) system encoded by a 200 kb virulence plasmid to enter epithelial cells and trigger apoptosis in macrophages. This TTS system comprises a secretion apparatus, translocators and effectors that transit through this apparatus, cytoplasmic chaperones and specific transcription regulators. The TTS apparatus assembled during growth of Shigella flexneri in broth is activated upon contact with epithelial cells. Transcription of approximately 15 genes encoding effectors, including IpaH proteins, is regulated by the TTS apparatus activity and controlled by MxiE, a transcription activator of the AraC family, and IpgC, the chaperone of the translocators IpaB and IpaC. We present evidence that MxiE is produced by a frameshift between a 59-codon open reading frame (ORF) (mxiEa) containing the translation start site and a 214-codon ORF (mxiEb) encoding the DNA binding domain of the protein. The mxiEa encoded N-terminal part of MxiE is required for MxiE function. Frameshifting efficiency was approximately 30% during growth in broth and was not modulated by the activity of secretion or the coactivator IpgC. Frameshifting involves slippage of RNA polymerase during transcription of mxiE, which results in the incorporation of one additional nucleotide in the mRNA and places mxiEa and mxiEb in the same reading frame. Frameshifting might represent an additional means of controlling gene expression under specific environmental conditions.

Molecular misreading: the frequency of dinucleotide deletions in neuronal mRNAs for beta-amyloid precursor protein and ubiquitin B

Human neuronal cells contain mutant beta-amyloid precursor protein (APP) and ubiquitin B (UBB) mRNAs, in which dinucleotide deletions ('Delta') are generated in/around GAGAG-motifs by an unknown mechanism referred to as 'Molecular Misreading.' The encoded frameshifted (+1) proteins accumulate in the neuropathological hallmarks of Alzheimer's disease (AD) and in other neurodegenerative and age-related diseases. To measure the concentration of Delta mRNAs, we developed a highly sensitive and specific assay, utilizing peptide nucleic acid-mediated PCR clamping, followed by cloning and colony hybridization with sequence-specific oligonucleotide probes. We found only a few molecules of Delta mRNA/microg of cellular RNA, at levels <10(-5) to 10(-6) x the concentration of WT mRNA, in RNA extracted from: (i) cultured human neuroblastoma cells grown under a variety of conditions, (ii) the frontal half of brains from wild type and XPA(-/-) DNA repair-deficient mice, and (iii) post-mortem temporal cortices from humans. Importantly, in RNA from the temporal cortices of AD and Down Syndrome patients that contain betaAPP+1 and UBB+1 immunoreactive cells, we found the same low levels of Delta mRNA. We infer that the accumulation of +1 proteins in neurons of these patients is not caused by an increase in the concentration of Delta mRNAs.

Transcriptional slippage in bacteria: distribution in sequenced genomes and utilization in IS element gene expression

To find a length of slippage-prone sequences at which selection against transcriptional slippage is evident, the transcription of repetitive runs of A and T of different lengths in 108 bacterial genomes was analyzed. IS element genes were found to exploit transcriptional slippage for regulation of gene expression.

Background

Transcription slippage occurs on certain patterns of repeat mononucleotides, resulting in synthesis of a heterogeneous population of mRNAs. Individual mRNA molecules within this population differ in the number of nucleotides they contain that are not specified by the template. When transcriptional slippage occurs in a coding sequence, translation of the resulting mRNAs yields more than one protein product. Except where the products of the resulting mRNAs have distinct functions, transcription slippage occurring in a coding region is expected to be disadvantageous. This probably leads to selection against most slippage-prone sequences in coding regions.

Results

To find a length at which such selection is evident, we analyzed the distribution of repetitive runs of A and T of different lengths in 108 bacterial genomes. This length varies significantly among different bacteria, but in a large proportion of available genomes corresponds to nine nucleotides. Comparative sequence analysis of these genomes was used to identify occurrences of 9A and 9T transcriptional slippage-prone sequences used for gene expression.

Conclusions

IS element genes are the largest group found to exploit this phenomenon. A number of genes with disrupted open reading frames (ORFs) have slippage-prone sequences at which transcriptional slippage would result in uninterrupted ORF restoration at the mRNA level. The ability of such genes to encode functional full-length protein products brings into question their annotation as pseudogenes and in these cases is pertinent to the significance of the term 'authentic frameshift' frequently assigned to such genes.

Genome-based peptide fingerprint scanning

We have implemented a method that identifies the genomic origins of sample proteins by scanning their peptide-mass fingerprint against the theoretical translation and proteolytic digest of an entire genome. Unlike previously reported techniques, this method requires no predefined ORF or protein annotations. Fixed-size windows along the genome sequence are scored by an equation accounting for the number of matching peptides, the number of missed enzymatic cleavages in each peptide, the number of in-frame stop codons within a window, the adjacency between peptides, and duplicate peptide matches. Statistical significance of matching regions is assessed by comparing their scores to scores from windows matching randomly generated mass data. Tests with samples from Saccharomyces cerevisiae mitochondria and Escherichia coli have demonstrated the ability to produce statistically significant identifications, agreeing with two commonly used programs, peptident and mascot, in 86% of samples analyzed. This genome fingerprint scanning method has the potential to aid in genome annotation, identify proteins for which annotation is incorrect or missing, and handle cases where sequencing errors have caused framing mistakes in the databases. It might also aid in the identification of proteins in which recoding events such as frameshifting or stop-codon read-through have occurred, elucidating alternative translation mechanisms. The prototype is implemented as a clientserver pair, allowing the distribution, among a set of cluster nodes, of a single or multiple genomes for concurrent analysis.

Factors affecting start site selection at the Escherichia coli fis promoter

Transcription initiation with CTP is an uncommon feature among Escherichia coli sigma(70) promoters. The fis promoter (fis P), which is subject to growth phase-dependent regulation, is among the few that predominantly initiate transcription with CTP. Mutations in this promoter that cause a switch from utilization of CTP to either ATP or GTP as the initiation nucleotide drastically alter its growth phase regulation pattern, suggesting that the choice of the primary initiating nucleotide can significantly affect its regulation. To better understand what factors influence this choice in fis P, we made use of a series of promoter mutations that altered the nucleotide or position used for initiation. Examination of these promoters indicates that start site selection is determined by a combination of factors that include preference for a nucleotide distance from the -10 region (8 > 7 > 9 >> 6 >> 10 > 11), initiation nucleotide preference (A = G >> CTP > or = UTP), the DNA sequence surrounding the initiation region, the position of the -35 region, and changes in the intracellular nucleoside triphosphate pools. We describe the effects that each of these factors has on start site selection in the fis P and discuss the interplay between position and nucleotide preference in this important process.

Transcription and nuclear transport of CAG/CTG trinucleotide repeats in yeast

Trinucleotide repeats are involved in several neurological disorders in humans. DNA sequences containing CAG/CTG repeats are prone to slippage during replication and double-strand break repair. The effects of trinucleotide repeats on transcription and on nuclear export were analyzed in vivo in yeast. Transcription of a CAG/CTG trinucleotide repeat in the 3'-untranslated region of a URA3 reporter gene leads to transcription of messenger RNAs several kilobases longer than the expected size. These long mRNAs form more readily when CAG rather than CTG repeats are transcribed. CAG- or CUG-containing transcripts show a non-homogeneous cellular localization. We propose that long mRNAs result from transcription slippage, and discuss the possible implications for human diseases.

Local variability of the phosphoglycerate kinase-triosephosphate isomerase fusion protein from Thermotoga maritima MSB 8

The pgk-tpi gene locus of Thermotoga maritima encodes both phosphoglycerate kinase (PGK) and a bienzyme complex consisting of a fusion protein of PGK with triosephosphate isomerase (TIM). No separate tpi gene for TIM is present in T. maritima. A frame-shift at the end of the pgk gene has been previously proposed as a mechanism to regulate the expression of the two protein variants [Schurig et al., EMBO J. 14 (1995), 442-451]. Surprisingly, the complete T. maritima genome was found to contain a pgk-tpi sequence not requiring the proposed frameshift mechanism. To clarify the apparent discrepancy, a variety of DNA sequencing techniques were applied, disclosing an anomalous local variability in the pgk-tpi fusion region. The comparison of different DNA samples and the mass spectrometric analysis of the amino acid sequence of the natural fusion protein from T. maritima MSB8 confirmed the local variability of the DNA variants. Since not all peptide masses could be assigned, further variations are conceivable, suggesting an even higher heterogeneity of the T. maritima MSB8 strain.

Construction of a multiframe vector to express coding sequences inEscherichia coli

Cloning of foreign DNA fragments for coding sequence analysis in Escherichia coli usually involves sets of three vectors. To simplify this, we constructed an expression vector named pMFV7 containing three ATG codons in different frames downstream of a Shine-Dalgarno sequence, assuming that the ribosome can use any of the three start codons in an alternative manner. Translation beginning at either of the start codons would drive the expression of any coding fragment cloned downstream. To test the feasibility of this proposal, we cloned DNA fragments of the lacZ gene in each of the possible reading frames downstream from pMFV7 start codons. Sequence analysis of the N-terminus regions around the fusion sites indicates that ribosomes indeed initiate translation at each of the three initiation codons. In one case, levels of β-galactosidase activity depended largely on the N-terminus of the translation products. We conclude that pMFV7 may be useful for expressing coding sequences regardless of their reading frame.Key words: translation initiation, in-frame gene cloning, expression vectors.

Transcription in Mycoplasma pneumoniae

Very little is understood of the structure of mycoplasma promoters, and this limits interpretation of genomic sequence data in these species. In this study the transcriptional start points of 22 genes of Mycoplasma pneumoniae were identified and the regions 5' to the start point compared. Although a strong consensus -10 region could be seen, there was only a weak consensus in the -35 region. A high proportion of transcripts had heterogeneous 5'-ends and characterisation of the sequence of the 5'-ends of two transcripts established that the heterogeneity was derived from initiation of transcription at reduced levels between 1 and 4 bases 5' to the major starting point. In addition to this apparently unique feature, a high proportion of transcripts lacked a 5' untranslated leader region that could contain a ribosomal binding site. Such leaderless transcripts are seen rarely in other bacterial species. Although the promoter regions for a number of members of lipoprotein multigene families were examined, no obvious explanation for regulation of expression was apparent. Using the data from this study an improved matrix for prediction of M.pneumoniae promoters was derived. Application of this matrix to the sequences immediately 3' and 5' to each predicted start codon in the genome suggested that most M. pneumoniae transcriptional start points were likely to occur between 5 and 30 bases 5' to the start codon.

Nonlinearity in genetic decoding: homologous DNA replicase genes use alternatives of transcriptional slippage or translational frameshifting

The tau and gamma subunits of DNA polymerase III are both encoded by a single gene in Escherichia coli and Thermus thermophilus. gamma is two-thirds the size of tau and shares virtually all its amino acid sequence with tau. E. coli and T. thermophilus have evolved very different mechanisms for setting the approximate 1:1 ratio between tau and gamma. Both mechanisms put ribosomes into alternate reading frames so that stop codons in the new frame serve to make the smaller gamma protein. In E. coli, approximately 50% of initiating ribosomes translate the dnaX mRNA conventionally to give tau, but the other 50% shift into the -1 reading frame at a specific site (A AAA AAG) in the mRNA to produce gamma. In T. thermophilus ribosomal frameshifting is not required: the dnaX mRNA is a heterogeneous population of molecules with different numbers of A residues arising from transcriptional slippage on a run of nine T residues in the DNA template. Translation of the subpopulation containing nine As (or +/- multiples of three As) yields tau. The rest of the population of mRNAs (containing nine +/- nonmultiples of three As) puts ribosomes into the alternate reading frames to produce the gamma protein(s). It is surprising that two rather similar dnaX sequences in E. coli and T. thermophilus lead to very different mechanisms of expression.

Transcriptional slippage of p53 gene enhanced by cellular damage in rat liver: monitoring the slippage by a yeast functional assay

The Long-Evans Cinnamon (LEC) rat is a mutant strain characterized by abnormal copper metabolism and a high incidence of hepatitis and hepatoma. Using a yeast-based assay which scores mutants in p53 gene transcripts as red colonies, we detected frequent mutations in the liver of LEC rats. The majority (50-60%) of these were frameshift mutations caused by the insertion of an extra adenine (A) in the regions containing six consecutive adenines. The rate of A insertion was calculated to be 6.9-9.0% of the total p53 cDNA. Insertions of an extra adenine were found almost exclusively in the mRNA (cDNA), especially in the (A)(6) tract located at the most 5'-side (exon 4) among the three (A)(6) tracts (exons 4, 7, and 8), but rarely in the corresponding sites of genomic DNA. Wild-type p53 cDNA was transcribed in vitro into mRNA with the use of SP6 RNA polymerase and tested by the yeast functional assay. Subsequent sequencing detected A insertions at an overall rate of 1.6% in exons 7 and 8 but none in exon 4. This indicates that the A insertion in the exon 4 (A)(6) tract was an in vivo phenomenon rather than an artifact in reverse transcription or polymerase chain reaction. The percentage of red colonies increased sharply to about 20% of the liver samples in the acute hepatitis stage, and returned to control level of those in the chronic hepatitis stage, and increased again slightly to those in the neoplastic stage. The percentage of red colonies correlated with the serum GOT level (r=0.96, p<0.001) but not with the contents of copper and 8-hydroxydeoxyguanosine in the liver of LEC rats. Ethanol treatment of hepatic cell lines also increased the rate of transcriptional slippage at the (A)(6) tract. These findings indicate that cellular damage is responsible for the increase in the rate of mutation at the transcriptional level, and suggest that cellular damage degrades transcriptional fidelity, thereby further impairing cellular functions.

A novel A-kinase anchoring protein in the heart interacts with G alpha 13

A cDNA of a tentative A-kinase anchoring protein, presumably coupled with heterotrimeric GTP binding protein alpha 13 subunit (G alpha 13), was cloned from a human heart cDNA library. It was approximately 650 bases and its mRNA was expressed in the heart. Homology search of DNA sequences revealed that it was a novel cDNA with 84% homology with the partial sequence of rabbit cDNA of AKAP 120 without a stop codon. 3'-Rapid Amplification of cDNA Ends (3'-RACE) and yeast functional assay were performed to determine the 3'-end of the cDNA and ribosomal frameshifting was suggested as a translational mechanism. Here we report that a protein encoded by the cDNA may be involved in intracellular signal transduction via the G alpha 13 and PKA in hearts.

Polymorphism and disorder of poly(dA).poly(dT) in fibers

The absence of proportionality between fiber length and interbase axial distance reveals a high degree of disorder in fibers of poly(dA).poly(dT) when the relative humidity (r.h.) varies. According to X-ray patterns, the increase of fiber length does not correspond to any conformational change. The B' structure is present alone and in a better-organized form at higher r.h. We noted a dependence of length values on the direction of variation of the r.h. An unusual minimum of fiber length is observed around 65% r.h. and explained as been due to an amorphous part of the fiber, essentially composed of single strands of poly(dA) and poly(dT). At a temperature higher than 30 degrees C, the transition B'(alpha)-B* takes place but is associated with an important decrease of fiber length in spite of constancy of the helical rise per base pair. It is shown that this behavior is also due to the amorphous part of the fiber and not to a conformational transition.

High-level expression of branching enzyme II from maize endosperm in Escherichia coli

The gene that encodes the mature branching enzyme II (BEII) protein from maize (Zea mays L.) endosperm was amplified by means of a polymerase chain reaction technique and inserted into a T7-based expression vector. Although this has been an efficient expression system of maize BEII in Escherichia coli, an example is presented in this report which allows a greater expression of mBEII protein from the bacterial system by changing only one codon. The key to the level of expression appears to be related to the conversion of the third thymine base in the 285 position codon of the mBEII cDNA to cytosine without altering the encoded mBEII protein product. The crude cell extracts of enzyme prepared from E.coli exhibited seven-fold higher expression of branching enzyme activity compared to expression of the native enzyme. The enzymes from wild-type and the silent mutation genes were purified. The proteins were indistinguishable kinetically and immunologically. Thus, we obtained a significantly improved expression of mBEII protein in the bacterial system.

tRNA imbalance promotes -1 frameshifting via near-cognate decoding

tRNAGly1 is the Escherichia coli glycine tRNA specific for GGG codons. A genetic selection for multicopy suppressors of a frameshift mutation has shown that increased levels of wild-type tRNAGly1 causes -1 frameshifting. Analysis of the suppression spectrum of this multicopy suppressor and peptide sequencing of the suppressed protein product showed that it promoted GG doublet decoding at the near-cognate GGA codons. It is proposed that increasing the concentration of the GGG-specific tRNAGly1 relative to the cognate GGA-decoding tRNAGly2 allows the near-cognate tRNA to read GGA codons. Near-cognate decoding of GGA codons by tRNAGly1 can occur by a two-out-of-three reading mechanism, in which only the first two bases of the GGA codon are paired with the anticodon, thus permitting doublet translocations. In mycoplasmas, a single tRNA typically decodes all four triplets of a codon family and introduction of a feature of the Mypoplasma mycoides tRNAGly responsible for non-discriminate decoding, a C at position 32, into the anticodon E. coli tRNAGly1, enhanced the efficiency of doublet decoding.

Growth phase dependent stop codon readthrough and shift of translation reading frame in Escherichia coli

Nonsense codon readthrough and changed translational reading frame were measured in different growth phases in E. coli. The strains used carry plasmid constructs with a translation assay reporter gene. This reporter gene contains an internal stop codon or a run of U-residues. Termination or frameshifting give rise to stable proteins that can be physically quantified on gels along with the complete protein products. Readthrough of the stop codon UGA by a nearcognate tRNA is several fold higher in active growth than in late exponential phase. In early exponential phase, about 7% of -1 frameshift at a U9 slippery sequence is detectable; upon entry to stationary phase this frameshifting increases to about 40% followed by a decrease in stationary phase. A similar increase is observed in the case of +1 reading frameshift at the U9 sequence, which increases from 13% in early exponential growth phase up to 38% at the beginning of stationary phase followed by a decrease. Thus, the levels of both stop codon readthrough and frameshifting are growth phase dependent, though not in an identical fashion.

Loss of intrinsic aminoglycoside resistance in Acinetobacter haemolyticus as a result of three distinct types of alterations in the aac(6')-Ig gene, including insertion of IS17

The distribution of the aac(6')-Ig gene, encoding aminoglycoside 6'-N-acetyltransferase-Ig [AAC(6')-Ig], was studied in 96 Acinetobacter haemolyticus strains and 12 proteolytic Acinetobacter strains, including Acinetobacter genomospecies 6, 13, and 14 and 3 unnamed species assigned to this genomic group by DNA-DNA hybridization. This gene was detected by DNA-DNA hybridization in all 96 A. haemolyticus strains and by PCR in 95 strains but was not detected in strains of other species, indicating that it may be used to identify A. haemolyticus. Three A. haemolyticus strains were susceptible to tobramycin and did not produce an aminoglycoside 6'-N-acetylating activity, although they contained aac(6')-Ig-related sequences. An analysis of three susceptible A. haemolyticus strains indicated that aminoglycoside resistance was abolished by the following three distinct mechanisms: (i) a point mutation in aac(6')-Ig that led to a Met56-->Arg substitution, which was shown by analysis of a revertant to be responsible for the loss of resistance; (ii) a polythymine insertion that altered the reading frame; and (iii) insertion of IS17, a new member of the IS903 family. These observations indicated that AAC(6')-Ig is not essential for the viability of A. haemolyticus, although the aac(6')-Ig gene was detected in all members of this species.

Regulation of upp expression in Escherichia coli by UTP-sensitive selection of transcriptional start sites coupled with UTP-dependent reiterative transcription

Expression of the upp gene of Escherichia coli, which encodes the pyrimidine salvage enzyme uracil phosphoribosyltransferase, is negatively regulated by pyrimidine availability. In this study, we demonstrate that this regulation occurs mainly by UTP-sensitive selection of alternative transcriptional start sites, which produces transcripts that differ in the ability to be productively elongated. The upp initially transcribed region contains the sequence GATTTTTTTTG (nontemplate strand). Transcription is initiated primarily at the first two bases in this sequence, designated G6 and A7 (counting from the promoter -10 region). High intracellular levels of UTP favor initiation at position A7; however, the resulting transcripts are subject to reiterative transcription (i.e., repetitive nucleotide addition) within the run of T residues in the initially transcribed region. The resulting AUUUUn (where n = 1 to >50) transcripts are not extended to include downstream upp sequences. In contrast, low intracellular levels of UTP strongly favor initiation at position G6, which results in transcripts that generally do not engage in reiterative transcription and thus can be normally elongated. This mechanism ensures that high levels of uracil phosphoribosyltransferase are produced only under conditions of pyrimidine limitation. The mechanisms that account for UTP-sensitive start site selection and different fates of upp transcripts, as well as the general use of UTP-dependent reiterative transcription in gene regulation, are discussed in detail.

Reading-frame restoration by transcriptional slippage at long stretches of adenine residues in mammalian cells

We previously characterized a mutant apoB allele (the apoB86 allele) that produces both a truncated apoB (apoB86) and a full-length apoB100. The mutant allele contained a deletion of a single cytosine in exon 26, creating a stretch of eight consecutive adenines in the -1 reading frame. The altered reading-frame allele was restored, with approximately 10% efficiency, by the transcriptional insertion of an extra adenine into the stretch of eight consecutive adenines, thereby accounting for the synthesis of the full-length apoB100. Here, we demonstrate that this reading-frame restoration does not occur when the long stretch of adenines is interrupted by a cytosine. To assess whether reading-frame restoration is unique to a single site in the apoB gene, the same mutation (eight consecutive adenines in the -1 reading frame) was inserted into another site within the apoB gene. Reading-frame restoration occurred at the second site and was abrogated when the stretch of adenines was interrupted by another base. Of note, a computerized analysis of human cDNA sequences revealed that long stretches of adenines in protein-coding sequences occur at a lower than predicted frequency, suggesting that evolution has selected against these sequences.

Attenuation of gene expression by a trinucleotide repeat-rich tract from the terminal exon of the rat hepatic polymeric immunoglobulin receptor gene

A 359 bp terminal exon fragment of the rat polymeric immunoglobulin receptor gene has been tested for biological effects. The fragment contains an S1 nuclease-sensitive microsatellite with d(GGA) and d(GAA) trinucleotide repeats that are expressed discordantly in the 3'UTRs of liver mRNAs encoded by the single copy gene. When human A293 cells are transfected with expression plasmids carrying this fragment in forward orientations, flanking or replacing poly(A) cassettes in the 3' ends of the transcription units, luciferase reporter gene expression is attenuated 47 to 59% or 98.5%, respectively. In contrast, when the fragment is tested similarly in reverse orientation, there is significantly less or no attenuation of gene expression. These observations, and computer models of partial triplet repeat DNA tertiary and RNA secondary structures, suggest that this fragment might regulate gene expression by orientation and position-dependent mechanisms at transcriptional and post-transcriptional levels.

Partial correction of a severe molecular defect in hemophilia A, because of errors during expression of the factor VIII gene

Although the molecular defect in patients in a Japanese family with mild to moderately severe hemophilia A was a deletion of a single nucleotide T within an A8TA2 sequence of exon 14 of the factor VIII gene, the severity of the clinical phenotype did not correspond to that expected of a frameshift mutation. A small amount of functional factor VIII protein was detected in the patient's plasma. Analysis of DNA and RNA molecules from normal and affected individuals and in vitro transcription/translation suggested a partial correction of the molecular defect, because of the following: (i) DNA replication/RNA transcription errors resulting in restoration of the reading frame and/or (ii) "ribosomal frameshifting" resulting in the production of normal factor VIII polypeptide and, thus, in a milder than expected hemophilia A. All of these mechanisms probably were promoted by the longer run of adenines, A10 instead of A8TA2, after the delT. Errors in the complex steps of gene expression therefore may partially correct a severe frameshift defect and ameliorate an expected severe phenotype.

Basic mechanisms of transcript elongation and its regulation

Ternary complexes of DNA-dependent RNA polymerase with its DNA template and nascent transcript are central intermediates in transcription. In recent years, several unusual biochemical reactions have been discovered that affect the progression of RNA polymerase in ternary complexes through various transcription units. These reactions can be signaled intrinsically, by nucleic acid sequences and the RNA polymerase, or extrinsically, by protein or other regulatory factors. These factors can affect any of these processes, including promoter proximal and promoter distal pausing in both prokaryotes and eukaryotes, and therefore play a central role in regulation of gene expression. In eukaryotic systems, at least two of these factors appear to be related to cellular transformation and human cancers. New models for the structure of ternary complexes, and for the mechanism by which they move along DNA, provide plausible explanations for novel biochemical reactions that have been observed. These models predict that RNA polymerase moves along DNA without the constant possibility of dissociation and consequent termination. A further prediction of these models is that the polymerase can move in a discontinuous or inchworm-like manner. Many direct predictions of these models have been confirmed. However, one feature of RNA chain elongation not predicted by the model is that the DNA sequence can determine whether the enzyme moves discontinuously or monotonically. In at least two cases, the encounter between the RNA polymerase and a DNA block to elongation appears to specifically induce a discontinuous mode of synthesis. These findings provide important new insights into the RNA chain elongation process and offer the prospect of understanding many significant biological regulatory systems at the molecular level.

In vitro transcription of a poly(dA) x poly(dT)-containing sequence is inhibited by interaction between the template and its transcripts

Transcription of poly(dA) x poly(dT)-containing sequences was investigated in vitro using plasmids carrying a (dA)34 x (dT)34 tract in the coding region of the lacZ gene. The efficiency of transcription of the (dT)34 sequence on the transcribing strand by Escherichia coli RNA polymerase was substantially lower (approximately 60%) than that of the (dA)34 sequence or of the control lacZ gene. Analysis of the transcription process of the (dT)34 sequence by T3 RNA polymerase showed that the transcription was frequently arrested or terminated at the middle as well as immediately proximal of the (dA)34 x (dT)34 tract, and it occurred more prominently following accumulation of transcription products. This inhibition was strongly enhanced by the addition of the oligonucleotide (dT)34 or poly(U) to the reaction mixture, while (dA)34 and the duplex (dA)34 x (dT)34 suppressed the inhibition. A similar transcriptional inhibition was also observed in transcription mediated by T7 RNA polymerase and eukaryotic RNA polymerase II. We also demonstrated RNA x DNA complex formation of the (dA)34 x (dT)34 tract with poly(U), but not with poly(A). These findings strongly suggest that poly(dT)-containing template sequences interact and form a complex with its transcription products, possibly an RNA x DNA triplex, which blocks further transcription. This would explain the instability of the plasmids transcribing mRNAs with poly(U) but not poly(A) tracts and the underrepresentation of poly(U) but not poly(A) tracts in mRNAs.

The yeast transcription terminator for RNA polymerase I is designed to prevent polymerase slippage

A transcription terminator for RNA polymerase I (polI) in the yeast, Saccharomyces cerevisiae, is composed of two essential elements, the 11bp binding site for Reb1p and an upstream T-rich element coding for the last 10-12 nucleotides of the terminated transcript. We now show that, if the upstream element is changed to homopolymer T residues, polI undergoes iterative slippage, long poly(U) tails are added to the transcript, and termination is impaired. Reinsertion of one or two non-T residues within a critical region prevents iterative slippage and reinstates termination. A survey of naturally occurring terminators reveals that many contain T-rich upstream regions with non-T residues situated appropriately to prevent slippage. We discuss the possibility that the first step in slippage, backward sliding of both the transcript and the catalytic center of the polymerase, may be an obligatory step in the normal termination process.

Phylogeny of the Mycoplasma mycoides cluster as determined by sequence analysis of the 16S rRNA genes from the two rRNA operons

The so-called Mycoplasma mycoides cluster consists of six species or subspecies of mycoplasmas (Mollicutes). These species are pathogenic for ruminants and some of them are of great concern in veterinary medicine. The members of the M. mycoides cluster have two rRNA operons (rrnA and rrnB). The nucleotide sequences of the 16S rRNA genes of 10 strains, representing all of the known species and subspecies of the M. mycoides cluster, were determined by direct automated solid-phase DNA sequencing. The sequences of both rRNA operons were determined by a novel strategy involving in vitro amplification by PCR with one operon-specific primer pair and one general primer pair. Interestingly, sequence differences (polymorphisms) between the two operons were observed for all strains. Two strains of M. capricolum subsp. capripneumoniae were sequenced, and 15 polymorphisms were found in the type strain (F38) and 17 polymorphisms were found in the other strain (4/2LC). Eight polymorphisms were found in the 16S rRNA genes of the M. mycoides subsp. mycoides small-colony type, and sequence length variations in a poly(A) region were observed in the 16S rRNA genes of the two operons of this species. Secondary-structure analysis showed that polymorphisms were present in both stem and loop regions. The nucleotide substitutions in the polymorphic sites of the stem regions often resulted in a change from a canonical to a noncanonical base pairing or vice versa. A compensatory mutation was never observed in the other nucleotide of the base pair. Phylogenetic analysis based on the 16S rRNA sequences indicated that Mycoplasma sp. strain PG50 should be included in the M. capricolum species group. Furthermore, the 16S rRNA sequences of M. mycoides subsp. capri and the M. mycoides subsp. mycoides large-colony type were 99.9% identical. We therefore suggest that these species be reclassified in a common species group (for instance, "Mycoplasma capri") distinct from the M. mycoides subsp. mycoides small-colony type, which formed an intermediate branch between the M. capricolum species group and the M. capri species group.