DNA sequences from the str operon of Escherichia coli

Fitness benefits of a synonymous substitution in an ancient EF-Tu gene depend on the genetic background

Synonymous mutations are changes to DNA sequence, which occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation elongation factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels as well as global polysome abundance on RNA transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.IMPORTANCEThis study explores the degree to which synonymous mutations in essential genes can influence adaptation in bacteria. An experimental system whereby an Escherichia coli strain harboring an engineered translation protein elongation factor-Tu (EF-Tu) was subjected to laboratory evolution. We find that a synonymous mutation acquired on the gene encoding for EF-Tu is conditionally beneficial for bacterial fitness. Our findings provide insight into the importance of the genetic background when a synonymous substitution is favored by natural selection and how such changes have the potential to impact evolution when critical cellular processes are involved.

No Evidence that Selection on Synonymous Codon Usage Affects Patterns of Protein Evolution in Bacteria

Bias in synonymous codon usage has been reported across all kingdoms of life. Evidence suggests that codon usage bias is often driven by selective pressures, typically for translational efficiency. These selective pressures have been shown to depress the rate at which synonymous sites evolve. We hypothesize that selection on synonymous codon use could also slow the rate of protein evolution if a non-synonymous mutation changes the codon from being preferred to unpreferred. We test this hypothesis by looking at patterns of protein evolution using polymorphism and substitution data in two bacterial species, Escherichia coli and Streptococcus pneumoniae. We find no evidence that non-synonymous mutations that change a codon from being unpreferred to preferred are more common than the opposite. Overall, selection on codon bias seems to have little influence over non-synonymous polymorphism or substitution patterns.

A beneficial synonymous substitution in EF-Tu is contingent on genetic background

Synonymous mutations are changes to DNA sequence that occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation Elongation Factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels, as well as the polysome abundance on global transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.

Enhanced effective codon numbers to understand codon usage bias

Codon usage bias is a well recognized phenomenon but the relative influence of its major causes: G+C content, mutational biases, and selection, are often difficult to disentangle. This paper presents methods to calculate modified effective codon numbers that allow the investigation of the sources of codon bias and how genes or organisms have their codon biases shaped. In particular, it demonstrates that variation in codon usage bias across organisms is likely driven more by likely mutational forces while the variation in codon usage bias within genomes is likely driven by codon selectional forces.

Effects of synonymous mutations beyond codon bias: The evidence for adaptive synonymous substitutions from microbial evolution experiments

Synonymous mutations are often assumed to be neutral with respect to fitness because they do not alter the encoded amino acid and so cannot be 'seen' by natural selection. Yet a growing body of evidence suggests that synonymous mutations can have fitness effects that drive adaptive evolution through their impacts on gene expression and protein folding. Here, we review what microbial experiments have taught us about the contribution of synonymous mutations to adaptation. A survey of site-directed mutagenesis experiments reveals the distributions of fitness effects for nonsynonymous and synonymous mutations are more similar, especially for beneficial mutations, than expected if all synonymous mutations were neutral, suggesting they should drive adaptive evolution more often than is typically observed. A review of experimental evolution studies where synonymous mutations have contributed to adaptation shows they can impact fitness through a range of mechanisms including the creation of illicit RNA polymerase binding sites impacting transcription and changes to mRNA folding stability that modulate translation. We suggest that clonal interference in evolving microbial populations may be the reason synonymous mutations play a smaller role in adaptive evolution than expected based on their observed fitness effects. We finish by discussing the impacts of falsely assuming synonymous mutations are neutral and discuss directions for future work exploring the role of synonymous mutations in adaptive evolution.

Dual-Seq reveals genome and transcriptome of Caedibacter taeniospiralis, obligate endosymbiont of Paramecium

Interest in host-symbiont interactions is continuously increasing, not only due to the growing recognition of the importance of microbiomes. Starting with the detection and description of novel symbionts, attention moves to the molecular consequences and innovations of symbioses. However, molecular analysis requires genomic data which is difficult to obtain from obligate intracellular and uncultivated bacteria. We report the identification of the Caedibacter genome, an obligate symbiont of the ciliate Paramecium. The infection does not only confer the host with the ability to kill other cells but also renders them immune against this effect. We obtained the C. taeniospiralis genome and transcriptome by dual-Seq of DNA and RNA from infected paramecia. Comparison of codon usage and expression level indicates that genes necessary for a specific trait of this symbiosis, i.e. the delivery of an unknown toxin, result from horizontal gene transfer hinting to the relevance of DNA transfer for acquiring new characters. Prediction of secreted proteins of Caedibacter as major agents of contact with the host implies, next to several toxin candidates, a rather uncharacterized secretome which appears to be highly adapted to this symbiosis. Our data provides new insights into the molecular establishment and evolution of this obligate symbiosis and for the pathway characterization of toxicity and immunity.

"Lost in translation: Seeing the forest by focusing on the trees"

A complex process translates messenger RNA (mRNA) base sequence into protein amino acid sequence. Transfer RNAs must recognize 3-base codons in the mRNA to insert the correct amino acids into the growing protein. Codon degeneracy makes decoding complicated in that multiple (synonymous) triplets can encode a single amino acid and multiple tRNAs can have the same anticodon. Over the last twenty years, new developments in structural biology, genome sequencing and bioinformatics has elucidated the intricacies of the ribosome structure and the details of the translation process. High throughput analyses of sequence information support the idea that mRNA folding has a major effect on expression for codons at the 5'-end of mRNA (N-terminal region of a polypeptide). Despite a forest of sequence data, significant details of the complex translation process can escape detection. However, a sensitive translation assay has allowed a single tree in this forest to be revealing.

Compositional Constraint Is the Key Force in Shaping Codon Usage Bias in Hemagglutinin Gene in H1N1 Subtype of Influenza A Virus

