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

Fluorescent reporters give new insights into antibiotics-induced nonsense and frameshift mistranslation

We developed a reporter system based on simultaneous expression of two fluorescent proteins: GFP as a reporter of the capacity of protein synthesis and mutated mScarlet-I as a reporter of translational errors. Because of the unique stop codons or frameshift mutations introduced into the mScarlet-I gene, red fluorescence was produced only after a mistranslation event. These reporters allowed us to estimate mistranslation at a single cell level using either flow cytometry or fluorescence microscopy. We found that laboratory strains of Escherichia coli are more prone to mistranslation compared to the clinical isolates. As relevant for uropathogenic E. coli, growth in human urine elevated translational frameshifting compared to standard laboratory media, whereas different standard media had a small effect on translational fidelity. Antibiotic-induced mistranslation was studied by using amikacin (aminoglycoside family) and azithromycin (macrolide family). Bactericidal amikacin induced preferably stop-codon readthrough at a moderate level. Bacteriostatic azithromycin on the other hand induced both frameshifting and stop-codon readthrough at much higher level. Single cell analysis revealed that fluorescent reporter-protein signal can be lost due to leakage from a fraction of bacteria in the presence of antibiotics, demonstrating the complexity of the antimicrobial activity.

A single-nucleotide resolution capillary gel electrophoresis workflow for poly(A) tail characterization in the development of mRNA therapeutics and vaccines

The 3' poly(A) tail is an important component of messenger RNA (mRNA). The length of the poly(A) tail has direct impact on the stability and translation efficiency of the mRNA molecule and is therefore considered to be a critical quality attribute (CQA) of mRNA-based therapeutics and vaccines. Various analytical methods have been developed to monitor this CQA. Methods like ion-pair reversed-phase liquid chromatography (IPRP-LC) can be used to quantify the percentage of mRNA with poly(A) tail but fail to provide further information on the actual length of poly(A). High-resolution methods such as liquid chromatography coupled with mass spectrometry (LC-MS) or next generation sequencing (NGS) can separate poly(A) tail length by one nucleotide (n/n + 1 resolution) but are complicated to implement for release testing of manufactured mRNA. In this study, a workflow utilizing capillary gel electrophoresis (CGE) for characterizing the poly(A) tail length of mRNA was developed. The CGE method demonstrated resolution comparable with the LC-MS method. With UV detection and the addition of poly(A) length markers, this method can provide poly(A) tail length information and can also provide quantitation of each poly(A) length, making it a suitable release method to monitor the CQA of poly(A) tail length.

Phenotypic mutations contribute to protein diversity and shape protein evolution

Errors in DNA replication generate genetic mutations, while errors in transcription and translation lead to phenotypic mutations. Phenotypic mutations are orders of magnitude more frequent than genetic ones, yet they are less understood. Here, we review the types of phenotypic mutations, their quantifications, and their role in protein evolution and disease. The diversity generated by phenotypic mutation can facilitate adaptive evolution. Indeed, phenotypic mutations, such as ribosomal frameshift and stop codon readthrough, sometimes serve to regulate protein expression and function. Phenotypic mutations have often been linked to fitness decrease and diseases. Thus, understanding the protein heterogeneity and phenotypic diversity caused by phenotypic mutations will advance our understanding of protein evolution and have implications on human health and diseases.

Structural and mechanistic basis of reiterative transcription initiation

Reiterative transcription initiation, observed at promoters that contain homopolymeric sequences at the transcription start site, generates RNA products having 5' sequences noncomplementary to the DNA template. Here, using crystallography and cryoelectron microscopy to define structures, protein-DNA photocrosslinking to map positions of RNAP leading and trailing edges relative to DNA, and single-molecule DNA nanomanipulation to assess RNA polymerase (RNAP)-dependent DNA unwinding, we show that RNA extension in reiterative transcription initiation 1) occurs without DNA scrunching; 2) involves a short, 2- to 3-bp, RNA-DNA hybrid; and 3) generates RNA that exits RNAP through the portal by which scrunched nontemplate-strand DNA exits RNAP in standard transcription initiation. The results establish that, whereas RNA extension in standard transcription initiation proceeds through a scrunching mechanism, RNA extension in reiterative transcription initiation proceeds through a slippage mechanism, with slipping of RNA relative to DNA within a short RNA-DNA hybrid, and with extrusion of RNA from RNAP through an alternative RNA exit.

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.

Chemical biology of non-canonical structures of nucleic acids for therapeutic applications

DNA forms not only the canonical duplex structure but also non-canonical structures. Most potential sequences that induce the formation of non-canonical structures are present in disease-related genes. Interestingly, biological reactions are inhibited or dysregulated by non-canonical structure formation in disease-related genes. To control biological reactions, methods for inducing the formation of non-canonical structures have been developed using small molecules and oligonucleotides. In this feature article, we review biological reactions such as replication, transcription, and reverse transcription controlled by non-canonical DNA structures formed by disease-related genes. Furthermore, we discuss recent studies aimed at developing methods for regulating these biological reactions using drugs targeting the DNA structure.

