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

Translation of Plant RNA Viruses

Plant RNA viruses encode essential viral proteins that depend on the host translation machinery for their expression. However, genomic RNAs of most plant RNA viruses lack the classical characteristics of eukaryotic cellular mRNAs, such as mono-cistron, 5′ cap structure, and 3′ polyadenylation. To adapt and utilize the eukaryotic translation machinery, plant RNA viruses have evolved a variety of translation strategies such as cap-independent translation, translation recoding on initiation and termination sites, and post-translation processes. This review focuses on advances in cap-independent translation and translation recoding in plant viruses.

Correction of frameshift mutations in the atpB gene by translational recoding in chloroplasts of Oenothera and tobacco

Translational recoding, also known as ribosomal frameshifting, is a process that causes ribosome slippage along the messenger RNA, thereby changing the amino acid sequence of the synthesized protein. Whether the chloroplast employs recoding is unknown. I-iota, a plastome mutant of Oenothera (evening primrose), carries a single adenine insertion in an oligoA stretch [11A] of the atpB coding region (encoding the β-subunit of the ATP synthase). The mutation is expected to cause synthesis of a truncated, nonfunctional protein. We report that a full-length AtpB protein is detectable in I-iota leaves, suggesting operation of a recoding mechanism. To characterize the phenomenon, we generated transplastomic tobacco lines in which the atpB reading frame was altered by insertions or deletions in the oligoA motif. We observed that insertion of two adenines was more efficiently corrected than insertion of a single adenine, or deletion of one or two adenines. We further show that homopolymeric composition of the oligoA stretch is essential for recoding, as an additional replacement of AAA lysine codon by AAG resulted in an albino phenotype. Our work provides evidence for the operation of translational recoding in chloroplasts. Recoding enables correction of frameshift mutations and can restore photoautotrophic growth in the presence of a mutation that otherwise would be lethal.

Integrated Proteomics and Transcriptomics Analyses Reveal the Transcriptional Slippage of P3N-PIPO in a Bymovirus

P3N-PIPO (P3 N-terminal fused with Pretty Interesting Potyviridae ORF), the movement protein of potyviruses, is expressed as a translational fusion with the N-terminus of P3 in potyviruses. As reported in previous studies, P3N-PIPO is expressed via transcriptional slippage at a conserved G2A6 slippery site in the genus Potyvirus. However, it is still unknown whether a similar expression mechanism of P3N-PIPO is used in the other genera of the family Potyviridae. Moreover, due to the extremely low expression level of P3N-PIPO in natural virus-infected plants, the peptides spanning the slippery site which provide direct evidence of the slippage at the protein level, have not been identified yet. In this study, a potato virus X (PVX)-based expression vector was utilized to investigate the expression mechanism of P3N-PIPO. A high expression level of the P3N-PIPO(WT) of turnip mosaic virus (TuMV, genus Potyvirus) was observed based on the PVX expression vector. For the first time, we successfully identified the peptides of P3N-PIPO spanning the slippery site by mass spectrometry. Likewise, the P3N-PIPO(WT) of wheat yellow mosaic virus (WYMV, genus Bymovirus) was also successfully expressed using the PVX expression vector. Integrated proteome and transcriptome analyses revealed that WYMV P3N-PIPO was expressed at the conserved G2A6 site through transcriptional slippage. Moreover, as revealed by mutagenesis analysis, Hexa-adenosine of the G2A6 site was important for the frameshift expression of P3N-PIPO in WYMV. According to our results, the PVX-based expression vector might be used as an excellent tool to study the expression mechanism of P3N-PIPO in Potyviridae. To the best of our knowledge, this is the first experimental evidence for the expression mechanism of P3N-PIPO in the genus Bymovirus, the only genus comprising bipartite virus species in the family Potyviridae.

Mycobacterium tuberculosis and M. bovis BCG Moreau Fumarate Reductase Operons Produce Different Polypeptides That May Be Related to Non-canonical Functions

Tuberculosis is a world widespread disease, caused by Mycobacterium tuberculosis (M.tb). Although considered an obligate aerobe, this organism can resist life-limiting conditions such as microaerophily mainly due to its set of enzymes responsible for energy production and coenzyme restoration under these conditions. One of these enzymes is fumarate reductase, an heterotetrameric complex composed of a catalytic (FrdA), an iron-sulfur cluster (FrdB) and two transmembrane (FrdC and FrdD) subunits involved in anaerobic respiration and important for the maintenance of membrane potential. In this work, aiming to further characterize this enzyme function in mycobacteria, we analyzed the expression of FrdB-containing proteins in M.tb and Mycobacterium bovis Bacillus Calmette–Guérin (BCG) Moreau, the Brazilian vaccine strain against tuberculosis. We identified three isoforms in both mycobacteria, two of them corresponding to the predicted encoded polypeptides of M.tb (27 kDa) and BCG Moreau (40 kDa) frd sequences, as due to an insertion on the latter’s operon a fused FrdBC protein is expected. The third 52 kDa band can be explained by a transcriptional slippage event, typically occurring when mutation arises in a repetitive region within a coding sequence, thought to reduce its impact allowing the production of both native and variant forms. Comparative modeling of the M.tb and BCG Moreau predicted protein complexes allowed the detection of subtle overall differences, showing a high degree of structure and maybe functional resemblance among them. Axenic growth and macrophage infection assays show that the frd locus is important for proper bacterial development in both scenarios, and that both M.tb’s and BCG Moreau’s alleles can partially revert the hampered phenotype of the knockout strain. Altogether, our results show that the frdABCD operon of Mycobacteria may have evolved to possess other yet non-described functions, such as those necessary during aerobic logarithmic growth and early stage steps of infection.

Immunohistochemical detection of the pro-apoptotic Bax∆2 protein in human tissues

The pro-apoptotic Bax isoform Bax∆2 was originally discovered in cancer patients with a microsatellite guanine deletion (G8 to G7). This deletion leads to an early stop codon; however, when combined with the alternative splicing of exon 2, the reading frame is restored allowing production of a full-length protein (Bax∆2). Unlike the parental Baxα, Bax∆2 triggers apoptosis through a non-mitochondrial pathway and the expression in human tissues was unknown. Here, we analyzed over 1000 tissue microarray samples from 13 different organs using immunohistochemistry. Bax∆2-positive cells were detected in all examined organs at low rates (1-5%) and mainly scattered throughout the connective tissues. Surprisingly, over 70% of normal colon samples scored high for BaxΔ2-positive staining. Only 7% of malignant colon samples scored high, with most high-grade tumors being negative. A similar pattern was observed in most organs examined. We also showed that both Baxα and Bax∆2 can co-exist in the same cells. Genotyping showed that the majority of Bax∆2-positive normal tissues contain no G7 mutation, but an unexpected high rate of G9 was observed. Although the underlying mechanism remains to be explored, the inverse correlation of Bax∆2 expression with tissue malignancy suggests that it may have a clinical implication in cancer development and treatment.