It is vital to unravel the codon usage bias in order to gain insights into the evolutionary forces dictating the viral evolution process. Influenza A virus has attracted attention of many investigators over the years due to high mutation rate and being cross-specific shift operational in the viral genome. Several authors have reported that the codon usage bias is low in influenza A viruses, citing mutational pressure as the decisive force shaping up the codon usage in these viruses. In this study, complete coding sequences of hemagglutinin genes for H1N1 subtype of influenza A virus have been explored for the possible codon usage bias acting upon these genes. The results indicate overall low bias with peaking ENC values. The GC content is found to be substantially low as against AT content in the silent codon sites. Significant correlations were observed in between the compositional parameters versus AT₃, implying the possible role of the latter in shaping codon usage profile in the viral hemagglutinin. The data showed conspicuously that the sequences were A redundant with most codons preferring nucleotide A over others in the third synonymous codon site. The results indicated the pivotal role of compositional pressure affecting codon usage in this virus.

The effects of codon context on in vivo translation speed

We developed a bacterial genetic system based on translation of the his operon leader peptide gene to determine the relative speed at which the ribosome reads single or multiple codons in vivo. Low frequency effects of so-called "silent" codon changes and codon neighbor (context) effects could be measured using this assay. An advantage of this system is that translation speed is unaffected by the primary sequence of the His leader peptide. We show that the apparent speed at which ribosomes translate synonymous codons can vary substantially even for synonymous codons read by the same tRNA species. Assaying translation through codon pairs for the 5'- and 3'- side positioning of the 64 codons relative to a specific codon revealed that the codon-pair orientation significantly affected in vivo translation speed. Codon pairs with rare arginine codons and successive proline codons were among the slowest codon pairs translated in vivo. This system allowed us to determine the effects of different factors on in vivo translation speed including Shine-Dalgarno sequence, rate of dipeptide bond formation, codon context, and charged tRNA levels.

Translation efficiency in yeasts correlates with nucleosome formation in promoters

Elongation efficiency index (EEI) was suggested earlier to estimate gene expression efficiency by nucleotide context of coding sequence in unicellular organisms. We have analyzed association between EEI and nucleosome formation potential (NFP) in 5' regulatory regions upstream translation initiation site (TIS) from two yeast species. Theoretical estimations of NFP based on DNA sequence were obtained by Recon method. Experimental estimation of nucleosome occupancy was obtained by high-throughput sequencing data of nucleosomal DNA in Saccharomyces cerevisiae . For the sample of all genes correlation coefficient was calculated between two vectors: vector of NFP values for fixed position relative to TIS and vector of EEI values. Profiles of correlation coefficients of NFP and EEI were counted in (-600; +600) regions relative to TIS for gene sequences extracted from GenBank. We found regions of strong negative dependence between NFP and EEI for all genes as well as for 10% highly expressed genes in Schizosaccharomyces pombe (10% of EEI-highest genes). At the same time, we found positive dependence between NFP and EEI for all genes and for low expressed genes in S. cerevisiae (10% of EEI-lowest genes). The association between NFP and EEI could be explained by evolutionary selection of context characteristics of nucleotide sequences for gene expression optimization.

Variations in Helicobacter pylori cytotoxin-associated genes and their influence in progression to gastric cancer: implications for prevention

Helicobacter pylori (HP) is a bacterium that colonizes the human stomach and can establish a long-term infection of the gastric mucosa. Persistent Hp infection often induces gastritis and is associated with the development of peptic ulcer disease, atrophic gastritis, and gastric adenocarcinoma. Virulent HP isolates harbor the cag (cytotoxin-associated genes) pathogenicity island (cagPAI), a 40 kb stretch of DNA that encodes components of a type IV secretion system (T4SS). This T4SS forms a pilus for the injection of virulence factors into host target cells, such as the CagA oncoprotein. We analyzed the genetic variability in cagA and other selected genes of the HP cagPAI (cagC, cagE, cagL, cagT, cagV and cag Gamma) using DNA extracted from frozen gastric biopsies or from clinical isolates. Study subjects were 95 cagA+ patients that were histologically diagnosed with chronic gastritis or gastric cancer in Venezuela and Mexico, areas with high prevalence of Hp infection. Sequencing reactions were carried out by both Sanger and next-generation pyrosequencing (454 Roche) methods. We found a total of 381 variants with unambiguous calls observed in at least 10% of the originally tested samples and reference strains. We compared the frequencies of these genetic variants between gastric cancer and chronic gastritis cases. Twenty-six SNPs (11 non-synonymous and 14 synonymous) showed statistically significant differences (P<0.05), and two SNPs, in position 1039 and 1041 of cagE, showed a highly significant association with cancer (p-value = 2.07×10⁻⁶), and the variant codon was located in the VirB3 homology domain of Agrobacterium. The results of this study may provide preliminary information to target antibiotic treatment to high-risk individuals, if effects of these variants are confirmed in further investigations.

The methylthiolation reaction mediated by the Radical-SAM enzymes

Over the past 10 years, considerable progress has been made in our understanding of the mechanistic enzymology of the Radical-SAM enzymes. It is now clear that these enzymes appear to be involved in a remarkably wide range of chemically challenging reactions. This review article highlights mechanistic and structural aspects of the methylthiotransferases (MTTases) sub-class of the Radical-SAM enzymes. The mechanism of methylthio insertion, now observed to be performed by three different enzymes is an exciting unsolved problem. This article is part of a Special Issue entitled: Radical SAM enzymes and Radical Enzymology.

Mapping codon usage of the translation initiation region in porcine reproductive and respiratory syndrome virus genome

BACKGROUND

Porcine reproductive and respitatory syndrome virus (PRRSV) is a recently emerged pathogen and severely affects swine populations worldwide. The replication of PRRSV is tightly controlled by viral gene expression and the codon usage of translation initiation region within each gene could potentially regulate the translation rate. Therefore, a better understanding of the codon usage pattern of the initiation translation region would shed light on the regulation of PRRSV gene expression.

RESULTS

In this study, the codon usage in the translation initiation region and in the whole coding sequence was compared in PRRSV ORF1a and ORFs2-7. To investigate the potential role of codon usage in affecting the translation initiation rate, we established a codon usage model for PRRSV translation initiation region. We observed that some non-preferential codons are preferentially used in the translation initiation region in particular ORFs. Although some positions vary with codons, they intend to use codons with negative CUB. Furthermore, our model of codon usage showed that the conserved pattern of CUB is not directly consensus with the conserved sequence, but shaped under the translation selection.

