Noncoding RNA Gene Detection Using Comparative Sequence Analysis

Overview

Journal BMC Bioinformatics

Publisher Biomed Central

Specialty Biology

Date 2002 Jan 22

PMID 11801179

Citations 208

Authors

E Rivas

S R Eddy

Affiliations

Soon will be listed here.

Abstract

Background: Noncoding RNA genes produce transcripts that exert their function without ever producing proteins. Noncoding RNA gene sequences do not have strong statistical signals, unlike protein coding genes. A reliable general purpose computational genefinder for noncoding RNA genes has been elusive.

Results: We describe a comparative sequence analysis algorithm for detecting novel structural RNA genes. The key idea is to test the pattern of substitutions observed in a pairwise alignment of two homologous sequences. A conserved coding region tends to show a pattern of synonymous substitutions, whereas a conserved structural RNA tends to show a pattern of compensatory mutations consistent with some base-paired secondary structure. We formalize this intuition using three probabilistic "pair-grammars": a pair stochastic context free grammar modeling alignments constrained by structural RNA evolution, a pair hidden Markov model modeling alignments constrained by coding sequence evolution, and a pair hidden Markov model modeling a null hypothesis of position-independent evolution. Given an input pairwise sequence alignment (e.g. from a BLASTN comparison of two related genomes) we classify the alignment into the coding, RNA, or null class according to the posterior probability of each class.

Conclusions: We have implemented this approach as a program, QRNA, which we consider to be a prototype structural noncoding RNA genefinder. Tests suggest that this approach detects noncoding RNA genes with a fair degree of reliability.

Citing Articles

TSS-Captur: a user-friendly pipeline for characterizing unclassified RNA transcripts.

Witte Paz M, Vogel T, Nieselt K NAR Genom Bioinform. 2024; 6(4):lqae168.

PMID: 39703424 PMC: 11655288. DOI: 10.1093/nargab/lqae168.

Comparative RNA Genomics.

Backofen R, Gorodkin J, Hofacker I, Stadler P Methods Mol Biol. 2024; 2802:347-393.

PMID: 38819565 DOI: 10.1007/978-1-0716-3838-5_12.

Whole-genome sequencing and evolutionary analysis of the wild edible mushroom, .

Li Y, Yang T, Qiao J, Liang J, Li Z, Sa W Front Microbiol. 2024; 14:1309703.

PMID: 38361578 PMC: 10868677. DOI: 10.3389/fmicb.2023.1309703.

Function and phylogeny support the independent evolution of an ASIC-like Deg/ENaC channel in the Placozoa.

Elkhatib W, Yanez-Guerra L, Mayorova T, Currie M, Singh A, Perera M Commun Biol. 2023; 6(1):951.

PMID: 37723223 PMC: 10507113. DOI: 10.1038/s42003-023-05312-0.

Quantification of the Diversity in Gene Structures Using the Principles of Polarization Mapping.

Zimnyakov D, Alonova M, Skripal A, Dobdin S, Feodorova V Curr Issues Mol Biol. 2023; 45(2):1720-1740.

PMID: 36826056 PMC: 9955201. DOI: 10.3390/cimb45020111.

References

Kelley R, Kuroda M . Noncoding RNA genes in dosage compensation and imprinting. Cell. 2000; 103(1):9-12. DOI: 10.1016/s0092-8674(00)00099-4. View

Huttenhofer A, Kiefmann M, OBrien J, Lehrach H, Bachellerie J, Brosius J . RNomics: an experimental approach that identifies 201 candidates for novel, small, non-messenger RNAs in mouse. EMBO J. 2001; 20(11):2943-53. PMC: 125495. DOI: 10.1093/emboj/20.11.2943. View

Selinger D, Cheung K, Mei R, Johansson E, Richmond C, Blattner F . RNA expression analysis using a 30 base pair resolution Escherichia coli genome array. Nat Biotechnol. 2000; 18(12):1262-8. DOI: 10.1038/82367. View

