Uncovering information on expression of natural antisense transcripts in Affymetrix MOE430 datasets
© Oeder et al; licensee BioMed Central Ltd. 2007
Received: 14 February 2007
Accepted: 28 June 2007
Published: 28 June 2007
The function and significance of the widespread expression of natural antisense transcripts (NATs) is largely unknown. The ability to quantitatively assess changes in NAT expression for many different transcripts in multiple samples would facilitate our understanding of this relatively new class of RNA molecules.
Here, we demonstrate that standard expression analysis Affymetrix MOE430 and HG-U133 GeneChips contain hundreds of probe sets that detect NATs. Probe sets carrying a "Negative Strand Matching Probes" annotation in NetAffx were validated using Ensembl by manual and automated approaches. More than 50 % of the 1,113 probe sets with "Negative Strand Matching Probes" on the MOE430 2.0 GeneChip were confirmed as detecting NATs. Expression of selected antisense transcripts as indicated by Affymetrix data was confirmed using strand-specific RT-PCR. Thus, Affymetrix datasets can be mined to reveal information about the regulated expression of a considerable number of NATs. In a correlation analysis of 179 sense-antisense (SAS) probe set pairs using publicly available data from 1637 MOE430 2.0 GeneChips a significant number of SAS transcript pairs were found to be positively correlated.
Standard expression analysis Affymetrix GeneChips can be used to measure many different NATs. The large amount of samples deposited in microarray databases represents a valuable resource for a quantitative analysis of NAT expression and regulation in different cells, tissues and biological conditions.
Natural antisense transcripts (NATs) are RNA molecules harboring sequences complementary to other transcripts. Expression of endogenous NATs was first observed in viruses  and prokaryotes , followed by its detection also in eukaryotes . The number of sense-antisense (SAS) transcripts slowly grew during the 1990's, as more examples were found in the study of the regulation of expression of individual genes. In some of these cases, a regulatory function for antisense RNA in controlling the mRNA or protein levels of the respective sense gene product was described (for reviews see ref. [4, 5]). The availability of large scale cDNA sequencing then led to the in silico discovery of significant numbers of SAS transcript pairs in different species, including yeast, plants and mammals [6–14]. Most recently, systematic genome-wide detection of transcription by strand-specific tiling oligonucleotide arrays or sequencing (CAGE) substantiated the notion that expression of antisense transcripts is a widespread phenomenon, accounting for up to 50 % of all transcripts [15–17]. To date, only few studies have analyzed quantitative changes in the expression of larger numbers of NATs under different experimental conditions. The function(s) of NATs and the regulation of their expression in relation to the corrresponding sense transcripts are largely unknown [4, 18]. One possible function of NATs is decreasing sense transcript abundance, e.g. by RNAi-mediated degradation . In this case, the levels of sense and antisense transcript would be expected to be negatively correlated. The effects of NATs are not necessarily mediated by RNA-RNA hybridization but can also be brought about by RNA-protein interactions [20, 21].
Oligonucleotide microarrays such as Affymetrix GeneChips detect labeled cRNA in a strand-specific manner. Standard labeling methods based on incorporation of a T7 binding site with an oligo-dT primer during first strand cDNA synthesis not only amplify selectively polyA mRNA, but also conserve strand-specificity of the resulting cRNA. This technology therefore in principle allows the measurement of transcripts in an orientation-specific manner. Tiling arrays have indeed been used for a detailed experimental investigation of transcriptional activity on human chromosomes [15, 22, 23]. These studies have greatly advanced our understanding of the extent of genome-wide transcriptional activity and the abundance of NATs and other non-coding RNA transcripts. However, the requirement of multiple arrays to perform a genome-wide tiling array analysis of a single RNA sample makes larger experiments prohibitively expensive, and therefore has to date precluded a comprehensive investigation of NATs expression and regulation in many different experimental conditions. In contrast, expression analysis GeneChips such as MOE430 for murine and U133 for human RNA detect and measure a comprehensive set of the majority of protein coding mRNAs on one single microarray, making global mRNA expression profiling experiments relatively affordable. Their extensive use in functional genomics studies has resulted in a wealth of data sets deposited in public data repositories. In this paper, we point to the presence of a considerable number of probe sets detecting NATs on Affymetrix GeneChip expression analysis arrays. Probe sets potentially detecting NAT were validated in silico, and microarray results showing expression of NATs were confirmed by strand-specific RT-PCR. The large number of Affymetrix datasets in public domain repositories like Gene Expression Omnibus (GEO) can therefore be mined for the expression and regulation of hundreds of NATs.
Results and Discussion
Affymetrix GeneChip probe sets with "Negative Strand Matching Probes": in silico validation
For comparison, we manually validated the orientation of a subset of 391 potential NAT-detecting probe sets by performing a Blast search using the probe set target sequence as query in Ensembl. Comparing the results of the automatic and manual classification approaches, we found a high concordance in the category "antisense" where 58/391 were called by the automatic versus 63/391 by the manual method (overlap of 51 probe sets). In the category "overlap", automatic classification gave more hits than manual inspection (120 versus 81 of 391 probe sets). Together, the manual and automatic classification called 653 probe sets on the MOE430 2.0 GeneChip as NAT-detecting probe sets, which are listed in Additional file 1 annotated with the gene symbol and the result of manual and automated validation. This annotation may be used to mine existing microarray data sets for expression and regulation of hundreds of NATs. Based on a rate of 57.5 % of confirmed NAT-detecting transcripts for MOE430 2.0 (640/1113 by the automatic classification), the human U113 Plus 2.0 GeneChip that has 3141 probe sets with a NetAffx annotation of "Negative Strand Matching Probes" can be expected to contain in the range of 1800 NAT-detecting probe sets.
While these large numbers of NAT-detecting probe sets should facilitate the quantitative study of expression and regulation of antisense transcripts, it has to be pointed out that only a fraction of all NATs can be detected by mining MOE430 or U133 expression analysis data sets: Because of the design of MOE430 microarrays, NAT-detecting probe sets are restricted to the 3' region of mRNAs and will usually not detect intronic non-coding RNAs. Since the labeling protocol used in most MOE430 data sets deposited in public databases is specific for poly-adenylated transcripts, all polyA- NATs are also not detected by the NAT-detecting probe sets confirmed in our study.
Confirmation of MOE430 2.0-detected antisense transcripts by strand-specific RT-PCR
Of the manually identified 144 probe sets detecting NATs on the MOE430 2.0 GeneChip, a significant fraction appears to be expressed based on present calls on at least 3 out of 12 arrays of 41 and 83 percent for "antisense" and "overlap" probe sets in a dataset obtained from macrophages stimulated in vitro (unpublished data). We next tested if the expression of NATs as indicated by these probe sets from the MOE430 2.0 GeneChip could be confirmed using strand-specific RT-PCR as an independent method.
Comparison of expression levels of sense-antisense pairs in mouse organs by RT-PCR
Correlation of SAS expression in multiple Affymetrix data sets
Taken together, we have shown here that Affymetrix MOE430 microarray data sets contain a considerable amount of quantitative information about the expression levels of hundreds of NATs. This information can be uncovered using the NetAffx annotation "Negative Strand Matching Probes" followed by manual or automated in silico validation in Ensembl, as we have done here for the MOE430 2.0 GeneChip, yielding the more than 600 NAT-detecting probe sets provided in Additional file 1. Experimental confirmation by strand-specific RT-PCR showed that most of the NATs indicated by microarray results were indeed expressed, interestingly at different levels between tissues.
Very recently, Werner et al. reported a similar approach to measure NATs by using the first version of the Affymetrix U74 mouse expression analysis GeneChips, that contained a large number of probe sets in reverse complementary orientation . They found a similar present call rate as reported by us for NAT-detecting probe sets on MOE430 arrays and also could confirm the expression of many SAS pairs by RT-PCR. While MOE430 and U133 GeneChips harbor a much lower percentage of "Negative Strand Matching Probes" probe sets than the U74 mouse set, the large and rapidly expanding number of publicly accessible datasets provides the unprecedented opportunity to extract hidden information about SAS transcript pairs expression as a by-product of gene expression profiling experiments.
Validation of the NetAffx probeset annotation "Negative Strand Match Probes" using the Ensembl genome browser
The manual and the automated probe set annotation are based on Ensembl  version 36 (released December 2005). Ensembl defines the probe set mapping to the genome and probe set association to genes as follows: each probe on the Affymetrix MOE430 2.0 array is mapped directly to the genome sequence. This mapping is not necessarily unique. Each probe set is then associated with one or more Ensembl genes by directly comparing the probe set to the set of cDNAs created from the Ensembl transcripts. Each transcript is associated with only one gene (although one gene may have many transcripts).
By manual inspection, the binding region was checked for the orientation of Ensembl known transcripts relative to the probe set. As the GeneChip microarray method analyzes cRNA, the orientations were evaluated as follows: If there is an annotated transcript in the same orientation and no annotated transcript in the reverse orientation, the probeset detects only sense transcripts (category "sense"). If there is an annotated transcript in the reverse orientation and no annotated transcript in the same orientation, the probe set detects only antisense transcripts (category "antisense"). If the probeset matches a genomic region where two annotated transcripts overlap, it detects transcripts that are sense and antisense (category "overlap").
In the automated classification, for a probe set to be considered "sense" to a gene, we required that the probe set be associated with only one gene and that the probes from the probe set are mapped to the genome on the same strand and overlapping the start of the gene by at least one base pair. "Antisense" probe sets map to the genome in a region where a gene is annotated on the opposite strand. Additionally, we require that "antisense" probe sets are not associated with any gene and have no mapped probes that overlap a gene on the same strand. Probes were classified as "cross-hybridizing" if they were associated with more than one gene and as "no transcript" if they are not associated with any gene and have no mapping that overlaps a gene by one base pair. Probes are classified as "overlap" if the genes they are associated with overlaps with a gene on the opposite strand by at least one base pair. This overlap is almost always associated with the untranslated regions.
RNA preparation and strand-specific reverse transcription PCR (RT-PCR)
Mouse organs were homogenized in peqGOLD TriFast (Peqlab) with an rotor-stator homogenizer, followed by extraction of RNA according to the manufacturer's instructions. Strand specific reverse transcription was performed using either the 5' or the 3' primer designed for a specific gene. The reaction was carried out in 10 μl on 100 ng RNA for one hour at 45°C using Superscript II (Invitrogen) and at 55°C using Superscript III. Subsequent PCR was performed with 5 μl of the RT reaction and the complemented primer (5' or 3') for 30 cycles for endpoint agarose gel analysis. Real-time qPCR using SYBRgreen was done for 40 cycles on an Applied Biosystems 7700 SDS.
natural antisense transcript
Work at the Microarray and Bioinformatics Core Unit in the Institute of Medical Microbiology, Immunology and Hygiene is funded by NGFN-2 grant FKZ 01GS0402 to RL, H. Wagner and R. Hoffmann. We thank Dr. Anne-Laure Boulesteix for statistical advice.
- Barrell BG, Air GM, Hutchison CA: Overlapping genes in bacteriophage phiX174. Nature. 1976, 264 (5581): 34-41. 10.1038/264034a0.PubMedView ArticleGoogle Scholar
- Tomizawa J, Itoh T, Selzer G, Som T: Inhibition of ColE1 RNA primer formation by a plasmid-specified small RNA. Proc Natl Acad Sci U S A. 1981, 78 (3): 1421-1425. 10.1073/pnas.78.3.1421.PubMed CentralPubMedView ArticleGoogle Scholar
- Williams T, Fried M: A mouse locus at which transcription from both DNA strands produces mRNAs complementary at their 3' ends. Nature. 1986, 322 (6076): 275-279. 10.1038/322275a0.PubMedView ArticleGoogle Scholar
- Lapidot M, Pilpel Y: Genome-wide natural antisense transcription: coupling its regulation to its different regulatory mechanisms. EMBO Rep. 2006, 7 (12): 1216-1222. 10.1038/sj.embor.7400857.PubMed CentralPubMedView ArticleGoogle Scholar
- Lavorgna G, Dahary D, Lehner B, Sorek R, Sanderson CM, Casari G: In search of antisense. Trends Biochem Sci. 2004, 29 (2): 88-94. 10.1016/j.tibs.2003.12.002.PubMedView ArticleGoogle Scholar
- Havilio M, Levanon EY, Lerman G, Kupiec M, Eisenberg E: Evidence for abundant transcription of non-coding regions in the Saccharomyces cerevisiae genome. BMC Genomics. 2005, 6 (1): 93-10.1186/1471-2164-6-93.PubMed CentralPubMedView ArticleGoogle Scholar
- Jen CH, Michalopoulos I, Westhead DR, Meyer P: Natural antisense transcripts with coding capacity in Arabidopsis may have a regulatory role that is not linked to double-stranded RNA degradation. Genome Biol. 2005, 6 (6): R51-10.1186/gb-2005-6-6-r51.PubMed CentralPubMedView ArticleGoogle Scholar
- Osato N, Yamada H, Satoh K, Ooka H, Yamamoto M, Suzuki K, Kawai J, Carninci P, Ohtomo Y, Murakami K, Matsubara K, Kikuchi S, Hayashizaki Y: Antisense transcripts with rice full-length cDNAs. Genome Biol. 2003, 5 (1): R5-10.1186/gb-2003-5-1-r5.PubMed CentralPubMedView ArticleGoogle Scholar
- Kiyosawa H, Yamanaka I, Osato N, Kondo S, Hayashizaki Y: Antisense transcripts with FANTOM2 clone set and their implications for gene regulation. Genome Res. 2003, 13 (6B): 1324-1334. 10.1101/gr.982903.PubMed CentralPubMedView ArticleGoogle Scholar
- Okazaki Y, Furuno M, Kasukawa T, Adachi J, Bono H, Kondo S, Nikaido I, Osato N, Saito R, Suzuki H, Yamanaka I, Kiyosawa H, Yagi K, Tomaru Y, Hasegawa Y, Nogami A, Schonbach C, Gojobori T, Baldarelli R, Hill DP, Bult C, Hume DA, Quackenbush J, Schriml LM, Kanapin A, Matsuda H, Batalov S, Beisel KW, Blake JA, Bradt D, Brusic V, Chothia C, Corbani LE, Cousins S, Dalla E, Dragani TA, Fletcher CF, Forrest A, Frazer KS, Gaasterland T, Gariboldi M, Gissi C, Godzik A, Gough J, Grimmond S, Gustincich S, Hirokawa N, Jackson IJ, Jarvis ED, Kanai A, Kawaji H, Kawasawa Y, Kedzierski RM, King BL, Konagaya A, Kurochkin IV, Lee Y, Lenhard B, Lyons PA, Maglott DR, Maltais L, Marchionni L, McKenzie L, Miki H, Nagashima T, Numata K, Okido T, Pavan WJ, Pertea G, Pesole G, Petrovsky N, Pillai R, Pontius JU, Qi D, Ramachandran S, Ravasi T, Reed JC, Reed DJ, Reid J, Ring BZ, Ringwald M, Sandelin A, Schneider C, Semple CA, Setou M, Shimada K, Sultana R, Takenaka Y, Taylor MS, Teasdale RD, Tomita M, Verardo R, Wagner L, Wahlestedt C, Wang Y, Watanabe Y, Wells C, Wilming LG, Wynshaw-Boris A, Yanagisawa M, Yang I, Yang L, Yuan Z, Zavolan M, Zhu Y, Zimmer A, Carninci P, Hayatsu N, Hirozane-Kishikawa T, Konno H, Nakamura M, Sakazume N, Sato K, Shiraki T, Waki K, Kawai J, Aizawa K, Arakawa T, Fukuda S, Hara A, Hashizume W, Imotani K, Ishii Y, Itoh M, Kagawa I, Miyazaki A, Sakai K, Sasaki D, Shibata K, Shinagawa A, Yasunishi A, Yoshino M, Waterston R, Lander ES, Rogers J, Birney E, Hayashizaki Y: Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs. Nature. 2002, 420 (6915): 563-573. 10.1038/nature01266.PubMedView ArticleGoogle Scholar
- Shendure J, Church GM: Computational discovery of sense-antisense transcription in the human and mouse genomes. Genome Biol. 2002, 3 (9): RESEARCH0044-10.1186/gb-2002-3-9-research0044.PubMed CentralPubMedView ArticleGoogle Scholar
- Yelin R, Dahary D, Sorek R, Levanon EY, Goldstein O, Shoshan A, Diber A, Biton S, Tamir Y, Khosravi R, Nemzer S, Pinner E, Walach S, Bernstein J, Savitsky K, Rotman G: Widespread occurrence of antisense transcription in the human genome. Nat Biotechnol. 2003, 21 (4): 379-386. 10.1038/nbt808.PubMedView ArticleGoogle Scholar
- Chen J, Sun M, Kent WJ, Huang X, Xie H, Wang W, Zhou G, Shi RZ, Rowley JD: Over 20% of human transcripts might form sense-antisense pairs. Nucleic Acids Res. 2004, 32 (16): 4812-4820. 10.1093/nar/gkh818.PubMed CentralPubMedView ArticleGoogle Scholar
- Sun M, Hurst LD, Carmichael GG, Chen J: Evidence for variation in abundance of antisense transcripts between multicellular animals but no relationship between antisense transcriptionand organismic complexity. Genome Res. 2006, 16 (7): 922-933. 10.1101/gr.5210006.PubMed CentralPubMedView ArticleGoogle Scholar
- Kampa D, Cheng J, Kapranov P, Yamanaka M, Brubaker S, Cawley S, Drenkow J, Piccolboni A, Bekiranov S, Helt G, Tammana H, Gingeras TR: Novel RNAs identified from an in-depth analysis of the transcriptome of human chromosomes 21 and 22. Genome Res. 2004, 14 (3): 331-342. 10.1101/gr.2094104.PubMed CentralPubMedView ArticleGoogle Scholar
- Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, Nishida H, Yap CC, Suzuki M, Kawai J, Suzuki H, Carninci P, Hayashizaki Y, Wells C, Frith M, Ravasi T, Pang KC, Hallinan J, Mattick J, Hume DA, Lipovich L, Batalov S, Engstrom PG, Mizuno Y, Faghihi MA, Sandelin A, Chalk AM, Mottagui-Tabar S, Liang Z, Lenhard B, Wahlestedt C: Antisense transcription in the mammalian transcriptome. Science. 2005, 309 (5740): 1564-1566. 10.1126/science.1112009.PubMedView ArticleGoogle Scholar
- Werner A: Natural antisense transcripts. RNA Biol. 2005, 2 (2): 53-62.PubMedView ArticleGoogle Scholar
- Mattick JS, Makunin IV: Non-coding RNA. Hum Mol Genet. 2006, 15 Spec No 1: R17-29. 10.1093/hmg/ddl046.PubMedView ArticleGoogle Scholar
- Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu JK: Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell. 2005, 123 (7): 1279-1291. 10.1016/j.cell.2005.11.035.PubMed CentralPubMedView ArticleGoogle Scholar
- Mattick JS: The functional genomics of noncoding RNA. Science. 2005, 309 (5740): 1527-1528. 10.1126/science.1117806.PubMedView ArticleGoogle Scholar
- Willingham AT, Orth AP, Batalov S, Peters EC, Wen BG, Aza-Blanc P, Hogenesch JB, Schultz PG: A strategy for probing the function of noncoding RNAs finds a repressor of NFAT. Science. 2005, 309 (5740): 1570-1573. 10.1126/science.1115901.PubMedView ArticleGoogle Scholar
- Bertone P, Stolc V, Royce TE, Rozowsky JS, Urban AE, Zhu X, Rinn JL, Tongprasit W, Samanta M, Weissman S, Gerstein M, Snyder M: Global identification of human transcribed sequences with genome tiling arrays. Science. 2004, 306 (5705): 2242-2246. 10.1126/science.1103388.PubMedView ArticleGoogle Scholar
- Cheng J, Kapranov P, Drenkow J, Dike S, Brubaker S, Patel S, Long J, Stern D, Tammana H, Helt G, Sementchenko V, Piccolboni A, Bekiranov S, Bailey DK, Ganesh M, Ghosh S, Bell I, Gerhard DS, Gingeras TR: Transcriptional maps of 10 human chromosomes at 5-nucleotide resolution. Science. 2005, 308 (5725): 1149-1154. 10.1126/science.1108625.PubMedView ArticleGoogle Scholar
- Affymetrix. [http://www.affymetrix.com]
- Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Fitzgerald S, Fernandez-Banet J, Graf S, Haider S, Hammond M, Herrero J, Holland R, Howe K, Howe K, Johnson N, Kahari A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Melsopp C, Megy K, Meidl P, Ouverdin B, Parker A, Prlic A, Rice S, Rios D, Schuster M, Sealy I, Severin J, Slater G, Smedley D, Spudich G, Trevanion S, Vilella A, Vogel J, White S, Wood M, Cox T, Curwen V, Durbin R, Fernandez-Suarez XM, Flicek P, Kasprzyk A, Proctor G, Searle S, Smith J, Ureta-Vidal A, Birney E: Ensembl 2007. Nucleic Acids Res. 2007, 35 (Database issue): D610-7. 10.1093/nar/gkl996.PubMed CentralPubMedView ArticleGoogle Scholar
- Lin L, Fevery J, Hiem Yap S: A novel strand-specific RT-PCR for detection of hepatitis C virus negative-strand RNA (replicative intermediate): evidence of absence or very low level of HCV replication in peripheral blood mononuclear cells. J Virol Methods. 2002, 100 (1-2): 97-105. 10.1016/S0166-0934(01)00399-8.PubMedView ArticleGoogle Scholar
- Werner A, Schmutzler G, Carlile M, Miles CG, Peters H: Expression profiling of antisense transcripts on DNA arrays. Physiol Genomics. 2006Google Scholar