Structural basis of reiterative transcription from the pyrG and pyrBI promoters by bacterial RNA polymerase

Reiterative transcription is a non-canonical form of RNA synthesis by RNA polymerase in which a ribonucleotide specified by a single base in the DNA template is repetitively added to the nascent RNA transcript. We previously determined the X-ray crystal structure of the bacterial RNA polymerase engaged in reiterative transcription from the pyrG promoter, which contains eight poly-G RNA bases synthesized using three C bases in the DNA as a template and extends RNA without displacement of the promoter recognition σ factor from the core enzyme. In this study, we determined a series of transcript initiation complex structures from the pyrG promoter using soak–trigger–freeze X-ray crystallography. We also performed biochemical assays to monitor template DNA translocation during RNA synthesis from the pyrG promoter and in vitro transcription assays to determine the length of poly-G RNA from the pyrG promoter variants. Our study revealed how RNA slips on template DNA and how RNA polymerase and template DNA determine length of reiterative RNA product. Lastly, we determined a structure of a transcript initiation complex at the pyrBI promoter and proposed an alternative mechanism of RNA slippage and extension requiring the σ dissociation from the core enzyme.

Transcription fidelity: New paradigms in epigenetic inheritance, genome instability and disease

RNA transcription errors are transient, yet frequent, events that do have consequences for the cell. However, until recently we lacked the tools to empirically measure and study these errors. Advances in RNA library preparation and next generation sequencing (NGS) have allowed the spectrum of transcription errors to be empirically measured over the entire transcriptome and in nascent transcripts. Combining these powerful methods with forward and reverse genetic strategies has refined our understanding of transcription factors known to enhance RNA accuracy and will enable the discovery of new candidates. Furthermore, these approaches will shed additional light on the complex interplay between transcription fidelity and other DNA transactions, such as replication and repair, and explore a role for transcription errors in cellular evolution and disease.

Propensity of a picornavirus polymerase to slip on potyvirus-derived transcriptional slippage sites

The substitution rates of viral polymerases have been studied extensively. However less is known about the tendency of these enzymes to 'slip' during RNA synthesis to produce progeny RNAs with nucleotide insertions or deletions. We recently described the functional utilization of programmed polymerase slippage in the family Potyviridae. This slippage results in either an insertion or a substitution, depending on whether the RNA duplex realigns following the insertion. In this study we investigated whether this phenomenon is a conserved feature of superfamily I viral RdRps, by inserting a range of potyvirus-derived slip-prone sequences into a picornavirus, Theiler's murine encephalomyelitis virus (TMEV). Deep-sequencing analysis of viral transcripts indicates that the TMEV polymerase 'slips' at the sequences U6-7 and A6-7 to insert additional nucleotides. Such sequences are under-represented within picornaviral genomes, suggesting that slip-prone sequences create a fitness cost. Nonetheless, the TMEV insertional and substitutional spectrum differed from that previously determined for the potyvirus polymerase.

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.

Frame-Insensitive Expression Cloning of Fluorescent Protein from Scolionema suvaense

Expression cloning from cDNA is an important technique for acquiring genes encoding novel fluorescent proteins. However, the probability of in-frame cDNA insertion following the first start codon of the vector is normally only 1/3, which is a cause of low cloning efficiency. To overcome this issue, we developed a new expression plasmid vector, pRSET-TriEX, in which transcriptional slippage was induced by introducing a DNA sequence of (dT)₁₄ next to the first start codon of pRSET. The effectiveness of frame-insensitive cloning was validated by inserting the gene encoding eGFP with all three possible frames to the vector. After transformation with one of these plasmids, E. coli cells expressed eGFP with no significant difference in the expression level. The pRSET-TriEX vector was then used for expression cloning of a novel fluorescent protein from Scolionema suvaense. We screened 3658 E. coli colonies transformed with pRSET-TriEX containing Scolionema suvaense cDNA, and found one colony expressing a novel green fluorescent protein, ScSuFP. The highest score in protein sequence similarity was 42% with the chain c of multi-domain green fluorescent protein like protein "ember" from Anthoathecata sp. Variations in the N- and/or C-terminal sequence of ScSuFP compared to other fluorescent proteins indicate that the expression cloning, rather than the sequence similarity-based methods, was crucial for acquiring the gene encoding ScSuFP. The absorption maximum was at 498 nm, with an extinction efficiency of 1.17 × 10⁵ M^-1·cm^-1. The emission maximum was at 511 nm and the fluorescence quantum yield was determined to be 0.6. Pseudo-native gel electrophoresis showed that the protein forms obligatory homodimers.

