Polyphasic characterization of a novel hot-spring cyanobacterium Thermocoleostomius sinensis gen et sp. nov. and genomic insights into its carbon concentration mechanism

Thermophilic cyanobacteria play a crucial role as primary producers in hot spring ecosystems, yet their microbiological, taxonomic, and ecological characteristics are not extensively studied. This study aimed to characterize a novel strain of thermophilic cyanobacteria, PKUAC-SCTA174 (A174), using a combination of traditional polyphasic methods and modern genomic-based approaches. The study included 16S rRNA-based phylogeny, ITS secondary structure prediction, morphological and habitat analyses, as well as high-quality genome sequencing with corresponding phylogenomic analyses. The results of the 16S rRNA, 16S-23S ITS secondary structure, morphology, and habitat analyses supported the classification of the strain as a member of a novel genus within the family Oculatellaceae, closely related to Albertania and Trichotorquatus. Genomic analysis revealed the presence of a sophisticated carbon-concentrating mechanism (CCM) in the strain, involving two CO2 uptake systems NDH-I3, and NDH-I4, three types of bicarbonate transporters (BCT1, bicA, sbtA,) and two distinct putative carboxysomal carbonic anhydrases (ccaA1 and ccaA2). The expression of CCM genes was investigated with a CO2 shift experiment, indicating varying transcript abundance among different carbon uptake systems. Based on the comprehensive characterization, the strain was delineated as Thermocoleostomius sinensis, based on the botanical code. The study of the complete genome of strain A174 contributes valuable insights into the genetic characteristics of the genus Thermocoleostomius and related organisms and provides a systematic understanding of thermophilic cyanobacteria. The findings presented here offer valuable data that can be utilized for future research in taxogenomics, ecogenomics, and geogenomics.

In vivo secondary structural analysis of Influenza A virus genomic RNA

Influenza A virus (IAV) is a respiratory virus that causes epidemics and pandemics. Knowledge of IAV RNA secondary structure in vivo is crucial for a better understanding of virus biology. Moreover, it is a fundament for the development of new RNA-targeting antivirals. Chemical RNA mapping using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) coupled with Mutational Profiling (MaP) allows for the thorough examination of secondary structures in low-abundance RNAs in their biological context. So far, the method has been used for analyzing the RNA secondary structures of several viruses including SARS-CoV-2 in virio and in cellulo. Here, we used SHAPE-MaP and dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) for genome-wide secondary structure analysis of viral RNA (vRNA) of the pandemic influenza A/California/04/2009 (H1N1) strain in both in virio and in cellulo environments. Experimental data allowed the prediction of the secondary structures of all eight vRNA segments in virio and, for the first time, the structures of vRNA5, 7, and 8 in cellulo. We conducted a comprehensive structural analysis of the proposed vRNA structures to reveal the motifs predicted with the highest accuracy. We also performed a base-pairs conservation analysis of the predicted vRNA structures and revealed many highly conserved vRNA motifs among the IAVs. The structural motifs presented herein are potential candidates for new IAV antiviral strategies.

SHAPE to Probe RNA Structure and RNA-Protein Interactions In Vitro

Selective 2' hydroxyl acylation analyzed by primer extension (SHAPE) is used to distinguish between the levels of flexibility of nucleotides regulated by base pairing or protein binding. In this method, a reagent reacts with the 2' hydroxyl group to form an adduct, which is then detected by reverse transcription reaction. The number of RNA molecules with an adduct at a specific nucleotide position indicates the SHAPE reactivity of that nucleotide. Here, we describe the method for probing the structure of an RNA in a protein-free or a protein-bound state by in vitro SHAPE.

Cooperative engagement and subsequent selective displacement of SR proteins define the pre-mRNA 3D structural scaffold for early spliceosome assembly

We recently reported that serine–arginine-rich (SR) protein-mediated pre-mRNA structural remodeling generates a pre-mRNA 3D structural scaffold that is stably recognized by the early spliceosomal components. However, the intermediate steps between the free pre-mRNA and the assembled early spliceosome are not yet characterized. By probing the early spliceosomal complexes in vitro and RNA-protein interactions in vivo, we show that the SR proteins bind the pre-mRNAs cooperatively generating a substrate that recruits U1 snRNP and U2AF65 in a splice signal-independent manner. Excess U1 snRNP selectively displaces some of the SR protein molecules from the pre-mRNA generating the substrate for splice signal-specific, sequential recognition by U1 snRNP, U2AF65 and U2AF35. Our work thus identifies a novel function of U1 snRNP in mammalian splicing substrate definition, explains the need for excess U1 snRNP compared to other U snRNPs in vivo, demonstrates how excess SR proteins could inhibit splicing, and provides a conceptual basis to examine if this mechanism of splicing substrate definition is employed by other splicing regulatory proteins.

Polyphasic Identification and Genomic Insights of Leptothermofonsia sichuanensis gen. sp. nov., a Novel Thermophilic Cyanobacteria Within Leptolyngbyaceae