Heterogeneous Sequences of Brain Cytoplasmic 200 RNA Formed by Multiple Adenine Nucleotide Insertions

Brain cytoplasmic 200 RNA (BC200 RNA), originally identified as a neuron-specific non-coding RNA, is also observed in various cancer cells that originate from non-neural cells. Studies have revealed diverse functions of BC200 RNA in cancer cells. Accordingly, we hypothesized that BC200 RNA might be modified in cancer cells to generate cancerous BC200 RNA responsible for its cancer-specific functions. Here, we report that BC200 RNA sequences are highly heterogeneous in cancer cells by virtue of multiple adenine nucleotide insertions in the internal A-rich region. The insertion of adenine nucleotides enhances BC200 RNAmediated translation inhibition, possibly by increasing the binding affinity of BC200 RNA for eIF4A (eukaryotic translation initiation factor 4A).

Localized Context-Dependent Effects of the "Ambush" Hypothesis: More Off-Frame Stop Codons Downstream of Shifty Codons

The ambush hypothesis speculates that off-frame stop codons increase translational efficiency after ribosomal frameshifts by stopping early frameshifted translation. Some evidences fit this hypothesis: (1) synonymous codon usages increase with their potential contribution to off-frame stops; (2) the genetic code assigns frequent amino acids to codon families contributing to off-frame stops; (3) positive biases for off-frame stops (AT rich) occur despite adverse nucleotide (GC) biases; and (4) mitochondrial off-frame stop codon densities increase with ribosomal structural instability, potential proxy of frameshift frequencies. In this study, analyses of vertebrate mitogenes and tRNA synthetase genes from all superkingdoms and viruses test a new prediction of the ambush hypothesis: sequences immediately downstream of frameshift-inducing homopolymer codons (AAA, CCC, GGG, and TTT) are off-frame stop rich. Codons immediately downstream of homopolymer codons form more than average off-frame stops, biases are stronger than for corresponding upstream distances and for any other group of synonymous codons. Sequences downstream of that high-density region are off-frame stop depleted. This decrease suggests that off-frame stops, combined with suppressor tRNAs regulate translation of overlapping coding sequences. Results show the predictive power of the ambush hypothesis, from macroevolutionary (genetic code structure) to detailed gene sequence anatomy.

Non-programmed transcriptional frameshifting is common and highly RNA polymerase type-dependent

Background

The viral or host systems for a gene expression assume repeatability of the process and high quality of the protein product. Since level and fidelity of transcription primarily determines the overall efficiency, all factors contributing to their decrease should be identified and optimized. Among many observed processes, non-programmed insertion/deletion (indel) of nucleotide during transcription (slippage) occurring at homopolymeric A/T sequences within a gene can considerably impact its expression. To date, no comparative study of the most utilized Escherichia coli and T7 bacteriophage RNA polymerases (RNAP) propensity for this type of erroneous mRNA synthesis has been reported. To address this issue we evaluated the influence of shift-prone A/T sequences by assessing indel-dependent phenotypic changes. RNAP-specific expression profile was examined using two of the most potent promoters, P_araBAD of E. coli and φ10 of phage T7.

Results

Here we report on the first systematic study on requirements for efficient transcriptional slippage by T7 phage and cellular RNAPs considering three parameters: homopolymer length, template type, and frameshift directionality preferences. Using a series of out-of-frame gfp reporter genes fused to a variety of A/T homopolymeric sequences we show that T7 RNAP has an exceptional potential for generating frameshifts and is capable of slipping on as few as three adenine or four thymidine residues in a row, in a flanking sequence-dependent manner. In contrast, bacterial RNAP exhibits a relatively low ability to baypass indel mutations and requires a run of at least 7 tymidine and even more adenine residues. This difference comes from involvement of various intrinsic proofreading properties. Our studies demonstrate distinct preference towards a specific homopolymer in slippage induction. Whereas insertion slippage performed by T7 RNAP (but not deletion) occurs tendentiously on poly(A) rather than on poly(T) runs, strong bias towards poly(T) for the host RNAP is observed.

Conclusions

Intrinsic RNAP slippage properties involve trade-offs between accuracy, speed and processivity of transcription. Viral T7 RNAP manifests far greater inclinations to the transcriptional slippage than E. coli RNAP. This possibly plays an important role in driving bacteriophage adaptation and therefore could be considered as beneficial. However, from biotechnological and experimental viewpoint, this might create some problems, and strongly argues for employing bacterial expression systems, stocked with proofreading mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1034-4) contains supplementary material, which is available to authorized users.

Firmicutes-enriched IS1447 represents a group of IS3-family insertion sequences exhibiting unique + 1 transcriptional slippage

BACKGROUND

Bacterial insertion sequences (ISs) are ubiquitous mobile genetic elements that play important roles in genome plasticity, cell adaptability, and function evolution. ISs of various families and subgroups contain significantly diverse molecular features and functional mechanisms that are not fully understood.

RESULTS

IS1447 is a member of the widespread IS3 family and was previously detected to have transposing activity in a typical thermophilic and cellulolytic microorganism Clostridium thermocellum. Phylogenetic analysis showed that IS1447-like elements are widely distributed in Firmicutes and possess unique features in the IS3 family. Therefore, IS1447 may represent a novel subgroup of the IS3 family. Unlike other well-known IS3 subgroups performing programmed - 1 translational frameshifting for the expression of the transposase, IS1447 exhibits transcriptional slippage in both the + 1 and - 1 directions, each with a frequency of ~ 16%, and only + 1 slippage results in full-length and functional transposase. The slippage-prone region of IS1447 contains a run of nine A nucleotides following a stem-loop structure in mRNA, but mutagenesis analysis indicated that seven of them are sufficient for the observed slippage. Western blot analysis indicated that IS1447 produces three types of transposases with alternative initiations. Furthermore, the IS1447-subgroup elements are abundant in the genomes of several cellulolytic bacteria.

CONCLUSION

Our result indicated that IS1447 represents a new Firmicutes-enriched subgroup of the IS3 family. The characterization of the novel IS3-family member will enrich our understanding of the transposition behavior of IS elements and may provide insight into developing IS-based mutagenesis tools for thermophiles.

Evaluation of GFP reporter utility for analysis of transcriptional slippage during gene expression

Background

Epimutations arising from transcriptional slippage seem to have more important role in regulating gene expression than earlier though. Since the level and the fidelity of transcription primarily determine the overall efficiency of gene expression, all factors contributing to their decrease should be identified and optimized.

Results

To examine the influence of A/T homopolymeric sequences on introduction of erroneous nucleotides by slippage mechanism green fluorescence protein (GFP) reporter was chosen. The in- or out-of-frame gfp gene was fused to upstream fragment with variable number of adenine or thymine stretches resulting in several hybrid GFP proteins with diverse amino acids at N-terminus. Here, by using T7 phage expression system we showed that the intensity of GFP fluorescence mainly depends on the number of the retained natural amino acids. While the lack of serine (S₂) residue results in negligible effects, the lack of serine and lysine (S₂K₃) contributed to a significant reduction in fluorescence by 2.7-fold for polyA-based in-frame controls and twofold for polyTs. What is more, N-terminal tails amino acid composition was rather of secondary importance, since the whole-cell fluorescence differed in a range of 9–18% between corresponding polyA- and polyT-based constructs.