Poly A tail length analysis of in vitro transcribed mRNA by LC-MS

The 3'-polyadenosine (poly A) tail of in vitro transcribed (IVT) mRNA was studied using liquid chromatography coupled to mass spectrometry (LC-MS). Poly A tails were cleaved from the mRNA using ribonuclease T1 followed by isolation with dT magnetic beads. Extracted tails were then analyzed by LC-MS which provided tail length information at single-nucleotide resolution. A 2100-nt mRNA with plasmid-encoded poly A tail lengths of either 27, 64, 100, or 117 nucleotides was used for these studies as enzymatically added poly A tails showed significant length heterogeneity. The number of As observed in the tails closely matched Sanger sequencing results of the DNA template, and even minor plasmid populations with sequence variations were detected. When the plasmid sequence contained a discreet number of poly As in the tail, analysis revealed a distribution that included tails longer than the encoded tail lengths. These observations were consistent with transcriptional slippage of T7 RNAP taking place within a poly A sequence. The type of RNAP did not alter the observed tail distribution, and comparison of T3, T7, and SP6 showed all three RNAPs produced equivalent tail length distributions. The addition of a sequence at the 3' end of the poly A tail did, however, produce narrower tail length distributions which supports a previously described model of slippage where the 3' end can be locked in place by having a G or C after the poly nucleotide region. Graphical abstract Determination of mRNA poly A tail length using magnetic beads and LC-MS.

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.

Genome-Wide Spectra of Transcription Insertions and Deletions Reveal That Slippage Depends on RNA:DNA Hybrid Complementarity

Advances in sequencing technologies have enabled direct quantification of genome-wide errors that occur during RNA transcription. These errors occur at rates that are orders of magnitude higher than rates during DNA replication, but due to technical difficulties such measurements have been limited to single-base substitutions and have not yet quantified the scope of transcription insertions and deletions. Previous reporter gene assay findings suggested that transcription indels are produced exclusively by elongation complex slippage at homopolymeric runs, so we enumerated indels across the protein-coding transcriptomes of Escherichia coli and Buchnera aphidicola, which differ widely in their genomic base compositions and incidence of repeat regions. As anticipated from prior assays, transcription insertions prevailed in homopolymeric runs of A and T; however, transcription deletions arose in much more complex sequences and were rarely associated with homopolymeric runs. By reconstructing the relocated positions of the elongation complex as inferred from the sequences inserted or deleted during transcription, we show that continuation of transcription after slippage hinges on the degree of nucleotide complementarity within the RNA:DNA hybrid at the new DNA template location.

The high level of mistakes generated during transcription can result in the accumulation of malfunctioning and misfolded proteins which can alter global gene regulation and in the expenditure of energy to degrade these nonfunctional proteins. The transcriptome-wide occurrence of base substitutions has been elucidated in bacteria, but information on transcription insertions and deletions—errors that potentially have more dire effects on protein function—is limited to reporter gene constructs. Here, we capture the transcriptome-wide spectrum of insertions and deletions in Escherichia coli and Buchnera aphidicola and show that they occur at rates approaching those of base substitutions. Knowledge of the full extent of sequences subject to transcription indels supports a new model of bacterial transcription slippage, one that relies on the number of complementary bases between the transcript and the DNA template to which it slipped.

X-ray crystal structure of a reiterative transcription complex reveals an atypical RNA extension pathway

Reiterative transcription is a noncanonical form of RNA synthesis in which a nucleotide specified by a single base in the DNA template is repetitively added to the nascent transcript. Here we determined the crystal structure of an RNA polymerase, the bacterial enzyme from Thermus thermophilus, engaged in reiterative transcription during transcription initiation at a promoter resembling the pyrG promoter of Bacillus subtilis The structure reveals that the reiterative transcript detours from the dedicated RNA exit channel and extends toward the main channel of the enzyme, thereby allowing RNA extension without displacement of the promoter recognition σ-factor. Nascent transcripts containing reiteratively added G residues are eventually extended by nonreiterative transcription, revealing an atypical pathway for the formation of a transcription elongation complex.

Recurrent emergence of structural variants of LTR retrotransposon CsRn1 evolving novel expression strategy and their selective expansion in a carcinogenic liver fluke, Clonorchis sinensis