CONCLUSIONS

The non-variation pattern with negative CUB in the PRRSV translation initiation region scanned by ribosomes is considered the rate-limiting step in the translation process.

Relative codon adaptation index, a sensitive measure of codon usage bias

We propose a simple, sensitive measure of synonymous codon usage bias, the Relative Codon Adaptation Index (rCAI), as a way to discriminate better between highly biased and unbiased regions, compared with the widely used Codon Adaptation Index (CAI). CAI is a geometric mean of the relative usage of codons in a gene, and is calculated using the codon usage table trained with a set of highly expressed genes. In contrast, rCAI is computed by subtracting the background codon usage trained with two noncoding frames of highly expressed genes from the codon usage in the coding frame. rCAI has higher signal-to-noise ratio than CAI, considering that noncoding frames would not show codon bias. Translation efficiency and protein abundance correlates comparably or better with rCAI than CAI or other measures such as ‘effective number of codons’ and ‘SCUMBLE offsets’. Within overlapping coding regions, one of the two coding frames dominates in codon usage bias according to rCAI. Presumably, rCAI could substitute CAI in diverse applications.

Forces that influence the evolution of codon bias

The frequencies of alternative synonymous codons vary both among species and among genes from the same genome. These patterns have been inferred to reflect the action of natural selection. Here we evaluate this in bacteria. While intragenomic variation in many species is consistent with selection favouring translationally optimal codons, much of the variation among species appears to be due to biased patterns of mutation. The strength of selection on codon usage can be estimated by two different approaches. First, the extent of bias in favour of translationally optimal codons in highly expressed genes, compared to that in genes where selection is weak, reveals the long-term effectiveness of selection. Here we show that the strength of selected codon usage bias is highly correlated with bacterial growth rate, suggesting that selection has favoured translational efficiency. Second, the pattern of bias towards optimal codons at polymorphic sites reveals the ongoing action of selection. Using this approach we obtained results that were completely consistent with the first method; importantly, the frequency spectra of optimal codons at polymorphic sites were similar to those predicted under an equilibrium model. Highly expressed genes in Escherichia coli appear to be under continuing strong selection, whereas selection is very weak in genes expressed at low levels.

Characterization of RimO, a new member of the methylthiotransferase subclass of the radical SAM superfamily

RimO, encoded by the yliG gene in Escherichia coli, has been recently identified in vivo as the enzyme responsible for the attachment of a methylthio group on the beta-carbon of Asp88 of the small ribosomal protein S12 [Anton, B. P., Saleh, L., Benner, J. S., Raleigh, E. A., Kasif, S., and Roberts, R. J. (2008) Proc. Natl. Acad. Sci. U.S.A. 105, 1826-1831]. To date, it is the only enzyme known to catalyze methylthiolation of a protein substrate; the four other naturally occurring methylthio modifications have been observed on tRNA. All members of the methylthiotransferase (MTTase) family, to which RimO belongs, have been shown to contain the canonical CxxxCxxC motif in their primary structures that is typical of the radical S-adenosylmethionine (SAM) family of proteins. MiaB, the only characterized MTTase, and the enzyme experimentally shown to be responsible for methylthiolation of N(6)-isopentenyladenosine of tRNA in E. coli and Thermotoga maritima, has been demonstrated to harbor two distinct [4Fe-4S] clusters. Herein, we report in vitro biochemical and spectroscopic characterization of RimO. We show by analytical and spectroscopic methods that RimO, overproduced in E. coli in the presence of iron-sulfur cluster biosynthesis proteins from Azotobacter vinelandii, contains one [4Fe-4S](2+) cluster. Reconstitution of this form of RimO (RimO(rcn)) with (57)Fe and sodium sulfide results in a protein that contains two [4Fe-4S](2+) clusters, similar to MiaB. We also show by mass spectrometry that RimO(rcn) catalyzes the attachment of a methylthio group to a peptide substrate analogue that mimics the loop structure bearing aspartyl 88 of the S12 ribosomal protein from E. coli. Kinetic analysis of this reaction shows that the activity of RimO(rcn) in the presence of the substrate analogue does not support a complete turnover. We discuss the possible requirement for an assembled ribosome for fully active RimO in vitro. Our findings are consistent with those of other enzymes that catalyze sulfur insertion, such as biotin synthase, lipoyl synthase, and MiaB.

Codon engineering for improved antibody expression in mammalian cells

While well established in bacterial hosts, the effect of coding sequence variation on protein expression in mammalian systems is poorly characterized outside of viral proteins or proteins from distant phylogenetic families. The potential impact is substantial given the extensive use of mammalian expression systems in research and manufacturing of protein biotherapeutics. We are studying the effect of codon engineering on expression of recombinant antibodies with an emphasis on developing manufacturing cell lines. CNTO 888, a human mAb specific for human MCP-1, was obtained by antibody phage display in collaboration with MorphoSys AG. The isolated DNA sequence of the antibody was biased towards bacterial codons, reflecting the engineering of the Fab library for phage display expression in Escherichia coli. We compared the expression of CNTO 888 containing the parental V-region sequences with two engineered coding variants. In the native codon exchanged (NCE) variant, the V-region codons were replaced with those used in naturally derived human antibody genes. In the human codon optimized (HCO) variant the V-region codons were those used at the highest frequency based on a human codon usage table. The antibody expression levels from stable transfections in mammalian host cells were measured. The HCO codon variant of CNTO 888 yielded the highest expressing cell lines and the highest average expression for the screened populations. This had a significant positive effect on the process to generate a CNTO 888 production cell line and indicates the potential to improve antibody expression in mammalian expression systems by codon engineering.

Molecular and functional properties of the psychrophilic elongation factor G from the Antarctic Eubacterium Pseudoalteromonas haloplanktis TAC 125