Conclusions

Here we present experimental evidence for utility of GFP reporter for accurate estimation of A/T homopolymeric sequence contribution in transcriptional slippage induction. We showed that the intensity of GFP hybrid fluorescence mainly depends on the number of retained natural amino acids, thus fluorescence raw data need to be referred to appropriate positive control. Moreover, only in case of GFP hybrids with relatively short N-terminal tags the fluorescence level solely reflects production yield, what further indicates the impact of an individual slippage sequence. Our results demonstrate that in contrast to the E. coli enzyme, T7 RNA polymerase exhibits extremely high propensity to slippage even on runs as short as 3 adenine or 4 thymine residues.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0999-3) contains supplementary material, which is available to authorized users.

Stimulation of reverse transcriptase generated cDNAs with specific indels by template RNA structure: retrotransposon, dNTP balance, RT-reagent usage

RNA dependent DNA-polymerases, reverse transcriptases, are key enzymes for retroviruses and retroelements. Their fidelity, including indel generation, is significant for their use as reagents including for deep sequencing. Here, we report that certain RNA template structures and G-rich sequences, ahead of diverse reverse transcriptases can be strong stimulators for slippage at slippage-prone template motif sequence 3′ of such ‘slippage-stimulatory’ structures. Where slippage is stimulated, the resulting products have one or more additional base(s) compared to the corresponding template motif. Such structures also inhibit slippage-mediated base omission which can be more frequent in the absence of a relevant stem–loop. Slippage directionality, base insertion and omission, is sensitive to the relative concentration ratio of dNTPs specified by the RNA template slippage-prone sequence and its 5′ adjacent base. The retrotransposon-derived enzyme TGIRT exhibits more slippage in vitro than the retroviral enzymes tested including that from HIV. Structure-mediated slippage may be exhibited by other polymerases and enrich gene expression. A cassette from Drosophila retrotransposon Dme1_chrX_2630566, a candidate for utilizing slippage for its GagPol synthesis, exhibits strong slippage in vitro. Given the widespread occurrence and importance of retrotransposons, systematic studies to reveal the extent of their functional utilization of RT slippage are merited.

Specific reverse transcriptase slippage at the HIV ribosomal frameshift sequence: potential implications for modulation of GagPol synthesis

Synthesis of HIV GagPol involves a proportion of ribosomes translating a U6A shift site at the distal end of the gag gene performing a programmed -1 ribosomal frameshift event to enter the overlapping pol gene. In vitro studies here show that at the same shift motif HIV reverse transcriptase generates -1 and +1 indels with their ratio being sensitive to the relative concentration ratio of dNTPs specified by the RNA template slippage-prone sequence and its 5' adjacent base. The GGG sequence 3' adjacent to the U6A shift/slippage site, which is important for ribosomal frameshifting, is shown here to limit reverse transcriptase base substitution and indel 'errors' in the run of A's in the product. The indels characterized here have either 1 more or less A, than the corresponding number of template U's. cDNA with 5 A's may yield novel Gag product(s), while cDNA with an extra base, 7 A's, may only be a minor contributor to GagPol polyprotein. Synthesis of a proportion of non-ribosomal frameshift derived GagPol would be relevant in efforts to identify therapeutically useful compounds that perturb the ratio of GagPol to Gag, and pertinent to the extent in which specific polymerase slippage is utilized in gene expression.

Mutational analysis of the Potyviridae transcriptional slippage site utilized for expression of the P3N-PIPO and P1N-PISPO proteins

The Potyviridae comprise the largest and most important family of RNA plant viruses. An essential overlapping ORF, termed pipo, resides in an internal region of the main polyprotein ORF. Recently, expression of pipo was shown to depend on programmed transcriptional slippage at a conserved GAAAAAA sequence, resulting in the insertion of an extra A into a proportion of viral transcripts, fusing the pipo ORF in frame with the 5' third of the polyprotein ORF. However, the sequence features that mediate slippage have not been characterized. Using a duplicate copy of the pipo slip site region fused into a different genomic location where it can be freely mutated, we investigated the sequence requirements for transcriptional slippage. We find that the leading G is not strictly required, but increased flanking sequence GC content correlates with higher insertion rates. A homopolymeric hexamer is optimal for producing mainly single-nucleotide insertions. We also identify an overabundance of G to A substitutions immediately 3'-adjacent to GAAAAAA in insertion-free transcripts, which we infer to result from a 'to-fro' form of slippage during positive-strand synthesis. Analysis of wild-type and reverse complement sequences suggests that slippage occurs preferentially during synthesis of poly(A) and therefore occurs mainly during positive-strand synthesis.

Common position of indels that cause deviations from canonical genome organization in different measles virus strains

BACKGROUND

The canonical genome organization of measles virus (MV) is characterized by total size of 15 894 nucleotides (nts) and defined length of every genomic region, both coding and non-coding. Only rarely have reports of strains possessing non-canonical genomic properties (possessing indels, with or without the change of total genome length) been published. The observed mutations are mutually compensatory in a sense that the total genome length remains polyhexameric. Although programmed and highly precise pseudo-templated nucleotide additions during transcription are inherent to polymerases of all viruses belonging to family Paramyxoviridae, a similar mechanism that would serve to non-randomly correct genome length, if an indel has occurred during replication, has so far not been described in the context of a complete virus genome.

METHODS

We compiled all complete MV genomic sequences (64 in total) available in open access sequence databases. Multiple sequence comparisons and phylogenetic analyses were performed with the aim of exploring whether non-recombinant and non-evolutionary linked measles strains that show deviations from canonical genome organization possess a common genetic characteristic.

RESULTS

In 11 MV sequences we detected deviations from canonical genome organization due to short indels located within homopolymeric stretches or next to them. In nine out of 11 identified non-canonical MV sequences, a common feature was observed: one mutation, either an insertion or a deletion, was located in a 28 nts long region in F gene 5' untranslated region (positions 5051-5078 in genomic cDNA of canonical strains). This segment is composed of five tandemly linked homopolymeric stretches, its consensus sequence is G6-7C7-8A6-7G1-3C5-6. Although none of the mononucleotide repeats within this segment has fixed length, the total number of nts in canonical strains is always 28. These nine non-canonical strains, as well as the tenth (not mutated in 5051-5078 segment), can be grouped in three clusters, based on their passage histories/epidemiological data/genetic similarities. There are no indications that the 3 clusters are evolutionary linked, other than the fact that they all belong to clade D.

CONCLUSIONS

A common narrow genomic region was found to be mutated in different, non-related, wild type strains suggesting that this region might have a function in non-random genome length corrections occurring during MV replication.

Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