Lease R, Belfort M . Riboregulation by DsrA RNA: trans-actions for global economy. Mol Microbiol. 2000; 38(4):667-72. DOI: 10.1046/j.1365-2958.2000.02162.x. View

Erdmann V, Barciszewska M, Szymanski M, Hochberg A, De Groot N, Barciszewski J . The non-coding RNAs as riboregulators. Nucleic Acids Res. 2000; 29(1):189-93. PMC: 29806. DOI: 10.1093/nar/29.1.189. View

Zuker M, Stiegler P . Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information. Nucleic Acids Res. 1981; 9(1):133-48. PMC: 326673. DOI: 10.1093/nar/9.1.133. View

Ikemura T . Codon usage and tRNA content in unicellular and multicellular organisms. Mol Biol Evol. 1985; 2(1):13-34. DOI: 10.1093/oxfordjournals.molbev.a040335. View

Le S, Chen J, Currey K, Maizel Jr J . A program for predicting significant RNA secondary structures. Comput Appl Biosci. 1988; 4(1):153-9. DOI: 10.1093/bioinformatics/4.1.153. View

Chen J, Le S, Shapiro B, Currey K, Maizel J . A computational procedure for assessing the significance of RNA secondary structure. Comput Appl Biosci. 1990; 6(1):7-18. DOI: 10.1093/bioinformatics/6.1.7. View

10.

Thomas J, Lea K, Blumenthal T . The spliceosomal snRNAs of Caenorhabditis elegans. Nucleic Acids Res. 1990; 18(9):2633-42. PMC: 330746. DOI: 10.1093/nar/18.9.2633. View

11.

Altschul S, Gish W, Miller W, Myers E, Lipman D . Basic local alignment search tool. J Mol Biol. 1990; 215(3):403-10. DOI: 10.1016/S0022-2836(05)80360-2. View

12.

Larsen N, Zwieb C . SRP-RNA sequence alignment and secondary structure. Nucleic Acids Res. 1991; 19(2):209-15. PMC: 333582. DOI: 10.1093/nar/19.2.209. View

13.

Altschul S . Amino acid substitution matrices from an information theoretic perspective. J Mol Biol. 1991; 219(3):555-65. PMC: 7130686. DOI: 10.1016/0022-2836(91)90193-a. View

14.

Henikoff S, Henikoff J . Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci U S A. 1992; 89(22):10915-9. PMC: 50453. DOI: 10.1073/pnas.89.22.10915. View

15.

Steinberg S, Misch A, Sprinzl M . Compilation of tRNA sequences and sequences of tRNA genes. Nucleic Acids Res. 1993; 21(13):3011-5. PMC: 309727. DOI: 10.1093/nar/21.13.3011. View

16.

Van de Peer Y, Van den Broeck I, De Rijk P, De Wachter R . Database on the structure of small ribosomal subunit RNA. Nucleic Acids Res. 1994; 22(17):3488-94. PMC: 308309. DOI: 10.1093/nar/22.17.3488. View

17.

Stormo G, Haussler D . Optimally parsing a sequence into different classes based on multiple types of evidence. Proc Int Conf Intell Syst Mol Biol. 1994; 2:369-75. View

18.

Bult C, White O, Olsen G, Zhou L, Fleischmann R, Sutton G . Complete genome sequence of the methanogenic archaeon, Methanococcus jannaschii. Science. 1996; 273(5278):1058-73. DOI: 10.1126/science.273.5278.1058. View

19.

Burge C, Karlin S . Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997; 268(1):78-94. DOI: 10.1006/jmbi.1997.0951. View

20.

Blattner F, Plunkett 3rd G, Bloch C, Perna N, Burland V, Riley M . The complete genome sequence of Escherichia coli K-12. Science. 1997; 277(5331):1453-62. DOI: 10.1126/science.277.5331.1453. View