RNA polymerase errors cause splicing defects and can be regulated by differential expression of RNA polymerase subunits

Characterization of Hepatitis B Virus Transcripts in Chronically HBV-Infected Chimpanzees and Patients Treated with ARC-520 siRNA Demonstrates Transcriptional Silencing of cccDNA

Full-length hepatitis B virus (HBV) transcripts of chimpanzees and patients treated with multidose (MD) HBV siRNA ARC-520 and entecavir (ETV) were characterized by single-molecule real-time (SMRT) sequencing, identifying multiple types of transcripts with the potential to encode HBx, HBsAg, HBeAg, core, and polymerase, as well as transcripts likely to be derived from dimers of dslDNA, and these differed between HBeAg-positive (HBeAg+) and HBeAg-negative (HBeAg-) individuals. HBV transcripts from the last follow-up ~30 months post-ARC-520 treatment were categorized from one HBeAg+ (one of two previously highly viremic patients that became HBeAg- upon treatment and had greatly reduced cccDNA products) and four HBeAg- patients. The previously HBeAg+ patient received a biopsy that revealed that he had 3.4 copies/cell cccDNA (two to three orders of magnitude more cccDNA than HBeAg- chimpanzees) but expressed primarily truncated X and HBsAg from iDNA, like two patients that were HBeAg- at the start of the study and had one copy/cell cccDNA. No HBV transcripts were detected in two other HBeAg- patients that had ~0.3 copies/cell cccDNA, one of which had seroconverted for HBsAg. The paucity of cccDNA-derived transcripts in the presence of high cccDNA demonstrates the transcriptional silencing of HBV following MD siRNA treatment with ETV.

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G·T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R1ρ experiments, we show that dG·rU and dT·rG mismatches in two RNA:DNA hybrids transiently form tautomeric (Genol·T/U $ \mathbin{\lower.3ex\hbox{$\buildrel\textstyle\rightarrow\over {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}}$}}$ G·Tenol/Uenol) and anionic (G·T-/U-) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG·rU- formed with a ten-fold higher propensity relative to dT-·rG and dG·dT- and this could be attributed to the lower pKa (ΔpKa ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G·T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G·U versus G·T, and indicate that anionic Watson-Crick-like G·U- could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Digital RNA sequencing using unique molecular identifiers enables ultrasensitive RNA mutation analysis

Mutation analysis is typically performed at the DNA level since most technical approaches are developed for DNA analysis. However, some applications, like transcriptional mutagenesis, RNA editing and gene expression analysis, require RNA analysis. Here, we combine reverse transcription and digital DNA sequencing to enable low error digital RNA sequencing. We evaluate yield, reproducibility, dynamic range and error correction rate for seven different reverse transcription conditions using multiplexed assays. The yield, reproducibility and error rate vary substantially between the specific conditions, where the yield differs 9.9-fold between the best and worst performing condition. Next, we show that error rates similar to DNA sequencing can be achieved for RNA using appropriate reverse transcription conditions, enabling detection of mutant allele frequencies <0.1% at RNA level. We also detect mutations at both DNA and RNA levels in tumor tissue using a breast cancer panel. Finally, we demonstrate that digital RNA sequencing can be applied to liquid biopsies, analyzing cell-free gene transcripts. In conclusion, we demonstrate that digital RNA sequencing is suitable for ultrasensitive RNA mutation analysis, enabling several basic research and clinical applications.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: Implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G•T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R 1ρ experiments, we show that dG•rU and dT•rG mismatches in two RNA:DNA hybrids transiently form tautomeric (G enol •T/U ⇄G•T enol /U enol ) and anionic (G•T - /U - ) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG•rU - formed with a ten-fold higher propensity relative to dT - •rG and dG•dT - and this could be attributed to the lower pK a (Δ pK a ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G•T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G•U versus G•T, and indicate that anionic Watson-Crick-like G•U - could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

When Size Really Matters: The Eccentricities of Dystrophin Transcription and the Hazards of Quantifying mRNA from Very Long Genes

At 2.3 megabases in length, the dystrophin gene is enormous: transcription of a single mRNA requires approximately 16 h. Principally expressed in skeletal muscle, the dystrophin protein product protects the muscle sarcolemma against contraction-induced injury, and dystrophin deficiency results in the fatal muscle-wasting disease, Duchenne muscular dystrophy. This gene is thus of key clinical interest, and therapeutic strategies aimed at eliciting dystrophin restoration require quantitative analysis of its expression. Approaches for quantifying dystrophin at the protein level are well-established, however study at the mRNA level warrants closer scrutiny: measured expression values differ in a sequence-dependent fashion, with significant consequences for data interpretation. In this manuscript, we discuss these nuances of expression and present evidence to support a transcriptional model whereby the long transcription time is coupled to a short mature mRNA half-life, with dystrophin transcripts being predominantly nascent as a consequence. We explore the effects of such a model on cellular transcriptional dynamics and then discuss key implications for the study of dystrophin gene expression, focusing on both conventional (qPCR) and next-gen (RNAseq) approaches.

Computational Analysis Predicts Correlations among Amino Acids in SARS-CoV-2 Proteomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a serious global challenge requiring urgent and permanent therapeutic solutions. These solutions can only be engineered if the patterns and rate of mutations of the virus can be elucidated. Predicting mutations and the structure of proteins based on these mutations have become necessary for early drug and vaccine design purposes in anticipation of future viral mutations. The amino acid composition (AAC) of proteomes and individual viral proteins provide avenues for exploitation since AACs have been previously used to predict structure, shape and evolutionary rates. Herein, the frequency of amino acid residues found in 1637 complete proteomes belonging to 11 SARS-CoV-2 variants/lineages were analyzed. Leucine is the most abundant amino acid residue in the SARS-CoV-2 with an average AAC of 9.658% while tryptophan had the least abundance of 1.11%. The AAC and ranking of lysine and glycine varied in the proteome. For some variants, glycine had higher frequency and AAC than lysine and vice versa in other variants. Tryptophan was also observed to be the most intolerant to mutation in the various proteomes for the variants used. A correlogram revealed a very strong correlation of 0.999992 between B.1.525 (Eta) and B.1.526 (Iota) variants. Furthermore, isoleucine and threonine were observed to have a very strong negative correlation of -0.912, while cysteine and isoleucine had a very strong positive correlation of 0.835 at p < 0.001. Shapiro-Wilk normality test revealed that AAC values for all the amino acid residues except methionine showed no evidence of non-normality at p < 0.05. Thus, AACs of SARS-CoV-2 variants can be predicted using probability and z-scores. AACs may be beneficial in classifying viral strains, predicting viral disease types, members of protein families, protein interactions and for diagnostic purposes. They may also be used as a feature along with other crucial factors in machine-learning based algorithms to predict viral mutations. These mutation-predicting algorithms may help in developing effective therapeutics and vaccines for SARS-CoV-2.

Identification of ADAR1 p150 and p110 Associated Edit Sites

Adenosine deaminase acting on RNA 1 (ADAR1) catalyzes adenosine-to-inosine editing on double-stranded RNA molecules and is involved in regulating cellular responses to endogenous and exogenous RNA. ADAR1 is the primary A-to-I editor of RNA in humans, and the majority of edit sites are found in a class of short interspersed nuclear elements called Alu elements, many of which are located in introns and 3' untranslated regions. Two ADAR1 protein isoforms, p110 (110 kDa) and p150 (150 kDa), are known to be coupled in expression, and decoupling the expression of these isoforms has revealed that the p150 isoform edits a broader range of targets compared to p110. Numerous methods for identification of ADAR1-associated edits have been developed, and we present here a specific method for identification of edit sites associated with individual ADAR1 isoforms.

Monoclonal antibody therapies in the management of SARS-CoV-2 infection

INTRODUCTION

Neutralizing antibodies (NAbs) that target key domains of the spike protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) may have therapeutic value because of their specificity. Depending on the targeted epitope, single agents may be effective, but combined treatment involving multiple NAbs may be necessary to prevent the emergence of resistant variants.

AREAS COVERED

This article highlights the accelerated regulatory processes established to facilitate the review and approval of potential therapies. An overview of treatment approaches for SARS-CoV-2 infection, with detailed examination of the preclinical and clinical evidence supporting the use of NAbs, is provided. Finally, insights are offered into the potential benefits and challenges associated with the use of these agents.

EXPERT OPINION

NAbs offer an effective, evidence-based therapeutic intervention during the early stages of SARS-CoV-2 infection when viral replication is the primary factor driving disease progression. As the pandemic progresses, appropriate use of NAbs will be important to minimize the risk of escape variants. Ultimately, the availability of effective treatments for COVID-19 will allow the establishment of treatment algorithms for minimizing the substantial rates of hospitalization, morbidity (including long COVID) and mortality currently associated with the disease.

Host-specific asymmetric accumulation of mutation types reveals that the origin of SARS-CoV-2 is consistent with a natural process

The capacity of RNA viruses to adapt to new hosts and rapidly escape the host immune system is largely attributable to de novo genetic diversity that emerges through mutations in RNA. Although the molecular spectrum of de novo mutations-the relative rates at which various base substitutions occur-are widely recognized as informative toward understanding the evolution of a viral genome, little attention has been paid to the possibility of using molecular spectra to infer the host origins of a virus. Here, we characterize the molecular spectrum of de novo mutations for SARS-CoV-2 from transcriptomic data obtained from virus-infected cell lines, enabled by the use of sporadic junctions formed during discontinuous transcription as molecular barcodes. We find that de novo mutations are generated in a replication-independent manner, typically on the genomic strand, and highly dependent on mutagenic mechanisms specific to the host cellular environment. De novo mutations will then strongly influence the types of base substitutions accumulated during SARS-CoV-2 evolution, in an asymmetric manner favoring specific mutation types. Consequently, similarities between the mutation spectra of SARS-CoV-2 and the bat coronavirus RaTG13, which have accumulated since their divergence strongly suggest that SARS-CoV-2 evolved in a host cellular environment highly similar to that of bats before its zoonotic transfer into humans. Collectively, our findings provide data-driven support for the natural origin of SARS-CoV-2.

A rapidly reversible mutation generates subclonal genetic diversity and unstable drug resistance

Most genetic changes have negligible reversion rates. As most mutations that confer resistance to an adverse condition (e.g., drug treatment) also confer a growth defect in its absence, it is challenging for cells to genetically adapt to transient environmental changes. Here, we identify a set of rapidly reversible drug-resistance mutations in Schizosaccharomyces pombe that are caused by microhomology-mediated tandem duplication (MTD) and reversion back to the wild-type sequence. Using 10,000× coverage whole-genome sequencing, we identify nearly 6,000 subclonal MTDs in a single clonal population and determine, using machine learning, how MTD frequency is encoded in the genome. We find that sequences with the highest-predicted MTD rates tend to generate insertions that maintain the correct reading frame, suggesting that MTD formation has shaped the evolution of coding sequences. Our study reveals a common mechanism of reversible genetic variation that is beneficial for adaptation to environmental fluctuations and facilitates evolutionary divergence.

A genome-scale yeast library with inducible expression of individual genes

The ability to switch a gene from off to on and monitor dynamic changes provides a powerful approach for probing gene function and elucidating causal regulatory relationships. Here, we developed and characterized YETI (Yeast Estradiol strains with Titratable Induction), a collection in which > 5,600 yeast genes are engineered for transcriptional inducibility with single-gene precision at their native loci and without plasmids. Each strain contains SGA screening markers and a unique barcode, enabling high-throughput genetics. We characterized YETI using growth phenotyping and BAR-seq screens, and we used a YETI allele to identify the regulon of Rof1, showing that it acts to repress transcription. We observed that strains with inducible essential genes that have low native expression can often grow without inducer. Analysis of data from eukaryotic and prokaryotic systems shows that native expression is a variable that can bias promoter-perturbing screens, including CRISPRi. We engineered a second expression system, Z3 EB42, that gives lower expression than Z3 EV, a feature enabling conditional activation and repression of lowly expressed essential genes that grow without inducer in the YETI library.

EmPC-seq: Accurate RNA-sequencing and Bioinformatics Platform to Map RNA Polymerases and Remove Background Error

Transcription errors can substantially affect metabolic processes in organisms by altering the epigenome and causing misincorporations in mRNA, which is translated into aberrant mutant proteins. Moreover, within eukaryotic genomes there are specific Transcription Error-Enriched genomic Loci (TEELs) which are transcribed by RNA polymerases with significantly higher error rates and hypothesized to have implications in cancer, aging, and diseases such as Down syndrome and Alzheimer's. Therefore, research into transcription errors is of growing importance within the field of genetics. Nevertheless, methodological barriers limit the progress in accurately identifying transcription errors. Pro-Seq and NET-Seq can purify nascent RNA and map RNA polymerases along the genome but cannot be used to identify transcriptional mutations. Here we present background Error Model-coupled Precision nuclear run-on Circular-sequencing (EmPC-seq), a method combining a nuclear run-on assay and circular sequencing with a background error model to precisely detect nascent transcription errors and effectively discern TEELs within the genome.

The evolution and clinical impact of hepatitis B virus genome diversity

The global burden of hepatitis B virus (HBV) is enormous, with 257 million persons chronically infected, resulting in more than 880,000 deaths per year worldwide. HBV exists as nine different genotypes, which differ in disease progression, natural history and response to therapy. HBV is an ancient virus, with the latest reports greatly expanding the host range of the Hepadnaviridae (to include fish and reptiles) and casting new light on the origins and evolution of this viral family. Although there is an effective preventive vaccine, there is no cure for chronic hepatitis B, largely owing to the persistence of a viral minichromosome that is not targeted by current therapies. HBV persistence is also facilitated through aberrant host immune responses, possibly due to the diverse intra-host viral populations that can respond to host-mounted and therapeutic selection pressures. This Review summarizes current knowledge on the influence of HBV diversity on disease progression and treatment response and the potential effect on new HBV therapies in the pipeline. The mechanisms by which HBV diversity can occur both within the individual host and at a population level are also discussed.

SVSR: A Program to Simulate Structural Variations and Generate Sequencing Reads for Multiple Platforms

Structural variation accounts for a major fraction of mutations in the human genome and confers susceptibility to complex diseases. Next generation sequencing along with the rapid development of computational methods provides a cost-effective procedure to detect such variations. Simulation of structural variations and sequencing reads with real characteristics is essential for benchmarking the computational methods. Here, we develop a new program, SVSR, to simulate five types of structural variations (indels, tandem duplication, CNVs, inversions, and translocations) and SNPs for the human genome and to generate sequencing reads with features from popular platforms (Illumina, SOLiD, 454, and Ion Torrent). We adopt a selection model trained from real data to predict copy number states, starting from the first site of a particular genome to the end. Furthermore, we utilize references of microbial genomes to produce insertion fragments and design probabilistic models to imitate inversions and translocations. Moreover, we create platform-specific errors and base quality profiles to generate normal, tumor, or normal-tumor mixture reads. Experimental results show that SVSR could capture more features that are realistic and generate datasets with satisfactory quality scores. SVSR is able to evaluate the performance of structural variation detection methods and guide the development of new computational methods.

The FKH domain in FOXP3 mRNA frequently contains mutations in hepatocellular carcinoma that influence the subcellular localization and functions of FOXP3

The transcription factor forkhead box P3 (FOXP3) is a biomarker for regulatory T cells and can also be expressed in cancer cells, but its function in cancer appears to be divergent. The role of hepatocyte-expressed FOXP3 in hepatocellular carcinoma (HCC) is unknown. Here, we collected tumor samples and clinical information from 115 HCC patients and used five human cancer cell lines. We examined FOXP3 mRNA sequences for mutations, used a luciferase assay to assess promoter activities of FOXP3's target genes, and employed mouse tumor models to confirm in vitro results. We detected mutations in the FKH domain of FOXP3 mRNAs in 33% of the HCC tumor tissues, but in none of the adjacent nontumor tissues. None of the mutations occurred at high frequency, indicating that they occurred randomly. Notably, the mutations were not detected in the corresponding regions of FOXP3 genomic DNA, and many of them resulted in amino acid substitutions in the FKH region, altering FOXP3's subcellular localization. FOXP3 delocalization from the nucleus to the cytoplasm caused loss of transcriptional regulation of its target genes, inactivated its tumor-inhibitory capability, and changed cellular responses to histone deacetylase (HDAC) inhibitors. More complex FKH mutations appeared to be associated with worse prognosis in HCC patients. We conclude that mutations in the FKH domain of FOXP3 mRNA frequently occur in HCC and that these mutations are caused by errors in transcription and are not derived from genomic DNA mutations. Our results suggest that transcriptional mutagenesis of FOXP3 plays a role in HCC.

Transcriptional Fidelity of Mitochondrial RNA Polymerase RpoTm from Arabidopsis thaliana

Fidelity of RNA synthesis is essential for the faithful transfer of information from DNA to RNA. A comprehensive analysis of the nucleotide selectivity by the mitochondrial RNA polymerase (RNAP) RpoTm, from Arabidopsis thaliana, has been carried out. The kinetic parameters for the incorporation of cognate, noncognate, and oxidized bases have been determined. The results establish high fidelity of mitochondrial transcription resembling those of replicative polymerases in the absence of repair. In addition, RpoTm incorporates oxidized nucleotides with similar efficiency compared with mismatches and is capable of extending the RNA beyond the insertion of the oxidized base. Furthermore, lesion bypass study on RpoTm demonstrates that the enzyme bypasses 8-oxo-guanine by insertion of adenine leading to C to A mutations in RNA. Homology modeling of RpoTm elongation complex allows delineation of the residues necessary for stabilizing the incoming NTP substrate and for posing the template nucleotide residue. Substitution of these residues leads to compromise in the activity of the enzyme corroborating their importance in RNA synthesis. Comparison of the data with T7 RNAPs indicates that low efficiency of misincorporation is a universal strategy used by single-subunit RNAPs for maintaining high fidelity in the absence of proofreading and repair activity in mitochondria.

An Omics Approach to Extracellular Vesicles from HIV-1 Infected Cells

Human Immunodeficiency Virus-1 (HIV-1) is the causative agent of Acquired Immunodeficiency Syndrome (AIDS), infecting nearly 37 million people worldwide. Currently, there is no definitive cure, mainly due to HIV-1's ability to enact latency. Our previous work has shown that exosomes, a small extracellular vesicle, from uninfected cells can activate HIV-1 in latent cells, leading to increased mostly short and some long HIV-1 RNA transcripts. This is consistent with the notion that none of the FDA-approved antiretroviral drugs used today in the clinic are transcription inhibitors. Furthermore, these HIV-1 transcripts can be packaged into exosomes and released from the infected cell. Here, we examined the differences in protein and nucleic acid content between exosomes from uninfected and HIV-1-infected cells. We found increased cyclin-dependent kinases, among other kinases, in exosomes from infected T-cells while other kinases were present in exosomes from infected monocytes. Additionally, we found a series of short antisense HIV-1 RNA from the 3' LTR that appears heavily mutated in exosomes from HIV-1-infected cells along with the presence of cellular noncoding RNAs and cellular miRNAs. Both physical and functional validations were performed on some of the key findings. Collectively, our data indicate distinct differences in protein and RNA content between exosomes from uninfected and HIV-1-infected cells, which can lead to different functional outcomes in recipient cells.

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.

Deoxyuridine in DNA has an inhibitory and promutagenic effect on RNA transcription by diverse RNA polymerases

dUTP is a close structural congener of dTTP and can be readily incorporated into DNA opposite to adenine during DNA replication leading to non-mutagenic dU/A base pairs ('uracilation'). We find that dU/A pairs located within DNA transcriptional templates optimized for either T7 RNA polymerase (T7 RNAP) or human RNA polymerase II (pol II) have inhibitory and mutagenic effects on transcription. The data for T7 RNAP establishes that even a single dU/A pair can inhibit promoter binding and transcription initiation up to 30-fold, and that inhibitory effects on transcription elongation are also possible. Sequencing of the mRNA transcribed from uniformly uracilated DNA templates by T7 RNAP indicated an increased frequency of transversion and insertion mutations compared to all T/A templates. Strong effects of dU/A pairs on cellular transcription activity and fidelity were also observed with RNA pol II using uracil base excision repair (UBER)-deficient human cells. At the highest levels of template uracilation, transcription by RNA pol II was completely blocked. We propose that these effects arise from the decreased thermodynamic stability and increased dynamics of dU/A pairs in DNA. The potential implications of these findings on gene regulation and disease are discussed.

RNA-Mediated Disease Mechanisms in Neurodegenerative Disorders

RNA is accurately entangled in virtually all pathways that maintain cellular homeostasis. To name but a few, RNA is the "messenger" between DNA encoded information and the resulting proteins. Furthermore, RNAs regulate diverse processes by forming DNA::RNA or RNA::RNA interactions. Finally, RNA itself can be the scaffold for ribonucleoprotein complexes, for example, ribosomes or cellular bodies. Consequently, disruption of any of these processes can lead to disease. This review describes known and emerging RNA-based disease mechanisms like interference with regular splicing, the anomalous appearance of RNA-protein complexes and uncommon RNA species, as well as non-canonical translation. Due to the complexity and entanglement of the above-mentioned pathways, only few drugs are available that target RNA-based disease mechanisms. However, advances in our understanding how RNA is involved in and modulates cellular homeostasis might pave the way to novel treatments.

Rationally designed perturbation factor drives evolution in Saccharomyces cerevisiae for industrial application

Saccharomyces cerevisiae strains with favorable characteristics are preferred for application in industries. However, the current ability to reprogram a yeast cell on the genome scale is limited due to the complexity of yeast ploids. In this study, a method named genome replication engineering-assisted continuous evolution (GREACE) was proved efficient in engineering S. cerevisiae with different ploids. Through iterative cycles of culture coupled with selection, GREACE could continuously improve the target traits of yeast by accumulating beneficial genetic modification in genome. The application of GREACE greatly improved the tolerance of yeast against acetic acid compared with their parent strain. This method could also be employed to improve yeast aroma profile and the phenotype could be stably inherited to the offspring. Therefore, GREACE method was efficient in S. cerevisiae engineering and it could be further used to evolve yeast with other specific characteristics.

Poor codon optimality as a signal to degrade transcripts with frameshifts

Frameshifting errors are common and mRNA quality control pathways, such as nonsense-mediated decay (NMD), exist to degrade these aberrant transcripts. Recent work has shown the existence of a genetic link between NMD and codon-usage mediated mRNA decay. Here we present computational evidence that these pathways are synergic for removing frameshifts.

O6-methylguanine-induced transcriptional mutagenesis reduces p53 tumor-suppressor function

The impact of DNA lesions on replication and mutagenesis is of high relevance for human health; however, the role of lesion-induced transcriptional mutagenesis (TM) in disease development is unknown. Here, the impact of O⁶-methylguanine–induced TM on p53 function as a tumor suppressor was investigated in human cells. Results showed that TM in 15% of the transcripts resulted in a reduced ability of p53 protein to transactivate genes that regulate cell-cycle arrest and induction of apoptosis. This resulted in the loss of functional cell-cycle checkpoints and in impaired activation of apoptosis, both canonical p53 tumor-suppressor functions. This work provides evidence that TM can induce phenotypic changes in mammalian cells that have important implications for its role in tumorigenesis.

Altered protein function due to mutagenesis plays an important role in disease development. This is perhaps most evident in tumorigenesis and the associated loss or gain of function of tumor-suppressor genes and oncogenes. The extent to which lesion-induced transcriptional mutagenesis (TM) influences protein function and its contribution to the development of disease is not well understood. In this study, the impact of O⁶-methylguanine on the transcription fidelity of p53 and the subsequent effects on the protein’s function as a regulator of cell death and cell-cycle arrest were examined in human cells. Levels of TM were determined by RNA-sequencing. In cells with active DNA repair, misincorporation of uridine opposite the lesion occurred in 0.14% of the transcripts and increased to 14.7% when repair by alkylguanine–DNA alkyltransferase was compromised. Expression of the dominant-negative p53 R248W mutant due to TM significantly reduced the transactivation of several established p53 target genes that mediate the tumor-suppressor function, including CDKN1A (p21) and BBC3 (PUMA). This resulted in deregulated signaling through the retinoblastoma protein and loss of G1/S cell-cycle checkpoint function. In addition, we observed impaired activation of apoptosis coupled to the reduction of the tumor-suppressor functions of p53. Taking these findings together, this work provides evidence that TM can induce phenotypic changes in mammalian cells that have important implications for the role of TM in tumorigenesis.

Transcriptional mutagenesis reduces splicing fidelity in mammalian cells

Splicing fidelity is essential to the maintenance of cellular functions and viability, and mutations or natural variations in pre-mRNA sequences and consequent alteration of splicing have been implicated in the etiology and progression of numerous diseases. The extent to which transcriptional errors or lesion-induced transcriptional mutagenesis (TM) influences splicing fidelity is not currently known. To investigate this, we employed site-specific DNA lesions on the transcribed strand of a minigene splicing reporter in normal mammalian cells. These were the common mutagenic lesions O⁶-methylguanine (O⁶-meG) and 8-oxoguanine (8-oxoG). The minigene splicing reporters were derived from lamin A (LMNA) and proteolipid protein 1 (PLP1), both with known links to human diseases that result from deregulated splicing. In cells with active DNA repair, 1–4% misincorporation occurred opposite the lesions, which increased to 20–40% when repair was compromised. Furthermore, our results reveal that TM at a splice site significantly reduces in vivo splicing fidelity, thereby changing the relative expression of alternative splicing forms in mammalian cells. These findings suggest that splicing defects caused by transcriptional errors can potentially lead to phenotypic cellular changes and increased susceptibility to the development of disease.

Misincorporation by RNA polymerase is a major source of transcription pausing in vivo

The transcription error rate estimated from mistakes in end product RNAs is 10−3–10−5. We analyzed the fidelity of nascent RNAs from all actively transcribing elongation complexes (ECs) in Escherichia coli and Saccharomyces cerevisiae and found that 1–3% of all ECs in wild-type cells, and 5–7% of all ECs in cells lacking proofreading factors are, in fact, misincorporated complexes. With the exception of a number of sequence-dependent hotspots, most misincorporations are distributed relatively randomly. Misincorporation at hotspots does not appear to be stimulated by pausing. Since misincorporation leads to a strong pause of transcription due to backtracking, our findings indicate that misincorporation could be a major source of transcriptional pausing and lead to conflicts with other RNA polymerases and replication in bacteria and eukaryotes. This observation implies that physical resolution of misincorporated complexes may be the main function of the proofreading factors Gre and TFIIS. Although misincorporation mechanisms between bacteria and eukaryotes appear to be conserved, the results suggest the existence of a bacteria-specific mechanism(s) for reducing misincorporation in protein-coding regions. The links between transcription fidelity, human disease, and phenotypic variability in genetically-identical cells can be explained by the accumulation of misincorporated complexes, rather than mistakes in mature RNA.

The landscape of transcription errors in eukaryotic cells

Accurate transcription is required for the faithful expression of genetic information. To understand the molecular mechanisms that control the fidelity of transcription, we used novel sequencing technology to provide the first comprehensive analysis of the fidelity of transcription in eukaryotic cells. Our results demonstrate that transcription errors can occur in any gene, at any location, and affect every aspect of protein structure and function. In addition, we show that multiple proteins safeguard the fidelity of transcription and provide evidence suggesting that errors that evade these layers of RNA quality control profoundly affect the physiology of living cells. Together, these observations demonstrate that there is an inherent limit to the faithful expression of the genome and suggest that the impact of mutagenesis on cellular health and fitness is substantially greater than currently appreciated.

Transcription fidelity and its roles in the cell

•

The Trigger Loop is one of the major determinants of transcription fidelity.

•

Intrinsic proofreading occurs via transcript-assisted cleavage.

•

Factor-assisted proofreading takes place via exchange of RNAP active centres.

•

Misincorporation is a major source of transcription pausing.

•

Another role of fidelity is the prevention of conflicts with other cellular processes.

Accuracy of transcription is essential for productive gene expression, and the past decade has brought new understanding of the mechanisms ensuring transcription fidelity. The discovery of a new catalytic domain, the Trigger Loop, revealed that RNA polymerase can actively choose the correct substrates. Also, the intrinsic proofreading activity was found to proceed via a ribozyme-like mechanism, whereby the erroneous nucleoside triphosphate (NTP) helps its own excision. Factor-assisted proofreading was shown to proceed through an exchange of active centres, a unique phenomenon among proteinaceous enzymes. Furthermore, most recent in vivo studies have revised the roles of transcription accuracy and proofreading factors, as not only required for production of errorless RNAs, but also for prevention of frequent misincorporation-induced pausing that may cause conflicts with fellow RNA polymerases and the replication machinery.

Transcriptional fidelities of human mitochondrial POLRMT, yeast mitochondrial Rpo41, and phage T7 single-subunit RNA polymerases

Single-subunit RNA polymerases (RNAPs) are present in phage T7 and in mitochondria of all eukaryotes. This RNAP class plays important roles in biotechnology and cellular energy production, but we know little about its fidelity and error rates. Herein, we report the error rates of three single-subunit RNAPs measured from the catalytic efficiencies of correct and all possible incorrect nucleotides. The average error rates of T7 RNAP (2 × 10^-6), yeast mitochondrial Rpo41 (6 × 10^-6), and human mitochondrial POLRMT (RNA polymerase mitochondrial) (2 × 10^-5) indicate high accuracy/fidelity of RNA synthesis resembling those of replicative DNA polymerases. All three RNAPs exhibit a distinctly high propensity for GTP misincorporation opposite dT, predicting frequent A→G errors in RNA with rates of ∼10^-4 The A→C, G→A, A→U, C→U, G→U, U→C, and U→G errors mostly due to pyrimidine-purine mismatches were relatively frequent (10^-5-10^-6), whereas C→G, U→A, G→C, and C→A errors from purine-purine and pyrimidine-pyrimidine mismatches were rare (10^-7-10^-10). POLRMT also shows a high C→A error rate on 8-oxo-dG templates (∼10^-4). Strikingly, POLRMT shows a high mutagenic bypass rate, which is exacerbated by TEFM (transcription elongation factor mitochondrial). The lifetime of POLRMT on terminally mismatched elongation substrate is increased in the presence of TEFM, which allows POLRMT to efficiently bypass the error and continue with transcription. This investigation of nucleotide selectivity on normal and oxidatively damaged DNA by three single-subunit RNAPs provides the basic information to understand the error rates in mitochondria and, in the case of T7 RNAP, to assess the quality of in vitro transcribed RNAs.

The histone variant H2A.Z promotes efficient cotranscriptional splicing in S. cerevisiae

In eukaryotes, a dynamic ribonucleic protein machine known as the spliceosome catalyzes the removal of introns from premessenger RNA (pre-mRNA). Recent studies show the processes of RNA synthesis and RNA processing to be spatio-temporally coordinated, indicating that RNA splicing takes place in the context of chromatin. H2A.Z is a highly conserved histone variant of the canonical histone H2A. In Saccharomyces cerevisiae, H2A.Z is deposited into chromatin by the SWR-C complex, is found near the 5' ends of protein-coding genes, and has been implicated in transcription regulation. Here we show that splicing of intron-containing genes in cells lacking H2A.Z is impaired, particularly under suboptimal splicing conditions. Cells lacking H2A.Z are especially dependent on a functional U2 snRNP (small nuclear RNA [snRNA] plus associated proteins), as H2A.Z shows extensive genetic interactions with U2 snRNP-associated proteins, and RNA sequencing (RNA-seq) reveals that introns with nonconsensus branch points are particularly sensitive to H2A.Z loss. Consistently, H2A.Z promotes efficient spliceosomal rearrangements involving the U2 snRNP, as H2A.Z loss results in persistent U2 snRNP association and decreased recruitment of downstream snRNPs to nascent RNA. H2A.Z impairs transcription elongation, suggesting that spliceosome rearrangements are tied to H2A.Z's role in elongation. Depletion of disassembly factor Prp43 suppresses H2A.Z-mediated splice defects, indicating that, in the absence of H2A.Z, stalled spliceosomes are disassembled, and unspliced RNAs are released. Together, these data demonstrate that H2A.Z is required for efficient pre-mRNA splicing and indicate a role for H2A.Z in coordinating the kinetics of transcription elongation and splicing.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

Viroid quasispecies revealed by deep sequencing

Viroids are non-coding single-stranded circular RNA molecules that replicate autonomously in infected host plants causing mild to lethal symptoms. Their genomes contain about 250-400 nucleotides, depending on viroid species. Members of the family Pospiviroidae, like the Potato spindle tuber viroid (PSTVd), replicate via an asymmetric rolling-circle mechanism using the host DNA-dependent RNA-Polymerase II in the nucleus, while members of Avsunviroidae are replicated in a symmetric rolling-circle mechanism probably by the nuclear-encoded polymerase in chloroplasts. Viroids induce the production of viroid-specific small RNAs (vsRNA) that can direct (post-)transcriptional gene silencing against host transcripts or genomic sequences. Here, we used deep-sequencing to analyze vsRNAs from plants infected with different PSTVd variants to elucidate the PSTVd quasipecies evolved during infection. We recovered several novel as well as previously known PSTVd variants that were obviously competent in replication and identified common strand-specific mutations. The calculated mean error rate per nucleotide position was less than [Formula: see text], quite comparable to the value of [Formula: see text] reported for a member of Avsunviroidae. The resulting error threshold allows the synthesis of longer-than-unit-length replication intermediates as required by the asymmetric rolling-circle mechanism of members of Pospiviroidae.

A link between transcription fidelity and pausing in vivo

Pausing by RNA polymerase is a major mechanism that regulates transcription elongation but can cause conflicts with fellow RNA polymerases and other cellular machineries. Here, we summarize our recent finding that misincorporation could be a major source of transcription pausing in vivo, and discuss the role of misincorporation-induced pausing.

My Cousin, My Enemy: quasispecies suppression of drug resistance

If a freshly minted genome contains a mutation that confers drug resistance, will it be selected in the presence of the drug? Not necessarily. During viral infections, newly synthesized viral genomes occupy the same cells as parent and other progeny genomes. If the antiviral target is chosen so that the drug-resistant progeny's growth is dominantly inhibited by the drug-susceptible members of its intracellular family, its outgrowth can be suppressed. Precedent for 'dominant drug targeting' as a deliberate approach to suppress the outgrowth of inhibitor-resistant viruses has been established for envelope variants of vesicular stomatitis virus and for capsid variants of poliovirus and dengue virus. Small molecules that stabilize oligomeric assemblages are a promising means to an unfit family to destroy the effectiveness of a newborn drug-resistant relative due to the co-assembly of drug-susceptible and drug-resistant monomers.

Genome stability versus transcript diversity

Our genome is protected from the introduction of mutations by high fidelity replication and an extensive network of DNA damage response and repair mechanisms. However, the expression of our genome, via RNA and protein synthesis, allows for more diversity in translating genetic information. In addition, the splicing process has become less stringent over evolutionary time allowing for a substantial increase in the diversity of transcripts generated. The result is a diverse transcriptome and proteome that harbor selective advantages over a more tightly regulated system. Here, we describe mechanisms in place that both safeguard the genome and promote translational diversity, with emphasis on post-transcriptional RNA processing.