Autonomous retrotransposons, in which replication and transcription are coupled, encode the essential gag and pol genes as a fusion or separate overlapping form(s) that are expressed in single transcripts regulated by a common upstream promoter. The element-specific expression strategies have driven development of relevant translational recoding mechanisms including ribosomal frameshifting to satisfy the protein stoichiometry critical for the assembly of infectious virus-like particles. Retrotransposons with different recoding strategies exhibit a mosaic distribution pattern across the diverse families of reverse transcribing elements, even though their respective distributions are substantially skewed towards certain family groups. However, only a few investigations to date have focused on the emergence of retrotransposons evolving novel expression strategy and causal genetic drivers of the structural variants. In this study, the bulk of genomic and transcribed sequences of a Ty3/gypsy-like CsRn1 retrotransposon in Clonorchis sinensis were analyzed for the comprehensive examination of its expression strategy. Our results demonstrated that structural variants with single open reading frame (ORF) have recurrently emerged from precedential CsRn1 copies encoding overlapping gag-pol ORFs by a single-nucleotide insertion in an upstream region of gag stop codon. In the parasite genome, some of the newly evolved variants appeared to undergo proliferative burst as active master lineages together with their ancestral copies. The genetic event was similarly observed in Opisthorchis viverrini, the closest neighbor of C. sinensis, whereas the resulting structural variants might have failed to overcome purifying selection and comprised minor remnant copies in the Opisthorchis genome.

Genetic Predisposition To Acquire a Polybasic Cleavage Site for Highly Pathogenic Avian Influenza Virus Hemagglutinin

Highly pathogenic avian influenza viruses with H5 and H7 hemagglutinin (HA) subtypes evolve from low-pathogenic precursors through the acquisition of multiple basic amino acid residues at the HA cleavage site. Although this mechanism has been observed to occur naturally only in these HA subtypes, little is known about the genetic basis for the acquisition of the polybasic HA cleavage site. Here we show that consecutive adenine residues and a stem-loop structure, which are frequently found in the viral RNA region encoding amino acids around the cleavage site of low-pathogenic H5 and H7 viruses isolated from waterfowl reservoirs, are important for nucleotide insertions into this RNA region. A reporter assay to detect nontemplated nucleotide insertions and deep-sequencing analysis of viral RNAs revealed that an increased number of adenine residues and enlarged stem-loop structure in the RNA region accelerated the multiple adenine and/or guanine insertions required to create codons for basic amino acids. Interestingly, nucleotide insertions associated with the HA cleavage site motif were not observed principally in the viral RNA of other subtypes tested (H1, H2, H3, and H4). Our findings suggest that the RNA editing-like activity is the key mechanism for nucleotide insertions, providing a clue as to why the acquisition of the polybasic HA cleavage site is restricted to the particular HA subtypes.

Influenza A viruses are divided into subtypes based on the antigenicity of the viral surface glycoproteins hemagglutinin (HA) and neuraminidase. Of the 16 HA subtypes (H1 to -16) maintained in waterfowl reservoirs of influenza A viruses, H5 and H7 viruses often become highly pathogenic through the acquisition of multiple basic amino acid residues at the HA cleavage site. Although this mechanism has been known since the 1980s, the genetic basis for nucleotide insertions has remained unclear. This study shows the potential role of the viral RNA secondary structure for nucleotide insertions and demonstrates a key mechanism explaining why the acquisition of the polybasic HA cleavage site is restricted to particular HA subtypes in nature. Our findings will contribute to better understanding of the ecology of influenza A viruses and will also be useful for the development of genetically modified vaccines against H5 and H7 influenza A viruses with increased stability.

In-silico prediction and deep-DNA sequencing validation indicate phase variation in 115 Neisseria meningitidis genes

BACKGROUND

The Neisseria meningitidis (Nm) chromosome shows a high abundance of simple sequence DNA repeats (SSRs) that undergo stochastic, reversible mutations at high frequency. This mechanism is reflected in an extensive phenotypic diversity that facilitates Nm adaptation to dynamic environmental changes. To date, phase-variable phenotypes mediated by SSRs variation have been experimentally confirmed for 26 Nm genes.

RESULTS

Here we present a population-scale comparative genomic analysis that identified 277 genes and classified them into 52 strong, 60 moderate and 165 weak candidates for phase variation. Deep-coverage DNA sequencing of single colonies grown overnight under non-selective conditions confirmed the presence of high-frequency, stochastic variation in 115 of them, providing circumstantial evidence for their phase variability. We confirmed previous observations of a predominance of variable SSRs within genes for components located on the cell surface or DNA metabolism. However, in addition we identified an unexpectedly broad spectrum of other metabolic functions, and most of the variable SSRs were predicted to induce phenotypic changes by modulating gene expression at a transcriptional level or by producing different protein isoforms rather than mediating on/off translational switching through frameshifts. Investigation of the evolutionary history of SSR contingency loci revealed that these loci were inherited from a Nm ancestor, evolved independently within Nm, or were acquired by Nm through lateral DNA exchange.

CONCLUSIONS

Overall, our results have identified a broader and qualitatively different phenotypic diversification of SSRs-mediated stochastic variation than previously documented, including its impact on central Nm metabolism.

Genome-wide In Silico Analysis, Characterization and Identification of Microsatellites in Spodoptera littoralis Multiple nucleopolyhedrovirus (SpliMNPV)