The molecular and functional properties of the elongation factor (EF) G from the psychrophilic Antarctic eubacterium Pseudoalteromonas haloplanktis (Ph) were studied. PhEF-G catalyzed protein synthesis in vitro that was inhibited by fusidic acid, an antibiotic specifically acting on EF-G. The EF interacted with GDP only in the presence of P. haloplanktis ribosome and fusidic acid with an affinity similar to that displayed by Escherichia coli EF-G. The psychrophilic translocase elicited a ribosome-dependent GTPase that was competitively inhibited by GDP, the slowly hydrolyzable GTP analog GppNHp, and the protein synthesis inhibitor ppGDP. The temperature dependence of the activity of PhEF-G reached its maximum at least 26 degrees C beyond the growth temperature of P. haloplanktis (4-20 degrees C). The heat inactivation profile of the ribosome-dependent GTPase of PhEF-G gave a temperature for half inactivation (46 degrees C), significantly lower than that for half denaturation measured by either UV- (57 degrees C) or fluorescence-melting (62 degrees C). This finding was attributed to a different effect of the temperature on the catalytic domain with respect to that elicited on the other domains constituting the EF, thus confirming the differential molecular flexibility present in psychrophilic enzymes. A molecular model, based on the 3D coordinates of a thermophilic EF-G, showed differences only in connecting loops.

Predicting gene expression level from codon usage bias

The "expression measure" of a gene, E(g), is a statistic devised to predict the level of gene expression from codon usage bias. E(g) has been used extensively to analyze prokaryotic genome sequences. We discuss 2 problems with this approach. First, the formulation of E(g) is such that genes with the strongest selected codon usage bias are not likely to have the highest predicted expression levels; indeed the correlation between E(g) and expression level is weak among moderate to highly expressed genes. Second, in some species, highly expressed genes do not have unusual codon usage, and so codon usage cannot be used to predict expression levels. We outline a simple approach, first to check whether a genome shows evidence of selected codon usage bias and then to assess the strength of bias in genes as a guide to their likely expression level; we illustrate this with an analysis of Shewanella oneidensis.

Codon usage domains over bacterial chromosomes

The geography of codon bias distributions over prokaryotic genomes and its impact upon chromosomal organization are analyzed. To this aim, we introduce a clustering method based on information theory, specifically designed to cluster genes according to their codon usage and apply it to the coding sequences of Escherichia coli and Bacillus subtilis. One of the clusters identified in each of the organisms is found to be related to expression levels, as expected, but other groups feature an over-representation of genes belonging to different functional groups, namely horizontally transferred genes, motility, and intermediary metabolism. Furthermore, we show that genes with a similar bias tend to be close to each other on the chromosome and organized in coherent domains, more extended than operons, demonstrating a role of translation in structuring bacterial chromosomes. It is argued that a sizeable contribution to this effect comes from the dynamical compartimentalization induced by the recycling of tRNAs, leading to gene expression rates dependent on their genomic and expression context.

Modulation of Streptomyces Leucine Aminopeptidase by Calcium

Streptomyces griseus leucine aminopeptidase (SGAP), which has two zinc atoms in its active site, is clinically important as a model for understanding the structure and mechanism of action of other metallopeptidases. SGAP is a calcium-activated and calcium-stabilized enzyme, and its activation by calcium correlates with substrate specificity. In our previous study, we found a non-calcium-modulated leucine aminopeptidase secreted by Streptomyces septatus, the primary structure of which showed 71% identity with SGAP. In this study, we constructed chimeras of SGAP and S. septatus aminopeptidase by using an in vivo DNA shuffling system and several mutant enzymes by site-directed mutagenesis to identify the key residues in this modulation by calcium. We identified the key residues Asp-173 and Asp-174 of SGAP associated with both SGAP activation and stabilization by calcium. We also showed that the known calcium-binding site, which is composed of Asp-3, Ile-4, Asp-262, and Asp-266 of SGAP, only contributes to SGAP stabilization by calcium. Furthermore, we identified an important residue, Glu-196, that functions in cooperation with Asp-173, Asp-174, and calcium to increase the catalytic activity of SGAP.

The relationship among gene expression, folding free energy and codon usage bias in Escherichia coli

Taking advantage of microarray data in Escherichia coli genome, the relationship among mRNA expression levels, folding free energy and codon usage bias are investigated. Our results indicate that mRNA expression is correlated to the stability of mRNA secondary structure and the codon usage bias. The decrease of the stability of mRNA structure contributes to the increase of mRNA expression. There is a negative correlation between codon adaptation index (CAI) and mRNA expression in genes with less stable structure. The relationship between the stability of mRNA structure and mRNA half-life indicates the stability of mRNA structure is different from mRNA half-life.

Repeat-length-independent broad-spectrum shuffling, a novel method of generating a random chimera library in vivo

We describe a novel method of random chimeragenesis based on highly frequent deletion formation in the Escherichia coli ssb-3 strain and a deletion-directed chimera selection system that uses the rpsL(+) gene as a reporter. It enables the selection of chimeras without target gene expression and can therefore be applied to cytotoxic targets. When this system was applied to phospholipase D genes from Streptomyces septatus TH-2 and Streptomyces halstedii subsp. scabies K6 (examples of cytotoxic targets), chimeragenesis occurred between short identical sequences at the corresponding position of the parental genes with large variations. Chimeragenesis was >1,000 times more frequent in the ssb-3 background than in the ssb(+) background. We called this system repeat-length-independent broad-spectrum shuffling. It enables the convenient chimeragenesis and functional study of chimeric proteins. In fact, we found two amino acid residues related to the thermostability of phospholipase D (Phe426 and Thr433) by comparing thermostability among the chimeric enzymes obtained.

Mistranslation induced by streptomycin provokes a RecABC/RuvABC-dependent mutator phenotype in Escherichia coli cells