Thermal environments are an important reservoir of thermophiles with significant ecological and biotechnological potentials. However, thermophilic isolates remain largely unrecovered from their habitats and are rarely systematically identified. In this study, we characterized using polyphasic approaches a thermophilic strain, PKUAC-SCTAE412 (E412 hereafter), recovered from Lotus Lake hot spring based in Ganzi prefecture, China. The results of 16S rRNA/16S-23S ITS phylogenies, secondary structure, and morphology comparison strongly supported that strain E412 represent a novel genus within Leptolyngbyaceae. This delineation was further confirmed by genome-based analyses [phylogenomic inference, average nucleotide/amino-acid identity, and the percentages of conserved proteins (POCP)]. Based on the botanical code, the isolate is herein delineated as Leptothermofonsia sichuanensis gen. sp. nov, a genus adjacent to recently delineated Kovacikia and Stenomitos. In addition, we successfully obtained the first complete genome of this new genus. Genomic analysis revealed its adaptations to the adverse hot spring environment and extensive molecular components related to mobile genetic elements, photosynthesis, and nitrogen metabolism. Moreover, the strain was capable of modifying the composition of its light-harvesting apparatus depending on the wavelength and photoperiod, showing chromatic adaptation capacity characteristic for T1 and T2 pigmentation types. Other physiological studies showed the strain’s ability to utilize sodium bicarbonate and various sulfur compounds. The strain was also shown to be diazotrophic. Interestingly, 24.6% of annotated protein-coding genes in the E412 genome were identified as putatively acquired, hypothesizing that a large number of genes acquired through HGT might contribute to the genome expansion and habitat adaptation of those thermophilic strains. Most the HGT candidates (69.4%) were categorized as metabolic functions as suggested by the KEGG analysis. Overall, the complete genome of strain E412 provides the first insight into the genomic feature of the genus Leptothermofonsia and lays the foundation for future global ecogenomic and geogenomic studies.

Nuclear magnetic resonance reveals a two hairpin equilibrium near the 3'-splice site of influenza A segment 7 mRNA that can be shifted by oligonucleotides

Influenza A kills hundreds of thousands of people globally every year and has the potential to generate more severe pandemics. Influenza A's RNA genome and transcriptome provide many potential therapeutic targets. Here, nuclear magnetic resonance (NMR) experiments suggest that one such target could be a hairpin loop of 8 nucleotides in a pseudoknot that sequesters a 3' splice site in canonical pairs until a conformational change releases it into a dynamic 2 × 2-nt internal loop. NMR experiments reveal that the hairpin loop is dynamic and able to bind oligonucleotides as short as pentamers. A 3D NMR structure of the complex contains 4 and likely 5 bp between pentamer and loop. Moreover, a hairpin sequence was discovered that mimics the equilibrium of the influenza hairpin between its structure in the pseudoknot and upon release of the splice site. Oligonucleotide binding shifts the equilibrium completely to the hairpin secondary structure required for pseudoknot folding. The results suggest this hairpin can be used to screen for compounds that stabilize the pseudoknot and potentially reduce splicing.

Secondary Structure of Influenza A Virus Genomic Segment 8 RNA Folded in a Cellular Environment

Influenza A virus (IAV) is a member of the single-stranded RNA (ssRNA) family of viruses. The most recent global pandemic caused by the SARS-CoV-2 virus has shown the major threat that RNA viruses can pose to humanity. In comparison, influenza has an even higher pandemic potential as a result of its high rate of mutations within its relatively short (<13 kbp) genome, as well as its capability to undergo genetic reassortment. In light of this threat, and the fact that RNA structure is connected to a broad range of known biological functions, deeper investigation of viral RNA (vRNA) structures is of high interest. Here, for the first time, we propose a secondary structure for segment 8 vRNA (vRNA8) of A/California/04/2009 (H1N1) formed in the presence of cellular and viral components. This structure shows similarities with prior in vitro experiments. Additionally, we determined the location of several well-defined, conserved structural motifs of vRNA8 within IAV strains with possible functionality. These RNA motifs appear to fold independently of regional nucleoprotein (NP)-binding affinity, but a low or uneven distribution of NP in each motif region is noted. This research also highlights several accessible sites for oligonucleotide tools and small molecules in vRNA8 in a cellular environment that might be a target for influenza A virus inhibition on the RNA level.

A Small Molecule Exploits Hidden Structural Features within the RNA Repeat Expansion That Causes c9ALS/FTD and Rescues Pathological Hallmarks

The hexanucleotide repeat expansion GGGGCC [r(G₄C₂)^exp] within intron 1 of C9orf72 causes genetically defined amyotrophic lateral sclerosis and frontotemporal dementia, collectively named c9ALS/FTD. , the repeat expansion causes neurodegeneration via deleterious phenotypes stemming from r(G₄C₂)^exp RNA gain- and loss-of-function mechanisms. The r(G₄C₂)^exp RNA folds into both a hairpin structure with repeating 1 × 1 nucleotide GG internal loops and a G-quadruplex structure. Here, we report the identification of a small molecule (CB253) that selectively binds the hairpin form of r(G₄C₂)^exp. Interestingly, the small molecule binds to a previously unobserved conformation in which the RNA forms 2 × 2 nucleotide GG internal loops, as revealed by a series of binding and structural studies. NMR and molecular dynamics simulations suggest that the r(G₄C₂)^exp hairpin interconverts between 1 × 1 and 2 × 2 internal loops through the process of strand slippage. We provide experimental evidence that CB253 binding indeed shifts the equilibrium toward the 2 × 2 GG internal loop conformation, inhibiting mechanisms that drive c9ALS/FTD pathobiology, such as repeat-associated non-ATG translation formation of stress granules and defective nucleocytoplasmic transport in various cellular models of c9ALS/FTD.