In this study, we undertook a survey to analyze the distribution and frequency of microsatellites or Simple Sequence Repeats (SSRs) in Spodoptera littoralis multiple nucleopolyhedrovirus (SpliMNPV) genome (isolate AN-1956). Out of the 55 microsatellite motifs, identified in the SpliMNPV-AN1956 genome using in silico analysis (inclusive of mono-, di-, tri- and hexa-nucleotide repeats), 39 were found to be distributed within coding regions (cSSRs), whereas 16 were observed to lie within intergenic or noncoding regions. Among the 39 motifs located in coding regions, 21 were located in annotated functional genes whilst 18 were identified in unknown functional genes (hypothetical proteins). Among the identified motifs, trinucleotide (80%) repeats were found to be the most abundant followed by dinucleotide (13%), mononucleotide (5%) and hexanucleotide (2%) repeats. The 39 motifs located within coding regions were further validated in vitro by using PCR analysis, while the 21 motifs located within known functional genes (15 genes) were characterized using nucleotide sequencing. A comparison of the sequence analysis data of the 21 sequenced cSSRs with the published sequences is presented. Finally, the developed SSR markers of the 39 motifs were further mapped/localized onto the SpliMNPV-AN1956 genome. In conclusion, the SSR markers specific to SpliMNPV, developed in this study, could be a useful tool for the identification of isolates and analysis of genetic diversity and viral evolutionary status.

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.

Truncated yet functional viral protein produced via RNA polymerase slippage implies underestimated coding capacity of RNA viruses

RNA viruses use various strategies to condense their genetic information into small genomes. Potyviruses not only use the polyprotein strategy, but also embed an open reading frame, pipo, in the P3 cistron in the -1 reading frame. PIPO is expressed as a fusion protein with the N-terminal half of P3 (P3N-PIPO) via transcriptional slippage of viral RNA-dependent RNA polymerase (RdRp). We herein show that clover yellow vein virus (ClYVV) produces a previously unidentified factor, P3N-ALT, in the +1 reading frame via transcriptional slippage at a conserved G(1-2)A(6-7) motif, as is the case for P3N-PIPO. The translation of P3N-ALT terminates soon, and it is considered to be a C-terminal truncated form of P3. In planta experiments indicate that P3N-ALT functions in cell-to-cell movement along with P3N-PIPO. Hence, all three reading frames are used to produce functional proteins. Deep sequencing of ClYVV RNA from infected plants endorses the slippage by viral RdRp. Our findings unveil a virus strategy that optimizes the coding capacity.

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.

DksA involvement in transcription fidelity buffers stochastic epigenetic change

DksA is an auxiliary transcription factor that interacts with RNA polymerase and influences gene expression. Depending on the promoter, DksA can be a positive or negative regulator of transcription initiation. Moreover, DksA has a substantial effect on transcription elongation where it prevents the collision of transcription and replication machineries, plays a key role in maintaining transcription elongation when translation and transcription are uncoupled and has been shown to be involved in transcription fidelity. Here, we assessed the role of DksA in transcription fidelity by monitoring stochastic epigenetic switching in the lac operon (with and without an error-prone transcription slippage sequence), partial phenotypic suppression of a lacZ nonsense allele, as well as monitoring the number of lacI mRNA transcripts produced in the presence and absence of DksA via an operon fusion and single molecule fluorescent in situ hybridization studies. We present data showing that DksA acts to maintain transcription fidelity in vivo and the role of DksA seems to be distinct from that of the GreA and GreB transcription fidelity factors.

The Evolutionary Potential of Phenotypic Mutations

Errors in protein synthesis, so-called phenotypic mutations, are orders-of-magnitude more frequent than genetic mutations. Here, we provide direct evidence that alternative protein forms and phenotypic variability derived from translational errors paved the path to genetic, evolutionary adaptations via gene duplication. We explored the evolutionary origins of Saccharomyces cerevisiae IDP3 - an NADP-dependent isocitrate dehydrogenase mediating fatty acids ß-oxidation in the peroxisome. Following the yeast whole genome duplication, IDP3 diverged from a cytosolic ancestral gene by acquisition of a C-terminal peroxisomal targeting signal. We discovered that the pre-duplicated cytosolic IDPs are partially localized to the peroxisome owing to +1 translational frameshifts that bypass the stop codon and unveil cryptic peroxisomal targeting signals within the 3'-UTR. Exploring putative cryptic signals in all 3'-UTRs of yeast genomes, we found that other enzymes related to NADPH production such as pyruvate carboxylase 1 (PYC1) might be prone to peroxisomal localization via cryptic signals. Using laboratory evolution we found that these translational frameshifts are rapidly imprinted via genetic single base deletions occurring within the very same gene location. Further, as exemplified here, the sequences that promote translational frameshifts are also more prone to genetic deletions. Thus, genotypes conferring higher phenotypic variability not only meet immediate challenges by unveiling cryptic 3'-UTR sequences, but also boost the potential for future genetic adaptations.

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.