Translational stress-induced mutagenesis (TSM) refers to the mutator phenotype observed in Escherichia coli cells expressing a mutant allele (mutA or mutC) of the glycine tRNA gene glyV (or glyW). Because of an anticodon mutation, expression of the mutA allele results in low levels of Asp-->Gly mistranslation. The mutA phenotype does not require lexA-regulated SOS mutagenesis functions, and appears to be suppressed in cells defective for RecABC-dependent homologous recombination functions. To test the hypothesis that the TSM response is mediated by non-specific mistranslation rather than specific Asp-->Gly misreading, we asked if streptomycin (Str), an aminoglycoside antibiotic known to promote mistranslation, can provoke a mutator phenotype. We report that Str induces a strong mutator phenotype in cells bearing certain alleles of rpsL, the gene encoding S12, an essential component of the ribosomal 30 S subunit. The phenotype is strikingly similar to that observed in mutA cells in its mutational specificity, as well as in its requirement for RecABC-mediated homologous recombination functions. Expression of Str-inducible mutator phenotype correlates with mistranslation efficiency in response to Str. Thus, mistranslation in general is able to induce the TSM response. The Str-inducible mutator phenotype described here defines a new functional class of rpsL alleles, and raises interesting questions on the mechanism of action of Str, and on bacterial response to antibiotic stress.

Streptomycin-resistant and streptomycin-dependent mutants of the extreme thermophile Thermus thermophilus

We have isolated spontaneous streptomycin-resistant, streptomycin-dependent and streptomycin-pseudo-dependent mutants of the thermophilic bacterium Thermus thermophilus IB-21. All mutant phenotypes were found to result from single amino acid substitutions located in the rpsL gene encoding ribosomal protein S12. Spontaneous suppressors of streptomycin dependence were also readily isolated. Thermus rpsL mutations were found to be very similar to rpsL mutations identified in mesophilic organisms. This similarity affords greater confidence in the utility of the crystal structures of Thermus ribosomes to interpret biochemical and genetic data obtained with Escherichia coli ribosomes. In the X-ray crystal structure of the T. thermophilus HB8 30 S subunit, the mutated residues are located in close proximity to one another and to helices 18, 27 and 44 of 16 S rRNA. X-ray crystallographic analysis of ribosomes from streptomycin-resistant, streptomycin-pseudo-dependent and streptomycin-dependent mutants described here is expected to reveal fundamental insights into the mechanism of tRNA selection, translocation, and conformational dynamics of the ribosome.

Absence of translationally selected synonymous codon usage bias in Helicobacter pylori

Synonymous codon usage in the complete genome of Helicobacter pylori was investigated. The moderate A+T-richness of the genome (G+C=39 mol%) is reflected in the overall synonymous codon usage but the frequencies of some codons cannot be explained by simple mutational biases. A low level of heterogeneity among genes was observed, but this does not appear to be due to varying mutational bias or translational selection. Some of the heterogeneity was due to amino acid composition variation among the encoded proteins, and some may be attributable to recent acquisition of genes from other species. Since Hel. pylori codon usage is not dominated by biased mutation patterns, the absence of evidence for translationally mediated selection among synonymous codons is striking. This has implications with regard to the life history of this species, and in particular suggests that Hel. pylori strains are not subject to periods of competitive exponential growth. Despite the lack of selected codon usage, base composition immediately after the translation initiation site is skewed, consistent with selection against secondary structure formation in this region.

Chromosome mapping of the sweet potato little leaf phytoplasma reveals genome heterogeneity within the phytoplasmas

To further understand the genomic diversity and genetic architecture of phytoplasmas, a physical and genetic map of the sweet potato little leaf (SPLL) strain V4 phytoplasma chromosome was determined. PFGE was used to determine the size of the SPLL-V4 genome, which was estimated to be 622 kb. A physical map was prepared by two-dimensional reciprocal digestions using the restriction endonucleases BssHII, Smal, Eagl and I-Ceul. Sixteen cleavage sites were located on the map. Southern hybridizations of digested SPLL-V4 chromosomal DNA were done using random clones and PCR-amplified genes as probes. This confirmed fragment positions and located the two rRNA operons and the linked fus/tuf genes encoding elongation factors G and Tu, respectively, on the physical map. An inversion of one of the rRNA operons was observed from hybridization data. Sequence analysis of one of the random clones identified a gid gene encoding a glucose-inhibited division protein. Digestions of the tomato big bud (TBB) phytoplasma chromosome with the same four enzymes revealed genome heterogeneity when compared to the closely related SPLL-V4, and a preliminary chromosome size for the TBB phytoplasma of 662 kb was estimated. This mapping information has revealed that significant genome diversity exists within the phytoplasmas.

Use of silent mutations in cDNA encoding human glutathione transferase M2-2 for optimized expression in Escherichia coli

Heterologous expression of human glutathione transferase M2-2 (GST M2-2) using Escherichia coli was improved 140-fold by mutating the cDNA expressing the enzyme. Expression of GST M2-2 from this cDNA clone, pKHXhGM2, generated approximately 190 mg protein per liter of bacterial culture, corresponding to approximately 12% of the total amount of soluble protein. The high-level-expressing cDNA was generated by oligonucleotide-directed mutagenesis introducing alternative silent mutations into the third nucleotide of codons 2, 4-7, and 10-14 in the 5' end of the cDNA coding region. The choice of alternative codons was restricted to those naturally occurring in highly biased genes in E. coli. Furthermore, the wild-type TAG stop codon at the 3' end was replaced with the two stop codons TAA and TGA in tandem to increase translation termination efficiency. The resulting partially randomized cDNA library was assayed for high-level expression using immunoscreening. Sequence similarities between the constructed high-level-expressing GST M2-2 cDNA and a similarly designed cDNA encoding the closely related human GST M1-1 suggest that the codons in the region immediately following the start codon are influential in achieving high-level expression. Pyrimidines seem to be more favorable than purines in the third position of codons in optimizing the expression of these enzymes in E. coli.

Modelling in Escherichia coli of mutations in mitoribosomal protein S12: novel mutant phenotypes of rpsL