Identification and characterization of a -1 reading frameshift in the heavy chain constant region of an IgG1 recombinant monoclonal antibody produced in CHO cells

Frameshifts lead to complete alteration of the intended amino acid sequences, and therefore may affect the biological activities of protein therapeutics and pose potential immunogenicity risks. We report here the identification and characterization of a novel -1 frameshift variant in a recombinant IgG1 therapeutic monoclonal antibody (mAb) produced in Chinese hamster ovary cells during the cell line selection studies. The variant was initially observed as an atypical post-monomer fragment peak in size exclusion chromatography. Characterization of the fragment peak using intact and reduced liquid chromatography-mass spectrometry (LC-MS) analyses determined that the fragment consisted of a normal light chain disulfide-linked to an aberrant 26 kDa fragment that could not be assigned to any HC fragment even after considering common modifications. Further analysis using LC-MS/MS peptide mapping revealed that the aberrant fragment contained the expected HC amino acid sequence (1-232) followed by a 20-mer novel sequence corresponding to expression of heavy chain DNA sequence in the -1 reading frame. Examination of the DNA sequence around the frameshift initiation site revealed that a mononucleotide repeat GGGGGG located in the IgG1 HC constant region was most likely the structural root cause of the frameshift. Rapid identification of the frameshift allowed us to avoid use of a problematic cell line containing the frameshift as the production cell line. The frameshift reported here may be observed in other mAb products and the hypothesis-driven analytical approaches employed here may be valuable for rapid identification and characterization of frameshift variants in other recombinant proteins.

Transcriptional slippage in the positive-sense RNA virus family Potyviridae

The family Potyviridae encompasses ~30% of plant viruses and is responsible for significant economic losses worldwide. Recently, a small overlapping coding sequence, termed pipo, was found to be conserved in the genomes of all potyvirids. PIPO is expressed as part of a frameshift protein, P3N-PIPO, which is essential for virus cell-to-cell movement. However, the frameshift expression mechanism has hitherto remained unknown. Here, we demonstrate that transcriptional slippage, specific to the viral RNA polymerase, results in a population of transcripts with an additional "A" inserted within a highly conserved GAAAAAA sequence, thus enabling expression of P3N-PIPO. The slippage efficiency is ~2% in Turnip mosaic virus and slippage is inhibited by mutations in the GAAAAAA sequence. While utilization of transcriptional slippage is well known in negative-sense RNA viruses such as Ebola, mumps and measles, to our knowledge this is the first report of its widespread utilization for gene expression in positive-sense RNA viruses.

Transcriptional Slippage and RNA Editing Increase the Diversity of Transcripts in Chloroplasts: Insight from Deep Sequencing of Vigna radiata Genome and Transcriptome

We performed deep sequencing of the nuclear and organellar genomes of three mungbean genotypes: Vigna radiata ssp. sublobata TC1966, V. radiata var. radiata NM92 and the recombinant inbred line RIL59 derived from a cross between TC1966 and NM92. Moreover, we performed deep sequencing of the RIL59 transcriptome to investigate transcript variability. The mungbean chloroplast genome has a quadripartite structure including a pair of inverted repeats separated by two single copy regions. A total of 213 simple sequence repeats were identified in the chloroplast genomes of NM92 and RIL59; 78 single nucleotide variants and nine indels were discovered in comparing the chloroplast genomes of TC1966 and NM92. Analysis of the mungbean chloroplast transcriptome revealed mRNAs that were affected by transcriptional slippage and RNA editing. Transcriptional slippage frequency was positively correlated with the length of simple sequence repeats of the mungbean chloroplast genome (R²=0.9911). In total, 41 C-to-U editing sites were found in 23 chloroplast genes and in one intergenic spacer. No editing site that swapped U to C was found. A combination of bioinformatics and experimental methods revealed that the plastid-encoded RNA polymerase-transcribed genes psbF and ndhA are affected by transcriptional slippage in mungbean and in main lineages of land plants, including three dicots (Glycine max, Brassica rapa, and Nicotiana tabacum), two monocots (Oryza sativa and Zea mays), two gymnosperms (Pinus taeda and Ginkgo biloba) and one moss (Physcomitrella patens). Transcript analysis of the rps2 gene showed that transcriptional slippage could affect transcripts at single sequence repeat regions with poly-A runs. It showed that transcriptional slippage together with incomplete RNA editing may cause sequence diversity of transcripts in chloroplasts of land plants.

Productive mRNA stem loop-mediated transcriptional slippage: Crucial features in common with intrinsic terminators

Escherichia coli and yeast DNA-dependent RNA polymerases are shown to mediate efficient nascent transcript stem loop formation-dependent RNA-DNA hybrid realignment. The realignment was discovered on the heteropolymeric sequence T5C5 and yields transcripts lacking a C residue within a corresponding U5C4. The sequence studied is derived from a Roseiflexus insertion sequence (IS) element where the resulting transcriptional slippage is required for transposase synthesis. The stability of the RNA structure, the proximity of the stem loop to the slippage site, the length and composition of the slippage site motif, and the identity of its 3' adjacent nucleotides (nt) are crucial for transcripts lacking a single C. In many respects, the RNA structure requirements for this slippage resemble those for hairpin-dependent transcription termination. In a purified in vitro system, the slippage efficiency ranges from 5% to 75% depending on the concentration ratios of the nucleotides specified by the slippage sequence and the 3' nt context. The only previous proposal of stem loop mediated slippage, which was in Ebola virus expression, was based on incorrect data interpretation. We propose a mechanical slippage model involving the RNAP translocation state as the main motor in slippage directionality and efficiency. It is distinct from previously described models, including the one proposed for paramyxovirus, where following random movement efficiency is mainly dependent on the stability of the new realigned hybrid. In broadening the scope for utilization of transcription slippage for gene expression, the stimulatory structure provides parallels with programmed ribosomal frameshifting at the translation level.

RNA editing by T7 RNA polymerase bypasses InDel mutations causing unexpected phenotypic changes

DNA-dependent T7 RNA polymerase (T7 RNAP) is the most powerful tool for both gene expression and in vitro transcription. By using a Next Generation Sequencing (NGS) approach we have analyzed the polymorphism of a T7 RNAP-generated mRNA pool of the mboIIM2 gene. We find that the enzyme displays a relatively high level of template-dependent transcriptional infidelity. The nucleotide misincorporations and multiple insertions in A/T-rich tracts of homopolymers in mRNA (0.20 and 0.089%, respectively) cause epigenetic effects with significant impact on gene expression that is disproportionally high to their frequency of appearance. The sequence-dependent rescue of single and even double InDel frameshifting mutants and wild-type phenotype recovery is observed as a result. As a consequence, a heterogeneous pool of functional and non-functional proteins of almost the same molecular mass is produced where the proteins are indistinguishable from each other upon ordinary analysis. We suggest that transcriptional infidelity as a general feature of the most effective RNAPs may serve to repair and/or modify a protein function, thus increasing the repertoire of phenotypic variants, which in turn has a high evolutionary potential.