Genetic evidence that both dNTP-stabilized and strand slippage mechanisms may dictate DNA polymerase errors within mononucleotide microsatellites

Mononucleotide microsatellites are tandem repeats of a single base pair, abundant within coding exons and frequent sites of mutation in the human genome. Because the repeated unit is one base pair, multiple mechanisms of insertion/deletion (indel) mutagenesis are possible, including strand-slippage, dNTP-stabilized, and misincorportion-misalignment. Here, we examine the effects of polymerase identity (mammalian Pols α, β, κ, and η), template sequence, dNTP pool size, and reaction temperature on indel errors during in vitro synthesis of mononucleotide microsatellites. We utilized the ratio of insertion to deletion errors as a genetic indicator of mechanism. Strikingly, we observed a statistically significant bias toward deletion errors within mononucleotide repeats for the majority of the 28 DNA template and polymerase combinations examined, with notable exceptions based on sequence and polymerase identity. Using mutator forms of Pol β did not substantially alter the error specificity, suggesting that mispairing-misalignment mechanism is not a primary mechanism. Based on our results for mammalian DNA polymerases representing three structurally distinct families, we suggest that dNTP-stabilized mutagenesis may be an alternative mechanism for mononucleotide microsatellite indel mutation. The change from a predominantly dNTP-stabilized mechanism to a strand-slippage mechanism with increasing microsatellite length may account for the differential rates of tandem repeat mutation that are observed genome-wide.

Lost in transcription: transient errors in information transfer

Errors in information transfer from DNA to RNA to protein are inevitable. Here, we focus on errors that occur in nascent transcripts during transcription, epimutations. Recent approaches using novel cDNA library preparation and next-generation sequencing begin to directly determine the rate of epimutation and allow analysis of the epimutational spectrum of transcription errors, the type and sequence context of the errors produced in a transcript by an RNA polymerase. The phenotypic consequences of transcription errors have been assessed using both forward and reverse epimutation systems. These studies reveal that transient transcription errors can produce a modification of cell phenotype, partial phenotypic suppression of a mutant allele, and a heritable change in cell phenotype, epigenetic switching in a bistable gene network.

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.

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.

New insights into transcription fidelity: thermal stability of non-canonical structures in template DNA regulates transcriptional arrest, pause, and slippage

The thermal stability and topology of non-canonical structures of G-quadruplexes and hairpins in template DNA were investigated, and the effect of non-canonical structures on transcription fidelity was evaluated quantitatively. We designed ten template DNAs: A linear sequence that does not have significant higher-order structure, three sequences that form hairpin structures, and six sequences that form G-quadruplex structures with different stabilities. Templates with non-canonical structures induced the production of an arrested, a slipped, and a full-length transcript, whereas the linear sequence produced only a full-length transcript. The efficiency of production for run-off transcripts (full-length and slipped transcripts) from templates that formed the non-canonical structures was lower than that from the linear. G-quadruplex structures were more effective inhibitors of full-length product formation than were hairpin structure even when the stability of the G-quadruplex in an aqueous solution was the same as that of the hairpin. We considered that intra-polymerase conditions may differentially affect the stability of non-canonical structures. The values of transcription efficiencies of run-off or arrest transcripts were correlated with stabilities of non-canonical structures in the intra-polymerase condition mimicked by 20 wt% polyethylene glycol (PEG). Transcriptional arrest was induced when the stability of the G-quadruplex structure (−ΔG^o₃₇) in the presence of 20 wt% PEG was more than 8.2 kcal mol⁻¹. Thus, values of stability in the presence of 20 wt% PEG are an important indicator of transcription perturbation. Our results further our understanding of the impact of template structure on the transcription process and may guide logical design of transcription-regulating drugs.

Disproportional plastome-wide increase of substitution rates and relaxed purifying selection in genes of carnivorous Lentibulariaceae

Carnivorous Lentibulariaceae exhibit the most sophisticated implementation of the carnivorous syndrome in plants. Their unusual lifestyle coincides with distinct genomic peculiarities such as the smallest angiosperm nuclear genomes and extremely high nucleotide substitution rates across all genomic compartments. Here, we report the complete plastid genomes from each of the three genera Pinguicula, Utricularia, and Genlisea, and investigate plastome-wide changes in their molecular evolution as the carnivorous syndrome unfolds. We observe a size reduction by up to 9% mostly due to the independent loss of genes for the plastid NAD(P)H dehydrogenase and altered proportions of plastid repeat DNA, as well as a significant plastome-wide increase of substitution rates and microstructural changes. Protein-coding genes across all gene classes show a disproportional elevation of nonsynonymous substitutions, particularly in Utricularia and Genlisea. Significant relaxation of purifying selection relative to noncarnivores occurs in the plastid-encoded fraction of the photosynthesis ATP synthase complex, the photosystem I, and in several other photosynthesis and metabolic genes. Shifts in selective regimes also affect housekeeping genes including the plastid-encoded polymerase, for which evidence for relaxed purifying selection was found once during the transition to carnivory, and a second time during the diversification of the family. Lentibulariaceae significantly exhibit enhanced rates of nucleotide substitution in most of the 130 noncoding regions. Various factors may underlie the observed patterns of relaxation of purifying selection and substitution rate increases, such as reduced net photosynthesis rates, alternative paths of nutrient uptake (including organic carbon), and impaired DNA repair mechanisms.