The rpsL gene of Escherichia coli encodes the highly conserved rps12 protein of the ribosomal accuracy centre. We have used the E. coli gene to model the phenotypic effects of specific substitutions found in the mitochondrial gene for rps12. Variants created by in vitro mutagenesis were tested in two different plasmid vector systems, in both streptomycin-sensitive and streptomycin-resistant hosts. A substitution with respect to eubacterial rps12 (K87-->Q), found in all metazoan and fungal mitochondrial orthologues thus far studied, is associated with low-level resistance to streptomycin and a modest (15%) drop in translational elongation rate, but without significant effects on translational accuracy. An amino-acid replacement at a highly conserved leucine residue (L56-->H), associated with the phenotype of sensitivity to mechanical vibration and hemizygous female lethality in Drosophila, creates a functionally inactive but structurally stable protein that is not assembled into ribosomes. The presence in the cell of the mutant, but not wild-type, rpsL greatly downregulates the level of a prominent polypeptide of approximately 50 kDa. These results indicate novel structure-function relationships in rps12 genes affecting translational function, ribosome assembly and drug sensitivity, and indicate a novel regulatory pathway that may influence ribosome biogenesis.

Chromosomal organization and nucleotide sequence of the genes coding for the elongation factors G and Tu of the apple proliferation phytoplasma

Genes coding for elongation factors G (fus) and Tu (tuf) of the non-culturable apple proliferation (AP) phytoplasma were cloned and sequenced. Arrangement of these genes and identification of the ribosomal protein gene rps7 upstream of the fus gene suggest a transcriptional organization similar to that of the streptomycin operon of Escherichia coli and other bacteria. The fus and tuf genes from other tested phytoplasmas were found to be similarly linked as in the AP agent. Thus, it is likely that they show a similar chromosomal arrangement. This organization would be in contrast to that of the phylogenetically distinctly different culturable mollicutes of the genus Mycoplasma in which the tuf and fus genes are separately transcribed.

Translational selection and molecular evolution

An interplay among experimental studies of protein synthesis, evolutionary theory, and comparisons of DNA sequence data has shed light on the roles of natural selection and genetic drift in 'silent' DNA evolution.

In vivo assay of protein-protein interactions in Hin-mediated DNA inversion

In order to form the catalytic nucleoprotein complex called the invertasome in the Hin-mediated DNA inversion reaction, interactions of the DNA-binding proteins Hin and Fis are required. Assays for these protein-protein interactions have been exploited with protein cross-linkers in vitro. In this study, an in vivo assay system that probes protein-protein interactions was developed. The formation of a DNA loop generated by protein interactions resulted in transcriptional repression of an artificially designed operon, which in turn increased the chance of survival of Escherichia coli host cells in a streptomycin-containing medium. Using this system, we were able to assay the Hin-Hin interaction that results in the pairing of the two recombination sites and protein interactions that result in the formation of the invertasome. This assay system also led us to find that an individual Hin dimer bound on a recombination site can form a stable complex with Fis bound on the recombinational enhancer; this finding has never been observed in in vitro studies. Possible pathways toward the formation of the invertasome are discussed based on the assay results for a previously reported Hin mutant.

Linkage map of Escherichia coli K-12, edition 10: the traditional map

This map is an update of the edition 9 map by Berlyn et al. (M. K. B. Berlyn, K. B. Low, and K. E. Rudd, p. 1715-1902, in F. C. Neidhardt et al., ed., Escherichia coli and Salmonella: cellular and molecular biology, 2nd ed., vol. 2, 1996). It uses coordinates established by the completed sequence, expressed as 100 minutes for the entire circular map, and adds new genes discovered and established since 1996 and eliminates those shown to correspond to other known genes. The latter are included as synonyms. An alphabetical list of genes showing map location, synonyms, the protein or RNA product of the gene, phenotypes of mutants, and reference citations is provided. In addition to genes known to correspond to gene sequences, other genes, often older, that are described by phenotype and older mapping techniques and that have not been correlated with sequences are included.

Metazoan nuclear genes for mitoribosomal protein S12

We have characterized nuclear genes for mitoribosomal protein S12 (mt-rps12) a major component of the ribosomal accuracy centre, in human, mouse and Drosophila melanogaster. In human and Drosophila, and probably also in mouse, there is a single intron within the coding region, located in the mitochondrial targeting pre-sequence. In humans, the mRNA structure is highly suggestive of translational regulation. In all three species, there is an amino-acid substitution with respect to eubacterial homologues in a residue implicated in aminoglycoside resistance. The only viable mutant allele of the Drosophila gene, associated with a bang-sensitive phenotype (paralysis upon mechanical vibration, arising from a mechanoreceptor cell defect) also has a novel substitution in a conserved region implicated in translational fidelity. Given the involvement of the mitoribosomal accuracy centre in human sensorineural deafness by virtue of rRNA mutations, our results indicate that this fly mutant may be a useful animal model of this disorder, and earmark the gene for mt-rps12 as a candidate in human hearing impairment.

A nucleotide polymorphism in ERCC1 in human ovarian cancer cell lines and tumor tissues

We studied the DNA sequence of the entire coding region of ERCC1 gene, in five cell lines established from human ovarian cancer (A2780, A2780/CP70, MCAS, OVCAR-3, SK-OV-3), 29 human ovarian cancer tumor tissue specimens, one human T-lymphocyte cell line (H9), and non-malignant human ovary tissue (NHO). Samples were assayed by PCR-SSCP and DNA sequence analyses. A silent mutation at codon 118 (site for restriction endonuclease MaeII) in exon 4 of the gene was detected in MCAS, OVCAR-3 and SK-OV-3 cells, and NHO. This mutation was a C-->T transition, that codes for the same amino acid: asparagine. This transition converts a common codon usage (AAC) to an infrequent codon usage (AAT), whereas frequency of use is reduced two-fold. This base change was associated with a detectable band shift on SSCP analysis. For the 29 ovarian cancer specimens, the same base change was observed in 15 tumor samples and was associated with the same band shift in exon 4. Cells and tumor tissue specimens that did not contain the C-->T transition, did not show the band shift in exon 4. Our data suggest that this alteration at codon 118 within the ERCC1 gene, may exist in platinum-sensitive and platinum-resistant ovarian cancer tissues.

Growth phase-dependent transcription of the Streptomyces ramocissimus tuf1 gene occurs from two promoters