Investigation of a Possible Link Between Vaccination and the 2010 Sheep Pox Epizootic in Morocco

Sheep pox is endemic in most parts of Northern Africa and has the potential to cause severe economic problems. Live attenuated vaccines are used in Morocco, and in many other countries, to control the disease. Sheep pox virus (SPPV) re-appeared in 2010 causing a nodular clinical form previously not observed in Morocco. The severe clinical signs observed during the course of this outbreak and initial reports citing similarity in nucleotide sequence between the Moroccan vaccine strain and field isolates warranted a more in depth analysis of this epizootic. In this study, sequence analysis showed that isolates obtained from four provinces of eastern Morocco were identical, demonstrating that a single SPPV strain was responsible for the 2010 epizootic. In addition, the genome fragments sequenced and phylogenetic analyses undertaken as part of this study showed significant differences between field isolates and the Moroccan vaccine strain. New PCR methods were developed to differentiate between wild-type isolates and vaccine strains of SPPV. Using these methods, no trace of wild-type SPPV was found in the vaccine and no evidence was found to suggest that the vaccine strain was causing clinical disease.

Transcriptional frameshifting rescues Citrobacter rodentium type VI secretion by the production of two length variants from the prematurely interrupted tssM gene

The Type VI secretion system (T6SS) mediates toxin delivery into both eukaryotic and prokaryotic cells. It is composed of a cytoplasmic structure resembling the tail of contractile bacteriophages anchored to the cell envelope through a membrane complex composed of the TssL and TssM inner membrane proteins and of the TssJ outer membrane lipoprotein. The C-terminal domain of TssM is required for its interaction with TssJ, and for the function of the T6SS. In Citrobacter rodentium, the tssM1 gene does not encode the C-terminal domain. However, the stop codon is preceded by a run of 11 consecutive adenosines. In this study, we demonstrate that this poly-A tract is a transcriptional slippery site that induces the incorporation of additional adenosines, leading to frameshifting, and hence the production of two TssM1 variants, including a full-length canonical protein. We show that both forms of TssM1, and the ratio between these two forms, are required for the function of the T6SS in C. rodentium. Finally, we demonstrate that the tssM gene associated with the Yersinia pseudotuberculosis T6SS-3 gene cluster is also subjected to transcriptional frameshifting.

Nonstandard decoding mechanisms lead to the synthesis of different protein variants from a single DNA sequence. These mechanisms are particularly important when the genome length has to be limited such as viral genomes, limited by the available space in the capsid, or to synthesize two different polypeptides that have distinct functional properties. Here, we report that tssM, a gene encoded within the Citrobacter rodentium Type VI secretion (T6S) gene cluster, is interrupted by a premature stop codon; however, the stop codon is preceded by a slippery site constituted by 11 consecutive adenosines. Reiterative transcription leads to the incorporation of additional nucleotides in the mRNA and therefore restores the original framing. As a consequence, two different TssM variants are created by transcriptional frameshifting, including a full-length 130-kDa protein and an 88-kDa truncated variant. We further show that both forms, and the ratio between these two forms, are required for the function of the transport apparatus. Interestingly, a similar mechanism regulates the synthesis of two TssM variants in Yersinia pseudotuberculosis.

The presence of an RNA:DNA hybrid that is prone to slippage promotes termination by T7 RNA polymerase

Intrinsic termination signals for multisubunit bacterial RNA polymerases (RNAPs) encode a GC-rich stem-loop structure followed by a polyuridine [poly(U)] tract, and it has been proposed that steric clash of the stem-loop with the exit pore of the RNAP imposes a shearing force on the RNA in the downstream RNA:DNA hybrid, resulting in misalignment of the active site. The structurally unrelated T7 RNAP terminates at a similar type of signal (TΦ), suggesting a common mechanism for termination. In the absence of a hairpin (passive conditions), T7 RNAP slips efficiently in both homopolymeric A and U tracts, and we have found that replacement of the U tract in TΦ with a slippage-prone A tract still allows efficient termination. Under passive conditions, incorporation of a single G residue following a poly(U) tract (which is the situation during termination at TΦ) results in a "locked" complex that is unable to extend the transcript. Our results support a model in which transmission of the shearing force generated by steric clash of the hairpin with the exit pore is promoted by the presence of a slippery tracts downstream, resulting in alterations in the active site and the formation of a locked complex that represents an early step in the termination pathway.

Analysis of tetra- and hepta-nucleotides motifs promoting -1 ribosomal frameshifting in Escherichia coli

Programmed ribosomal -1 frameshifting is a non-standard decoding process occurring when ribosomes encounter a signal embedded in the mRNA of certain eukaryotic and prokaryotic genes. This signal has a mandatory component, the frameshift motif: it is either a Z_ZZN tetramer or a X_XXZ_ZZN heptamer (where ZZZ and XXX are three identical nucleotides) allowing cognate or near-cognate repairing to the -1 frame of the A site or A and P sites tRNAs. Depending on the signal, the frameshifting frequency can vary over a wide range, from less than 1% to more than 50%. The present study combines experimental and bioinformatics approaches to carry out (i) a systematic analysis of the frameshift propensity of all possible motifs (16 Z_ZZN tetramers and 64 X_XXZ_ZZN heptamers) in Escherichia coli and (ii) the identification of genes potentially using this mode of expression amongst 36 Enterobacteriaceae genomes. While motif efficiency varies widely, a major distinctive rule of bacterial -1 frameshifting is that the most efficient motifs are those allowing cognate re-pairing of the A site tRNA from ZZN to ZZZ. The outcome of the genomic search is a set of 69 gene clusters, 59 of which constitute new candidates for functional utilization of -1 frameshifting.

Bacteriophage λ N protein inhibits transcription slippage by Escherichia coli RNA polymerase

Transcriptional slippage is a class of error in which ribonucleic acid (RNA) polymerase incorporates nucleotides out of register, with respect to the deoxyribonucleic acid (DNA) template. This phenomenon is involved in gene regulation mechanisms and in the development of diverse diseases. The bacteriophage λ N protein reduces transcriptional slippage within actively growing cells and in vitro. N appears to stabilize the RNA/DNA hybrid, particularly at the 5′ end, preventing loss of register between transcript and template. This report provides the first evidence of a protein that directly influences transcriptional slippage, and provides a clue about the molecular mechanism of transcription termination and N-mediated antitermination.

Bacterial insertion sequences: their genomic impact and diversity

