RNA Secondary Structure Prediction

Impact of Genomic Deletion RD16 on the Expression of the Mycobacterium bovis BCG Moreau VapBC47 Toxin-Antitoxin System

Mycobacterium bovis BCG is the only vaccine against tuberculosis. The variable forms of cultivation throughout the years, before seed-lots were developed, allowed in vitro evolution of the original strain, generating a family of vaccines with different phenotypic and genotypic characteristics. Molecular studies revealed regions of difference (RDs) in the genomes of the various BCG strains. This work aims to characterize the gene pair rv3407-rv3408 (vapB47-vapC47), coding for a toxin-antitoxin system of the VapBC family, and to evaluate possible transcriptional effects due to the adjacent BCG Moreau-specific genomic deletion RD16. We show that these genes are co-transcribed in BCG strains Moreau and Pasteur, and that the inactivation of an upstream transcriptional repressor (Rv3405c) due to RD16 has a polar effect, leading to increased vapBC47 expression. Furthermore, we detect VapB47 DNA binding in vitro, dependent on a 5' vapB47 sequence that contributes to a palindrome, spanning the promoter and coding region. Our data shed light on the regulation of VapBC systems and on the impact of the BCG Moreau RD16 deletion in the expression of adjacent genes, contributing to a better understanding of BCG Moreau physiology.

Antisense oligonucleotide is a promising intervention for liver diseases

As the body's critical metabolic organ, the liver plays an essential role in maintaining proper body homeostasis. However, as people's living standards have improved and the number of unhealthy lifestyles has increased, the liver has become overburdened. These have made liver disease one of the leading causes of death worldwide. Under the influence of adverse factors, liver disease progresses from simple steatosis to hepatitis, to liver fibrosis, and finally to cirrhosis and cancer, followed by increased mortality. Until now, there has been a lack of accepted effective treatments for liver disease. Based on current research, antisense oligonucleotide (ASO), as an alternative intervention for liver diseases, is expected to be an effective treatment due to its high efficiency, low toxicity, low dosage, strong specificity, and additional positive characteristics. In this review, we will first introduce the design, modification, delivery, and the mechanisms of ASO, and then summarize the application of ASO in liver disease treatment, including in non-alcoholic fatty liver disease (NAFLD), hepatitis, liver fibrosis, and liver cancer. Finally, we discuss challenges and perspectives on the transfer of ASO drugs into clinical use. This review provides a current and comprehensive understanding of the integrative and systematic functions of ASO for its use in liver disease.

Deep learning models for RNA secondary structure prediction (probably) do not generalize across families

MOTIVATION

The secondary structure of RNA is of importance to its function. Over the last few years, several papers attempted to use machine learning to improve de novo RNA secondary structure prediction. Many of these papers report impressive results for intra-family predictions but seldom address the much more difficult (and practical) inter-family problem.

RESULTS

We demonstrate that it is nearly trivial with convolutional neural networks to generate pseudo-free energy changes, modelled after structure mapping data that improve the accuracy of structure prediction for intra-family cases. We propose a more rigorous method for inter-family cross-validation that can be used to assess the performance of learning-based models. Using this method, we further demonstrate that intra-family performance is insufficient proof of generalization despite the widespread assumption in the literature and provide strong evidence that many existing learning-based models have not generalized inter-family.

AVAILABILITY AND IMPLEMENTATION

Source code and data are available at https://github.com/marcellszi/dl-rna.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Epithelial cells-enriched lncRNA SNHG8 regulates chromatin condensation by binding to Histone H1s

Linker histone H1 proteins contain many variants in mammalian and can stabilize the condensed state of chromatin by binding to nucleosomes and promoting a more inaccessible structure of DNA. However, it is poorly understood how the binding of histone H1s to chromatin DNA is regulated. Screened as one of a collection of epithelial cells-enriched long non-coding RNAs (lncRNAs), here we found that small nucleolar RNA host gene 8 (SNHG8) is a chromatin-localized lncRNA and presents strong interaction and phase separation with histone H1 variants. Moreover, SNHG8 presents stronger ability to bind H1s than linker DNA, and outcompetes linker DNA for H1 binding. Consequently, loss of SNHG8 increases the amount of H1s that bind to chromatin, promotes chromatin condensation, and induces an epithelial differentiation-associated gene expression pattern. Collectively, our results propose that the highly abundant SNHG8 in epithelial cells keeps histone H1 variants out of nucleosome and its loss contributes to epithelial cell differentiation.

Association of ABCG2 polymorphisms with susceptibility to anti-tuberculosis drug-induced hepatotoxicity in the Chinese population