The str operon of Streptomyces ramocissimus contains the genes for ribosomal proteins S12 (rpsL) and S7 (rpsG) and for the polypeptide chain elongation factors G (EF-G) (fus) and Tu (EF-Tu) (tuf). This kirromycin producer contains three tuf or tuf-like genes; tuf1 encodes the regular EF-Tu and is located immediately downstream of fus. In vivo and in vitro transcription analysis revealed a transcription start site directly upstream of S. ramocissimus tuf1, in addition to the operon promoter rpsLp. Transcription from these promoters appeared to be growth phase dependent, diminishing drastically upon entry into stationary phase and at the onset of production of the EF-Tu-targeted antibiotic kirromycin. In surface-grown cultures, a second round of tuf1 transcription, coinciding with aerial mycelium formation and kirromycin production, was observed. The tuf1-specific promoter (tuf1p) was located in the intercistronic region between fus and tuf1 by high-resolution S1 mapping, in vitro transcription, and in vivo promoter probing. During logarithmic growth, the tuf1p and rpsLp transcripts are present at comparable levels. In contrast to Escherichia coli, which has two almost identical tuf genes, the gram-positive S. ramocissimus contains only tuf1 for its regular EF-Tu. High levels of EF-Tu may therefore be achieved by the compensatory activity of tuf1p.

Use of rpsL for dominance selection and gene replacement in Streptomyces roseosporus

We developed a gene replacement system using the rpsL gene of Streptomyces roseosporus and demonstrated its utility by constructing a deletion in the S. roseosporus glnA gene. A 1.3-kb BamHI fragment that hybridized to the Mycobacterium smegmatis rpsL gene was subcloned from an S. roseosporus cosmid library and sequenced. Plasmid pRHB514 containing the rpsL gene conferred streptomycin sensitivity (Sm(S)) to the Sm(r) S. roseosporus TH149. The temperature-sensitive plasmid pRHB543 containing rpsL and the S. roseosporus glnA gene disrupted with a hygromycin resistance (Hm(r)) gene was introduced into S. roseosporus TH149, and recombinants containing single and double crossovers were obtained after a temperature increase. Southern hybridization analysis revealed that single crossovers occurred in the glnA or rpsL genes and that double crossovers resulted in replacement of the chromosomal glnA gene with the disrupted glnA. Glutamine synthetase activity was undetectable in the recombinant containing the disrupted glnA gene.

Induction of actinorhodin production by rpsL (encoding ribosomal protein S12) mutations that confer streptomycin resistance in Streptomyces lividans and Streptomyces coelicolor A3(2)

A strain of Streptomyces lividans, TK24, was found to produce a pigmented antibiotic, actinorhodin, although S. lividans normally does not produce this antibiotic. Genetic analyses revealed that a streptomycin-resistant mutation str-6 in strain TK24 is responsible for induction of antibiotic synthesis. DNA sequencing showed that str-6 is a point mutation in the rpsL gene encoding ribosomal protein S12, changing Lys-88 to Glu. Gene replacement experiments with the Lys88-->Glu str allele demonstrated unambiguously that the str mutation is alone responsible for the activation of actinorhodin production observed. In contrast, the strA1 mutation, a genetic marker frequently used for crosses, did not restore actinorhodin production and was found to result in an amino acid alteration of Lys-43 to Asn. Induction of actinorhodin production was also detected in strain TK21, which does not harbor the str-6 mutation, when cells were incubated with sufficient streptomycin or tetracycline to reduce the cell's growth rate, and 40 and 3% of streptomycin- or tetracycline-resistant mutants, respectively, derived from strain TK21 produced actinorhodin. Streptomycin-resistant mutations also blocked the inhibitory effects of relA and brgA mutations on antibiotic production, aerial mycelium formation or both. These str mutations changed Lys-88 to Glu or Arg and Arg-86 to His in ribosomal protein S12. The decrease in streptomycin production in relC mutants in Streptomyces griseus could also be abolished completely by introducing streptomycin-resistant mutations, although the impairment in antibiotic production due to bldA (in Streptomyces coelicolor) or afs mutations (in S. griseus) was not eliminated. These results indicate that the onset and extent of secondary metabolism in Streptomyces spp. is significantly controlled by the translational machinery.

A chimeric disposition of the elongation factor genes in Rickettsia prowazekii

An exceptional disposition of the elongation factor genes is observed in Rickettsia prowazekii, in which there is only one tuf gene, which is distant from the lone fus gene. In contrast, the closely related bacterium Agrobacterium tumefaciens has the normal bacterial arrangement of two tuf genes, of which one is tightly linked to the fus gene. Analysis of the flanking sequences of the single tuf gene in R. prowazekii shows that it is preceded by two of the four tRNA genes located in the 5' region of the Escherichia coli tufB gene and that it is followed by rpsJ as well as associated ribosomal protein genes, which in E. coli are located downstream of the tufA gene. The fus gene is located within the str operon and is followed by one tRNA gene as well as by the genes secE and nusG, which are located in the 3' region of tufB in E. coli. This atypical disposition of genes suggests that intrachromosomal recombination between duplicated tuf genes has contributed to the evolution of the unique genomic architecture of R. prowazekii.

Ribosome-mediated translational pause and protein domain organization

Because regions on the messenger ribonucleic acid differ in the rate at which they are translated by the ribosome and because proteins can fold cotranslationally on the ribosome, a question arises as to whether the kinetics of translation influence the folding events in the growing nascent polypeptide chain. Translationally slow regions were identified on mRNAs for a set of 37 multidomain proteins from Escherichia coli with known three-dimensional structures. The frequencies of individual codons in mRNAs of highly expressed genes from E. coli were taken as a measure of codon translation speed. Analysis of codon usage in slow regions showed a consistency with the experimentally determined translation rates of codons; abundant codons that are translated with faster speeds compared with their synonymous codons were found to be avoided; rare codons that are translated at an unexpectedly higher rate were also found to be avoided in slow regions. The statistical significance of the occurrence of such slow regions on mRNA spans corresponding to the oligopeptide domain termini and linking regions on the encoded proteins was assessed. The amino acid type and the solvent accessibility of the residues coded by such slow regions were also examined. The results indicated that protein domain boundaries that mark higher-order structural organization are largely coded by translationally slow regions on the RNA and are composed of such amino acids that are stickier to the ribosome channel through which the synthesized polypeptide chain emerges into the cytoplasm. The translationally slow nucleotide regions on mRNA possess the potential to form hairpin secondary structures and such structures could further slow the movement of ribosome. The results point to an intriguing correlation between protein synthesis machinery and in vivo protein folding. Examination of available mutagenic data indicated that the effects of some of the reported mutations were consistent with our hypothesis.