Insertion sequences (ISs), arguably the smallest and most numerous autonomous transposable elements (TEs), are important players in shaping their host genomes. This review focuses on prokaryotic ISs. We discuss IS distribution and impact on genome evolution. We also examine their effects on gene expression, especially their role in activating neighbouring genes, a phenomenon of particular importance in the recent upsurge of bacterial antibiotic resistance. We explain how ISs are identified and classified into families by a combination of characteristics including their transposases (Tpases), their overall genetic organisation and the accessory genes which some ISs carry. We then describe the organisation of autonomous and nonautonomous IS‐related elements. This is used to illustrate the growing recognition that the boundaries between different types of mobile element are becoming increasingly difficult to define as more are being identified. We review the known Tpase types, their different catalytic activities used in cleaving and rejoining DNA strands during transposition, their organisation into functional domains and the role of this in regulation. Finally, we consider examples of prokaryotic IS domestication. In a more speculative section, we discuss the necessity of constructing more quantitative dynamic models to fully appreciate the continuing impact of TEs on prokaryotic populations.

Direct assessment of transcription fidelity by high-resolution RNA sequencing

Cancerous and aging cells have long been thought to be impacted by transcription errors that cause genetic and epigenetic changes. Until now, a lack of methodology for directly assessing such errors hindered evaluation of their impact to the cells. We report a high-resolution Illumina RNA-seq method that can assess noncoded base substitutions in mRNA at 10⁻⁴–10⁻⁵ per base frequencies in vitro and in vivo. Statistically reliable detection of changes in transcription fidelity through ∼10³ nt DNA sites assures that the RNA-seq can analyze the fidelity in a large number of the sites where errors occur. A combination of the RNA-seq and biochemical analyses of the positions for the errors revealed two sequence-specific mechanisms that increase transcription fidelity by Escherichia coli RNA polymerase: (i) enhanced suppression of nucleotide misincorporation that improves selectivity for the cognate substrate, and (ii) increased backtracking of the RNA polymerase that decreases a chance of error propagation to the full-length transcript after misincorporation and provides an opportunity to proofread the error. This method is adoptable to a genome-wide assessment of transcription fidelity.

Identification of the nature of reading frame transitions observed in prokaryotic genomes

Our goal was to identify evolutionary conserved frame transitions in protein coding regions and to uncover an underlying functional role of these structural aberrations. We used the ab initio frameshift prediction program, GeneTack, to detect reading frame transitions in 206 991 genes (fs-genes) from 1106 complete prokaryotic genomes. We grouped 102 731 fs-genes into 19 430 clusters based on sequence similarity between protein products (fs-proteins) as well as conservation of predicted position of the frameshift and its direction. We identified 4010 pseudogene clusters and 146 clusters of fs-genes apparently using recoding (local deviation from using standard genetic code) due to possessing specific sequence motifs near frameshift positions. Particularly interesting was finding of a novel type of organization of the dnaX gene, where recoding is required for synthesis of the longer subunit, τ. We selected 20 clusters of predicted recoding candidates and designed a series of genetic constructs with a reporter gene or affinity tag whose expression would require a frameshift event. Expression of the constructs in Escherichia coli demonstrated enrichment of the set of candidates with sequences that trigger genuine programmed ribosomal frameshifting; we have experimentally confirmed four new families of programmed frameshifts.

Isolation and characterization of RNA polymerase rpoB mutations that alter transcription slippage during elongation in Escherichia coli

Transcription fidelity is critical for maintaining the accurate flow of genetic information. The study of transcription fidelity has been limited because the intrinsic error rate of transcription is obscured by the higher error rate of translation, making identification of phenotypes associated with transcription infidelity challenging. Slippage of elongating RNA polymerase (RNAP) on homopolymeric A/T tracts in DNA represents a special type of transcription error leading to disruption of open reading frames in Escherichia coli mRNA. However, the regions in RNAP involved in elongation slippage and its molecular mechanism are unknown. We constructed an A/T tract that is out of frame relative to a downstream lacZ gene on the chromosome to examine transcriptional slippage during elongation. Further, we developed a genetic system that enabled us for the first time to isolate and characterize E. coli RNAP mutants with altered transcriptional slippage in vivo. We identified several amino acid residues in the β subunit of RNAP that affect slippage in vivo and in vitro. Interestingly, these highly clustered residues are located near the RNA strand of the RNA-DNA hybrid in the elongation complex. Our E. coli study complements an accompanying study of slippage by yeast RNAP II and provides the basis for future studies on the mechanism of transcription fidelity.

The fidelity of transcription: RPB1 (RPO21) mutations that increase transcriptional slippage in S. cerevisiae

The fidelity of RNA synthesis depends on both accurate template-mediated nucleotide selection and proper maintenance of register between template and RNA. Loss of register, or transcriptional slippage, is particularly likely on homopolymeric runs in the template. Transcriptional slippage can alter the coding capacity of mRNAs and is used as a regulatory mechanism. Here we describe mutations in the largest subunit of Saccharomyces cerevisiae RNA polymerase II that substantially increase the level of transcriptional slippage. Alleles of RPB1 (RPO21) with elevated slippage rates were identified among 6-azauracil-sensitive mutants and were also isolated using a slippage-dependent reporter gene. Biochemical characterization of polymerase II isolated from these mutants confirms elevated levels of transcriptional slippage.

Hypomorphic glycosyltransferase alleles and recoding at contingency loci influence glycan microheterogeneity in the protein glycosylation system of Neisseria species

As more bacterial protein glycosylation systems are identified and characterized, a central question that arises is, what governs the prevalence of particular glycans associated with them? In addition, accumulating evidence shows that bacterial protein glycans can be subject to the phenomenon of microheterogeneity, in which variant glycan structures are found at specific attachment sites of a given glycoprotein. Although factors underlying microheterogeneity in reconstituted expression systems have been identified and modeled, those impacting natural systems largely remain enigmatic. On the basis of a sensitive and specific glycan serotyping system, microheterogeneity has been reported for the broad-spectrum, O-linked protein glycosylation system in species within the genus Neisseria. To elucidate the mechanisms involved, a genetic approach was used to identify a hypomorphic allele of pglA (encoding the PglA galactosyltransferase) as a significant contributor to simultaneous expression of multiple glycoforms. Moreover, this phenotype was mapped to a single amino acid polymorphism in PglA. Further analyses revealed that many pglA phase-off variants (containing out-of-frame configurations in simple nucleotide repeats within the open reading frame) were associated with disproportionally high levels of the N,N'-diacetylbacillosamine-Gal disaccharide glycoform generated by PglA. This phenotype is emblematic of nonstandard decoding involving programmed ribosomal frameshifting and/or programmed transcriptional realignment. Together, these findings provide new information regarding the mechanisms of neisserial protein glycan microheterogeneity and the anticipatory nature of contingency loci.

A pilot study of bacterial genes with disrupted ORFs reveals a surprising profusion of protein sequence recoding mediated by ribosomal frameshifting and transcriptional realignment

