ABSTRACT
End-to-end RNA sequencing methods that capture 5'-sequence content without cumbersome library manipulations are of great interest, particularly for analysis of long RNAs. While template-switching methods have been developed for RNA sequencing by distributive short-read RTs, such as the MMLV RT enzymes used in SMART-Seq methods, they have not been adapted to leverage the power of ultraprocessive RTs, such as those that derive from group II self-splicing introns. To facilitate this transition, we dissected the individual processes that guide the enzymatic specificity and efficiency of the multi-step template switching reaction carried out by RT enzymes, in this case, by a well-characterized enzyme known as MarathonRT. Remarkably, this is the first study of its kind, for any RT. First, we characterized and optimized the enzymatic nontemplated addition (NTA) reaction that occurs when the RT enzyme extends past the RNA 5'-terminus, and we determined the nucleotide specificity of the NTA reaction. We then evaluated the binding specificity of specialized template-switching oligonucleotides, optimizing their sequences and chemical properties to guide efficient template switching reaction. Having dissected and optimized these individual steps, we then unified them into a procedure for performing RNA sequencing with MarathonRT enzymes, using a well-characterized RNA reference set. The resulting reads span a six-log range in transcript concentration and accurately represent the input RNA identities in both length and composition. We also performed RNA-seq starting from total human RNA and poly(A)-enriched RNA, with short and long-read sequencing demonstrating that MarathonRT enhances the discovery of unseen RNA molecules by conventional RT. Altogether, by employing mechanistic enzymology on RT enzymes and using them to modify RNA-seq technologies, we have generated a new pipeline for rapid, accurate sequencing of complex RNA libraries containing mixtures of long RNA transcripts.
ABSTRACT
Circulating bat coronaviruses represent a pandemic threat. However, our understanding of bat coronavirus pathogenesis and transmission potential is limited by the lack of phenotypically characterized strains. We created molecular clones for the two closest known relatives of SARS-CoV-2, BANAL-52 and BANAL-236. We demonstrated that BANAL-CoVs and SARS-CoV-2 have similar replication kinetics in human bronchial epithelial cells. However, BANAL-CoVs have impaired replication in human nasal epithelial cells and in the upper airway of mice. We also observed reduced pathogenesis in mice and diminished transmission in hamsters. Further, we observed that diverse bat coronaviruses evade interferon and downregulate major histocompatibility complex class I. Collectively, our study demonstrates that despite high genetic similarity across bat coronaviruses, prediction of pandemic potential of a virus necessitates functional characterization. Finally, the restriction of bat coronavirus replication in the upper airway highlights that transmission potential and innate immune restriction can be uncoupled in this high-risk family of emerging viruses.
Subject(s)
COVID-19 , Chiroptera , Immunity, Innate , SARS-CoV-2 , Virus Replication , Animals , Humans , SARS-CoV-2/immunology , SARS-CoV-2/genetics , SARS-CoV-2/physiology , Chiroptera/virology , Chiroptera/immunology , COVID-19/transmission , COVID-19/virology , COVID-19/immunology , Mice , Cricetinae , Immune Evasion , Epithelial Cells/virology , Epithelial Cells/immunology , Coronavirus Infections/transmission , Coronavirus Infections/immunology , Coronavirus Infections/virology , Coronavirus/immunology , Coronavirus/genetics , Coronavirus/classification , Coronavirus/physiology , Coronavirus/pathogenicity , Cell Line , FemaleABSTRACT
West Nile virus (WNV) is an arthropod-borne, positive-sense RNA virus that poses an increasing global threat due to warming climates and lack of effective therapeutics. Like other enzootic viruses, little is known about how host context affects the structure of the full-length RNA genome. Here, we report a complete secondary structure of the entire WNV genome within infected mammalian and arthropod cell lines. Our analysis affords structural insights into multiple, conserved aspects of flaviviral biology. We show that the WNV genome folds with minimal host dependence, and we prioritize well-folded regions for functional validation using structural homology between hosts as a guide. Using structure-disrupting, antisense locked nucleic acids, we then demonstrate that the WNV genome contains riboregulatory structures with conserved and host-specific functional roles. These results reveal promising RNA drug targets within flaviviral genomes, and they highlight the therapeutic potential of ASO-LNAs as both WNV-specific and pan-flaviviral therapeutic agents.
Subject(s)
Genome, Viral , RNA, Viral , West Nile virus , West Nile virus/genetics , Animals , RNA, Viral/genetics , RNA, Viral/metabolism , Humans , Cell Line , Nucleic Acid Conformation , West Nile Fever/virology , Host Specificity/genetics , Host-Pathogen Interactions/geneticsABSTRACT
Pharmacological activation of the retinoic acid-inducible gene I (RIG-I) pathway holds promise for increasing tumor immunogenicity and improving the response to immune checkpoint inhibitors (ICIs). However, the potency and clinical efficacy of 5'-triphosphate RNA (3pRNA) agonists of RIG-I are hindered by multiple pharmacological barriers, including poor pharmacokinetics, nuclease degradation, and inefficient delivery to the cytosol where RIG-I is localized. Here, we address these challenges through the design and evaluation of ionizable lipid nanoparticles (LNPs) for the delivery of 3p-modified stem-loop RNAs (SLRs). Packaging of SLRs into LNPs (SLR-LNPs) yielded surface charge-neutral nanoparticles with a size of â¼100 nm that activated RIG-I signaling in vitro and in vivo. SLR-LNPs were safely administered to mice via both intratumoral and intravenous routes, resulting in RIG-I activation in the tumor microenvironment (TME) and the inhibition of tumor growth in mouse models of poorly immunogenic melanoma and breast cancer. Significantly, we found that systemic administration of SLR-LNPs reprogrammed the breast TME to enhance the infiltration of CD8+ and CD4+ T cells with antitumor function, resulting in enhanced response to αPD-1 ICI in an orthotopic EO771 model of triple-negative breast cancer. Therapeutic efficacy was further demonstrated in a metastatic B16.F10 melanoma model, with systemically administered SLR-LNPs significantly reducing lung metastatic burden compared to combined αPD-1 + αCTLA-4 ICI. Collectively, these studies have established SLR-LNPs as a translationally promising immunotherapeutic nanomedicine for potent and selective activation of RIG-I with the potential to enhance response to ICIs and other immunotherapeutic modalities.
Subject(s)
Immunotherapy , Nanoparticles , Animals , Female , Humans , Mice , Cell Line, Tumor , Lipids/chemistry , Mice, Inbred C57BL , Nanoparticles/chemistry , Tumor Microenvironment/drug effectsABSTRACT
Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.
Subject(s)
Antifungal Agents , Candida albicans , Introns , Humans , Antifungal Agents/pharmacology , Candida albicans/drug effects , Candida albicans/genetics , Introns/genetics , RNA Splicing/genetics , RNA, Fungal/metabolismABSTRACT
Hsp70 (70 kDa heat shock protein) performs molecular chaperone functions by assisting the folding of newly synthesized and misfolded proteins, thereby counteracting various cell stresses and preventing multiple diseases, including neurodegenerative disorders and cancers. It is well established that, immediately after heat shock, Hsp70 gene expression is mediated by a canonical mechanism of cap-dependent translation. However, the molecular mechanism of Hsp70 expression during heat shock remains elusive. Intriguingly, the 5' end of Hsp70 messenger RNA (mRNA) appears to form a compact structure with the potential to regulate protein expression in a cap-independent manner. Here, we determined the minimal length of the mHsp70 5'-terminal mRNA sequence that is required for RNA folding into a highly compact structure. This span of this RNA element was mapped and the secondary structure characterized by chemical probing, resulting in a secondary structural model that includes multiple stable stems, including one containing the canonical start codon. All of these components, including a short stretch of the 5' open reading frame (ORF), were shown to be vital for RNA folding. This work provides a structural basis for future investigations on the role of translational regulatory structures in the 5' untranslated region and ORF sequences of Hsp70 during heat shock.
ABSTRACT
Although next-generation sequencing (NGS) technologies have revolutionized our ability to sequence DNA with high-throughput, the chain termination-based Sanger sequencing method remains a widely used approach for DNA sequence analysis due to its simplicity, low cost and high accuracy. In particular, high accuracy makes Sanger sequencing the "gold standard" for sequence validation in basic research and clinical applications. During the early days of Sanger sequencing development, reverse transcriptase (RT)-based RNA sequencing was also explored and showed great promise, but the approach did not acquire popularity over time due to the limited processivity and low template unwinding capability of Avian Myeloblastosis Virus (AMV) RT, and other RT enzymes available at the time. RNA molecules have complex features, often containing repetitive sequences and stable secondary or tertiary structures. While these features are required for RNA biological function, they represent strong obstacles for retroviral RTs. Repetitive sequences and stable structures cause reverse transcription errors and premature primer extension stops, making chain termination-based methods unfeasible. MarathonRT is an ultra-processive RT encoded group II intron that can copy RNA molecules of any sequence and structure in a single cycle, making it an ideal RT enzyme for Sanger RNA sequencing. In this chapter, we upgrade the Sanger RNA sequencing method by replacing AMV RT with MarathonRT, providing a simple, robust method for direct RNA sequence analysis. The guidance for troubleshooting and further optimization are also provided.
Subject(s)
RNA-Directed DNA Polymerase , RNA , RNA/genetics , RNA-Directed DNA Polymerase/genetics , RNA-Directed DNA Polymerase/metabolism , Repetitive Sequences, Nucleic Acid , High-Throughput Nucleotide Sequencing , Sequence Analysis, RNAABSTRACT
RNA is playing an ever-growing role in molecular biology and biomedicine due to the many ways it influences gene expression and its increasing use in modern therapeutics. Hence, production of RNA molecules in large quantity and high purity has become essential for advancing basic scientific research and for developing next-generation therapeutics. T7 RNA polymerase (RNAP) is a DNA-dependent RNA polymerase of bacteriophage origin and it is the most widely-utilized tool enzyme for producing RNA. Here we describe a set of robust methods for in vitro transcribing RNA molecules from DNA templates using T7 RNAP, along with a set of subsequent RNA purification schemes. In the first part of this chapter, we provide the general method for T7 RNAP-based in vitro transcription and technical notes for troubleshooting failed or inefficient transcription. We also provide modified protocols for preparing specialized RNA transcripts. In the second part, we provide two purification methods using either gel-based denaturing purification or size exclusion column-based non-denaturing purification for isolating high-purity RNA products from transcription reaction mixtures and preparing them for downstream applications. This chapter is designed to provide researchers with versatile ways to efficiently generate RNA molecules of interest and a troubleshooting guide should they encounter problems while working with in vitro transcription using T7 RNAP.
Subject(s)
RNA , Transcription, Genetic , DNA-Directed RNA Polymerases/genetics , DNA-Directed RNA Polymerases/metabolism , DNA , Bacteriophage T7/genetics , Bacteriophage T7/metabolismABSTRACT
RNA molecules play important roles in numerous normal cellular processes and disease states, from protein coding to gene regulation. RT-PCR, applying the power of polymerase chain reaction (PCR) to RNA by coupling reverse transcription with PCR, is one of the most important techniques to characterize RNA transcripts and monitor gene expression. The ability to analyze full-length RNA transcripts and detect their expression is critical to decipher their biological functions. However, due to the low processivity of retroviral reverse transcriptases (RTs), we can only monitor a small fraction of long RNA transcripts, especially those containing stable secondary and tertiary structures. The full-length sequences can only be deduced by computational analysis, which is often misleading. Group II intron-encoded RTs are a new type of RT enzymes. They have evolved specialized structural elements that unwind template structures and maintain close contact with the RNA template. Therefore, group II intron-encoded RTs are more processive than the retroviral RTs. The discovery, optimization and deployment of processive group II intron RTs provide us the opportunity to analyze RNA transcripts with single molecule resolution. MarathonRT, the most processive group II intron RT, has been extensively optimized for processive reverse transcription. In this chapter, we use MarathonRT to deliver a general protocol for long amplicon generation by RT-PCR, and also provide guidance for troubleshooting and further optimization.
Subject(s)
RNA-Directed DNA Polymerase , RNA , RNA/chemistry , Reverse Transcriptase Polymerase Chain Reaction , RNA-Directed DNA Polymerase/genetics , RNA-Directed DNA Polymerase/chemistry , RNA-Directed DNA Polymerase/metabolism , Polymerase Chain Reaction , IntronsABSTRACT
The group II intron ribonucleoprotein is an archetypal splicing system with numerous mechanistic parallels to the spliceosome, including excision of lariat introns1,2. Despite the importance of branching in RNA metabolism, structural understanding of this process has remained elusive. Here we present a comprehensive analysis of three single-particle cryogenic electron microscopy structures captured along the splicing pathway. They reveal the network of molecular interactions that specifies the branchpoint adenosine and positions key functional groups to catalyse lariat formation and coordinate exon ligation. The structures also reveal conformational rearrangements of the branch helix and the mechanism of splice site exchange that facilitate the transition from branching to ligation. These findings shed light on the evolution of splicing and highlight the conservation of structural components, catalytic mechanism and dynamical strategies retained through time in premessenger RNA splicing machines.
Subject(s)
Biocatalysis , Introns , Nucleic Acid Conformation , RNA Splicing , Adenosine/metabolism , Cryoelectron Microscopy , Exons , RNA Precursors/chemistry , RNA Precursors/metabolism , RNA Precursors/ultrastructure , RNA Splice SitesABSTRACT
RIG-I is an essential innate immune receptor that responds to infection by RNA viruses. The RIG-I signaling cascade is mediated by a series of post-translational modifications, the most important of which is ubiquitination of the RIG-I Caspase Recruitment Domains (CARDs) by E3 ligase Riplet. This is required for interaction between RIG-I and its downstream adapter protein MAVS, but the mechanism of action remains unclear. Here we show that Riplet is required for RIG-I signaling in the presence of both short and long dsRNAs, establishing that Riplet activation does not depend upon RIG-I filament formation on long dsRNAs. Likewise, quantitative Riplet-RIG-I affinity measurements establish that Riplet interacts with RIG-I regardless of whether the receptor is bound to RNA. To understand this, we solved high-resolution cryo-EM structures of RIG-I/RNA/Riplet complexes, revealing molecular interfaces that control Riplet-mediated activation and enabling the formulation of a unified model for the role of Riplet in signaling.
Subject(s)
Immunity, Innate , Ubiquitin-Protein Ligases , Ubiquitin-Protein Ligases/metabolism , DEAD Box Protein 58/metabolism , Signal Transduction , Ubiquitination , RNA, Double-StrandedABSTRACT
Spliceosomal snRNPs are multicomponent particles that undergo a complex maturation pathway. Human Sm-class snRNAs are generated as 3'-end extended precursors, which are exported to the cytoplasm and assembled together with Sm proteins into core RNPs by the SMN complex. Here, we provide evidence that these pre-snRNA substrates contain compact, evolutionarily conserved secondary structures that overlap with the Sm binding site. These structural motifs in pre-snRNAs are predicted to interfere with Sm core assembly. We model structural rearrangements that lead to an open pre-snRNA conformation compatible with Sm protein interaction. The predicted rearrangement pathway is conserved in Metazoa and requires an external factor that initiates snRNA remodeling. We show that the essential helicase Gemin3, which is a component of the SMN complex, is crucial for snRNA structural rearrangements during snRNP maturation. The SMN complex thus facilitates ATP-driven structural changes in snRNAs that expose the Sm site and enable Sm protein binding.
Subject(s)
RNA Precursors , RNA, Small Nuclear , Humans , RNA, Small Nuclear/metabolism , SMN Complex Proteins/metabolism , RNA Precursors/metabolism , HeLa Cells , Ribonucleoproteins, Small Nuclear/metabolism , snRNP Core Proteins/geneticsABSTRACT
Chemical modifications are essential regulatory elements that modulate the behavior and function of cellular RNAs. Despite recent advances in sequencing-based RNA modification mapping, methods combining accuracy and speed are still lacking. Here, we introduce MRT-ModSeq for rapid, simultaneous detection of multiple RNA modifications using MarathonRT. MRT-ModSeq employs distinct divalent cofactors to generate 2-D mutational profiles that are highly dependent on nucleotide identity and modification type. As a proof of concept, we use the MRT fingerprints of well-studied rRNAs to implement a general workflow for detecting RNA modifications. MRT-ModSeq rapidly detects positions of diverse modifications across a RNA transcript, enabling assignment of m1acp3Y, m1A, m3U, m7G and 2'-OMe locations through mutation-rate filtering and machine learning. m1A sites in sparsely modified targets, such as MALAT1 and PRUNE1 could also be detected. MRT-ModSeq can be trained on natural and synthetic transcripts to expedite detection of diverse RNA modification subtypes across targets of interest.
Subject(s)
RNA Processing, Post-Transcriptional , RNA, Ribosomal , Mutation , RNA, Ribosomal/chemistry , RNA, Ribosomal/genetics , Sequence Analysis, RNA/methods , HumansABSTRACT
Here, we describe the discovery of compounds that inhibit self-splicing in group II introns. Using docking calculations, we targeted the catalytic active site within the Oceanobacillus iheyensis group IIC intron and virtually screened a library of lead-like compounds. From this initial virtual screen, we identified three unique scaffolds that inhibit splicing in vitro. Additional tests revealed that an analog of the lead scaffold inhibits splicing in an intron-dependent manner. Furthermore, this analog exhibited activity against the group II intron from a different class: the yeast ai5γ IIB intron. The splicing inhibitors we identified could serve as chemical tools for developing group II intron-targeted antifungals, and, more broadly, our results highlight the potential of in silico techniques for identifying bioactive hits against structured and functionally complex RNAs.
Subject(s)
Antifungal Agents , RNA , Introns , Catalysis , Catalytic Domain , RNA Splicing , Saccharomyces cerevisiaeABSTRACT
RNA tertiary structures from experiments or computational predictions often contain missing atoms, which prevent analyses requiring full atomic structures. Current programs for RNA reconstruction can be slow, inaccurate, and/or require specific atoms to be present in the input. We present Arena (Atomic Reconstruction of RNA), which reconstructs a full atomic RNA structure from residues that can have as few as one atom. Arena first fills in missing atoms and then iteratively refines their placement to reduce nonideal geometries. We benchmarked Arena on a dataset of 361 RNA structures, where Arena achieves high accuracy and speed compared to other structure reconstruction programs. For example, Arena was used to reconstruct full atomic structures from a single phosphorus atom per nucleotide to, on average, within 3.63 Å RMSD of the experimental structure, while virtually removing all clashes and running in <3 s, which is 353× and 46× faster than state-of-the-art programs PDBFixer and C2A, respectively. The Arena source code is available at https://github.com/pylelab/Arena and the webserver at https://zhanggroup.org/Arena/.
Subject(s)
RNA , Software , RNA/chemistry , Models, Molecular , Nucleotides , Nucleic Acid ConformationABSTRACT
Chemical modifications are essential regulatory elements that modulate the behavior and function of cellular RNAs. Despite recent advances in sequencing-based RNA modification mapping, methods combining accuracy and speed are still lacking. Here, we introduce MRT- ModSeq for rapid, simultaneous detection of multiple RNA modifications using MarathonRT. MRT-ModSeq employs distinct divalent cofactors to generate 2-D mutational profiles that are highly dependent on nucleotide identity and modification type. As a proof of concept, we use the MRT fingerprints of well-studied rRNAs to implement a general workflow for detecting RNA modifications. MRT-ModSeq rapidly detects positions of diverse modifications across a RNA transcript, enabling assignment of m1acp3Y, m1A, m3U, m7G and 2'-OMe locations through mutation-rate filtering and machine learning. m1A sites in sparsely modified targets, such as MALAT1 and PRUNE1 could also be detected. MRT-ModSeq can be trained on natural and synthetic transcripts to expedite detection of diverse RNA modification subtypes across targets of interest.
ABSTRACT
The multiple sequence alignment (MSA) is the entry point of many RNA structure modeling tasks, such as prediction of RNA secondary structure (rSS) and contacts. However, there are few automated programs for generating high quality MSAs of target RNA molecules. We have developed rMSA, a hierarchical pipeline for sensitive search and accurate alignment of RNA homologs for a target RNA. On a diverse set of 365 non-redundant RNA structures, rMSA significantly outperforms an existing MSA generation method (RNAcmap) by approximately 20% and 5% higher F1-scores for rSS and long-range contact prediction, respectively. rMSA is available at https://zhanggroup.org/rMSA/ and https://github.com/pylelab/rMSA.
Subject(s)
Nucleic Acid Conformation , RNA , Software , Algorithms , RNA/genetics , RNA/chemistry , Sequence AlignmentABSTRACT
Compact RNA structural motifs control many aspects of gene expression, but we lack methods for finding these structures in the vast expanse of multi-kilobase RNAs. To adopt specific 3-D shapes, many RNA modules must compress their RNA backbones together, bringing negatively charged phosphates into close proximity. This is often accomplished by recruiting multivalent cations (usually Mg2+), which stabilize these sites and neutralize regions of local negative charge. Coordinated lanthanide ions, such as terbium (III) (Tb3+), can also be recruited to these sites, where they induce efficient RNA cleavage, thereby revealing compact RNA 3-D modules. Until now, Tb3+ cleavage sites were monitored via low-throughput biochemical methods only applicable to small RNAs. Here we present Tb-seq, a high-throughput sequencing method for detecting compact tertiary structures in large RNAs. Tb-seq detects sharp backbone turns found in RNA tertiary structures and RNP interfaces, providing a way to scan transcriptomes for stable structural modules and potential riboregulatory motifs.
Subject(s)
RNA , Terbium , Nucleic Acid Conformation , RNA/metabolism , Terbium/metabolism , Terbium/pharmacology , Nucleotide Motifs , CationsABSTRACT
Recent improvements in computational and experimental techniques for obtaining protein structures have resulted in an explosion of 3D coordinate data. To cope with the ever-increasing sizes of structure databases, this work proposes the Protein Data Compression (PDC) format, which compresses coordinates and temperature factors of full-atomic and Cα-only protein structures. Without loss of precision, PDC results in 69% to 78% smaller file sizes than Protein Data Bank (PDB) and macromolecular Crystallographic Information File (mmCIF) files with standard GZIP compression. It uses â¼60% less space than existing compression algorithms specific to macromolecular structures. PDC optionally performs lossy compression with minimal sacrifice of precision, which allows reduction of file sizes by another 79%. Conversion between PDC, mmCIF and PDB formats is typically achieved within 0.02 s. The compactness and fast reading/writing speed of PDC make it valuable for storage and analysis of large quantity of tertiary structural data. Database URL https://github.com/kad-ecoli/pdc.
Subject(s)
Data Compression , Software , Algorithms , Proteins/chemistry , Databases, ProteinABSTRACT
Capturing different conformations of receptor proteins that are complexed with ligands by single-particle cryo-EM facilitates our understanding toward the mechanisms of ligand recognition and receptor activation cascades. Here, we present a protocol for capturing RNA-sensing innate immune receptors, such as RIG-I, in multiple conformations by single-particle cryo-EM. We describe steps for protein-ligand sample preparation, data acquisition, and image processing covering focused three-dimensional classification. This protocol can be adapted to capture the dynamic behavior of other receptors that can be stabilized. For complete details on the use and execution of this protocol, please refer to Wang and Pyle (2022).1.