Discovery of a pre-mRNA structural scaffold as a contributor to the mammalian splicing code

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML-a model pre-mRNA substrate-guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5' and 3' splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

A map of the SARS-CoV-2 RNA structurome

SARS-CoV-2 has exploded throughout the human population. To facilitate efforts to gain insights into SARS-CoV-2 biology and to target the virus therapeutically, it is essential to have a roadmap of likely functional regions embedded in its RNA genome. In this report, we used a bioinformatics approach, ScanFold, to deduce the local RNA structural landscape of the SARS-CoV-2 genome with the highest likelihood of being functional. We recapitulate previously-known elements of RNA structure and provide a model for the folding of an essential frameshift signal. Our results find that SARS-CoV-2 is greatly enriched in unusually stable and likely evolutionarily ordered RNA structure, which provides a large reservoir of potential drug targets for RNA-binding small molecules. Results are enhanced via the re-analyses of publicly-available genome-wide biochemical structure probing datasets that are broadly in agreement with our models. Additionally, ScanFold was updated to incorporate experimental data as constraints in the analysis to facilitate comparisons between ScanFold and other RNA modelling approaches. Ultimately, ScanFold was able to identify eight highly structured/conserved motifs in SARS-CoV-2 that agree with experimental data, without explicitly using these data. All results are made available via a public database (the RNAStructuromeDB: https://structurome.bb.iastate.edu/sars-cov-2) and model comparisons are readily viewable at https://structurome.bb.iastate.edu/sars-cov-2-global-model-comparisons.

Scale-up fermentation of Escherichia coli for the production of recombinant endoglucanase from Clostridium thermocellum

Endoglucanase (EC 3.2.1.4) catalysing the hydrolysis of β-1.4-glycosidic linkage of cellulose molecules is an enzyme of tremendous industrial importance. The present study describes a response surface methodology based predicted model to deduce a set of fermentation conditions for optimum growth and activity of recombinant endoglucanase in E. coli BL21 (DE3). Numerous significant parameters including fermentation media composition, temperature (Celsius), pH and agitation rate (rpm) were analysed systemically by employing central composite design. This effort reports highly efficient recombinant endoglucanase overproduction (6.9 gl^-1 of biomass) with 30% expression by E. coli in modified M9NG media incubated at 37 °C and pH 7 agitated at 200 rpm. Addition of 3 mM glucose and 24 mM glycerol in the M9NG media has shown positive effect on the enzyme yield and activity. The CMCase activity experimentally estimated was found to be 1185 U/mg with the optimized parameters. The outcomes of both the responses by the predicted quadratic model were found in consensus with the obtained values. Our results well depicted the favourable conditions to further scale-up the volumetric yield of other relevant recombinant enzymes and proteins.

Relative Information Gain: Shannon entropy-based measure of the relative structural conservation in RNA alignments

Structural characterization of RNAs is a dynamic field, offering many modelling possibilities. RNA secondary structure models are usually characterized by an encoding that depicts structural information of the molecule through string representations or graphs. In this work, we provide a generalization of the BEAR encoding (a context-aware structural encoding we previously developed) by expanding the set of alignments used for the construction of substitution matrices and then applying it to secondary structure encodings ranging from fine-grained to more coarse-grained representations. We also introduce a re-interpretation of the Shannon Information applied on RNA alignments, proposing a new scoring metric, the Relative Information Gain (RIG). The RIG score is available for any position in an alignment, showing how different levels of detail encoded in the RNA representation can contribute differently to convey structural information. The approaches presented in this study can be used alongside state-of-the-art tools to synergistically gain insights into the structural elements that RNAs and RNA families are composed of. This additional information could potentially contribute to their improvement or increase the degree of confidence in the secondary structure of families and any set of aligned RNAs.

Structural Features of Small Molecules Targeting the RNA Repeat Expansion That Causes Genetically Defined ALS/FTD

Genetically defined amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), collectively named c9ALS/FTD, are triggered by hexanucleotide GGGGCC repeat expansions [r(G₄C₂)^exp] within the C9orf72 gene. In these diseases, neuronal loss occurs through an interplay of deleterious phenotypes, including r(G₄C₂)^exp RNA gain-of-function mechanisms. Herein, we identified a benzimidazole derivative, CB096, that specifically binds to a repeating 1 × 1 GG internal loop structure, 5'CGG/3'GGC, that is formed when r(G₄C₂)^exp folds. Structure-activity relationship (SAR) studies and molecular dynamics (MD) simulations were used to define the molecular interactions formed between CB096 and r(G₄C₂)^exp that results in the rescue of disease-associated pathways. Overall, this study reveals a unique structural feature within r(G₄C₂)^exp that can be exploited for the development of lead medicines and chemical probes.

Structural disruption of exonic stem-loops immediately upstream of the intron regulates mammalian splicing

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Intronic TP53 Polymorphisms Are Associated with Increased Δ133TP53 Transcript, Immune Infiltration and Cancer Risk