Background The accumulation of endogenous hepatotoxin protoporphyrin IX (PPIX) in the liver was proposed to be a novel mechanism of anti-tuberculosis drug-induced hepatotoxicity (ATDH). ATP-binding cassette transporter G2 (ABCG2) plays an important role in modulating PPIX concentrations. This study aimed to explore the role of ABCG2 genetic polymorphisms in the risk of ATDH in Chinese patients.Methods A 1:4 matched case-control study was performed among 202 ATDH cases and 808 controls. Conditional logistic regression model was used to estimate the association between genotypes and the risk of ATDH by odds ratios (ORs) with 95% confidence intervals (CIs).Results Male patients with CC genotype of rs2622605 had an increased risk of ATDH (adjusted OR =1.615, 95% CI: 1.119-2.332, P = 0.011). The peak value of alkaline phosphatase was significantly higher in male patients with CC genotype of rs2622605 than in those with TT + TC genotype during antituberculosis treatment (102.0 U/L vs. 98.0 U/L, P = 0.029).Conclusions This is the first attempt to evaluate the association between ABCG2 genetic variants and the risk of ATDH. Based on the 1:4 matched case-control study, the polymorphism at rs2622605 in the ABCG2 gene may be associated with the susceptibility to ATDH in Chinese male patients.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Mapping the conformational landscape of the neutral network of RNA sequences that connect two functional distinctly different ribozymes

During evolution of an RNA world, the development of enzymatic function was essential. Such enzymatic function was linked to RNA sequences capable of adopting specific RNA folds that possess catalytic pockets to promote catalysis. Within this primordial RNA world, initially evolved self-replicating ribozymes presumably mutated to ribozymes with new functions. In 2000, Schultes and Bartel investigated such conversion from one ribozyme to a new ribozyme with distinctly different catalytic functions. Within a neutral network that linked these two prototype ribozymes, a single RNA chain could be identified that exhibited both enzymatic functions. This system serves as a paradigm for an evolutionary system that allows neutral drifts by stepwise mutation from one ribozyme into a different ribozyme without loss of intermittent function. Here, we investigated this complex functional diversification by analyzing several RNA sequences within this neutral network between two ribozymes with class III ligase activity and with self-cleavage reactivity. We utilized rapid RNA sample preparation for NMR spectroscopic studies together with SHAPE analysis and in-line probing to characterize secondary structure changes within the neutral network. Our investigations allowed delineation of the 2 nd structure space and by comparison with the previously determined catalytic function allowed correlation of the structure-function relation of ribozyme function in this neutral network.

Predicting the Structure of a Viroid : Structure, Structure Distribution, Consensus Structure, and Structure Drawing

Viroids are small non-coding RNAs that require a special sequence and structure to be replicated and transported by the host machinery. Many of these features can be predicted and later experimentally verified. Here, we will present workflows to predict viroid structures and draw the predicted structures in a pleasing and descriptive way using recently developed software.

Translation inhibitory elements from Hoxa3 and Hoxa11 mRNAs use uORFs for translation inhibition

During embryogenesis, Hox mRNA translation is tightly regulated by a sophisticated molecular mechanism that combines two RNA regulons located in their 5'UTR. First, an internal ribosome entry site (IRES) enables cap-independent translation. The second regulon is a translation inhibitory element or TIE, which ensures concomitant cap-dependent translation inhibition. In this study, we deciphered the molecular mechanisms of mouse Hoxa3 and Hoxa11 TIEs. Both TIEs possess an upstream open reading frame (uORF) that is critical to inhibit cap-dependent translation. However, the molecular mechanisms used are different. In Hoxa3 TIE, we identify an uORF which inhibits cap-dependent translation and we show the requirement of the non-canonical initiation factor eIF2D for this process. The mode of action of Hoxa11 TIE is different, it also contains an uORF but it is a minimal uORF formed by an uAUG followed immediately by a stop codon, namely a 'start-stop'. The 'start-stop' sequence is species-specific and in mice, is located upstream of a highly stable stem loop structure which stalls the 80S ribosome and thereby inhibits cap-dependent translation of Hoxa11 main ORF.

The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

Light-Dependent Translation Change of Arabidopsis psbA Correlates with RNA Structure Alterations at the Translation Initiation Region