Bacterial genome annotations contain a number of coding sequences (CDSs) that, in spite of reading frame disruptions, encode a single continuous polypeptide. Such disruptions have different origins: sequencing errors, frameshift, or stop codon mutations, as well as instances of utilization of nontriplet decoding. We have extracted over 1,000 CDSs with annotated disruptions and found that about 75% of them can be clustered into 64 groups based on sequence similarity. Analysis of the clusters revealed deep phylogenetic conservation of open reading frame organization as well as the presence of conserved sequence patterns that indicate likely utilization of the nonstandard decoding mechanisms: programmed ribosomal frameshifting (PRF) and programmed transcriptional realignment (PTR). Further enrichment of these clusters with additional homologous nucleotide sequences revealed over 6,000 candidate genes utilizing PRF or PTR. Analysis of the patterns of conservation apparently associated with nontriplet decoding revealed the presence of both previously characterized frameshift-prone sequences and a few novel ones. Since the starting point of our analysis was a set of genes with already annotated disruptions, it is highly plausible that in this study, we have identified only a fraction of all bacterial genes that utilize PRF or PTR. In addition to the identification of a large number of recoded genes, a surprising observation is that nearly half of them are expressed via PTR-a mechanism that, in contrast to PRF, has not yet received substantial attention.

Regulation of gene expression by reiterative transcription

Gene regulation involves many different types of transcription control mechanisms, including mechanisms based on reiterative transcription in which nucleotides are repetitively added to the 3' end of a nascent transcript due to upstream transcript slippage. In these mechanisms, reiterative transcription is typically modulated by interactions between RNA polymerase and its nucleoside triphosphate substrates without the involvement of regulatory proteins. This review describes the current state of knowledge of gene regulation involving reiterative transcription. It focuses on the methods by which reiterative transcription is controlled and emphasizes the different fates of transcripts produced by this reaction. The review also includes a discussion of possible new and fundamentally different mechanisms of gene regulation that rely on conditional reiterative transcription.

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.

DARNED: a DAtabase of RNa EDiting in humans

MOTIVATION

RNA editing is a phenomenon, which is responsible for the alteration of particular nucleotides in RNA sequences relative to their genomic templates. Recently, a large number of RNA editing instances in humans have been identified using bioinformatic screens and high-throughput experimental investigations utilizing next-generation sequencing technologies. However, the available data on RNA editing are not uniform and difficult to access.

RESULTS

Here, we describe a new database DARNED (DAtabase of RNa EDiting) that provides centralized access to available published data related to RNA editing. RNA editing locations are mapped on the reference human genome. The current release of the database contains information on approximately 42 000 human genome coordinates corresponding to RNA locations that undergo RNA editing, mostly involving adenosine-to-inosine (A-to-I) substitutions. The data can be queried using a range of genomic coordinates, their corresponding functional localization in RNA molecules [Exons, Introns, CoDing Sequence (CDS) and UnTranslated Regions (UTRs)] and information regarding tissue/organ/cell sources where RNA editing has been observed. It is also possible to obtain RNA editing information for a specific gene or an RNA molecule using corresponding accession numbers. Search results provide information on the number of expressed sequence tags (ESTs) supporting edited and genomic bases, functional localization of RNA editing and existence of known single nucleotide polymorphisms (SNPs). Editing data can be explored in UCSC and Ensembl genome browsers, in conjunction with additional data provided by these popular genome browsers. DARNED has been designed for researchers seeking information on RNA editing and for the developers of novel algorithms for its prediction.

AVAILABILITY

DARNED is accessible at http://darned.ucc.ie

Simple sequence repeats and genome plasticity in Streptococcus agalactiae

Simple sequence repeats (SSRs) and their role in phase variation have been extensively studied in Gram-negative organisms, where they have been associated with antigenic variation and other adaptation strategies. In this study, we apply comparative genomics in order to find evidence of slipped-strand mispairing in the human Gram-positive pathogen Streptococcus agalactiae. In two consecutive screenings, 2,233 (650 + 1,583) SSRs were identified in our reference genome 2603V/R, and these loci were examined in seven other S. agalactiae genomes. A total of 56 SSR loci were found to exhibit variation, where gain or loss of repeat units was observed in at least one other genome, resulting in aberrant genotypes. Homopolymeric adenine tracts predominated among the repeats that varied. Positional analysis revealed that long polyadenine tracts were overrepresented in the 5' ends of open reading frames (ORFs) and underrepresented in the 3' ends. Repeat clustering in ORFs was also examined, and the highest degree of clustering was observed for a capsule biosynthesis gene and a pilus sortase. A statistical analysis of observed over expected ratios suggested a selective pressure against long homopolymeric tracts. Altered phenotypes were verified for three genes encoding surface-attached proteins, in which frameshifts or fusions led to truncation of proteins and/or affected surface localization through loss or gain of the cell wall sorting signal. The data suggest that SSRs contributes to genome plasticity in S. agalactiae but that the bet-hedging strategy is different from Gram-negative organisms.

Genomic, proteomic, and transcriptomic analysis of virulent and avirulent Rickettsia prowazekii reveals its adaptive mutation capabilities

Rickettsia prowazekii, the agent of epidemic typhus, is an obligate intracellular bacterium that is transmitted to human beings by the body louse. Several strains that differ considerably in virulence are recognized, but the genetic basis for these variations has remained unknown since the initial description of the avirulent vaccine strain nearly 70 yr ago. We use a recently developed murine model of epidemic typhus and transcriptomic, proteomic, and genetic techniques to identify the factors associated with virulence. We identified four phenotypes of R. prowazekii that differed in virulence, associated with the up-regulation of antiapoptotic genes or the interferon I pathway in the host cells. Transcriptional and proteomic analyses of R. prowazekii surface protein expression and protein methylation varied with virulence. By sequencing a virulent strain and using comparative genomics, we found hotspots of mutations in homopolymeric tracts of poly(A) and poly(T) in eight genes in an avirulent strain that split and inactivated these genes. These included recO, putative methyltransferase, and exported protein. Passage of the avirulent Madrid E strain in cells or in experimental animals was associated with a cascade of gene reactivations, beginning with recO, that restored the virulent phenotype. An area of genomic plasticity appears to determine virulence in R. prowazekii and represents an example of adaptive mutation for this pathogen.

Slip into something more functional: selection maintains ancient frameshifts in homopolymeric sequences

Mutational hotspots offer significant sources of genetic variability upon which selection can act. However, with a few notable exceptions, we know little about the dynamics and fitness consequences of mutations in these regions. Here, we explore evolutionary forces shaping homopolymeric tracts that are especially vulnerable to slippage errors during replication and transcription. Such tracts are typically eliminated by selection from most bacterial sequences, yet persist in genomes of endosymbionts with small effective population sizes (N(e)) and biased base compositions. Focusing on Blochmannia, a bacterial endosymbiont of ants, we track the divergence of genes that contain frameshift mutations within long (9-11 bp) polyA or polyT tracts. Earlier experimental work documented that transcriptional slippage restores the reading frame in a fraction of messenger RNA molecules and thereby rescues the function of frameshifted genes. In this study, we demonstrate a surprising persistence of these frameshifts and associated tracts for millions of years. Across the genome of this ant mutualist, rates of indel mutation within homopolymeric tracts far exceed the synonymous mutation rate, indicating that long-term conservation of frameshifts within these tracts is inconsistent with neutrality. In addition, the homopolymeric tracts themselves are more conserved than expected by chance, given extensive neutral substitutions that occur elsewhere in the genes sampled. These data suggest an unexpected role for slippage-prone DNA tracts and highlight a new mechanism for their persistence. That is, when such tracts contain a frameshift, transcriptional slippage plays a critical role in rescuing gene function. In such cases, selection will purge nucleotide changes interrupting the slippery tract so that otherwise volatile sequences become frozen in evolutionary time. Although the advantage of the frameshift itself is less clear, it may offer a mechanism to lower effective gene expression by reducing but not eliminating transcripts that encode full-length proteins.