We investigated the influence of selected TP53 SNPs in exon 4 and intron 4 on cancer risk, clinicopathological features and expression of TP53 isoforms. The intron 4 SNPs were significantly over-represented in cohorts of mixed cancers compared to three ethnically matched controls, suggesting they confer increased cancer risk. Further analysis showed that heterozygosity at rs1042522(GC) and either of the two intronic SNPs rs9895829(TC) and rs2909430(AG) confer a 2.34-5.35-fold greater risk of developing cancer. These SNP combinations were found to be associated with shorter patient survival for glioblastoma and prostate cancer. Additionally, these SNPs were associated with tumor-promoting inflammation as evidenced by high levels of infiltrating immune cells and expression of the Δ133TP53 and TP53β transcripts. We propose that these SNP combinations allow increased expression of the Δ133p53 isoforms to promote the recruitment of immune cells that create an immunosuppressive environment leading to cancer progression.

Genome and sequence determinants governing the expression of horizontally acquired DNA in bacteria

While horizontal gene transfer is prevalent across the biosphere, the regulatory features that enable expression and functionalization of foreign DNA remain poorly understood. Here, we combine high-throughput promoter activity measurements and large-scale genomic analysis of regulatory regions to investigate the cross-compatibility of regulatory elements (REs) in bacteria. Functional characterization of thousands of natural REs in three distinct bacterial species revealed distinct expression patterns according to RE and recipient phylogeny. Host capacity to activate foreign promoters was proportional to their genomic GC content, while many low GC regulatory elements were both broadly active and had more transcription start sites across hosts. The difference in expression capabilities could be explained by the influence of the host GC content on the stringency of the AT-rich canonical σ70 motif necessary for transcription initiation. We further confirm the generalizability of this model and find widespread GC content adaptation of the σ70 motif in a set of 1,545 genomes from all major bacterial phyla. Our analysis identifies a key mechanism by which the strength of the AT-rich σ70 motif relative to a host's genomic GC content governs the capacity for expression of acquired DNA. These findings shed light on regulatory adaptation in the context of evolving genomic composition.

A domain-level DNA strand displacement reaction enumerator allowing arbitrary non-pseudoknotted secondary structures

Information technologies enable programmers and engineers to design and synthesize systems of startling complexity that nonetheless behave as intended. This mastery of complexity is made possible by a hierarchy of formal abstractions that span from high-level programming languages down to low-level implementation specifications, with rigorous connections between the levels. DNA nanotechnology presents us with a new molecular information technology whose potential has not yet been fully unlocked in this way. Developing an effective hierarchy of abstractions may be critical for increasing the complexity of programmable DNA systems. Here, we build on prior practice to provide a new formalization of 'domain-level' representations of DNA strand displacement systems that has a natural connection to nucleic acid biophysics while still being suitable for formal analysis. Enumeration of unimolecular and bimolecular reactions provides a semantics for programmable molecular interactions, with kinetics given by an approximate biophysical model. Reaction condensation provides a tractable simplification of the detailed reactions that respects overall kinetic properties. The applicability and accuracy of the model is evaluated across a wide range of engineered DNA strand displacement systems. Thus, our work can serve as an interface between lower-level DNA models that operate at the nucleotide sequence level, and high-level chemical reaction network models that operate at the level of interactions between abstract species.

Structural transition of replicable RNAs during in vitro evolution with Qβ replicase

Single-stranded RNAs (ssRNAs) are utilized as genomes in some viruses and also in experimental models of ancient life-forms, owing to their simplicity. One of the largest problems for ssRNA replication is the formation of double-stranded RNA (dsRNA), a dead-end product for ssRNA replication. A possible strategy to avoid dsRNA formation is to create strong intramolecular secondary structures of ssRNA. To design ssRNAs that efficiently replicate by Qβ replicase with minimum dsRNA formation, we previously proposed the "fewer unpaired GC rule." According to this rule, ssRNAs that have fewer unpaired G and C bases in the secondary structure should efficiently replicate with less dsRNA formation. However, the validity of this rule still needs to be examined, especially for longer ssRNAs. Here, we analyze nine long ssRNAs that successively appeared during an in vitro evolution of replicable ssRNA by Qβ replicase and examine whether this rule can explain the structural transitions of the RNAs. We found that these ssRNAs improved their template abilities step-by-step with decreasing dsRNA formation as mutations accumulated. We then examine the secondary structures of all the RNAs by a chemical modification method. The analysis of the structures revealed that the probabilities of unpaired G and C bases tended to decrease gradually in the course of evolution. The decreases were caused by the local structural changes around the mutation sites in most of the cases. These results support the validity of the "fewer unpaired GC rule" to efficiently design replicable ssRNAs by Qβ replicase, useful for more complex ssRNA replication systems.

Targeted Delivery of C/EBPα-saRNA by RNA Aptamers Shows Anti-tumor Effects in a Mouse Model of Advanced PDAC

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive malignancies; it preferentially metastasizes to the liver and is the main cause of death from this disease. In previous studies, small activating RNA against CCAAT/enhancer-binding protein-α (C/EBPα-saRNA) demonstrated efficacy of PDAC in a local subcutaneous tumor model. In this study, we focused on the efficacy of C/EBPα-saRNA in advanced stage PDAC. For targeted delivery, we selected a new anti-transferrin receptor aptamer (TR14), which demonstrated a high binding affinity to target proteins. The TR14 aptamer was internalized with clathrin-mediated endocytosis, distributed in early endosome, late endosome, and lysosome subcellularly. To investigate its anti-tumor effects to advanced PDAC, we conjugated C/EBPα-saRNA to TR14. Treatment of pancreatic cancer cells with the conjugates upregulated expression of C/EBPα and its downstream target p21, and inhibited cell proliferation. For in vivo assays, we established an advanced PDAC mouse model by engrafting luciferase reporter-PANC-1 cells directly into the livers of non-obese diabetic/severe combined immunodeficiency (NOD/SCID) mice. After treatment of aptamer-C/EBPα conjugates, we observed significant reduction of tumor growth in this advanced PDAC mouse model. Combinational treatment of the conjugates with gemcitabine also demonstrated enhanced anti-tumor effects in advanced PDAC. This suggests that aptamer-C/EBPα conjugates could be used as an adjuvant, along with other conventional anti-cancer drugs in advanced PDAC. In conclusion, targeted delivery of C/EBPα-saRNAs by aptamers might have potential therapeutic effects in advanced PDAC.