Correlated occurrence and bypass of frame-shifting insertion-deletions (InDels) to give functional proteins

Short insertions and deletions (InDels) comprise an important part of the natural mutational repertoire. InDels are, however, highly deleterious, primarily because two-thirds result in frame-shifts. Bypass through slippage over homonucleotide repeats by transcriptional and/or translational infidelity is known to occur sporadically. However, the overall frequency of bypass and its relation to sequence composition remain unclear. Intriguingly, the occurrence of InDels and the bypass of frame-shifts are mechanistically related - occurring through slippage over repeats by DNA or RNA polymerases, or by the ribosome, respectively. Here, we show that the frequency of frame-shifting InDels, and the frequency by which they are bypassed to give full-length, functional proteins, are indeed highly correlated. Using a laboratory genetic drift, we have exhaustively mapped all InDels that occurred within a single gene. We thus compared the naive InDel repertoire that results from DNA polymerase slippage to the frame-shifting InDels tolerated following selection to maintain protein function. We found that InDels repeatedly occurred, and were bypassed, within homonucleotide repeats of 3–8 bases. The longer the repeat, the higher was the frequency of InDels formation, and the more frequent was their bypass. Besides an expected 8A repeat, other types of repeats, including short ones, and G and C repeats, were bypassed. Although obtained in vitro, our results indicate a direct link between the genetic occurrence of InDels and their phenotypic rescue, thus suggesting a potential role for frame-shifting InDels as bridging evolutionary intermediates.

Homonucleotide repeats are exceptionally prone to insertions and/or deletions of bases (InDels). However, unless they occur in a multiplicity of 3 bases, InDels disrupt the reading frame and are thus expected to be purged from coding regions. Homonucleotide repeats, however, are also vulnerable to slippage by RNA polymerases and the ribosome. Using laboratory evolution techniques, we systematically mapped the occurrence of InDels within a given gene, before and after selection. Our data indicate that frame-shifting InDels were frequently bypassed to give functional proteins at surprisingly high frequencies. Further, we found a strict correlation between the repeat length, the frequency of occurrence of InDels at the DNA level, and the likelihood of bypass by transcriptional/translational slippage. Our results suggest that frame-shifting InDels might comprise functional evolutionary intermediates, and an effective mean of sequence divergence (e.g. when an adjacent InDel restores the frame, resulting in altered sequence and, potentially, in an altered protein structure).

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.

Heritable change caused by transient transcription errors

Transmission of cellular identity relies on the faithful transfer of information from the mother to the daughter cell. This process includes accurate replication of the DNA, but also the correct propagation of regulatory programs responsible for cellular identity. Errors in DNA replication (mutations) and protein conformation (prions) can trigger stable phenotypic changes and cause human disease, yet the ability of transient transcriptional errors to produce heritable phenotypic change ('epimutations') remains an open question. Here, we demonstrate that transcriptional errors made specifically in the mRNA encoding a transcription factor can promote heritable phenotypic change by reprogramming a transcriptional network, without altering DNA. We have harnessed the classical bistable switch in the lac operon, a memory-module, to capture the consequences of transient transcription errors in living Escherichia coli cells. We engineered an error-prone transcription sequence (A9 run) in the gene encoding the lac repressor and show that this 'slippery' sequence directly increases epigenetic switching, not mutation in the cell population. Therefore, one altered transcript within a multi-generational series of many error-free transcripts can cause long-term phenotypic consequences. Thus, like DNA mutations, transcriptional epimutations can instigate heritable changes that increase phenotypic diversity, which drives both evolution and disease.

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.

Multiple mechanisms contribute to leakiness of a frameshift mutation in canine cone-rod dystrophy

Mutations in RPGRIP1 are associated with early onset retinal degenerations in humans and dogs. Dogs homozygous for a 44 bp insertion including a polyA₂₉ tract potentially leading to premature truncation of the protein, show cone rod degeneration. This is rapid and blinding in a colony of dogs in which the mutation was characterised but in dogs with the same mutation in the pet population there is very variable disease severity and rate of progression.

Objective

We hypothesized that this variability must be associated with leakiness of the RPGRIP1 mutation, allowing continued RPGRIP1 production. The study was designed to discover mechanisms that might allow such leakiness.