Beta-methylthio-aspartic acid: identification of a novel posttranslational modification in ribosomal protein S12 from Escherichia coli

Utilizing microscale chemical derivatization reactions and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry, we have identified a novel posttranslational modification of aspartic acid, beta-methylthio-aspartic acid. The modified residue is located at position 88 in ribosomal protein S12 from Escherichia coli, a phylogenetically conserved protein that has been implicated in maintaining translational accuracy of the ribosome.

Conservation of the organization of the mitochondrial nad3 and rps12 genes in evolutionarily distant angiosperms

The organization of the genes nad3 and rps12 has been investigated in the mitochondrial genome of two dicotyledonous plants - Helianthus and Magnolia - and one monocotyledonous plant (Allium). These plants all contain a complete rps12 gene downstream of the nad3 gene. This arrangement is thus highly conserved within angiosperms. The two genes are co-transcribed and the transcript is modified at several positions by RNA editing of the C to U-type, thus confirming that both genes encode functional proteins. Some 26, 35 and 27 editing events have been identified in the PCR-derived nad3-rps12 cDNA population from sunflower, Magnolia and onion, respectively. Editing of the nad3-rps12 transcript is thus more extensive in Magnolia than in the other angiosperms so far investigated and radically changes the genomically encoded polypeptide sequence. A novel species-specific codon modification was observed in Magnolia. Several homologous sites show differences in editing pattern among plant species. A C-to-U alteration is also found in the non-coding region separating the nad3 and rps12 genes in sunflower. The PCR-derived cDNA populations from the nad3-rps12 loci analysed were found to be differently edited. In addition the plant species show marked variations in the completeness of RNA editing, with only the Magnolia nad3 mRNA being edited fully.

Impaired splicing of the rps12 transcript in ribosome-deficient plastids

Analysis of RNA maturation in ribosome-deficient plastids of four non-allelic barley mutants revealed an increased accumulation and altered processing of transcripts of the ribosomal protein gene CS12 (rps12) compared to normal chloroplasts. The three exons of rps12 are part of two different polycistronic transcription units. Generation of mature rps12-mRNA involves both cis- and trans-splicing. In ribosome-deficient plastids, the cis-intron separating exons 2 and 3 remains entirely unspliced whereas the splicing of the bipartite rps12 trans-intron between exon 1 and exon 2 occurs, but at a reduced level. A comparison of the 3' and 5' ends of the two RNAs that are generally assumed to interact during trans-splicing showed a difference in the processing pathways of 3' rps12 transcripts between mutated and normal chloroplasts. Nonetheless, the final products were identical.

Protease-dependent streptomycin sensitivity in E. coli--a system for protease inhibitor selection

We have developed a bacterial cell system in which the activity of an expressed heterologous protease confers a dominant streptomycin-sensitive (strs) phenotype on the cell. This phenotype owes its high selectivity to the fact that streptomycin (strep) resistance, which is conferred on E. coli by mutants of ribosomal protein S12, is highly recessive to strep sensitivity. Thus, when strep-resistant (strr) strains of E. coli are transformed to co-express the wild-type allele of S12 in addition to the mutant allele, their sensitivity to strep increases by a factor of 100-1000. Similarly, we found that when the same strr strains were transformed to co-express a heterologous protease and an inactive fusion of S12 with a substrate of the protease, the strep sensitivity of the cells increased approximately 100-fold. This effect was strictly dependent on correct cleavage of the S12 precursor, required only modest levels of expression of protease and substrate, and could be competitively inhibited by co-expression of an alternative substrate gene. This system thus appears to be well-suited to the identification of protease inhibitors, either by selection from libraries of endogenously expressed random peptide-encoding genes, or by screening synthetic or natural products libraries. Protease-dependent dominant phenotypes may be more sensitive and appropriate than the more commonly used recessive phenotypes for proteases which are activating enzymes.

Chromosomal organization and nucleotide sequence of the genes for elongation factors EF-1 alpha and EF-2 and ribosomal proteins S7 and S10 of the hyperthermophilic archaeum Desulfurococcus mobilis

The Desulfurococcus mobilis genes fus (encoding EF-2) and tuf (for EF-1 alpha) were cloned and sequenced together with genes for ribosomal proteins S10 (rps10) and S7 (rps7). Unlike Methanococcus, which displays the bacterial-like fus and tuf gene context 5'-rps12-rps7-fus-tuf-3', and similar to Sulfolobus and Pyrococcus, the Desulfurococcus fus gene (734 codons) has a distinct chromosomal location. Moreover, tuf (441 codons) is the promoter-proximal unit of a three-gene cluster comprising the genes rps10 (98 codons) and tRNA(Ser); the arrangement of the cluster is 5'-tuf-91 bp spacer -rps10-138 bp spacer -tRNA(Ser)-3' and the tuf gene is preceded by a canonical archaeal promoter. The D. mobilis gene rps7 (198 codons) is located further upstream from tuf (535 bp 'silent' intergenic spacing) and no rps12 homolog occurs in its immediate vicinity. Also, judging from putative promoter and transcription termination sequences, rps7 appears to be separately transcribed. Analysis of the predicted fus and tuf gene products revealed the three consensus motifs characteristic of GTP-binding proteins, and the fus-encoded EF-2 protein also displayed the consensus sequence required for ADP-ribosylation by Diphtheria toxin. Both EF sequences were definitely crenarchaeal by comparison with available homologs from other Archaea. Outgroup-rooted phylogenies derived from the sequences of ribosomal proteins S10 and S7 yielded the Sulfolobus-Desulfurococcus association at a high bootstrap confidence level.