Massively parallel RNA device engineering in mammalian cells with RNA-Seq

Synthetic RNA-based genetic devices dynamically control a wide range of gene-regulatory processes across diverse cell types. However, the limited throughput of quantitative assays in mammalian cells has hindered fast iteration and interrogation of sequence space needed to identify new RNA devices. Here we report developing a quantitative, rapid and high-throughput mammalian cell-based RNA-Seq assay to efficiently engineer RNA devices. We identify new ribozyme-based RNA devices that respond to theophylline, hypoxanthine, cyclic-di-GMP, and folinic acid from libraries of ~22,700 sequences in total. The small molecule responsive devices exhibit low basal expression and high activation ratios, significantly expanding our toolset of highly functional ribozyme switches. The large datasets obtained further provide conserved sequence and structure motifs that may be used for rationally guided design. The RNA-Seq approach offers a generally applicable strategy for developing broad classes of RNA devices, thereby advancing the engineering of genetic devices for mammalian systems.

Identifying and validating small molecules interacting with RNA (SMIRNAs)

High throughput sequencing has revolutionized our ability to identify aberrant RNA expression and mutations that cause or contribute to disease. These data can be used directly to design oligonucleotide-based modalities using Watson-Crick pairing to target unstructured regions in an RNA. A complementary, although more difficult, strategy to deactivate a malfunctioning RNA is to target highly structured regions with small molecules. Indeed, RNA structures are directly causative of disease. Herein, we discuss emerging strategies to design high affinity, selective, bioactive ligands targeting RNA, or small molecules interacting with RNA (SMIRNAs), and target validation and profiling methods. An experimental foundation is required for a lead identification strategy for RNA structures, constructed from a library-vs.-library screen that probes vast libraries of small molecules for binding RNA three dimensional folds. Dubbed 2-dimensional combinatorial screening (2DCS), the resulting data can be mined against transcriptomes or the composite of RNAs that are produced in an organism to define folded RNA structures that can be targeted. By applying SMIRNAs to cells and using target validation tools such as Chemical Cross-Linking and Isolation by Pull-down (Chem-CLIP) and Small Molecule Nucleic Acid Profiling by Cleavage Applied to RNA (RiboSNAP), all targets engaged in cells can be defined, along with rules for molecular recognition to affect RNA biology. This chapter will describe lessons learned in applying these approaches in vitro, in cells, and in pre-clinical animal models of disease, enabling SMIRNAs to capture opportunities in chemical biology.

In Vivo Use of a Multi-DNA Aptamer-Based Payload/Targeting System To Study Dopamine Dysregulation in the Central Nervous System

The delivery of therapeutics across the blood-brain barrier remains a considerable challenge in investigating central nervous system related processes. In this work, a liposome vehicle was surface-modified with an aptamer that binds to the transferrin receptor and was loaded with two different dopamine-binding aptamer payloads. This system was effectively used to promote the delivery of the aptamer cargo from the peripheral injection site into the brain. The effect of these delivered aptamers on behavior was investigated in vivo in a locomotor task. The first dopamine binding aptamer assessed was a DNA aptamer, the binding of which had been previously validated through the aptamer-based biosensor development reported by several independent research groups. The second aptamer investigated was the result of a novel in vitro selection experiment described herein. Our data suggest that systemic administration of the modified liposomes led to delivery of the dopamine aptamers into the brain. Fluorescence microscopy revealed differential distribution of fluorescence based on the presence or absence of the transferrin receptor aptamer on the surface of fluorescently modified liposomes. In a behavioral experiment using cocaine administration to induce elevated concentrations of neural dopamine, systemic pretreatment with the dopamine aptamer-loaded liposomes reduced cocaine-induced hyperlocomotion. Multiple controls including a transferrin-negative liposome control and transferrin-positive liposomes loaded with either a nonbinding, base-substituted dopamine aptamer or a random oligonucleotide were investigated. None of these controls altered cocaine-induced hyperlocomotion. Chronic systemic administration of the modified liposomes produced no deleterious neurobehavioral or neural degenerative effects. Importantly, this work is one example of an application for this versatile multiaptamer payload/targeting system. Its general application is limited only by the availability of aptamers for specific neural targets.

Human T-cell leukemia virus type 1 Gag domains have distinct RNA-binding specificities with implications for RNA packaging and dimerization

