1
|
Devi C, Ranjan P, Raj S, Das P. Computational exploration of protein structure dynamics and RNA structural consequences of PKD1 missense variants: implications in ADPKD pathogenesis. 3 Biotech 2024; 14:211. [PMID: 39188533 PMCID: PMC11344749 DOI: 10.1007/s13205-024-04057-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2024] [Accepted: 08/14/2024] [Indexed: 08/28/2024] Open
Abstract
We analyzed the impact of nine previously identified missense PKD1 variants from our studies, including c.6928G > A p.G2310R, c.8809G > A p.E2937K, c.2899 T > C p.W967R, c.6284A > G p.D2095G, c.6644G > A p.R2215Q, c.7810G > A p.D2604N, c.11249G > C p.R3750P, c.1001C > T p.T334M, and c.3101A > G p.N1034S on RNA structures and PC1 protein structure dynamics utilizing computational tools. RNA structure analysis was done using short RNA snippets of 41 nucleotides with the variant position at the 21st nucleotide, ensuring 20 bases on both sides. The secondary structures of these RNA snippets were predicted using RNAstructure. Structural changes of the mutants compared to the wild type were analyzed using the MutaRNA webserver. Molecular dynamics (MD) simulation of PC1 wild-type and mutant protein regions were performed using GROMACS 2018 (GROMOS96 54a7 force field). Findings revealed that five variants including c.8809G > A (p.E2937K), c.11249G > C (p.R3750P), c.3101A > G (p.N1034S), c.6928G > A (p.G2310R), c.6644G > A (p.R2215Q) exhibited major alterations in RNA structures and thereby their interactions with other proteins or RNAs affecting protein structure dynamics. While certain variants have minimal impact on RNA conformations, their observed alterations in MD simulations indicate impact on protein structure dynamics highlighting the importance of evaluating the functional consequences of genetic variants by considering both RNA and protein levels. The study also emphasizes that each missense variant exerts a unique impact on RNA stability, and protein structure dynamics, potentially contributing to the heterogeneous clinical manifestations and progression observed in Autosomal Dominant Polycystic Kidney Disease (ADPKD) patients offering a novel perspective in this direction. Thus, the utility of studying the structure dynamics through computational tools can help in prioritizing the variants for their functional implications, understanding the molecular mechanisms underlying variability in ADPKD presentation and developing targeted therapeutic interventions. Supplementary Information The online version contains supplementary material available at 10.1007/s13205-024-04057-9.
Collapse
Affiliation(s)
- Chandra Devi
- Centre for Genetic Disorders, Institute of Science, Banaras Hindu University, Varanasi, Uttar Pradesh 221005 India
| | - Prashant Ranjan
- Centre for Genetic Disorders, Institute of Science, Banaras Hindu University, Varanasi, Uttar Pradesh 221005 India
| | - Sonam Raj
- National Cancer Institute, National Institutes of Health, Bethesda, MD 20892 USA
| | - Parimal Das
- Centre for Genetic Disorders, Institute of Science, Banaras Hindu University, Varanasi, Uttar Pradesh 221005 India
| |
Collapse
|
2
|
Hertz LM, White EN, Kuznedelov K, Cheng L, Yu AM, Kakkaramadam R, Severinov K, Chen A, Lucks J. The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch. Nucleic Acids Res 2024; 52:4466-4482. [PMID: 38567721 PMCID: PMC11077080 DOI: 10.1093/nar/gkae231] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Revised: 02/17/2024] [Accepted: 03/20/2024] [Indexed: 04/16/2024] Open
Abstract
A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.
Collapse
Affiliation(s)
- Laura M Hertz
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL 60208, USA
| | - Elise N White
- Department of Chemistry and the RNA Institute, University at Albany, Albany, NY 12222, USA
| | | | - Luyi Cheng
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL 60208, USA
| | - Angela M Yu
- Department of Electrical and Computer Engineering, University of Washington, Seattle, WA 98195, USA
| | - Rivaan Kakkaramadam
- Department of Chemistry and the RNA Institute, University at Albany, Albany, NY 12222, USA
| | - Konstantin Severinov
- Waksman Institute of Microbiology, Rutgers University, Piscataway, NJ 08854, USA
| | - Alan Chen
- Department of Chemistry and the RNA Institute, University at Albany, Albany, NY 12222, USA
| | - Julius B Lucks
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL 60208, USA
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, IL 60208, USA
- Center for Synthetic Biology, Northwestern University, Evanston, IL 60208, USA
| |
Collapse
|
3
|
Zhukova M, Schedl P, Shidlovskii YV. The role of secondary structures in the functioning of 3' untranslated regions of mRNA: A review of functions of 3' UTRs' secondary structures and hypothetical involvement of secondary structures in cytoplasmic polyadenylation in Drosophila. Bioessays 2024; 46:e2300099. [PMID: 38161240 PMCID: PMC11337203 DOI: 10.1002/bies.202300099] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Revised: 12/11/2023] [Accepted: 12/12/2023] [Indexed: 01/03/2024]
Abstract
3' untranslated regions (3' UTRs) of mRNAs have many functions, including mRNA processing and transport, translational regulation, and mRNA degradation and stability. These different functions require cis-elements in 3' UTRs that can be either sequence motifs or RNA structures. Here we review the role of secondary structures in the functioning of 3' UTRs and discuss some of the trans-acting factors that interact with these secondary structures in eukaryotic organisms. We propose potential participation of 3'-UTR secondary structures in cytoplasmic polyadenylation in the model organism Drosophila melanogaster. Because the secondary structures of 3' UTRs are essential for post-transcriptional regulation of gene expression, their disruption leads to a wide range of disorders, including cancer and cardiovascular diseases. Trans-acting factors, such as STAU1 and nucleolin, which interact with 3'-UTR secondary structures of target transcripts, influence the pathogenesis of neurodegenerative diseases and tumor metastasis, suggesting that they are possible therapeutic targets.
Collapse
Affiliation(s)
- Mariya Zhukova
- Laboratory of Gene Expression Regulation in Development, Russian Academy of Sciences, Institute of Gene Biology, Moscow, Russia
| | - Paul Schedl
- Laboratory of Gene Expression Regulation in Development, Russian Academy of Sciences, Institute of Gene Biology, Moscow, Russia
- Department of Molecular Biology, Princeton University, Princeton, New Jersey, USA
| | - Yulii V. Shidlovskii
- Laboratory of Gene Expression Regulation in Development, Russian Academy of Sciences, Institute of Gene Biology, Moscow, Russia
- Department of Biology and General Genetics, Sechenov First Moscow State Medical University (Sechenov University), Moscow, Russia
| |
Collapse
|
4
|
Hertz LM, White EN, Kuznedelov K, Cheng L, Yu AM, Kakkaramadam R, Severinov K, Chen A, Lucks JB. The Effect of Pseudoknot Base Pairing on Cotranscriptional Structural Switching of the Fluoride Riboswitch. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.05.570056. [PMID: 38106011 PMCID: PMC10723315 DOI: 10.1101/2023.12.05.570056] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of E. coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37 °C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65 °C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.
Collapse
Affiliation(s)
- Laura M Hertz
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL 60208, USA
| | - Elise N White
- Department of Chemistry and the RNA Institute, University at Albany, Albany, NY 12222, USA
| | | | - Luyi Cheng
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL 60208, USA
| | - Angela M Yu
- Department of Electrical and Computer Engineering, University of Washington, Seattle, WA 98195, USA
| | - Rivaan Kakkaramadam
- Department of Chemistry and the RNA Institute, University at Albany, Albany, NY 12222, USA
| | - Konstantin Severinov
- Waksman Institute of Microbiology, Rutgers University, Piscataway, NJ 08854, USA
| | - Alan Chen
- Department of Chemistry and the RNA Institute, University at Albany, Albany, NY 12222, USA
| | - Julius B Lucks
- Interdisciplinary Biological Sciences Graduate Program, Northwestern University, Evanston, IL 60208, USA
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, IL 60208, USA
- Center for Synthetic Biology, Northwestern University, Evanston, IL 60208, USA
| |
Collapse
|
5
|
Behairy MY, Eid RA, Otifi HM, Mohammed HM, Alshehri MA, Asiri A, Aldehri M, Zaki MSA, Darwish KM, Elhady SS, El-Shaer NH, Eldeen MA. Unraveling Extremely Damaging IRAK4 Variants and Their Potential Implications for IRAK4 Inhibitor Efficacy. J Pers Med 2023; 13:1648. [PMID: 38138875 PMCID: PMC10744719 DOI: 10.3390/jpm13121648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 11/02/2023] [Accepted: 11/07/2023] [Indexed: 12/24/2023] Open
Abstract
Interleukin-1-receptor-associated kinase 4 (IRAK4) possesses a crucial function in the toll-like receptor (TLR) signaling pathway, and the dysfunction of this molecule could lead to various infectious and immune-related diseases in addition to cancers. IRAK4 genetic variants have been linked to various types of diseases. Therefore, we conducted a comprehensive analysis to recognize the missense variants with the most damaging impacts on IRAK4 with the employment of diverse bioinformatics tools to study single-nucleotide polymorphisms' effects on function, stability, secondary structures, and 3D structure. The residues' location on the protein domain and their conservation status were investigated as well. Moreover, docking tools along with structural biology were engaged in analyzing the SNPs' effects on one of the developed IRAK4 inhibitors. By analyzing IRAK4 gene SNPs, the analysis distinguished ten variants as the most detrimental missense variants. All variants were situated in highly conserved positions on an important protein domain. L318S and L318F mutations were linked to changes in IRAK4 secondary structures. Eight SNPs were revealed to have a decreasing effect on the stability of IRAK4 via both I-Mutant 2.0 and Mu-Pro tools, while Mu-Pro tool identified a decreasing effect for the G198E SNP. In addition, detrimental effects on the 3D structure of IRAK4 were also discovered for the selected variants. Molecular modeling studies highlighted the detrimental impact of these identified SNP mutant residues on the druggability of the IRAK4 ATP-binding site towards the known target inhibitor, HG-12-6, as compared to the native protein. The loss of important ligand residue-wise contacts, altered protein global flexibility, increased steric clashes, and even electronic penalties at the ligand-binding site interfaces were all suggested to be associated with SNP models for hampering the HG-12-6 affinity towards IRAK4 target protein. This given model lays the foundation for the better prediction of various disorders relevant to IRAK4 malfunction and sheds light on the impact of deleterious IRAK4 variants on IRAK4 inhibitor efficacy.
Collapse
Affiliation(s)
- Mohammed Y. Behairy
- Department of Microbiology and Immunology, Faculty of Pharmacy, University of Sadat City, Sadat City 32897, Egypt;
| | - Refaat A. Eid
- Department of Pathology, College of Medicine, King Khalid University, Abha P.O. Box 61421, Saudi Arabia; (R.A.E.); (H.M.O.)
| | - Hassan M. Otifi
- Department of Pathology, College of Medicine, King Khalid University, Abha P.O. Box 61421, Saudi Arabia; (R.A.E.); (H.M.O.)
| | - Heitham M. Mohammed
- Department of Anatomy, College of Medicine, King Khalid University, Abha P.O. Box 61421, Saudi Arabia; (H.M.M.); (M.A.); (M.S.A.Z.)
| | - Mohammed A. Alshehri
- Department of Child Health, College of Medicine, King Khalid University, Abha P.O. Box 62529, Saudi Arabia; (M.A.A.)
| | - Ashwag Asiri
- Department of Child Health, College of Medicine, King Khalid University, Abha P.O. Box 62529, Saudi Arabia; (M.A.A.)
| | - Majed Aldehri
- Department of Anatomy, College of Medicine, King Khalid University, Abha P.O. Box 61421, Saudi Arabia; (H.M.M.); (M.A.); (M.S.A.Z.)
| | - Mohamed Samir A. Zaki
- Department of Anatomy, College of Medicine, King Khalid University, Abha P.O. Box 61421, Saudi Arabia; (H.M.M.); (M.A.); (M.S.A.Z.)
| | - Khaled M. Darwish
- Department of Medicinal Chemistry, Faculty of Pharmacy, Suez Canal University, Ismailia 41522, Egypt;
| | - Sameh S. Elhady
- Department of Natural Products, Faculty of Pharmacy, King Abdulaziz University, Jeddah 21589, Saudi Arabia;
| | - Nahla H. El-Shaer
- Department of Zoology, Faculty of Science, Zagazig University, Zagazig 44511, Egypt;
| | - Muhammad Alaa Eldeen
- Department of Zoology, Faculty of Science, Zagazig University, Zagazig 44511, Egypt;
| |
Collapse
|
6
|
Zhou S, Van Bortle K. The Pol III transcriptome: Basic features, recurrent patterns, and emerging roles in cancer. WILEY INTERDISCIPLINARY REVIEWS. RNA 2023; 14:e1782. [PMID: 36754845 PMCID: PMC10498592 DOI: 10.1002/wrna.1782] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/06/2022] [Revised: 01/13/2023] [Accepted: 01/18/2023] [Indexed: 02/10/2023]
Abstract
The RNA polymerase III (Pol III) transcriptome is universally comprised of short, highly structured noncoding RNA (ncRNA). Through RNA-protein interactions, the Pol III transcriptome actuates functional activities ranging from nuclear gene regulation (7SK), splicing (U6, U6atac), and RNA maturation and stability (RMRP, RPPH1, Y RNA), to cytoplasmic protein targeting (7SL) and translation (tRNA, 5S rRNA). In higher eukaryotes, the Pol III transcriptome has expanded to include additional, recently evolved ncRNA species that effectively broaden the footprint of Pol III transcription to additional cellular activities. Newly evolved ncRNAs function as riboregulators of autophagy (vault), immune signaling cascades (nc886), and translation (Alu, BC200, snaR). Notably, upregulation of Pol III transcription is frequently observed in cancer, and multiple ncRNA species are linked to both cancer progression and poor survival outcomes among cancer patients. In this review, we outline the basic features and functions of the Pol III transcriptome, and the evidence for dysregulation and dysfunction for each ncRNA in cancer. When taken together, recurrent patterns emerge, ranging from shared functional motifs that include molecular scaffolding and protein sequestration, overlapping protein interactions, and immunostimulatory activities, to the biogenesis of analogous small RNA fragments and noncanonical miRNAs, augmenting the function of the Pol III transcriptome and further broadening its role in cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Processing of Small RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.
Collapse
Affiliation(s)
- Sihang Zhou
- Department of Cell and Developmental Biology, University of Illinois Urbana-Champaign, Urbana, Illinois, USA
| | - Kevin Van Bortle
- Department of Cell and Developmental Biology, University of Illinois Urbana-Champaign, Urbana, Illinois, USA
- Cancer Center at Illinois, University of Illinois Urbana-Champaign, Urbana, Illinois, USA
| |
Collapse
|
7
|
How does precursor RNA structure influence RNA processing and gene expression? Biosci Rep 2023; 43:232489. [PMID: 36689327 PMCID: PMC9977717 DOI: 10.1042/bsr20220149] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Revised: 01/17/2023] [Accepted: 01/23/2023] [Indexed: 01/24/2023] Open
Abstract
RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.
Collapse
|
8
|
Behairy MY, Soltan MA, Eldeen MA, Abdulhakim JA, Alnoman MM, Abdel-Daim MM, Otifi H, Al-Qahtani SM, Zaki MSA, Alsharif G, Albogami S, Jafri I, Fayad E, Darwish KM, Elhady SS, Eid RA. HBD-2 variants and SARS-CoV-2: New insights into inter-individual susceptibility. Front Immunol 2022; 13:1008463. [PMID: 36569842 PMCID: PMC9780532 DOI: 10.3389/fimmu.2022.1008463] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2022] [Accepted: 11/17/2022] [Indexed: 12/13/2022] Open
Abstract
Background A deep understanding of the causes of liability to SARS-CoV-2 is essential to develop new diagnostic tests and therapeutics against this serious virus in order to overcome this pandemic completely. In the light of the discovered role of antimicrobial peptides [such as human b-defensin-2 (hBD-2) and cathelicidin LL-37] in the defense against SARS-CoV-2, it became important to identify the damaging missense mutations in the genes of these molecules and study their role in the pathogenesis of COVID-19. Methods We conducted a comprehensive analysis with multiple in silico approaches to identify the damaging missense SNPs for hBD-2 and LL-37; moreover, we applied docking methods and molecular dynamics analysis to study the impact of the filtered mutations. Results The comprehensive analysis reveals the presence of three damaging SNPs in hBD-2; these SNPs were predicted to decrease the stability of hBD-2 with a damaging impact on hBD-2 structure as well. G51D and C53G mutations were located in highly conserved positions and were associated with differences in the secondary structures of hBD-2. Docking-coupled molecular dynamics simulation analysis revealed compromised binding affinity for hBD-2 SNPs towards the SARS-CoV-2 spike domain. Different protein-protein binding profiles for hBD-2 SNPs, in relation to their native form, were guided through residue-wise levels and differential adopted conformation/orientation. Conclusions The presented model paves the way for identifying patients prone to COVID-19 in a way that would guide the personalization of both the diagnostic and management protocols for this serious disease.
Collapse
Affiliation(s)
- Mohammed Y. Behairy
- Department of Microbiology and Immunology, Faculty of Pharmacy, University of Sadat City, Sadat City, Egypt,*Correspondence: Mohamed A Soltan, ; Mohammed Y. Behairy,
| | - Mohamed A. Soltan
- Department of Microbiology and immunology, Faculty of Pharmacy, Sinai University – Kantara Branch, Ismailia, Egypt,*Correspondence: Mohamed A Soltan, ; Mohammed Y. Behairy,
| | - Muhammad Alaa Eldeen
- Cell Biology, Histology & Genetics Division, Biology Department, Faculty of Science, Zagazig University, Zagazig, Egypt
| | - Jawaher A. Abdulhakim
- Medical Laboratory Department, College of Applied Medical Sciences, Taibah University, Yanbu, Saudi Arabia
| | - Maryam M. Alnoman
- Biology Department, Faculty of Science, Taibah University, Yanbu, Saudi Arabia
| | - Mohamed M. Abdel-Daim
- Department of Pharmaceutical Sciences, Pharmacy Program, Batterjee Medical College, Jeddah, Saudi Arabia,Pharmacology Department, Faculty of Veterinary Medicine, Suez Canal University, Ismailia, Egypt
| | - Hassan Otifi
- Pathology Department, College of Medicine, King Khalid University, Abha, Saudi Arabia
| | - Saleh M. Al-Qahtani
- Department of Child Health, College of Medicine, King Khalid University, Abha, Saudi Arabia
| | - Mohamed Samir A. Zaki
- Anatomy Department, College of Medicine, King Khalid University, Abha, Saudi Arabia,Department of Histology and Cell Biology, College of Medicine, Zagazig University, Zagazig, Egypt
| | - Ghadi Alsharif
- College of Clinical Laboratory Sciences, King Saud bin Abdulaziz University for Health Sciences, Jeddah, Saudi Arabia
| | - Sarah Albogami
- Department of Biotechnology, College of Sciences, Taif University, Taif, Saudi Arabia
| | - Ibrahim Jafri
- Department of Biotechnology, College of Sciences, Taif University, Taif, Saudi Arabia
| | - Eman Fayad
- Department of Biotechnology, College of Sciences, Taif University, Taif, Saudi Arabia
| | - Khaled M. Darwish
- Department of Medicinal Chemistry, Faculty of Pharmacy, Suez Canal University, Ismailia, Egypt
| | - Sameh S. Elhady
- Department of Natural Products, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Refaat A. Eid
- Pathology Department, College of Medicine, King Khalid University, Abha, Saudi Arabia
| |
Collapse
|
9
|
Waldern JM, Kumar J, Laederach A. Disease-associated human genetic variation through the lens of precursor and mature RNA structure. Hum Genet 2022; 141:1659-1672. [PMID: 34741198 PMCID: PMC9072596 DOI: 10.1007/s00439-021-02395-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 10/26/2021] [Indexed: 12/14/2022]
Abstract
Disease-associated variants (DAVs) are commonly considered either through a genomic lens that describes variant function at the DNA level, or at the protein function level if the variant is translated. Although the genomic and proteomic effects of variation are well-characterized, genetic variants disrupting post-transcriptional regulation is another mechanism of disease that remains understudied. Specific RNA sequence motifs mediate post-transcriptional regulation both in the nucleus and cytoplasm of eukaryotic cells, often by binding to RNA-binding proteins or other RNAs. However, many DAVs map far from these motifs, which suggests deeper layers of post-transcriptional mechanistic control. Here, we consider a transcriptomic framework to outline the importance of post-transcriptional regulation as a mechanism of disease-causing single-nucleotide variation in the human genome. We first describe the composition of the human transcriptome and the importance of abundant yet overlooked components such as introns and untranslated regions (UTRs) of messenger RNAs (mRNAs). We present an analysis of Human Gene Mutation Database variants mapping to mRNAs and examine the distribution of causative disease-associated variation across the transcriptome. Although our analysis confirms the importance of post-transcriptional regulatory motifs, a majority of DAVs do not directly map to known regulatory motifs. Therefore, we review evidence that regions outside these well-characterized motifs can regulate function by RNA structure-mediated mechanisms in all four elements of an mRNA: exons, introns, 5' and 3' UTRs. To this end, we review published examples of riboSNitches, which are single-nucleotide variants that result in a change in RNA structure that is causative of the disease phenotype. In this review, we present the current state of knowledge of how DAVs act at the transcriptome level, both through altering post-transcriptional regulatory motifs and by the effects of RNA structure.
Collapse
Affiliation(s)
- Justin M Waldern
- Department of Biology, University of North Carolina, Chapel Hill, NC, 27599, USA
| | - Jayashree Kumar
- Department of Biology, University of North Carolina, Chapel Hill, NC, 27599, USA
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA
| | - Alain Laederach
- Department of Biology, University of North Carolina, Chapel Hill, NC, 27599, USA.
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA.
| |
Collapse
|
10
|
Sun W, Ding L, Zhang H. The Potential Role of RNA Structure in Crop Molecular Breeding. FRONTIERS IN PLANT SCIENCE 2022; 13:868771. [PMID: 35586218 PMCID: PMC9108716 DOI: 10.3389/fpls.2022.868771] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/03/2022] [Accepted: 03/21/2022] [Indexed: 06/15/2023]
Abstract
The continually growing human population creates a concomitantly increasing demand for nutritious crops with high yields. Advances in high throughput sequencing technologies have revealed the genetic architecture of major crops. This includes extensive information enabling comprehensive genetic markers for breeding selection, new gene discoveries, and novel gene regulatory strategies for crop editing. RNA structure is an important type of genetic feature, essential for post-transcriptional regulation of gene expression. Here, we summarize recent advances in genome-wide RNA structure studies in crops and review the associated RNA structure-mediated regulation of gene expression. We also discuss the functional importance of those single nucleotide variations that induce large RNA structure disparities. Lastly, we discuss the potential role of RNA structure in crop molecular breeding.
Collapse
|
11
|
Zeng Z, Bromberg Y. Inferring Potential Cancer Driving Synonymous Variants. Genes (Basel) 2022; 13:778. [PMID: 35627162 PMCID: PMC9140830 DOI: 10.3390/genes13050778] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Revised: 04/25/2022] [Accepted: 04/26/2022] [Indexed: 02/01/2023] Open
Abstract
Synonymous single nucleotide variants (sSNVs) are often considered functionally silent, but a few cases of cancer-causing sSNVs have been reported. From available databases, we collected four categories of sSNVs: germline, somatic in normal tissues, somatic in cancerous tissues, and putative cancer drivers. We found that screening sSNVs for recurrence among patients, conservation of the affected genomic position, and synVep prediction (synVep is a machine learning-based sSNV effect predictor) recovers cancer driver variants (termed proposed drivers) and previously unknown putative cancer genes. Of the 2.9 million somatic sSNVs found in the COSMIC database, we identified 2111 proposed cancer driver sSNVs. Of these, 326 sSNVs could be further tagged for possible RNA splicing effects, RNA structural changes, and affected RBP motifs. This list of proposed cancer driver sSNVs provides computational guidance in prioritizing the experimental evaluation of synonymous mutations found in cancers. Furthermore, our list of novel potential cancer genes, galvanized by synonymous mutations, may highlight yet unexplored cancer mechanisms.
Collapse
Affiliation(s)
- Zishuo Zeng
- Department of Biochemistry and Microbiology, Rutgers University, New Brunswick, NJ 08873, USA
| | - Yana Bromberg
- Department of Biochemistry and Microbiology, Rutgers University, New Brunswick, NJ 08873, USA
- Department of Genetics, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
12
|
Ferrero-Serrano Á, Sylvia MM, Forstmeier PC, Olson AJ, Ware D, Bevilacqua PC, Assmann SM. Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis. Genome Biol 2022; 23:101. [PMID: 35440059 PMCID: PMC9017077 DOI: 10.1186/s13059-022-02656-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Accepted: 03/20/2022] [Indexed: 11/23/2022] Open
Abstract
Background Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.” Results We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation. Conclusion We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation. Supplementary Information The online version contains supplementary material available at 10.1186/s13059-022-02656-4.
Collapse
Affiliation(s)
- Ángel Ferrero-Serrano
- Department of Biology, Pennsylvania State University, University Park, State College, PA, 16802, USA.
| | - Megan M Sylvia
- Department of Biology, Pennsylvania State University, University Park, State College, PA, 16802, USA
| | - Peter C Forstmeier
- Department of Biochemistry, Microbiology, and Molecular Biology, Pennsylvania State University, University Park, State College, PA, 16802, USA
| | - Andrew J Olson
- Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, 11724, USA
| | - Doreen Ware
- Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, 11724, USA.,USDA ARS NAA Robert W. Holley Center for Agriculture and Health, Ithaca, NY, 14853, USA
| | - Philip C Bevilacqua
- Department of Biochemistry, Microbiology, and Molecular Biology, Pennsylvania State University, University Park, State College, PA, 16802, USA.,Department of Chemistry, Pennsylvania State University, University Park, State College, PA, 16802, USA.,Center for RNA Molecular Biology, Pennsylvania State University, University Park, State College, PA, 16802, USA
| | - Sarah M Assmann
- Department of Biology, Pennsylvania State University, University Park, State College, PA, 16802, USA. .,Center for RNA Molecular Biology, Pennsylvania State University, University Park, State College, PA, 16802, USA.
| |
Collapse
|
13
|
Yang X, Yu H, Sun W, Ding L, Li J, Cheema J, Ramirez-Gonzalez R, Zhao X, Martín AC, Lu F, Liu B, Uauy C, Ding Y, Zhang H. Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry. Genome Biol 2021; 22:326. [PMID: 34847934 PMCID: PMC8638558 DOI: 10.1186/s13059-021-02549-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Accepted: 11/19/2021] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry. RESULTS Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat. CONCLUSIONS We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.
Collapse
Affiliation(s)
- Xiaofei Yang
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China
- Department of Cell and Developmental Biology, John Innes Centre, Norwich Research Park, Norwich, NR4 7UH, UK
| | - Haopeng Yu
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China
- Department of Cell and Developmental Biology, John Innes Centre, Norwich Research Park, Norwich, NR4 7UH, UK
| | - Wenqing Sun
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China
| | - Ling Ding
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China
| | - Ji Li
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China
| | - Jitender Cheema
- Department of Cell and Developmental Biology, John Innes Centre, Norwich Research Park, Norwich, NR4 7UH, UK
| | | | - Xuebo Zhao
- State Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, The Innovative Academy of Seed Design, Chinese Academy of Sciences, Beijing, China
- University of Chinese Academy of Sciences, Beijing, China
- CAS-JIC Centre of Excellence for Plant and Microbial Science (CEPAMS), Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, China
| | - Azahara C Martín
- Department of Crop Genetics, John Innes Centre, Norwich Research Park, Norwich, NR4 7UH, UK
| | - Fei Lu
- State Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, The Innovative Academy of Seed Design, Chinese Academy of Sciences, Beijing, China
- University of Chinese Academy of Sciences, Beijing, China
- CAS-JIC Centre of Excellence for Plant and Microbial Science (CEPAMS), Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, China
| | - Bao Liu
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China
| | - Cristobal Uauy
- Department of Crop Genetics, John Innes Centre, Norwich Research Park, Norwich, NR4 7UH, UK
| | - Yiliang Ding
- Department of Cell and Developmental Biology, John Innes Centre, Norwich Research Park, Norwich, NR4 7UH, UK.
| | - Huakun Zhang
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China.
| |
Collapse
|
14
|
Lackey L, Coria A, Ghosh AJ, Grayeski P, Hatfield A, Shankar V, Platig J, Xu Z, Ramos SBV, Silverman EK, Ortega VE, Cho MH, Hersh CP, Hobbs BD, Castaldi P, Laederach A. Alternative poly-adenylation modulates α1-antitrypsin expression in chronic obstructive pulmonary disease. PLoS Genet 2021; 17:e1009912. [PMID: 34784346 PMCID: PMC8631626 DOI: 10.1371/journal.pgen.1009912] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Revised: 11/30/2021] [Accepted: 10/25/2021] [Indexed: 01/07/2023] Open
Abstract
α1-anti-trypsin (A1AT), encoded by SERPINA1, is a neutrophil elastase inhibitor that controls the inflammatory response in the lung. Severe A1AT deficiency increases risk for Chronic Obstructive Pulmonary Disease (COPD), however, the role of A1AT in COPD in non-deficient individuals is not well known. We identify a 2.1-fold increase (p = 2.5x10-6) in the use of a distal poly-adenylation site in primary lung tissue RNA-seq in 82 COPD cases when compared to 64 controls and replicate this in an independent study of 376 COPD and 267 controls. This alternative polyadenylation event involves two sites, a proximal and distal site, 61 and 1683 nucleotides downstream of the A1AT stop codon. To characterize this event, we measured the distal ratio in human primary tissue short read RNA-seq data and corroborated our results with long read RNA-seq data. Integrating these results with 3' end RNA-seq and nanoluciferase reporter assay experiments we show that use of the distal site yields mRNA transcripts with over 50-fold decreased translation efficiency and A1AT expression. We identified seven RNA binding proteins using enhanced CrossLinking and ImmunoPrecipitation precipitation (eCLIP) with one or more binding sites in the SERPINA1 3' UTR. We combined these data with measurements of the distal ratio in shRNA knockdown experiments, nuclear and cytoplasmic fractionation, and chemical RNA structure probing. We identify Quaking Homolog (QKI) as a modulator of SERPINA1 mRNA translation and confirm the role of QKI in SERPINA1 translation with luciferase reporter assays. Analysis of single-cell RNA-seq showed differences in the distribution of the SERPINA1 distal ratio among hepatocytes, macrophages, αβ-Tcells and plasma cells in the liver. Alveolar Type 1,2, dendritic cells and macrophages also vary in their distal ratio in the lung. Our work reveals a complex post-transcriptional mechanism that regulates alternative polyadenylation and A1AT expression in COPD.
Collapse
Affiliation(s)
- Lela Lackey
- Department of Genetics and Biochemistry, Center for Human Genetics, Clemson University, Greenwood, South Carolina, United States of America
| | - Aaztli Coria
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina, United States of America
- Department of Biochemistry and Biophysics, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Auyon J. Ghosh
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
- Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Phil Grayeski
- Curriculum in Genetics and Molecular Biology, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Abigail Hatfield
- Department of Genetics and Biochemistry, Center for Human Genetics, Clemson University, Greenwood, South Carolina, United States of America
| | - Vijay Shankar
- Department of Genetics and Biochemistry, Center for Human Genetics, Clemson University, Greenwood, South Carolina, United States of America
| | - John Platig
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Zhonghui Xu
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Silvia B. V. Ramos
- Department of Biochemistry and Biophysics, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Edwin K. Silverman
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
- Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Victor E. Ortega
- Department of Internal Medicine, Division of Respiratory Medicine, Center for Individualized Medicine, Mayo Clinic, Scottsdale, Arizona, United States of America
| | - Michael H. Cho
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
- Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Craig P. Hersh
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
- Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Brian D. Hobbs
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
- Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Peter Castaldi
- Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
- Division of Internal Medicine and Primary Care, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Alain Laederach
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina, United States of America
| |
Collapse
|
15
|
Martens L, Rühle F, Witten A, Meder B, Katus HA, Arbustini E, Hasenfuß G, Sinner MF, Kääb S, Pankuweit S, Angermann C, Bornberg-Bauer E, Stoll M. A genetic variant alters the secondary structure of the lncRNA H19 and is associated with dilated cardiomyopathy. RNA Biol 2021; 18:409-415. [PMID: 34313541 DOI: 10.1080/15476286.2021.1952756] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022] Open
Abstract
lncRNAs are at the core of many regulatory processes and have also been recognized to be involved in various complex diseases. They affect gene regulation through direct interactions with RNA, DNA or proteins. Accordingly, lncRNA structure is likely to be essential for their regulatory function. Point mutations, which manifest as SNPs (single nucleotide polymorphisms) in genome screens, can substantially alter their function and, subsequently, the expression of their downstream regulated genes. To test the effect of SNPs on structure, we investigated lncRNAs associated with dilated cardiomyopathy. Among 322 human candidate lncRNAs, we demonstrate first the significant association of an SNP located in lncRNA H19 using data from 1084 diseased and 751 control patients. H19 is generally highly expressed in the heart, with a complex expression pattern during heart development. Next, we used MFE (minimum free energy) folding to demonstrate a significant refolding in the secondary structure of this 861 nt long lncRNA. Since MFE folding may overlook the importance of sub-optimal structures, we showed that this refolding also manifests in the overall Boltzmann structure ensemble. There, the composition of structures is tremendously affected in their thermodynamic probabilities through the genetic variant. Finally, we confirmed these results experimentally, using SHAPE-Seq, corroborating that SNPs affecting such structures may explain hidden genetic variance not accounted for through genome wide association studies. Our results suggest that structural changes in lncRNAs, and lncRNA H19 in particular, affect regulatory processes and represent optimal targets for further in-depth studies probing their molecular interactions.
Collapse
Affiliation(s)
- Leonie Martens
- Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
| | - Frank Rühle
- Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
| | - Anika Witten
- Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
| | - Benjamin Meder
- German Centre for Cardiovascular Research (DZHK), Partner Site Heidelberg/Mannheim, Heidelberg, Germany.,Department of Cardiology, Heidelberg University, Heidelberg, Germany.,Genome Technology Center Stanford, Department of Genetics, Stanford University, Stanford, United States
| | - Hugo A Katus
- German Centre for Cardiovascular Research (DZHK), Partner Site Heidelberg/Mannheim, Heidelberg, Germany.,Department of Cardiology, Heidelberg University, Heidelberg, Germany
| | - Eloisa Arbustini
- Centre for Inherited Cardiovascular Diseases, IRCCS Foundation, University Hospital Policlinico San Matteo, Pavia, Italy
| | - Gerd Hasenfuß
- Department of Cardiology and Pneumology, University Medical Center Göttingen, Göttingen, Germany.,German Center for Cardiovascular Research (DZHK), Partner Site, Göttingen, Germany
| | - Moritz F Sinner
- Department of Cardiology, University Hospital, LMU Munich, Munich, Germany.,German Centre for Cardiovascular Research (DZHK), Partner Site: Munich Heart Alliance, Munich, Germany
| | - Stefan Kääb
- Department of Cardiology, University Hospital, LMU Munich, Munich, Germany.,German Centre for Cardiovascular Research (DZHK), Partner Site: Munich Heart Alliance, Munich, Germany
| | - Sabine Pankuweit
- Department of Cardiology, University Hospital Giessen and Marburg, Marburg, Germany
| | - Christiane Angermann
- Comprehensive Heart Failure Center, University Hospital and University of Würzburg, Würzburg, Germany
| | - Erich Bornberg-Bauer
- Institute for Evolution and Biodiversity, University of Münster, Münster, Germany
| | - Monika Stoll
- Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany.,Department of Biochemistry, Genetic Epidemiology and Statistical Genetics, CARIM School for Cardiovascular Diseases, Maastricht Center for Systems Biology (MaCSBio), Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
16
|
Miladi M, Raden M, Diederichs S, Backofen R. MutaRNA: analysis and visualization of mutation-induced changes in RNA structure. Nucleic Acids Res 2020; 48:W287-W291. [PMID: 32392303 PMCID: PMC7319544 DOI: 10.1093/nar/gkaa331] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2020] [Revised: 04/14/2020] [Accepted: 04/22/2020] [Indexed: 12/25/2022] Open
Abstract
RNA molecules fold into complex structures as a result of intramolecular interactions between their nucleotides. The function of many non-coding RNAs and some cis-regulatory elements of messenger RNAs highly depends on their fold. Single-nucleotide variants (SNVs) and other types of mutations can disrupt the native function of an RNA element by altering its base pairing pattern. Identifying the effect of a mutation on an RNA’s structure is, therefore, a crucial step in evaluating the impact of mutations on the post-transcriptional regulation and function of RNAs within the cell. Even though a single nucleotide variation can have striking impacts on the structure formation, interpreting and comparing the impact usually needs expertise and meticulous efforts. Here, we present MutaRNA, a web server for visualization and interpretation of mutation-induced changes on the RNA structure in an intuitive and integrative fashion. To this end, probabilities of base pairing and position-wise unpaired probabilities of wildtype and mutated RNA sequences are computed and compared. Differential heatmap-like dot plot representations in combination with circular plots and arc diagrams help to identify local structure abberations, which are otherwise hidden in standard outputs. Eventually, MutaRNA provides a comprehensive and comparative overview of the mutation-induced changes in base pairing potentials and accessibility. The MutaRNA web server is freely available at http://rna.informatik.uni-freiburg.de/MutaRNA.
Collapse
Affiliation(s)
- Milad Miladi
- Bioinformatics Group, Department of Computer Science, University of Freiburg, Georges-Koehler-Allee 106, 79110 Freiburg, Germany
| | - Martin Raden
- Bioinformatics Group, Department of Computer Science, University of Freiburg, Georges-Koehler-Allee 106, 79110 Freiburg, Germany
| | - Sven Diederichs
- Division of Cancer Research, Department of Thoracic Surgery, Faculty of Medicine, German Cancer Consortium (DKTK), University of Freiburg, 79085 Freiburg, Germany.,Division of RNA Biology and Cancer, German Cancer Research Center (DKFZ) and National Center for Tumor Diseases (NCT), 69120 Heidelberg, Germany
| | - Rolf Backofen
- Bioinformatics Group, Department of Computer Science, University of Freiburg, Georges-Koehler-Allee 106, 79110 Freiburg, Germany.,Signalling Research Centres BIOSS and CIBSS, University of Freiburg, Schaenzlestr. 18, 79104 Freiburg, Germany
| |
Collapse
|
17
|
Lin J, Chen Y, Zhang Y, Ouyang Z. Identification and analysis of RNA structural disruptions induced by single nucleotide variants using Riprap and RiboSNitchDB. NAR Genom Bioinform 2020; 2:lqaa057. [PMID: 33575608 PMCID: PMC7671322 DOI: 10.1093/nargab/lqaa057] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2019] [Revised: 06/22/2020] [Accepted: 08/11/2020] [Indexed: 11/27/2022] Open
Abstract
RNA conformational alteration has significant impacts on cellular processes and phenotypic variations. An emerging genetic factor of RNA conformational alteration is a new class of single nucleotide variant (SNV) named riboSNitch. RiboSNitches have been demonstrated to be involved in many genetic diseases. However, identifying riboSNitches is notably difficult as the signals of RNA structural disruption are often subtle. Here, we introduce a novel computational framework–RIboSNitch Predictor based on Robust Analysis of Pairing probabilities (Riprap). Riprap identifies structurally disrupted regions around any given SNVs based on robust analysis of local structural configurations between wild-type and mutant RNA sequences. Compared to previous approaches, Riprap shows higher accuracy when assessed on hundreds of known riboSNitches captured by various experimental RNA structure probing methods including the parallel analysis of RNA structure (PARS) and the selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE). Further, Riprap detects the experimentally validated riboSNitch that regulates human catechol-O-methyltransferase haplotypes and outputs structurally disrupted regions precisely at base resolution. Riprap provides a new approach to interpreting disease-related genetic variants. In addition, we construct a database (RiboSNitchDB) that includes the annotation and visualization of all presented riboSNitches in this study as well as 24 629 predicted riboSNitches from human expression quantitative trait loci.
Collapse
Affiliation(s)
- Jianan Lin
- Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, MA 01003, USA
| | - Yang Chen
- Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, MA 01003, USA
| | - Yuping Zhang
- Department of Statistics, University of Connecticut, Storrs, CT 06269, USA
| | - Zhengqing Ouyang
- Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, MA 01003, USA
| |
Collapse
|
18
|
Matyjasik MM, Batey RT. Structural basis for 2'-deoxyguanosine recognition by the 2'-dG-II class of riboswitches. Nucleic Acids Res 2020; 47:10931-10941. [PMID: 31598729 PMCID: PMC6847200 DOI: 10.1093/nar/gkz839] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Revised: 09/12/2019] [Accepted: 10/05/2019] [Indexed: 12/25/2022] Open
Abstract
A recent bioinformatic analysis of well-characterized classes of riboswitches uncovered subgroups unable to bind to the regulatory molecule of the parental class. Within the guanine/adenine class, seven groups of RNAs were identified that deviate from the consensus sequence at one or more of three positions directly involved purine nucleobase recognition, one of which was validated as a second class of 2'-deoxyguanosine riboswitch (called 2'-dG-II). To understand how 2'-dG-II riboswitches recognize their cognate ligand and how they differ from a previously identified class of 2'-deoxyguanosine binding riboswitches, we have solved the crystal structure of a 2'-dG-II aptamer domain bound to 2'-deoxyguanosine. This structure reveals a global architecture similar to other members of the purine riboswitch family, but contains key differences within the ligand binding core. Defining the 2'-dG-II riboswitches is a two-nucleotide insertion in the three-way junction that promotes novel base-base interactions. Unlike 2'-dG-I riboswitches, the 2'-dG-II class only requires local changes to the ligand binding pocket of the guanine/adenine class to achieve a change in ligand preference. Notably, members of the 2'-dG-II family have variable ability to discriminate between 2'-deoxyguanosine and riboguanosine, suggesting that a subset of 2'-dG-II riboswitches may bind either molecule to regulate gene expression.
Collapse
Affiliation(s)
- Michal M Matyjasik
- Department of Biochemistry, University of Colorado at Boulder, Campus Box 596, Boulder, CO 80309-0596, USA
| | - Robert T Batey
- Department of Biochemistry, University of Colorado at Boulder, Campus Box 596, Boulder, CO 80309-0596, USA
| |
Collapse
|
19
|
Developing parallel ant colonies filtered by deep learned constrains for predicting RNA secondary structure with pseudo-knots. Neurocomputing 2020. [DOI: 10.1016/j.neucom.2019.12.041] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
20
|
Miladi M, Sokhoyan E, Houwaart T, Heyne S, Costa F, Grüning B, Backofen R. GraphClust2: Annotation and discovery of structured RNAs with scalable and accessible integrative clustering. Gigascience 2019; 8:giz150. [PMID: 31808801 PMCID: PMC6897289 DOI: 10.1093/gigascience/giz150] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2019] [Revised: 08/23/2019] [Accepted: 11/20/2019] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND RNA plays essential roles in all known forms of life. Clustering RNA sequences with common sequence and structure is an essential step towards studying RNA function. With the advent of high-throughput sequencing techniques, experimental and genomic data are expanding to complement the predictive methods. However, the existing methods do not effectively utilize and cope with the immense amount of data becoming available. RESULTS Hundreds of thousands of non-coding RNAs have been detected; however, their annotation is lagging behind. Here we present GraphClust2, a comprehensive approach for scalable clustering of RNAs based on sequence and structural similarities. GraphClust2 bridges the gap between high-throughput sequencing and structural RNA analysis and provides an integrative solution by incorporating diverse experimental and genomic data in an accessible manner via the Galaxy framework. GraphClust2 can efficiently cluster and annotate large datasets of RNAs and supports structure-probing data. We demonstrate that the annotation performance of clustering functional RNAs can be considerably improved. Furthermore, an off-the-shelf procedure is introduced for identifying locally conserved structure candidates in long RNAs. We suggest the presence and the sparseness of phylogenetically conserved local structures for a collection of long non-coding RNAs. CONCLUSIONS By clustering data from 2 cross-linking immunoprecipitation experiments, we demonstrate the benefits of GraphClust2 for motif discovery under the presence of biological and methodological biases. Finally, we uncover prominent targets of double-stranded RNA binding protein Roquin-1, such as BCOR's 3' untranslated region that contains multiple binding stem-loops that are evolutionary conserved.
Collapse
Affiliation(s)
- Milad Miladi
- Bioinformatics Group, Department of Computer Science, University of Freiburg, Georges-Koehler-Allee 106, 79110 Freiburg, Germany
| | - Eteri Sokhoyan
- Bioinformatics Group, Department of Computer Science, University of Freiburg, Georges-Koehler-Allee 106, 79110 Freiburg, Germany
| | - Torsten Houwaart
- Institute of Medical Microbiology and Hospital Hygiene, University of Dusseldorf, Universitaetsstr. 1, 40225 Dusseldorf, Germany
| | - Steffen Heyne
- Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Stuebeweg 51, 79108 Freiburg, Germany
| | - Fabrizio Costa
- Department of Computer Science, University of Exeter, North Park Road, EX4 4QF Exeter, UK
| | - Björn Grüning
- Bioinformatics Group, Department of Computer Science, University of Freiburg, Georges-Koehler-Allee 106, 79110 Freiburg, Germany
- ZBSA Centre for Biological Systems Analysis, University of Freiburg, Hauptstr. 1, 79104 Freiburg, Germany
| | - Rolf Backofen
- Bioinformatics Group, Department of Computer Science, University of Freiburg, Georges-Koehler-Allee 106, 79110 Freiburg, Germany
- ZBSA Centre for Biological Systems Analysis, University of Freiburg, Hauptstr. 1, 79104 Freiburg, Germany
- Signalling Research Centres BIOSS and CIBSS, University of Freiburg, Schaenzlestr. 18, 79104 Freiburg, Germany
| |
Collapse
|
21
|
Venkata Subbaiah KC, Hedaya O, Wu J, Jiang F, Yao P. Mammalian RNA switches: Molecular rheostats in gene regulation, disease, and medicine. Comput Struct Biotechnol J 2019; 17:1326-1338. [PMID: 31741723 PMCID: PMC6849081 DOI: 10.1016/j.csbj.2019.10.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2019] [Revised: 09/30/2019] [Accepted: 10/07/2019] [Indexed: 01/12/2023] Open
Abstract
Alteration of RNA structure by environmental signals is a fundamental mechanism of gene regulation. For example, the riboswitch is a noncoding RNA regulatory element that binds a small molecule and causes a structural change in the RNA, thereby regulating transcription, splicing, or translation of an mRNA. The role of riboswitches in metabolite sensing and gene regulation in bacteria and other lower species was reported almost two decades ago, but riboswitches have not yet been discovered in mammals. An analog of the riboswitch, the protein-directed RNA switch (PDRS), has been identified as an important regulatory mechanism of gene expression in mammalian cells. RNA-binding proteins and microRNAs are two major executors of PDRS via their interaction with target transcripts in mammals. These protein-RNA interactions influence cellular functions by integrating environmental signals and intracellular pathways from disparate stimuli to modulate stability or translation of specific mRNAs. The discovery of a riboswitch in eukaryotes that is composed of a single class of thiamine pyrophosphate (TPP) suggests that additional ligand-sensing RNAs may be present to control eukaryotic or mammalian gene expression. In this review, we focus on protein-directed RNA switch mechanisms in mammals. We offer perspectives on the potential discovery of mammalian protein-directed and compound-dependent RNA switches that are related to human disease and medicine.
Collapse
Affiliation(s)
- Kadiam C Venkata Subbaiah
- Aab Cardiovascular Research Institute, Department of Medicine, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States
| | - Omar Hedaya
- Aab Cardiovascular Research Institute, Department of Medicine, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States.,Department of Biochemistry & Biophysics, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States
| | - Jiangbin Wu
- Aab Cardiovascular Research Institute, Department of Medicine, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States
| | - Feng Jiang
- Aab Cardiovascular Research Institute, Department of Medicine, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States.,Department of Biochemistry & Biophysics, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States
| | - Peng Yao
- Aab Cardiovascular Research Institute, Department of Medicine, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States.,Department of Biochemistry & Biophysics, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States.,The Center for RNA Biology, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States.,The Center for Biomedical Informatics, University of Rochester School of Medicine & Dentistry, Rochester, NY 14586, United States
| |
Collapse
|
22
|
Spasic A, Assmann SM, Bevilacqua PC, Mathews DH. Modeling RNA secondary structure folding ensembles using SHAPE mapping data. Nucleic Acids Res 2019; 46:314-323. [PMID: 29177466 PMCID: PMC5758915 DOI: 10.1093/nar/gkx1057] [Citation(s) in RCA: 61] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2017] [Accepted: 10/30/2017] [Indexed: 12/22/2022] Open
Abstract
RNA secondary structure prediction is widely used for developing hypotheses about the structures of RNA sequences, and structure can provide insight about RNA function. The accuracy of structure prediction is known to be improved using experimental mapping data that provide information about the pairing status of single nucleotides, and these data can now be acquired for whole transcriptomes using high-throughput sequencing. Prior methods for using these experimental data focused on predicting structures for sequences assuming that they populate a single structure. Most RNAs populate multiple structures, however, where the ensemble of strands populates structures with different sets of canonical base pairs. The focus on modeling single structures has been a bottleneck for accurately modeling RNA structure. In this work, we introduce Rsample, an algorithm for using experimental data to predict more than one RNA structure for sequences that populate multiple structures at equilibrium. We demonstrate, using SHAPE mapping data, that we can accurately model RNA sequences that populate multiple structures, including the relative probabilities of those structures. This program is freely available as part of the RNAstructure software package.
Collapse
Affiliation(s)
- Aleksandar Spasic
- Department of Biochemistry & Biophysics, University of Rochester Medical Center, Rochester, NY 14642, USA.,Center for RNA Biology, University of Rochester Medical Center, Rochester, NY 14642, USA
| | - Sarah M Assmann
- Department of Biology, Pennsylvania State University, University Park, PA 16802, USA
| | - Philip C Bevilacqua
- Department of Chemistry, Department of Biochemistry & Molecular Biology, Center for RNA Molecular Biology, Pennsylvania State University, University Park, PA 16802, USA
| | - David H Mathews
- Department of Biochemistry & Biophysics, University of Rochester Medical Center, Rochester, NY 14642, USA.,Center for RNA Biology, University of Rochester Medical Center, Rochester, NY 14642, USA.,Department of Biostatistics & Computational Biology, University of Rochester Medical Center, Rochester, NY 14642, USA
| |
Collapse
|
23
|
Genome-wide RNA structurome reprogramming by acute heat shock globally regulates mRNA abundance. Proc Natl Acad Sci U S A 2018; 115:12170-12175. [PMID: 30413617 PMCID: PMC6275526 DOI: 10.1073/pnas.1807988115] [Citation(s) in RCA: 70] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The heat shock response is crucial for organism survival in natural environments. RNA structure is known to influence numerous processes related to gene expression, but there have been few studies on the global RNA structurome as it prevails in vivo. Moreover, how heat shock rapidly affects RNA structure genome-wide in living systems remains unknown. We report here in vivo heat-regulated RNA structuromes. We applied Structure-seq chemical [dimethyl sulfate (DMS)] structure probing to rice (Oryza sativa L.) seedlings with and without 10 min of 42 °C heat shock and obtained structural data on >14,000 mRNAs. We show that RNA secondary structure broadly regulates gene expression in response to heat shock in this essential crop species. Our results indicate significant heat-induced elevation of DMS reactivity in the global transcriptome, revealing RNA unfolding over this biological temperature range. Our parallel Ribo-seq analysis provides no evidence for a correlation between RNA unfolding and heat-induced changes in translation, in contrast to the paradigm established in prokaryotes, wherein melting of RNA thermometers promotes translation. Instead, we find that heat-induced DMS reactivity increases correlate with significant decreases in transcript abundance, as quantified from an RNA-seq time course, indicating that mRNA unfolding promotes transcript degradation. The mechanistic basis for this outcome appears to be mRNA unfolding at both 5' and 3'-UTRs that facilitates access to the RNA degradation machinery. Our results thus reveal unexpected paradigms governing RNA structural changes and the eukaryotic RNA life cycle.
Collapse
|
24
|
Wang H, Lu X, Chen F, Ding Y, Zheng H, Wang L, Zhang G, Yang J, Bai Y, Li J, Wu J, Zhou M, Xu L. Landscape of SNPs-mediated lncRNA structural variations and their implication in human complex diseases. Brief Bioinform 2018; 21:85-95. [PMID: 30379995 DOI: 10.1093/bib/bby102] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2018] [Revised: 09/14/2018] [Accepted: 09/23/2018] [Indexed: 12/17/2022] Open
Abstract
An increasing number of functional studies shows that long noncoding RNAs (lncRNAs) are involved in many aspects of cellular physiology and fulfills a wide variety of regulatory roles at almost every stage of gene expression. A major feature of lncRNAs is the highly folded modular domains in transcripts. With improved modeling and definition, it is now feasible to explore and gain novel insights into the structural-functional relationship of lncRNAs and their association with complex human diseases. In this study, we utilized an automatic computational pipeline to scan lncRNA architecture at the genome-wide scale and to obtain a landscape of functional domains. An accurate alignment algorithm was performed to identify 40 triple pairs between single-nucleotide polymorphisms (SNPs), lncRNAs and diseases. In order to detect the potential contribution of a lncRNA's modular character, we estimated and evaluated structural rearrangements, which were derived from disease-associated SNPs. In addition, we focused on annotating and comparing the global and local heterogeneity of the wild-type and mutant lncRNAs. Assessing lncRNA architecture has yielded how variations in structured regions impact the molecular mechanisms of lncRNAs and how SNPs disturb binding and recruiting ability. These observations are the first glimpse of the 'lncRNA structurome' and make it possible to robustly explore and assemble intricate space conformation and their stress variation. This result also successfully demonstrates that lncRNA transcripts contain a complex structural landscape and highlights the proposed contribution of lncRNA structure in controlling RNA functions and disease mechanisms.
Collapse
Affiliation(s)
- Hong Wang
- Harbin Medical University.,Wenzhou Medical University
| | - Xiaoyan Lu
- Harbin Medical University.,Wenzhou Medical University
| | - Fukun Chen
- Harbin Medical University.,Wenzhou Medical University
| | - Yu Ding
- Harbin Medical University.,Wenzhou Medical University
| | - Hewei Zheng
- Harbin Medical University.,Wenzhou Medical University
| | - Lianzong Wang
- Harbin Medical University.,Wenzhou Medical University
| | - Guosi Zhang
- Harbin Medical University.,Wenzhou Medical University
| | - Jiaxin Yang
- Harbin Medical University.,Wenzhou Medical University
| | - Yu Bai
- Harbin Medical University.,Wenzhou Medical University
| | - Jing Li
- Harbin Medical University.,Wenzhou Medical University
| | - Jingqi Wu
- Harbin Medical University.,Wenzhou Medical University
| | | | | |
Collapse
|
25
|
A Functional riboSNitch in the 3' Untranslated Region of FKBP5 Alters MicroRNA-320a Binding Efficiency and Mediates Vulnerability to Chronic Post-Traumatic Pain. J Neurosci 2018; 38:8407-8420. [PMID: 30150364 DOI: 10.1523/jneurosci.3458-17.2018] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Revised: 07/11/2018] [Accepted: 07/13/2018] [Indexed: 01/30/2023] Open
Abstract
Previous studies have shown that common variants of the gene coding for FK506-binding protein 51 (FKBP5), a critical regulator of glucocorticoid sensitivity, affect vulnerability to stress-related disorders. In a previous report, FKBP5 rs1360780 was identified as a functional variant because of its effect on gene methylation. Here we report evidence for a novel functional FKBP5 allele, rs3800373. This study assessed the association between rs3800373 and post-traumatic chronic pain in 1607 women and men from two ethnically diverse human cohorts. The molecular mechanism through which rs3800373 affects adverse outcomes was established via in silico, in vivo, and in vitro analyses. The rs3800373 minor allele predicted worse adverse outcomes after trauma exposure, such that individuals with the minor (risk) allele developed more severe post-traumatic chronic musculoskeletal pain. Among these individuals, peritraumatic circulating FKBP5 expression levels increased as cortisol and glucocorticoid receptor (NR3C1) mRNA levels increased, consistent with increased glucocorticoid resistance. Bioinformatic, in vitro, and mutational analyses indicate that the rs3800373 minor allele reduces the binding of a stress- and pain-associated microRNA, miR-320a, to FKBP5 via altering the FKBP5 mRNA 3'UTR secondary structure (i.e., is a riboSNitch). This results in relatively greater FKBP5 translation, unchecked by miR-320a. Overall, these results identify an important gene-miRNA interaction influencing chronic pain risk in vulnerable individuals and suggest that exogenous methods to achieve targeted reduction in poststress FKBP5 mRNA expression may constitute useful therapeutic strategies.SIGNIFICANCE STATEMENT FKBP5 is a critical regulator of the stress response. Previous studies have shown that dysregulation of the expression of this gene plays a role in the pathogenesis of chronic pain development as well as a number of comorbid neuropsychiatric disorders. In the current study, we identified a functional allele (rs3800373) in the 3'UTR of FKBP5 that influences vulnerability to chronic post-traumatic pain in two ethnic cohorts. Using multiple complementary experimental approaches, we show that the FKBP5 rs3800373 minor allele alters the secondary structure of FKBP5 mRNA, decreasing the binding of a stress- and pain-associated microRNA, miR-320a. This results in relatively greater FKBP5 translation, unchecked by miR-320a, increasing glucocorticoid resistance and increasing vulnerability to post-traumatic pain.
Collapse
|
26
|
Lackey L, Coria A, Woods C, McArthur E, Laederach A. Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure. RNA (NEW YORK, N.Y.) 2018; 24:513-528. [PMID: 29317542 PMCID: PMC5855952 DOI: 10.1261/rna.064469.117] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2017] [Accepted: 01/04/2018] [Indexed: 05/22/2023]
Abstract
The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1 Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.
Collapse
Affiliation(s)
- Lela Lackey
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Aaztli Coria
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Chanin Woods
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Evonne McArthur
- School of Medicine, Vanderbilt University, Nashville, Tennessee 37232, USA
| | - Alain Laederach
- Department of Biology, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| |
Collapse
|
27
|
Woods CT, Laederach A. Classification of RNA structure change by 'gazing' at experimental data. Bioinformatics 2018; 33:1647-1655. [PMID: 28130241 PMCID: PMC5447233 DOI: 10.1093/bioinformatics/btx041] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2016] [Accepted: 01/20/2017] [Indexed: 11/12/2022] Open
Abstract
Motivation Mutations (or Single Nucleotide Variants) in folded RiboNucleic Acid structures that cause local or global conformational change are riboSNitches. Predicting riboSNitches is challenging, as it requires making two, albeit related, structure predictions. The data most often used to experimentally validate riboSNitch predictions is Selective 2' Hydroxyl Acylation by Primer Extension, or SHAPE. Experimentally establishing a riboSNitch requires the quantitative comparison of two SHAPE traces: wild-type (WT) and mutant. Historically, SHAPE data was collected on electropherograms and change in structure was evaluated by 'gel gazing.' SHAPE data is now routinely collected with next generation sequencing and/or capillary sequencers. We aim to establish a classifier capable of simulating human 'gazing' by identifying features of the SHAPE profile that human experts agree 'looks' like a riboSNitch. Results We find strong quantitative agreement between experts when RNA scientists 'gaze' at SHAPE data and identify riboSNitches. We identify dynamic time warping and seven other features predictive of the human consensus. The classSNitch classifier reported here accurately reproduces human consensus for 167 mutant/WT comparisons with an Area Under the Curve (AUC) above 0.8. When we analyze 2019 mutant traces for 17 different RNAs, we find that features of the WT SHAPE reactivity allow us to improve thermodynamic structure predictions of riboSNitches. This is significant, as accurate RNA structural analysis and prediction is likely to become an important aspect of precision medicine. Availability and Implementation The classSNitch R package is freely available at http://classsnitch.r-forge.r-project.org . Contact alain@email.unc.edu. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Chanin Tolson Woods
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Alain Laederach
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| |
Collapse
|
28
|
Lin L, McKerrow WH, Richards B, Phonsom C, Lawrence CE. Characterization and visualization of RNA secondary structure Boltzmann ensemble via information theory. BMC Bioinformatics 2018; 19:82. [PMID: 29506466 PMCID: PMC5836418 DOI: 10.1186/s12859-018-2078-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2017] [Accepted: 02/20/2018] [Indexed: 12/26/2022] Open
Abstract
Background The nearest neighbor model and associated dynamic programming algorithms allow for the efficient estimation of the RNA secondary structure Boltzmann ensemble. However because a given RNA secondary structure only contains a fraction of the possible helices that could form from a given sequence, the Boltzmann ensemble is multimodal. Several methods exist for clustering structures and finding those modes. However less focus is given to exploring the underlying reasons for this multimodality: the presence of conflicting basepairs. Information theory, or more specifically mutual information, provides a method to identify those basepairs that are key to the secondary structure. Results To this end we find most informative basepairs and visualize the effect of these basepairs on the secondary structure. Knowing whether a most informative basepair is present tells us not only the status of the particular pair but also provides a large amount of information about which other pairs are present or not present. We find that a few basepairs account for a large amount of the structural uncertainty. The identification of these pairs indicates small changes to sequence or stability that will have a large effect on structure. Conclusion We provide a novel algorithm that uses mutual information to identify the key basepairs that lead to a multimodal Boltzmann distribution. We then visualize the effect of these pairs on the overall Boltzmann ensemble.
Collapse
Affiliation(s)
- Luan Lin
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, 20993, MD, USA
| | - Wilson H McKerrow
- Division of Applied Mathematics, Brown University, Providence, 02912, RI, USA
| | | | - Chukiat Phonsom
- Department of Mathematics, University of Southern California, Los Angeles, 90089, CA, USA
| | - Charles E Lawrence
- Division of Applied Mathematics, Brown University, Providence, 02912, RI, USA.
| |
Collapse
|
29
|
An RNA structure-mediated, posttranscriptional model of human α-1-antitrypsin expression. Proc Natl Acad Sci U S A 2017; 114:E10244-E10253. [PMID: 29109288 PMCID: PMC5703279 DOI: 10.1073/pnas.1706539114] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
Protein and mRNA expression are in most cases poorly correlated, which suggests that the posttranscriptional regulatory program of a cell is an important component of gene expression. This regulatory network is still poorly understood, including how RNA structure quantitatively contributes to translational control. We present here a series of structural and functional experiments that together allow us to derive a quantitative, structure-dependent model of translation that accurately predicts translation efficiency in reporter assays and primary human tissue for a complex and medically important protein, α-1-antitrypsin. Our model demonstrates the importance of accurate, experimentally derived RNA structural models partnered with Kozak sequence information to explain protein expression and suggests a strategy by which α-1-antitrypsin expression may be increased in diseased individuals. Chronic obstructive pulmonary disease (COPD) affects over 65 million individuals worldwide, where α-1-antitrypsin deficiency is a major genetic cause of the disease. The α-1-antitrypsin gene, SERPINA1, expresses an exceptional number of mRNA isoforms generated entirely by alternative splicing in the 5′-untranslated region (5′-UTR). Although all SERPINA1 mRNAs encode exactly the same protein, expression levels of the individual mRNAs vary substantially in different human tissues. We hypothesize that these transcripts behave unequally due to a posttranscriptional regulatory program governed by their distinct 5′-UTRs and that this regulation ultimately determines α-1-antitrypsin expression. Using whole-transcript selective 2′-hydroxyl acylation by primer extension (SHAPE) chemical probing, we show that splicing yields distinct local 5′-UTR secondary structures in SERPINA1 transcripts. Splicing in the 5′-UTR also changes the inclusion of long upstream ORFs (uORFs). We demonstrate that disrupting the uORFs results in markedly increased translation efficiencies in luciferase reporter assays. These uORF-dependent changes suggest that α-1-antitrypsin protein expression levels are controlled at the posttranscriptional level. A leaky-scanning model of translation based on Kozak translation initiation sequences alone does not adequately explain our quantitative expression data. However, when we incorporate the experimentally derived RNA structure data, the model accurately predicts translation efficiencies in reporter assays and improves α-1-antitrypsin expression prediction in primary human tissues. Our results reveal that RNA structure governs a complex posttranscriptional regulatory program of α-1-antitrypsin expression. Crucially, these findings describe a mechanism by which genetic alterations in noncoding gene regions may result in α-1-antitrypsin deficiency.
Collapse
|
30
|
Martens L, Rühle F, Stoll M. LncRNA secondary structure in the cardiovascular system. Noncoding RNA Res 2017; 2:137-142. [PMID: 30159432 PMCID: PMC6084829 DOI: 10.1016/j.ncrna.2017.12.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2017] [Revised: 11/28/2017] [Accepted: 12/08/2017] [Indexed: 01/27/2023] Open
Abstract
Long non-coding RNAs (lncRNAs) have been increasingly studied during the past decade. This led to an immense number of annotated transcripts, out of which many were linked to a diverse range of biological mechanisms and diseases. Due to the variety of their regulatory potential, they are seen as an important link in understanding complex epigenetic mechanisms. Prominent examples of lncRNAs in the cardiovascular system are ANRIL, Braveheart, MALAT1 and HOTAIR which have been excessively studied. But despite the impressive number of described transcripts, only a few examples are characterized functionally. One way to do this is to identify accessible structural domains in the RNA secondary structure which have the ability to bind to DNA, RNA or proteins. Through recent improvements in computational as well as experimental methods, this exploration of secondary structure became not only more efficient than traditional methods like crystallization, but also feasible to investigate whole genome RNA structures. The purpose of this review is to highlight the recent advances in secondary structure probing methods and how these can be applied in order to investigate the functional roles of lncRNAs in the cardiovascular system.
Collapse
Affiliation(s)
- Leonie Martens
- Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
| | - Frank Rühle
- Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
| | - Monika Stoll
- Department of Genetic Epidemiology, Institute of Human Genetics, University of Münster, Münster, Germany
- Department of Biochemistry, Genetic Epidemiology and Statistical Genetics, CARIM School for Cardiovascular Diseases, Maastricht Center for Systems Biology (MaCSBio), Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
31
|
Woods CT, Lackey L, Williams B, Dokholyan NV, Gotz D, Laederach A. Comparative Visualization of the RNA Suboptimal Conformational Ensemble In Vivo. Biophys J 2017. [PMID: 28625696 PMCID: PMC5529173 DOI: 10.1016/j.bpj.2017.05.031] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
When a ribonucleic acid (RNA) molecule folds, it often does not adopt a single, well-defined conformation. The folding energy landscape of an RNA is highly dependent on its nucleotide sequence and molecular environment. Cellular molecules sometimes alter the energy landscape, thereby changing the ensemble of likely low-energy conformations. The effects of these energy landscape changes on the conformational ensemble are particularly challenging to visualize for large RNAs. We have created a robust approach for visualizing the conformational ensemble of RNAs that is well suited for in vitro versus in vivo comparisons. Our method creates a stable map of conformational space for a given RNA sequence. We first identify single point mutations in the RNA that maximally sample suboptimal conformational space based on the ensemble’s partition function. Then, we cluster these diverse ensembles to identify the most diverse partition functions for Boltzmann stochastic sampling. By using, to our knowledge, a novel nestedness distance metric, we iteratively add mutant suboptimal ensembles to converge on a stable 2D map of conformational space. We then compute the selective 2′ hydroxyl acylation by primer extension (SHAPE)-directed ensemble for the RNA folding under different conditions, and we project these ensembles on the map to visualize. To validate our approach, we established a conformational map of the Vibrio vulnificus add adenine riboswitch that reveals five classes of structures. In the presence of adenine, projection of the SHAPE-directed sampling correctly identified the on-conformation; without the ligand, only off-conformations were visualized. We also collected the whole-transcript in vitro and in vivo SHAPE-MaP for human β-actin messenger RNA that revealed similar global folds in both conditions. Nonetheless, a comparison of in vitro and in vivo data revealed that specific regions exhibited significantly different SHAPE-MaP profiles indicative of structural rearrangements, including rearrangement consistent with binding of the zipcode protein in a region distal to the stop codon.
Collapse
Affiliation(s)
- Chanin T Woods
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina; Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - Lela Lackey
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - Benfeard Williams
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - Nikolay V Dokholyan
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - David Gotz
- Carolina Health Informatics Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina; School of Information and Library Science, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - Alain Laederach
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina; Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina.
| |
Collapse
|
32
|
Choudhary K, Deng F, Aviran S. Comparative and integrative analysis of RNA structural profiling data: current practices and emerging questions. QUANTITATIVE BIOLOGY 2017; 5:3-24. [PMID: 28717530 PMCID: PMC5510538 DOI: 10.1007/s40484-017-0093-6] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2016] [Revised: 12/08/2016] [Accepted: 12/15/2016] [Indexed: 12/30/2022]
Abstract
BACKGROUND Structure profiling experiments provide single-nucleotide information on RNA structure. Recent advances in chemistry combined with application of high-throughput sequencing have enabled structure profiling at transcriptome scale and in living cells, creating unprecedented opportunities for RNA biology. Propelled by these experimental advances, massive data with ever-increasing diversity and complexity have been generated, which give rise to new challenges in interpreting and analyzing these data. RESULTS We review current practices in analysis of structure profiling data with emphasis on comparative and integrative analysis as well as highlight emerging questions. Comparative analysis has revealed structural patterns across transcriptomes and has become an integral component of recent profiling studies. Additionally, profiling data can be integrated into traditional structure prediction algorithms to improve prediction accuracy. CONCLUSIONS To keep pace with experimental developments, methods to facilitate, enhance and refine such analyses are needed. Parallel advances in analysis methodology will complement profiling technologies and help them reach their full potential.
Collapse
Affiliation(s)
| | | | - Sharon Aviran
- Department of Biomedical Engineering and Genome Center, University of California at Davis, Davis, CA 95616, USA
| |
Collapse
|
33
|
|
34
|
Abstract
Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.
Collapse
Affiliation(s)
- Philip C Bevilacqua
- Department of Chemistry.,Department of Biochemistry and Molecular Biology.,Center for RNA Molecular Biology
| | | | - Zhao Su
- Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802;
| | - Sarah M Assmann
- Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802;
| |
Collapse
|
35
|
Kutchko KM, Laederach A. Transcending the prediction paradigm: novel applications of SHAPE to RNA function and evolution. WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [PMID: 27396578 PMCID: PMC5179297 DOI: 10.1002/wrna.1374] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/04/2016] [Revised: 04/29/2016] [Accepted: 05/23/2016] [Indexed: 12/31/2022]
Abstract
Selective 2′‐hydroxyl acylation analyzed by primer extension (SHAPE) provides information on RNA structure at single‐nucleotide resolution. It is most often used in conjunction with RNA secondary structure prediction algorithms as a probabilistic or thermodynamic restraint. With the recent advent of ultra‐high‐throughput approaches for collecting SHAPE data, the applications of this technology are extending beyond structure prediction. In this review, we discuss recent applications of SHAPE data in the transcriptomic context and how this new experimental paradigm is changing our understanding of these experiments and RNA folding in general. SHAPE experiments probe both the secondary and tertiary structure of an RNA, suggesting that model‐free approaches for within and comparative RNA structure analysis can provide significant structural insight without the need for a full structural model. New methods incorporating SHAPE at different nucleotide resolutions are required to parse these transcriptomic data sets to transcend secondary structure modeling with global structural metrics. These ‘multiscale’ approaches provide deeper insights into RNA global structure, evolution, and function in the cell. WIREs RNA 2017, 8:e1374. doi: 10.1002/wrna.1374 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Katrina M Kutchko
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Alain Laederach
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.,Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| |
Collapse
|
36
|
Hamada M, Ono Y, Kiryu H, Sato K, Kato Y, Fukunaga T, Mori R, Asai K. Rtools: a web server for various secondary structural analyses on single RNA sequences. Nucleic Acids Res 2016; 44:W302-7. [PMID: 27131356 PMCID: PMC4987903 DOI: 10.1093/nar/gkw337] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2016] [Accepted: 04/15/2016] [Indexed: 11/12/2022] Open
Abstract
The secondary structures, as well as the nucleotide sequences, are the important features of RNA molecules to characterize their functions. According to the thermodynamic model, however, the probability of any secondary structure is very small. As a consequence, any tool to predict the secondary structures of RNAs has limited accuracy. On the other hand, there are a few tools to compensate the imperfect predictions by calculating and visualizing the secondary structural information from RNA sequences. It is desirable to obtain the rich information from those tools through a friendly interface. We implemented a web server of the tools to predict secondary structures and to calculate various structural features based on the energy models of secondary structures. By just giving an RNA sequence to the web server, the user can get the different types of solutions of the secondary structures, the marginal probabilities such as base-paring probabilities, loop probabilities and accessibilities of the local bases, the energy changes by arbitrary base mutations as well as the measures for validations of the predicted secondary structures. The web server is available at http://rtools.cbrc.jp, which integrates software tools, CentroidFold, CentroidHomfold, IPKnot, CapR, Raccess, Rchange and RintD.
Collapse
Affiliation(s)
- Michiaki Hamada
- Department of Electrical Engineering and Bioscience, Faculty of Science and Engineering, Waseda University, 55N-06-10, 3-4-1, Okubo Shinjuku-ku, Tokyo 169-8555, Japan Artificial Intelligence Research Center (AIRC), National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, 135-0064 Tokyo, Japan
| | - Yukiteru Ono
- IMSBIO Co., Ltd, 4-21-1-601 Higashi-Ikebukuro, Toshima-ku, Tokyo 170-0013, Japan
| | - Hisanori Kiryu
- Graduate School of Frontier Sciences, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8562, Japan
| | - Kengo Sato
- Department of Biosciences and Informatics, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Yuki Kato
- Center for iPS Cell Research and Application (CiRA), Kyoto University, Shogoin, Sakyo-ku, Kyoto 606-8507, Japan
| | - Tsukasa Fukunaga
- Graduate School of Frontier Sciences, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8562, Japan
| | - Ryota Mori
- Graduate School of Frontier Sciences, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8562, Japan
| | - Kiyoshi Asai
- Artificial Intelligence Research Center (AIRC), National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, 135-0064 Tokyo, Japan Graduate School of Frontier Sciences, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8562, Japan
| |
Collapse
|
37
|
Long non-coding RNA Databases in Cardiovascular Research. GENOMICS PROTEOMICS & BIOINFORMATICS 2016; 14:191-9. [PMID: 27049585 PMCID: PMC4996844 DOI: 10.1016/j.gpb.2016.03.001] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/18/2016] [Revised: 03/16/2016] [Accepted: 03/17/2016] [Indexed: 12/05/2022]
Abstract
With the rising interest in the regulatory functions of long non-coding RNAs (lncRNAs) in complex human diseases such as cardiovascular diseases, there is an increasing need in public databases offering comprehensive and integrative data for all aspects of these versatile molecules. Recently, a variety of public data repositories that specialized in lncRNAs have been developed, which make use of huge high-throughput data particularly from next-generation sequencing (NGS) approaches. Here, we provide an overview of current lncRNA databases covering basic and functional annotation, lncRNA expression and regulation, interactions with other biomolecules, and genomic variants influencing the structure and function of lncRNAs. The prominent lncRNA antisense noncoding RNA in the INK4 locus (ANRIL), which has been unequivocally associated with coronary artery disease through genome-wide association studies (GWAS), serves as an example to demonstrate the features of each individual database.
Collapse
|
38
|
Rogers E, Heitsch C. New insights from cluster analysis methods for RNA secondary structure prediction. WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 7:278-94. [PMID: 26971529 DOI: 10.1002/wrna.1334] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/08/2015] [Revised: 12/03/2015] [Accepted: 12/17/2015] [Indexed: 01/12/2023]
Abstract
A widening gap exists between the best practices for RNA secondary structure prediction developed by computational researchers and the methods used in practice by experimentalists. Minimum free energy predictions, although broadly used, are outperformed by methods which sample from the Boltzmann distribution and data mine the results. In particular, moving beyond the single structure prediction paradigm yields substantial gains in accuracy. Furthermore, the largest improvements in accuracy and precision come from viewing secondary structures not at the base pair level but at lower granularity/higher abstraction. This suggests that random errors affecting precision and systematic ones affecting accuracy are both reduced by this 'fuzzier' view of secondary structures. Thus experimentalists who are willing to adopt a more rigorous, multilayered approach to secondary structure prediction by iterating through these levels of granularity will be much better able to capture fundamental aspects of RNA base pairing. WIREs RNA 2016, 7:278-294. doi: 10.1002/wrna.1334 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Emily Rogers
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0765, USA
| | - Christine Heitsch
- School of Mathematics, Georgia Institute of Technology, Atlanta, GA 30332-0160, USA
| |
Collapse
|
39
|
Hu X, Wu Y, Lu ZJ, Yip KY. Analysis of sequencing data for probing RNA secondary structures and protein–RNA binding in studying posttranscriptional regulations. Brief Bioinform 2015; 17:1032-1043. [DOI: 10.1093/bib/bbv106] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2015] [Revised: 10/11/2015] [Indexed: 11/12/2022] Open
|