mRNA secondary structure influences translation. Proteins that modulate the mRNA secondary structure around the translation initiation region may regulate translation in plastids. To test this hypothesis, we exposed Arabidopsis thaliana to high light, which induces translation of psbA mRNA encoding the D1 subunit of photosystem II. We assayed translation by ribosome profiling and applied two complementary methods to analyze in vivo RNA secondary structure: DMS-MaPseq and SHAPE-seq. We detected increased accessibility of the translation initiation region of psbA after high light treatment, likely contributing to the observed increase in translation by facilitating translation initiation. Furthermore, we identified the footprint of a putative regulatory protein in the 5' UTR of psbA at a position where occlusion of the nucleotide sequence would cause the structure of the translation initiation region to open up, thereby facilitating ribosome access. Moreover, we show that other plastid genes with weak Shine-Dalgarno sequences (SD) are likely to exhibit psbA-like regulation, while those with strong SDs do not. This supports the idea that changes in mRNA secondary structure might represent a general mechanism for translational regulation of psbA and other plastid genes.

Advances in the Bioinformatics Knowledge of mRNA Polyadenylation in Baculovirus Genes

Baculoviruses are a group of insect viruses with large circular dsDNA genomes exploited in numerous biotechnological applications, such as the biological control of agricultural pests, the expression of recombinant proteins or the gene delivery of therapeutic sequences in mammals, among others. Their genomes encode between 80 and 200 proteins, of which 38 are shared by all reported species. Thanks to multi-omic studies, there is remarkable information about the baculoviral proteome and the temporality in the virus gene expression. This allows some functional elements of the genome to be very well described, such as promoters and open reading frames. However, less information is available about the transcription termination signals and, consequently, there are still imprecisions about what are the limits of the transcriptional units present in the baculovirus genomes and how is the processing of the 3′ end of viral mRNA. Regarding to this, in this review we provide an update about the characteristics of DNA signals involved in this process and we contribute to their correct prediction through an exhaustive analysis that involves bibliography information, data mining, RNA structure and a comprehensive study of the core gene 3′ ends from 180 baculovirus genomes.

Expression of four new ferritins from grooved carpet shell clam Ruditapes decussatus challenged with Perkinsus olseni and metals (Cd, Cu and Zn)

Iron has a fundamental role in life and in its biochemical reactions but, when in excess, it can promote the formation of free radicals which can lead to cell death. Therefore, managing the levels of iron is essential to regulate the production of oxidative stress related to iron, and ferritins are one of the main protein families involved in this process. Ferritins are ≈480 kDa multimeric proteins composed by 24 subunits, each with 19-26 kDa, which can accumulate up to 4500 iron atoms. Besides their role in managing iron bioavailability, they have also developed a role in organism immunity and defence present throughout evolution. In this work, we identified and characterized, for the first time, four different ferritin subunits in the clam Ruditapes decussatus, a bivalve commercially and ecologically important along the south Atlantic coast and in the Mediterranean basin, which is a major target of the parasitic protozoa Perkinsus olseni, considered one of the main causes of high levels of clam mortality. Following phylogenetic annotation, the four ferritins subunits identified were subdivided into two cytosolic and two secreted forms. All four subunits maintain the canonical ferritin structure with four main helices α (A-D) and a small helix (E), but the secreted ferritins present an additional helix in their N-terminal region (F), located after the signal peptide and with possible antimicrobial properties. Additionally, we identified in ferritin 4 an extra helix α (G) located between helices B and C. These alpha helix domains revealed high degree of similarity with antimicrobial peptides associated with antibacterial and antifungal activities. Analysis of the expression of these subunits showed that ferritins 1 and 2 are ubiquitously expressed while ferritins 3 and 4 are present mainly in visceral mass. Ferritin 1 lacked a putative functional iron response element (IRE) and appeared to be under a tight regulation. Ferritins 2 and 3 showed a strong response to infection by parasite Perkinsus olseni in contrast to ferritin 4, whose main response was related to exposure to a combination of metals. The synergistic effect between metals and infection promoted a general upregulation of the four ferritins. In conclusion, our results suggest that ferritins, besides their function in iron and metals detoxification, may play a determinant role in clam immune response.

Assembly of infectious enteroviruses depends on multiple, conserved genomic RNA-coat protein contacts

Picornaviruses are important viral pathogens, but despite extensive study, the assembly process of their infectious virions is still incompletely understood, preventing the development of anti-viral strategies targeting this essential part of the life cycle. We report the identification, via RNA SELEX and bioinformatics, of multiple RNA sites across the genome of a typical enterovirus, enterovirus-E (EV-E), that each have affinity for the cognate viral capsid protein (CP) capsomer. Many of these sites are evolutionarily conserved across known EV-E variants, suggesting they play essential functional roles. Cryo-electron microscopy was used to reconstruct the EV-E particle at ~2.2 Å resolution, revealing extensive density for the genomic RNA. Relaxing the imposed symmetry within the reconstructed particles reveals multiple RNA-CP contacts, a first for any picornavirus. Conservative mutagenesis of the individual RNA-contacting amino acid side chains in EV-E, many of which are conserved across the enterovirus family including poliovirus, is lethal but does not interfere with replication or translation. Anti-EV-E and anti-poliovirus aptamers share sequence similarities with sites distributed across the poliovirus genome. These data are consistent with the hypothesis that these RNA-CP contacts are RNA Packaging Signals (PSs) that play vital roles in assembly and suggest that the RNA PSs are evolutionarily conserved between pathogens within the family, augmenting the current protein-only assembly paradigm for this family of viruses.