Human T-cell leukemia virus type 1 (HTLV-1) is the first retrovirus that has conclusively been shown to cause human diseases. In HIV-1, specific interactions between the nucleocapsid (NC) domain of the Gag protein and genomic RNA (gRNA) mediate gRNA dimerization and selective packaging; however, the mechanism for gRNA packaging in HTLV-1, a deltaretrovirus, is unclear. In other deltaretroviruses, the matrix (MA) and NC domains of Gag are both involved in gRNA packaging, but MA binds nucleic acids with higher affinity and has more robust chaperone activity, suggesting that this domain may play a primary role. Here, we show that the MA domain of HTLV-1, but not the NC domain, binds short hairpin RNAs derived from the putative gRNA packaging signal. RNA probing of the HTLV-1 5' leader and cross-linking studies revealed that the primer-binding site and a region within the putative packaging signal form stable hairpins that interact with MA. In addition to a previously identified palindromic dimerization initiation site (DIS), we identified a new DIS in HTLV-1 gRNA and found that both palindromic sequences bind specifically the NC domain. Surprisingly, a mutant partially defective in dimer formation in vitro exhibited a significant increase in RNA packaging into HTLV-1-like particles, suggesting that efficient RNA dimerization may not be strictly required for RNA packaging in HTLV-1. Moreover, the lifecycle of HTLV-1 and other deltaretroviruses may be characterized by NC and MA functions that are distinct from those of the corresponding HIV-1 proteins, but together provide the functions required for viral replication.

Disruption of a stem-loop structure located upstream of pseudoknot domain in Tobacco mosaic virus enhanced its infectivity and viral RNA accumulation

A predicted stem-loop structure of 25 nucleotides, located in the coat protein (CP) gene and 3'-UTR sequences of Tobacco mosaic virus (TMV), was validated previously (Guo et al., 2015). In this study, both disrupted stem-loop and nucleotide deletion mutants of TMV replicated more rapidly in Nicotiana benthamiana protoplasts. The TMV mutant with a complete mirrored stem-loop structure showed similar level of viral RNA accumulation as TMV. Recovering the stem-loop structure also resulted in a similar replication level as TMV. All these mutants induced necrosis in N. benthamiana and assembled into typical rigid rod-shaped virions. TMV mutant without the stem-loop structure induced more local lesions in Chenopodium quinoa. When the putative stem-loop structure in Tomato mosaic virus (ToMV) was disrupted, the mutant also showed an enhanced virus replication. This suggests that the stem-loop structure of TMV is a new cis-acting element with a role in virus replication.

Metagenomic mining of regulatory elements enables programmable species-selective gene expression

Robust and predictably performing synthetic circuits rely on the use of well-characterized regulatory parts across different genetic backgrounds and environmental contexts. Here we report the large-scale metagenomic mining of thousands of natural 5' regulatory sequences from diverse bacteria, and their multiplexed gene expression characterization in industrially relevant microbes. We identified sequences with broad and host-specific expression properties that are robust in various growth conditions. We also observed substantial differences between species in terms of their capacity to utilize exogenous regulatory sequences. Finally, we demonstrate programmable species-selective gene expression that produces distinct and diverse output patterns in different microbes. Together, these findings provide a rich resource of characterized natural regulatory sequences and a framework that can be used to engineer synthetic gene circuits with unique and tunable cross-species functionality and properties, and also suggest the prospect of ultimately engineering complex behaviors at the community level.

Evaluation of global HIV/SIV envelope gp120 RNA structure and evolution within and among infected hosts

Lentiviral RNA genomes contain structural elements that play critical roles in viral replication. Although structural features of 5'-untranslated regions have been well characterized, attempts to identify important structures in other genomic regions by Selective 2'-Hydroxyl Acylation analyzed by Primer Extension (SHAPE) have led to conflicting structural and mechanistic conclusions. Previous approaches accounted neither for sequence heterogeneity that is ubiquitous in viral populations, nor for selective constraints operating at the protein level. We developed an approach that augments SHAPE with phylogenetic analyses and applied it to investigate structure in coding regions (cRNA) within the HIV and SIV envelope genes. Analysis of single-genome SHAPE data with phylogenetic information from diverse lentiviral sequences argues against the conservation of a putative global gp120 RNA structure but points to the existence of core RNA sub-structures. Our findings establish a framework for considering sequence heterogeneity and protein function in de novo RNA structure inference approaches.

JNSViewer-A JavaScript-based Nucleotide Sequence Viewer for DNA/RNA secondary structures

Many tools are available for visualizing RNA or DNA secondary structures, but there is scarce implementation in JavaScript that provides seamless integration with the increasingly popular web computational platforms. We have developed JNSViewer, a highly interactive web service, which is bundled with several popular tools for DNA/RNA secondary structure prediction and can provide precise and interactive correspondence among nucleotides, dot-bracket data, secondary structure graphs, and genic annotations. In JNSViewer, users can perform RNA secondary structure predictions with different programs and settings, add customized genic annotations in GFF format to structure graphs, search for specific linear motifs, and extract relevant structure graphs of sub-sequences. JNSViewer also allows users to choose a transcript or specific segment of Arabidopsis thaliana genome sequences and predict the corresponding secondary structure. Popular genome browsers (i.e., JBrowse and BrowserGenome) were integrated into JNSViewer to provide powerful visualizations of chromosomal locations, genic annotations, and secondary structures. In addition, we used StructureFold with default settings to predict some RNA structures for Arabidopsis by incorporating in vivo high-throughput RNA structure profiling data and stored the results in our web server, which might be a useful resource for RNA secondary structure studies in plants. JNSViewer is available at http://bioinfolab.miamioh.edu/jnsviewer/index.html.