Recode-2: new design, new search tools, and many more genes

‘Recoding’ is a term used to describe non-standard read-out of the genetic code, and encompasses such phenomena as programmed ribosomal frameshifting, stop codon readthrough, selenocysteine insertion and translational bypassing. Although only a small proportion of genes utilize recoding in protein synthesis, accurate annotation of ‘recoded’ genes lags far behind annotation of ‘standard’ genes. In order to address this issue, provide a service to researchers in the field, and offer training data for developers of gene-annotation software, we have gathered together known cases of recoding within the Recode database. Recode-2 is an improved and updated version of the database. It provides access to detailed information on genes known to utilize translational recoding and allows complex search queries, browsing of recoding data and enhanced visualization of annotated sequence elements. At present, the Recode-2 database stores information on approximately 1500 genes that are known to utilize recoding in their expression—a factor of approximately three increase over the previous version of the database. Recode-2 is available at http://recode.ucc.ie

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.

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.

CapA, an autotransporter protein of Campylobacter jejuni, mediates association with human epithelial cells and colonization of the chicken gut

Two putative autotransporter proteins, CapA and CapB, were identified in silico from the genome sequence of Campylobacter jejuni NCTC11168. The genes encoding each protein contain homopolymeric tracts, suggestive of phase variation mediated by a slipped-strand mispairing mechanism; in each case the gene sequence contained frameshifts at these positions. The C-terminal two-thirds of the two genes, as well as a portion of the predicted signal peptides, were identical; the remaining N-terminal portions were gene specific. Both genes were cloned and expressed; recombinant polypeptides were purified and used to raise rabbit polyclonal monospecific antisera. Using immunoblotting, expression of the ca.116-kDa CapA protein was demonstrated for in vitro-grown cells of strain NCTC11168, for 4 out of 11 recent human fecal isolates, and for 2 out of 8 sequence-typed strains examined. Expression of CapB was not detected for any of the strains tested. Surface localization of CapA was demonstrated by subcellular fractionation and immunogold electron microscopy. Export of CapA was inhibited by globomycin, reinforcing the bioinformatic prediction that the protein is a lipoprotein. A capA insertion mutant had a significantly reduced capacity for association with and invasion of Caco-2 cells and failed to colonize and persist in chickens, indicating that CapA plays a role in host association and colonization by Campylobacter. In view of this demonstrated role, we propose that CapA stands for Campylobacter adhesion protein A.

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.

Mutagenesis via IS transposition in Deinococcus radiodurans

Analysis of the complete genome indicates that insertion sequences (ISs) are abundant in the radio-resistant bacterium Deinococcus radiodurans. By developing a forward mutagenesis assay to detect any inactivation events in D. radiodurans, we found that in the presence of an active mismatch repair system 75% of the mutations to trimethoprim-resistance (Tmp(R)) resulted from an IS insertion into the thyA coding region. Analysis of their distribution among the spontaneous Tmp(R) mutants indicated that five different ISs were transpositionally active. A type II Miniature Inverted-repeat Transposable Element (MITE), related to one of the deinococcal ISs, was also discovered as an insertion into thyA. Seven additional genomic copies of this MITE element were identified by BLASTN. Gamma-ray irradiation of D. radiodurans led to an increase of up to 10-fold in the frequency of Tmp(R) mutants. Analysis of the induced mutations in cells exposed to 10 kGy indicated that gamma-irradiation induced transposition of ISDra2 approximately 100-fold. A 50-fold induction of ISDra2 transposition was also observed in cells exposed to 600 J m(-2) UV-irradiation. Point mutations to rifampicin resistance (Rif(R)) were also induced by gamma-irradiation to reach a plateau at 2 kGy. The plateau value represented a 16-fold increase in the mutant frequency over the background. Although error-free repair strategies predominate in D. radiodurans, an upregulation of transposition, as well as induction of point mutations in cells recovering from DNA damage, provide a genetic variability that may have long-term evolutionary consequences on the fitness of this organism in its habitat.

Recoding in bacteriophages and bacterial IS elements

Dynamic shifts between open reading frames and the redefinition of codon meaning at specific sites, programmed by signals in mRNA, permits versatility of gene expression. Such alterations are characteristic of organisms in all domains of life and serve a variety of functional purposes. In this article, we concentrate on programmed ribosomal frameshifting, stop codon read-through and transcriptional slippage in the decoding of phage genes and bacterial mobile elements. Together with their eukaryotic counterparts, the genes encoding these elements are the richest known source of nonstandard decoding. Recent analyses revealed several novel sequences encoding programmed alterations in gene decoding and provide a glimpse of the emerging picture.

Identification of Hepta- and Octo-Uridine stretches as sole signals for programmed +1 and -1 ribosomal frameshifting during translation of SARS-CoV ORF 3a variants

Programmed frameshifting is one of the translational recoding mechanisms that read the genetic code in alternative ways. This process is generally programmed by signals at defined locations in a specific mRNA. In this study, we report the identification of hepta- and octo-uridine stretches as sole signals for programmed +1 and -1 ribosomal frameshifting during translation of severe acute respiratory syndrome coronavirus (SARS-CoV) ORF 3a variants. SARS-CoV ORF 3a encodes a minor structural protein of 274 amino acids. Over the course of cloning and expression of the gene, a mixed population of clones with six, seven, eight and nine T stretches located 14 nt downstream of the initiation codon was found. In vitro and in vivo expression of clones with six, seven and eight Ts, respectively, showed the detection of the full-length 3a protein. Mutagenesis studies led to the identification of the hepta- and octo-uridine stretches as slippery sequences for efficient frameshifting. Interestingly, no stimulatory elements were found in the sequences upstream or downstream of the slippage site. When the hepta- and octo-uridine stretches were used to replace the original slippery sequence of the SARS-CoV ORF 1a and 1b, efficient frameshift events were observed. Furthermore, the efficiencies of frameshifting mediated by the hepta- and octo-uridine stretches were not affected by mutations introduced into a downstream stem-loop structure that totally abolish the frameshift event mediated by the original slippery sequence of ORF 1a and 1b. Taken together, this study identifies the hepta- and octo-uridine stretches that function as sole elements for efficient +1 and -1 ribosomal frameshift events.