Abstract
Human and mouse genome sequences contain roughly 100,000 regions that are unalignable in primary sequence and neighbor corresponding alignable regions between both organisms. These pairs are generally assumed to be nonconserved, although the level of structural conservation between these has never been investigated. Owing to the limitations in computational methods, comparative genomics has been lacking the ability to compare such nonconserved sequence regions for conserved structural RNA elements. We have investigated the presence of structural RNA elements by conducting a local structural alignment, using FOLDALIGN, on a subset of these 100,000 corresponding regions and estimate that 1800 contain common RNA structures. Comparing our results with the recent mapping of transcribed fragments (transfrags) in human, we find that high-scoring candidates are twice as likely to be found in regions overlapped by transfrags than regions that are not overlapped by transfrags. To verify the coexpression between predicted candidates in human and mouse, we conducted expression studies by RT-PCR and Northern blotting on mouse candidates, which overlap with transfrags on human chromosome 20. RT-PCR results confirmed expression of 32 out of 36 candidates, whereas Northern blots confirmed four out of 12 candidates. Furthermore, many RT-PCR results indicate differential expression in different tissues. Hence, our findings suggest that there are corresponding regions between human and mouse, which contain expressed non-coding RNA sequences not alignable in primary sequence.
MeSH Terms
Animals
Base Pairing
Base Sequence
Chickens/genetics
Chromosome Mapping
Chromosomes, Human, Pair 20
Conserved Sequence
Dogs
Genome
Genome, Human
Humans
Mice/genetics
Nucleic Acid Conformation
RNA/chemistry
Rats
Sequence Analysis, RNA/statistics & numerical data
Sequence Homology, Nucleic Acid
Software
Transcription, Genetic
Authors & Affiliations
5 authors, click to expand affiliations / ORCID
Torarinsson Elfar
Division of Genetics and Bioinformatics, IBHV, The Royal Veterinary and Agricultural University, 1870 Frederiksberg C, Denmark.
Sawera Milena
Havgaard Jakob H
Fredholm Merete
Gorodkin Jan
References (17)
17 references, click to expand
-
Primer3 on the WWW for general users and for biologist programmers.
Methods Mol Biol. 2000;132:365-86
PMID: 10547847
-
Accumulation of unstable promoter-associated transcripts upon loss of the nuclear exosome subunit Rrp6p in Saccharomyces cerevisiae.
Proc Natl Acad Sci U S A. 2006 Feb 28;103(9):3262-7
PMID: 16484372
-
No evidence that mRNAs have lower folding free energies than random sequences with the same dinucleotide distribution.
Nucleic Acids Res. 1999 Dec 15;27(24):4816-22
PMID: 10572183
-
BLAT--the BLAST-like alignment tool.
Genome Res. 2002 Apr;12(4):656-64
PMID: 11932250
-
Transcript abundance in yeast varies over six orders of magnitude.
J Biol Chem. 2002 Apr 26;277(17):14363-6
PMID: 11882647
-
Rfam: an RNA family database.
Nucleic Acids Res. 2003 Jan 1;31(1):439-41
PMID: 12520045
-
Noncoding RNA gene detection using comparative sequence analysis.
BMC Bioinformatics. 2001;2:8
PMID: 11801179
-
Exploiting conserved structure for faster annotation of non-coding RNAs without loss of accuracy.
Bioinformatics. 2004 Aug 4;20 Suppl 1:i334-41
PMID: 15262817
-
Significance of nucleotide sequence alignments: a method for random sequence permutation that preserves dinucleotide and codon usage.
Mol Biol Evol. 1985 Nov;2(6):526-38
PMID: 3870875
-
Basic local alignment search tool.
J Mol Biol. 1990 Oct 5;215(3):403-10
PMID: 2231712
-
Pairwise local structural alignment of RNA sequences with sequence similarity less than 40%.
Bioinformatics. 2005 May 1;21(9):1815-24
PMID: 15657094
-
Transcriptional maps of 10 human chromosomes at 5-nucleotide resolution.
Science. 2005 May 20;308(5725):1149-54
PMID: 15790807
-
Cryptic pol II transcripts are degraded by a nuclear quality control pathway involving a new poly(A) polymerase.
Cell. 2005 Jun 3;121(5):725-37
PMID: 15935759
-
The transcriptional landscape of the mammalian genome.
Science. 2005 Sep 2;309(5740):1559-63
PMID: 16141072
-
A systematic search for new mammalian noncoding RNAs indicates little conserved intergenic transcription.
BMC Genomics. 2005;6:104
PMID: 16083503
-
Mapping of conserved RNA secondary structures predicts thousands of functional noncoding RNAs in the human genome.
Nat Biotechnol. 2005 Nov;23(11):1383-90
PMID: 16273071
-
Experimental validation of the regulated expression of large numbers of non-coding RNAs from the mouse genome.
Genome Res. 2006 Jan;16(1):11-9
PMID: 16344565