PROBer Provides a General Toolkit for Analyzing Sequencing-Based Toeprinting Assays

A number of sequencing-based transcriptase drop-off assays have recently been developed to probe post-transcriptional dynamics of RNA-protein interaction, RNA structure, and RNA modification. Although these assays survey a diverse set of epitranscriptomic marks, we use the term toeprinting assays since they share methodological similarities. Their interpretation is predicated on addressing a similar computational challenge: how to learn isoform-specific chemical modification profiles in the face of complex read multi-mapping. We introduce PROBer, a statistical model and associated software, that addresses this challenge for the analysis of toeprinting assays. PROBer takes sequencing data as input and outputs estimated transcript abundances and isoform-specific modification profiles. Results on both simulated and biological data demonstrate that PROBer significantly outperforms individual methods tailored for specific toeprinting assays. Since the space of toeprinting assays is ever expanding and these assays are likely to be performed and analyzed together, we believe PROBer's unified data analysis solution will be valuable to the RNA community.

Structure of Ty1 Internally Initiated RNA Influences Restriction Factor Expression

The long-terminal repeat retrotransposon Ty1 is the most abundant mobile genetic element in many Saccharomyces cerevisiae isolates. Ty1 retrotransposons contribute to the genetic diversity of host cells, but they can also act as an insertional mutagen and cause genetic instability. Interestingly, retrotransposition occurs at a low level despite a high level of Ty1 RNA, even though S. cerevisiae lacks the intrinsic defense mechanisms that other eukaryotes use to prevent transposon movement. p22 is a recently discovered Ty1 protein that inhibits retrotransposition in a dose-dependent manner. p22 is a truncated form of Gag encoded by internally initiated Ty1i RNA that contains two closely-spaced AUG codons. Mutations of either AUG codon compromise p22 translation. We found that both AUG codons were utilized and that translation efficiency depended on the Ty1i RNA structure. Structural features that stimulated p22 translation were context dependent and present only in Ty1i RNA. Destabilization of the 5′ untranslated region (5′ UTR) of Ty1i RNA decreased the p22 level, both in vitro and in vivo. Our data suggest that protein factors such as Gag could contribute to the stability and translational activity of Ty1i RNA through specific interactions with structural motifs in the RNA.

Variants in mitochondrial tRNA gene may not be associated with thyroid carcinoma

Thyroid cancer is a very common form of endocrine system malignancy. To date, the molecular mechanism underlying thyroid cancer remains poorly understood. Studies of oncocytic tumors have led to a hypothesis which proposes that defects in oxidative phosphorylation (OX- PHOS) may result in a compensatory increase in mitochondrial replication and gene expression. As a result, mitochondrial DNA (mtDNA) mutation analysis has become a useful tool to explore the molecular basis of this disease. Among these mutations, mitochondrial transfer RNAs (mttRNAs) are the hot spots for pathogenic mutations associated with thyroid cancer. However, due to its high mutation rate, the role of mt-tRNA variants in thyroid cancer is still controversial. To address this problem, in this study, we reassessed seven reported mt-tRNA variants: tRNA^Asp G7521A, tRNA^Arg T10411C and T10463C, tRNA^Leu(CUN) A12308G, tRNA^Ile G4292C and C4312T, and tRNA^Ala T5655C, in clinical manifestations of thyroid cancer. We first performed the phylogenetic conservation analysis for these variants; moreover, we used a bioinformatic tool to compare the minimum free energy (G) of mt-tRNA with and without mutations. Most strikingly, none of these variants caused the significant change of the G between the wild-type and the mutant form, suggesting that they may not play an important roles in thyroid cancer. In addition, we screened the frequency of the “pathogenic” A12308G alternation in 300 patients with thyroid cancer and 200 healthy controls. We found that there were five patients and three control subjects carrying this variant. It seemed that the A12308G variant may be a common polymorphism in the human population. Taken together, our study indicated that variants in mt-tRNA genes may not play active roles in patients with thyroid cancer.

Inforna 2.0: A Platform for the Sequence-Based Design of Small Molecules Targeting Structured RNAs

The development of small molecules that target RNA is challenging yet, if successful, could advance the development of chemical probes to study RNA function or precision therapeutics to treat RNA-mediated disease. Previously, we described Inforna, an approach that can mine motifs (secondary structures) within target RNAs, which is deduced from the RNA sequence, and compare them to a database of known RNA motif-small molecule binding partners. Output generated by Inforna includes the motif found in both the database and the desired RNA target, lead small molecules for that target, and other related meta-data. Lead small molecules can then be tested for binding and affecting cellular (dys)function. Herein, we describe Inforna 2.0, which incorporates all known RNA motif-small molecule binding partners reported in the scientific literature, a chemical similarity searching feature, and an improved user interface and is freely available via an online web server. By incorporation of interactions identified by other laboratories, the database has been doubled, containing 1936 RNA motif-small molecule interactions, including 244 unique small molecules and 1331 motifs. Interestingly, chemotype analysis of the compounds that bind RNA in the database reveals features in small molecule chemotypes that are privileged for binding. Further, this updated database expanded the number of cellular RNAs to which lead compounds can be identified.

The Conservation and Function of RNA Secondary Structure in Plants