Methods

We analysed alternate start sites and splicing of RPGRIP1 transcripts; variability of polyA_n length in the insertion and slippage at polyA_n during transcription/translation.

Results and Significance

We observed a low rate of use of alternative start codons having potential to allow forms of transcript not including the insertion, with the possibility of encoding truncated functional RPGRIP1 protein isoforms. Complex alternative splicing was observed, but did not increase this potential. Variable polyA_n length was confirmed in DNA from different RPGRIP1^−/− dogs, yet polyA_n variability did not correspond with the clinical phenotypes and no individual was found that carried a polyA_n tract capable of encoding an in-frame variant. Remarkably though, in luciferase reporter gene assays, out-of-frame inserts still allowed downstream reporter gene expression at some 40% of the efficiency of in-frame controls. This indicates a major role of transcriptional or translational frameshifting in RPGRIP1 expression. The known slippage of reverse transcriptases as well as RNA polymerases and thermostable DNA polymerases on oligoA homopolymers meant that we could not distinguish whether the majority of slippage was transcriptional or translational. This leakiness at the mutation site may allow escape from severe effects of the mutation for some dogs.

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.

Promoter swapping unveils the role of the Citrobacter rodentium CTS1 type VI secretion system in interbacterial competition

The type VI secretion system (T6SS) is a versatile secretion machine dedicated to various functions in Gram-negative bacteria, including virulence toward eukaryotic cells and antibacterial activity. Activity of T6SS might be followed in vitro by the release of two proteins, Hcp and VgrG, in the culture supernatant. Citrobacter rodentium, a rodent pathogen, harbors two T6SS gene clusters, cts1 and cts2. Reporter fusion and Hcp release assays suggested that the CTS1 T6SS was not produced or not active. The cts1 locus is composed of two divergent operons. We therefore developed a new vector allowing us to swap the two divergent endogenous promoters by P(tac) and P(BAD) using the λ red recombination technology. Artificial induction of both promoters demonstrated that the CTS1 T6SS is functional as shown by the Hcp release assay and confers on C. rodentium a growth advantage in antibacterial competition experiments with Escherichia coli.

Isolation and characterization of transcription fidelity mutants

Accurate transcription is an essential step in maintaining genetic information. Error-prone transcription has been proposed to contribute to cancer, aging, adaptive mutagenesis, and mutagenic evolution of retroviruses and retrotransposons. The mechanisms controlling transcription fidelity and the biological consequences of transcription errors are poorly understood. Because of the transient nature of mRNAs and the lack of reliable experimental systems, the identification and characterization of defects that increase transcription errors have been particularly challenging. In this review we describe novel genetic screens for the isolation of fidelity mutants in both Saccharomyces cerevisiae and Escherichia coli RNA polymerases. We obtained and characterized two distinct classes of mutants altering NTP misincorporation and transcription slippage both in vivo and in vitro. Our study not only validates the genetic schemes for the isolation of RNA polymerase mutants that alter fidelity, but also sheds light on the mechanism of transcription accuracy. This article is part of a Special Issue entitled: Chromatin in time and space.

Error prevention and mitigation as forces in the evolution of genes and genomes

Why are short introns rarely a multiple of three nucleotides long? Why do essential genes cluster? Why are genes in operons often lined up in the order in which they are needed in the encoded pathway? In this Opinion article, we argue that these and many other - ostensibly disparate - observations are all pieces of an emerging picture in which multiple aspects of gene anatomy and genome architecture have evolved in response to error-prone gene expression.

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.

Poly(T) variation in heteroderid nematode mitochondrial genomes is predominantly an artefact of amplification

We assessed the rate of in vitro polymerase errors at polythymidine [poly(T)] tracts in the mitochondrial DNA (mtDNA) of a heteroderid nematode (Heterodera cajani). The mtDNA of these nematodes contain unusually high numbers of poly(T) tracts, and have previously been suggested to contain biological poly(T) length variation. However, using a cloned molecule, we observed that poly(T) variation was generated in vitro at regions containing more than six consecutive Ts. This artefactual error rate was estimated at 7.3 × 10(-5) indels/poly(T) tract >6 Ts/cycle. This rate was then compared to the rate of poly(T) variation detected after the amplification of a biological sample, in order to estimate the 'biological + artefactual' rate of poly(T) variation. There was no significant difference between the artefactual and the artefactual + biological rates, suggesting that the majority of poly(T) variation in the biological sample was artefactual. We then examined the generation of poly(T) variation in a range of templates with tracts up to 16 Ts long, utilizing a range of Heteroderidae species. We observed that T deletions occurred five times more frequently than insertions, and a trend towards increasing error rates with increasing poly(T) tract length. These findings have significant implications for studies involving genomes with many homopolymer tracts.