Assessing the accuracy of direct-coupling analysis for RNA contact prediction

Many noncoding RNAs are known to play a role in the cell directly linked to their structure. Structure prediction based on the sole sequence is, however, a challenging task. On the other hand, thanks to the low cost of sequencing technologies, a very large number of homologous sequences are becoming available for many RNA families. In the protein community, the idea of exploiting the covariance of mutations within a family to predict the protein structure using the direct-coupling-analysis (DCA) method has emerged in the last decade. The application of DCA to RNA systems has been limited so far. We here perform an assessment of the DCA method on 17 riboswitch families, comparing it with the commonly used mutual information analysis and with state-of-the-art R-scape covariance method. We also compare different flavors of DCA, including mean-field, pseudolikelihood, and a proposed stochastic procedure (Boltzmann learning) for solving exactly the DCA inverse problem. Boltzmann learning outperforms the other methods in predicting contacts observed in high-resolution crystal structures.

Secondary Structure of Chloroplast mRNAs In Vivo and In Vitro

mRNA secondary structure can influence gene expression, e.g., by influencing translation initiation. The probing of in vivo mRNA secondary structures is therefore necessary to understand what determines the efficiency and regulation of gene expression. Here, in vivo mRNA secondary structure was analyzed using dimethyl sulfate (DMS)-MaPseq and compared to in vitro-folded RNA. We used an approach to analyze specific, full-length transcripts. To test this approach, we chose low, medium, and high abundant mRNAs. We included both monocistronic and multicistronic transcripts. Because of the slightly alkaline pH of the chloroplast stroma, we could probe all four nucleotides with DMS. The structural information gained was evaluated using the known structure of the plastid 16S rRNA. This demonstrated that the results obtained for adenosines and cytidines were more reliable than for guanosines and uridines. The majority of mRNAs analyzed were less structured in vivo than in vitro. The in vivo secondary structure of the translation initiation region of most tested genes appears to be optimized for high translation efficiency.

Novel gRNA design pipeline to develop broad-spectrum CRISPR/Cas9 gRNAs for safe targeting of the HIV-1 quasispecies in patients

The CRISPR/Cas9 system has been proposed as a cure strategy for HIV. However, few published guide RNAs (gRNAs) are predicted to cleave the majority of HIV-1 viral quasispecies (vQS) observed within and among patients. We report the design of a novel pipeline to identify gRNAs that target HIV across a large number of infected individuals. Next generation sequencing (NGS) of LTRs from 269 HIV-1-infected samples in the Drexel CARES Cohort was used to select gRNAs with predicted broad-spectrum activity. In silico, D-LTR-P4-227913 (package of the top 4 gRNAs) accounted for all detectable genetic variation within the vQS of the 269 samples and the Los Alamos National Laboratory HIV database. In silico secondary structure analyses from NGS indicated extensive TAR stem-loop malformations predicted to inactivate proviral transcription, which was confirmed by reduced viral gene expression in TZM-bl or P4R5 cells. Similarly, a high sensitivity in vitro CRISPR/Cas9 cleavage assay showed that the top-ranked gRNA was the most effective at cleaving patient-derived HIV-1 LTRs from five patients. Furthermore, the D-LTR-P4-227913 was predicted to cleave a median of 96.1% of patient-derived sequences from other HIV subtypes. These results demonstrate that the gRNAs possess broad-spectrum cutting activity and could contribute to an HIV cure.

Selecting Fully-Modified XNA Aptamers Using Synthetic Genetics

This unit describes the application of "synthetic genetics," i.e., the replication of xeno nucleic acids (XNAs), artificial analogs of DNA and RNA bearing alternative backbone or sugar congeners, to the directed evolution of synthetic oligonucleotide ligands (XNA aptamers) specific for target proteins or nucleic acid motifs, using a cross-chemistry selective exponential enrichment (X-SELEX) approach. Protocols are described for synthesis of diverse-sequence XNA repertoires (typically 10¹⁴ molecules) using DNA templates, isolation and panning for functional XNA sequences using targets immobilized on solid phase or gel shift induced by target binding in solution, and XNA reverse transcription to allow cDNA amplification or sequencing. The method may be generally applied to select fully-modified XNA aptamers specific for a wide range of target molecules. © 2018 by John Wiley & Sons, Inc.

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.