RNA transcripts fold into secondary structures via intricate patterns of base pairing. These secondary structures impart catalytic, ligand binding, and scaffolding functions to a wide array of RNAs, forming a critical node of biological regulation. Among their many functions, RNA structural elements modulate epigenetic marks, alter mRNA stability and translation, regulate alternative splicing, transduce signals, and scaffold large macromolecular complexes. Thus, the study of RNA secondary structure is critical to understanding the function and regulation of RNA transcripts. Here, we review the origins, form, and function of RNA secondary structure, focusing on plants. We then provide an overview of methods for probing secondary structure, from physical methods such as X-ray crystallography and nuclear magnetic resonance (NMR) imaging to chemical and nuclease probing methods. Combining these latter methods with high-throughput sequencing has enabled them to scale across whole transcriptomes, yielding tremendous new insights into the form and function of RNA secondary structure.

Genome-Wide Approaches for RNA Structure Probing

RNA molecules of all types fold into complex secondary and tertiary structures that are important for their function and regulation. Structural and catalytic RNAs such as ribosomal RNA (rRNA) and transfer RNA (tRNA) are central players in protein synthesis, and only function through their proper folding into intricate three-dimensional structures. Studies of messenger RNA (mRNA) regulation have also revealed that structural elements embedded within these RNA species are important for the proper regulation of their total level in the transcriptome. More recently, the discovery of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) has shed light on the importance of RNA structure to genome, transcriptome, and proteome regulation. Due to the relatively small number, high conservation, and importance of structural and catalytic RNAs to all life, much early work in RNA structure analysis mapped out a detailed view of these molecules. Computational and physical methods were used in concert with enzymatic and chemical structure probing to create high-resolution models of these fundamental biological molecules. However, the recent expansion in our knowledge of the importance of RNA structure to coding and regulatory RNAs has left the field in need of faster and scalable methods for high-throughput structural analysis. To address this, nuclease and chemical RNA structure probing methodologies have been adapted for genome-wide analysis. These methods have been deployed to globally characterize thousands of RNA structures in a single experiment. Here, we review these experimental methodologies for high-throughput RNA structure determination and discuss the insights gained from each approach.

Automated RNA 3D Structure Prediction with RNAComposer

RNAs adopt specific structures to perform their activities and these are critical to virtually all RNA-mediated processes. Because of difficulties in experimentally assessing structures of large RNAs using NMR, X-ray crystallography, or cryo-microscopy, there is currently great demand for new high-resolution 3D structure prediction methods. Recently we reported on RNAComposer, a knowledge-based method for the fully automated RNA 3D structure prediction from a user-defined secondary structure. RNAComposer method is especially suited for structural biology users. Since our initial report in 2012, both servers, freely available at http://rnacomposer.ibch.poznan.pl and http://rnacomposer.cs.put.poznan.pl have been often visited. Therefore this chapter provides guidance for using RNAComposer and discusses points that should be considered when predicting 3D RNA structure. An application example presents current scope and limitations of RNAComposer.

Native Purification and Analysis of Long RNAs

The purification and analysis of long noncoding RNAs (lncRNAs) in vitro is a challenge, particularly if one wants to preserve elements of functional structure. Here, we describe a method for purifying lncRNAs that preserves the cotranscriptionally derived structure. The protocol avoids the misfolding that can occur during denaturation-renaturation protocols, thus facilitating the folding of long RNAs to a native-like state. This method is simple and does not require addition of tags to the RNA or the use of affinity columns. LncRNAs purified using this type of native purification protocol are amenable to biochemical and biophysical analysis. Here, we describe how to study lncRNA global compaction in the presence of divalent ions at equilibrium using sedimentation velocity analytical ultracentrifugation and analytical size-exclusion chromatography as well as how to use these uniform RNA species to determine robust lncRNA secondary structure maps by chemical probing techniques like selective 2'-hydroxyl acylation analyzed by primer extension and dimethyl sulfate probing.

Thermodynamic and kinetic folding of riboswitches

Riboswitches are structured RNA regulatory elements located in the 5'-UTRs of mRNAs. Ligand-binding induces a structural rearrangement in these RNA elements, effecting events in downstream located coding sequences. Since they do not require proteins for their functions, they are ideally suited for computational analysis using the toolbox of RNA structure prediction methods. By their very definition riboswitch function depends on structural change. Methods that consider only the thermodynamic equilibrium of an RNA are therefore of limited use. Instead, one needs to employ computationally more expensive methods that consider the energy landscape and the folding dynamics on that landscape. Moreover, for the important class of kinetic riboswitches, the mechanism of riboswitch function can only be understood in the context of co-transcriptional folding. We present a computational approach to simulate the dynamic behavior of riboswitches during co-transcriptional folding in the presence and absence of a ligand. Our investigations show that the abstraction level of RNA secondary structure in combination with a dynamic folding landscape approach is expressive enough to understand how riboswitches perform their function. We apply our approach to a experimentally validated theophylline-binding riboswitch.

Improving RNA secondary structure prediction with structure mapping data

Methods to probe RNA secondary structure, such as small molecule modifying agents, secondary structure-specific nucleases, inline probing, and SHAPE chemistry, are widely used to study the structure of functional RNA. Computational secondary structure prediction programs can incorporate probing data to predict structure with high accuracy. In this chapter, an overview of current methods for probing RNA secondary structure is provided, including modern high-throughput methods. Methods for guiding secondary structure prediction algorithms using these data are explained, and best practices for using these data are provided. This chapter concludes by listing a number of open questions about how to best use probing data, and what these data can provide.