Computational discovery of sense-antisense transcription in the human and mouse genomes

RNA silencing in mammalian oocytes and early embryos

RNA silencing is a common term for homology-dependent silencing phenomena found in the majority of eukaryotic species. RNA silencing pathways share several conserved components. The common denominator of these pathways is the presence of specific, short (21-25 nt) RNA molecules generated from different double-stranded RNA substrates by a specific RNase III activity. Short RNA molecules serve as a template for sequence-specific effects including transcriptional silencing, mRNA degradation, and inhibition of translation. This review will discuss possible roles of RNA silencing pathways in mouse oocytes and early embryos as well as the use of RNA silencing for experimental inhibition of gene expression in this model system.

Antisense RNA and RNAi in protozoan parasites: working hard or hardly working?

The complex life cycles of many protozoan parasites require the ability to respond to environmental and developmental cues through regulated gene expression. Traditionally, parasitologists have investigated these mechanisms by identifying and characterizing proteins that are necessary for the regulated expression of the genetic material. Although often successful, it is clear that protein-mediated gene regulation is only part of a complex story in which RNA itself is endowed with regulatory functions. Herein, we review both the known and potential regulatory roles of two types of RNA pathways within protozoan parasites: the RNA interference pathway and natural antisense transcripts. A better understanding of the native role of these pathways will not only enhance our understanding of the biology of these organisms but also aid in the development of more robust tools for reverse genetic analysis in this post-genomic era.

Non-coding antisense transcription detected by conventional and single-stranded cDNA microarray

BACKGROUND

Recent studies revealed that many mammalian protein-coding genes also transcribe their complementary strands. This phenomenon raises questions regarding the validity of data obtained from double-stranded cDNA microarrays since hybridization to both strands may occur. Here, we wanted to analyze experimentally the incidence of antisense transcription in human cells and to estimate their influence on protein coding expression patterns obtained by double-stranded microarrays. Therefore, we profiled transcription of sense and antisense independently by using strand-specific cDNA microarrays.

RESULTS

Up to 88% of expressed protein coding loci displayed concurrent expression from the complementary strand. Antisense transcription is cell specific and showed a strong tendency to be positively correlated to the expression of the sense counterparts. Even if their expression is wide-spread, detected antisense signals seem to have a limited distorting effect on sense profiles obtained with double-stranded probes.

CONCLUSION

Antisense transcription in humans can be far more common than previously estimated. However, it has limited influence on expression profiles obtained with conventional cDNA probes. This can be explained by a biological phenomena and a bias of the technique: a) a co-ordinate sense and antisense expression variation and b) a bias for sense-hybridization to occur with more efficiency, presumably due to variable exonic overlap between antisense transcripts.

Antisense transcription in the human cytomegalovirus transcriptome

Human cytomegalovirus (HCMV) infections are prevalent in human populations and can cause serious diseases, especially in those with compromised or immature immune systems. The HCMV genome of 230 kb is among the largest of the herpesvirus genomes. Although the entire sequence of the laboratory-adapted AD169 strain of HCMV has been available for 18 years, the precise number of viral genes is still in question. We undertook an analysis of the HCMV transcriptome as an approach to enumerate and analyze the gene products of HCMV. Transcripts of HCMV-infected fibroblasts were isolated at different times after infection and used to generate cDNA libraries representing different temporal classes of viral genes. cDNA clones harboring viral sequences were selected and subjected to sequence analysis. Of the 604 clones analyzed, 45% were derived from genomic regions predicted to be noncoding. Additionally, at least 55% of the cDNA clones in this study were completely or partially antisense to known or predicted HCMV genes. The remarkable accumulation of antisense transcripts during infection suggests that currently available genomic maps based on open-reading-frame and other in silico analyses may drastically underestimate the true complexity of viral gene products. These findings also raise the possibility that aspects of both the HCMV life cycle and genome organization are influenced by antisense transcription. Correspondingly, virus-derived noncoding and antisense transcripts may shed light on HCMV pathogenesis and may represent a new class of targets for antiviral therapies.

Origin of phenotypes: genes and transcripts

While the concept of a gene has been helpful in defining the relationship of a portion of a genome to a phenotype, this traditional term may not be as useful as it once was. Currently, "gene" has come to refer principally to a genomic region producing a polyadenylated mRNA that encodes a protein. However, the recent emergence of a large collection of unannotated transcripts with apparently little protein coding capacity, collectively called transcripts of unknown function (TUFs), has begun to blur the physical boundaries and genomic organization of genic regions with noncoding transcripts often overlapping protein-coding genes on the same (sense) and opposite strand (antisense). Moreover, they are often located in intergenic regions, making the genic portions of the human genome an interleaved network of both annotated polyadenylated and nonpolyadenylated transcripts, including splice variants with novel 5' ends extending hundreds of kilobases. This complex transcriptional organization and other recently observed features of genomes argue for the reconsideration of the term "gene" and suggests that transcripts may be used to define the operational unit of a genome.

Distinct expression patterns of natural antisense transcripts in Arabidopsis

It has been shown that overlapping cis-natural antisense transcripts (cis-NATs) can form a regulatory circuit in which small RNAs derived from one transcript regulate stability of the other transcript, which manifests itself as anticorrelated expression. However, little is known about how widespread antagonistic expression of cis-NATs is. We have determined how frequently cis-NAT pairs, which make up 7.4% of annotated transcription units in the Arabidopsis (Arabidopsis thaliana) genome, show anticorrelated expression patterns. Indeed, global expression profiles of pairs of cis-NATs on average have significantly lower pairwise Pearson correlation coefficients than other pairs of neighboring genes whose transcripts do not overlap. However, anticorrelated expression that is greater than expected by chance is found in only a small number of cis-NAT pairs. The degree of anticorrelation does not depend on the length of the overlap or on the distance of the 5' ends of the transcripts. Consistent with earlier findings, cis-NATs do not exhibit an increased likelihood to give rise to small RNAs, as determined from available small RNA sequences and massively parallel signature sequencing tags. However, the overlapping regions of cis-NATs appeared to be enriched for small RNA loci compared to nonoverlapping regions. Furthermore, expression of cis-NATs was not disproportionately affected in various RNA-silencing mutants. Our results demonstrate that there is a trend toward anticorrelated expression of cis-NAT pairs in Arabidopsis, but currently available data do not produce a strong signature of small RNA-mediated silencing for this process.

Uncovering information on expression of natural antisense transcripts in Affymetrix MOE430 datasets

Background

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.

Results

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.

Conclusion

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.

Sense-antisense pairs in mammals: functional and evolutionary considerations

Analysis of a catalog of S-AS pairs in the human and mouse genomes revealed several putative roles for natural antisense transcripts and showed that some are artifacts of cDNA library construction.

Background

A significant number of genes in mammalian genomes are being found to have natural antisense transcripts (NATs). These sense-antisense (S-AS) pairs are believed to be involved in several cellular phenomena.

Results

Here, we generated a catalog of S-AS pairs occurring in the human and mouse genomes by analyzing different sources of expressed sequences available in the public domain plus 122 massively parallel signature sequencing (MPSS) libraries from a variety of human and mouse tissues. Using this dataset of almost 20,000 S-AS pairs in both genomes we investigated, in a computational and experimental way, several putative roles that have been assigned to NATs, including gene expression regulation. Furthermore, these global analyses allowed us to better dissect and propose new roles for NATs. Surprisingly, we found that a significant fraction of NATs are artifacts produced by genomic priming during cDNA library construction.

Conclusion

We propose an evolutionary and functional model in which alternative polyadenylation and retroposition account for the origin of a significant number of functional S-AS pairs in mammalian genomes.

Computational analysis reveals 43% antisense transcription in 1182 transcripts in mouse muscle

It is increasingly evident that there is a widespread antisense transcription in the human and other eukaryotic genomes. However, the concept of general antisense expression is rarely investigated. We retrieved and correlated the expression of sense and antisense sequences of 1182 mouse transcripts to assess the prevalence of antisense transcription. We contrasted 20 Affymetrix MGU74A version 1 mouse genome chips to 12 MGU74A version 2 chips. For these 1182 transcripts, the version 1 chips contained the antisense sequences of the transcripts presented on the version 2 chips. The original data was taken from the GEO database. As the Affymetrix data is semi quantitative, the relative levels of antisense partners were analysed. We detected antisense transcription with an overall magnitude of 43% relative to sense transcription in the investigated transcripts. The average MGU74A version 1 expression is shifted towards smaller expression values (MGU74A version 1: 525.1; version 2: 1219.1; t-test: p < 0.001). A direct correlation between sense and antisense expression values could not be observed. Genes with high inverse expression values may be correlated to the investigated condition: genes where sense/control and control/antisense ratios were above two may be included in the pathogenetic pathways associated with dystrophin deficiency. The ratio of sense to antisense transcription varied between different chromosomes. We conclude that antisense transcription is a common phenomenon in the mouse genome and may have indirect regulatory functions.

Potential pitfalls in the accuracy of analysis of natural sense-antisense RNA pairs by reverse transcription-PCR

Background

The ability to accurately measure patterns of gene expression is essential in studying gene function. The reverse transcription polymerase chain reaction (RT-PCR) has become the method of choice for the detection and measurement of RNA expression patterns in both cells and small quantities of tissue. Our previous results show that there is a significant production of primer-independent cDNA synthesis using a popular RNase H^- RT enzyme. A PCR product was amplified from RT reactions that were carried out without addition of RT-primer. This finding jeopardizes the accuracy of RT-PCR when analyzing RNA that is expressed in both orientations. Current literature findings suggest that naturally occurring antisense expression is widespread in the mammalian transcriptome and consists of both coding and non-coding regulatory RNA. The primary purpose of this present study was to investigate the occurrence of primer-independent cDNA synthesis and how it may influence the accuracy of detection of sense-antisense RNA pairs.

Results

Our findings on cellular RNA and in vitro synthesized RNA suggest that these products are likely the results of RNA self-priming to generate random cDNA products, which contributes to the loss of strand specificity. The use of RNase H⁺ RT enzyme and carrying the RT reaction at high temperature (50°C) greatly improved the strand specificity of the RT-PCR detection.

Conclusion

While RT PCR is a basic method used for the detection and quantification of RNA expression in cells, primer-independent cDNA synthesis can interfere with RT specificity, and may lead to misinterpretation of the results, especially when both sense and antisense RNA are expressed. For accurate interpretation of the results, it is essential to carry out the appropriate negative controls.

RNAs everywhere: genome-wide annotation of structured RNAs

Starting with the discovery of microRNAs and the advent of genome-wide transcriptomics, non-protein-coding transcripts have moved from a fringe topic to a central field research in molecular biology. In this contribution we review the state of the art of "computational RNomics", i.e., the bioinformatics approaches to genome-wide RNA annotation. Instead of rehashing results from recently published surveys in detail, we focus here on the open problem in the field, namely (functional) annotation of the plethora of putative RNAs. A series of exploratory studies are used to provide non-trivial examples for the discussion of some of the difficulties.

Transcriptional interferences in cis natural antisense transcripts of humans and mice

For a significant fraction of mRNAs, their expression is regulated by other RNAs, including cis natural antisense transcripts (cis-NATs) that are complementary mRNAs transcribed from opposite strands of DNA at the same genomic locus. The regulatory mechanism of mRNA expression by cis-NATs is unknown, although a few possible explanations have been proposed. To understand this regulatory mechanism, we conducted a large-scale analysis of the currently available data and examined how the overlapping arrangements of cis-NATs affect their expression level. Here, we show that for both human and mouse the expression level of cis-NATs decreases as the length of the overlapping region increases. In particular, the proportions of the highly expressed cis-NATs in all cis-NATs examined were approximately 36 and 47% for human and mouse, respectively, when the overlapping region was <200 bp. However, both proportions decreased to virtually zero when the overlapping regions were >2000 bp in length. Moreover, the distribution of the expression level of cis-NATs changes according to different types of the overlapping pattern of cis-NATs in the genome. These results are consistent with the transcriptional collision model for the regulatory mechanism of gene expression by cis-NATs.

Serial analysis of gene expression in sugarcane (Saccharum spp.) leaves revealed alternative C4 metabolism and putative antisense transcripts

Sugarcane (Saccharum spp.) is a highly efficient biomass and sugar producing crop. Leaf reactions have been considered as potential rate-limiting step for sucrose accumulation in sugarcane stalks. To characterize the sugarcane leaf transcriptome, field-grown mature leaves from cultivar "SP80-3280" were analyzed using Serial Analysis of Gene Expression (SAGE). From 480 sequenced clones, 9,482 valid tags were extracted, with 5,227 unique sequences, from which 3,659 (70%) matched at least a sugarcane assembled sequence (SAS) with putative function; while 872 tags (16.7%) matched SAS with unknown function; 523 (10%) matched SAS without a putative annotation; and only 173 (3.3%) did not match any sugarcane ESTs. Based on gene ontology (GO), photosystem (PS) I reaction center was identified as the most frequent gene product location, followed by the remaining sites of PS I, PS II and thylakoid complexes. For metabolic processes, photosynthesis light harvesting complexes; carbon fixation; and chlorophyll biosynthesis were the most enriched GO-terms. Considering the alternative photosynthetic C(4) cycles, tag frequencies related to phosphoenolpyruvate carboxykinase (PEPCK) and aspartate aminotransferase compared to those for NADP(+)-malic enzyme (NADP-ME) and NADP-malate dehydrogenase, suggested that PEPCK-type decarboxylation appeared to predominate over NADP-ME in mature leaves, although both may occur, opposite to currently assumed in sugarcane. From the unique tag set, 894 tags (17.1%) were assigned as potentially derived from antisense transcripts, while 73 tags (1.4%) were assigned to more than one SAS, suggesting the occurrence of alternative processing. The occurrence of antisense was validated by quantitative reverse transcription amplification. Sugarcane leaf transcriptome provided new insights for functional studies associated with sucrose synthesis and accumulation.

Localizing hotspots of antisense transcription

Analysis of the transcriptome by computational and experimental methods has established that sense-antisense transcriptional units are a common phenomenon. Although the regulatory potential of antisense transcripts has been experimentally verified in a number of studies, the biological importance of sense-antisense regulation of gene expression is still a matter of debate. Here, we report the identification of sequence features that are associated with antisense transcription. We show that the sequence composition of the first exon and the 5'end of the first intron of many human genes is similar to the sequence composition observed in promoter regions as measured by the density of known transcription regulatory motifs. Cloned intron-derived fragments were found to possess bidirectional promoter activity. In agreement with the reported abundance of antisense transcripts overlapping the 5'UTR, mapping of the 5'ends of antisense transcripts to the corresponding sense transcripts revealed that the first exon and the 5'end of the first intron are hotspots of antisense transcription as measured by the number of antisense transcription start sites per unit sequence. CpG dinucleotide suppression that is typically weak in non-methylated promoter regions is similarly weakened upstream as well as downstream of the first exon. In support of antisense transcripts playing a regulatory role, we find that 5'UTRs and first exons of genes with overlapping antisense transcripts are significantly longer than the genomic average. Interestingly, a similar size distribution of 5'UTRs and first exons is observed for genes silenced by CpG island methylation in human cancer.

Androgen responsive intronic non-coding RNAs

BACKGROUND

Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line.

RESULTS

Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5-5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation.

CONCLUSION

Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.

Eukaryotic regulatory RNAs: an answer to the 'genome complexity' conundrum

A large portion of the eukaryotic genome is transcribed as noncoding RNAs (ncRNAs). While once thought of primarily as "junk," recent studies indicate that a large number of these RNAs play central roles in regulating gene expression at multiple levels. The increasing diversity of ncRNAs identified in the eukaryotic genome suggests a critical nexus between the regulatory potential of ncRNAs and the complexity of genome organization. We provide an overview of recent advances in the identification and function of eukaryotic ncRNAs and the roles played by these RNAs in chromatin organization, gene expression, and disease etiology.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

BACKGROUND

RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known.

RESULTS

A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages.

CONCLUSION

We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

Expression profiling of antisense transcripts on DNA arrays

The majority of mouse genes are estimated to undergo bidirectional transcription; however, their tissue-specific distribution patterns and physiological significance are largely unknown. This is in part due to the lack of methodology to routinely assess the expression of natural antisense transcripts (NATs) on a large scale. Here we tested whether commercial DNA arrays can be used to monitor antisense transcription in mouse kidney and brain. We took advantage of the reversely annotated oligonucleotides on the U74 mouse genome array from Affymetrix that hybridize to NATs overlapping with the sense transcript in the area of the probe set. In RNA samples from mouse kidney and brain, 11.9% and 10.1%, respectively, of 5,652 potential NATs returned positive and about half of the antisense RNAs were detected in both tissues, which was similar to the fraction of sense transcripts expressed in both tissues. Notably, we found that the majority of NATs are related to the sense transcriptome since corresponding sense transcripts were detected for 92.5% (kidney) and 74.5% (brain) of the detected antisense RNAs. Antisense RNA transcription was confirmed by real-time PCR and included additional RNA samples from heart, thymus, and liver. The randomly selected transcripts showed tissue specific expression patterns and varying sense/antisense ratios. The results indicate that antisense transcriptomes are tissue specific, and although pairing of sense/antisense transcripts are known to result in rapid degradation, our data provide proof of principle that the sensitivity of commercial DNA arrays is sufficient to assess NATs in total RNA of whole organs.

Retinoic acid regulation of eye and testis-specific transcripts within a complex locus

We previously used a yeast-based enhancer trap to identify a strong, retinoic acid response element (RARE). We have now characterized testis and eye transcripts that are adjacent to this regulatory element. Bioinformatics analysis of expressed sequence tag (EST) clones and RNase protection, reverse transcription-PCR, and Northern blot assays indicate that these two RNAs are transcribed from the same locus on opposite template strands. This positions the RARE upstream of the testis transcript and downstream of the eye transcript. Additionally, these two RNAs are embedded within the third intron of the 329kbp gene that encodes the Zinc Finger and BTB domain containing 7C protein (Zbtb7C). We present evidence indicating that the testis transcript is expressed primarily in spermatocytes and/or early round spermatids. Furthermore, our analyses of transcript levels in eyes and testes isolated from vitamin A deficient mice or from mice with defects in retinoid storage or signaling indicate that retinoids are required for expression in vivo.

Natural antisense and noncoding RNA transcripts as potential drug targets

Information on the complexity of mammalian RNA transcription has increased greatly in the past few years. Notably, thousands of sense transcripts (conventional protein-coding genes) have antisense transcript partners, most of which are noncoding. Interestingly, a number of antisense transcripts regulate the expression of their sense partners, either in a discordant (antisense knockdown results in sense-transcript elevation) or concordant (antisense knockdown results in concomitant sense-transcript reduction) manner. Two new pharmacological strategies based on the knockdown of antisense RNA transcripts by siRNA (or another RNA targeting principle) are proposed in this review. In the case of discordant regulation, knockdown of antisense transcript elevates the expression of the conventional (sense) gene, thereby conceivably mimicking agonist-activator action. In the case of concordant regulation, knockdown of antisense transcript, or concomitant knockdown of antisense and sense transcripts, results in an additive or even synergistic reduction of the conventional gene expression. Although both strategies have been demonstrated to be valid in cell culture, it remains to be seen whether they provide advantages in other contexts.

Prediction of trans-antisense transcripts in Arabidopsis thaliana

BACKGROUND

Natural antisense transcripts (NATs) are coding or non-coding RNAs with sequence complementarity to other transcripts (sense transcripts). These RNAs could potentially regulate the expression of their sense partner(s) at either the transcriptional or post-transcriptional level. Experimental and computational methods have demonstrated the widespread occurrence of NATs in eukaryotes. However, most previous studies only focused on cis-NATs with little attention being paid to NATs that originate in trans.

RESULTS

We have performed a genome-wide screen of trans-NATs in Arabidopsis thaliana and identified 1,320 putative trans-NAT pairs. An RNA annealing program predicted that most trans-NATs could form extended double-stranded RNA duplexes with their sense partners. Among trans-NATs with available expression data, more than 85% were found in the same tissue as their sense partners; of these, 67% were found in the same cell as their sense partners at comparable expression levels. For about 60% of Arabidopsis trans-NATs, orthologs of at least one transcript of the pair also had trans-NAT partners in either Populus trichocarpa or Oryza sativa. The observation that 430 transcripts had both putative cis- and trans-NATs implicates multiple regulations by antisense transcripts. The potential roles of trans-NATs in inducing post-transcriptional gene silencing and in regulating alternative splicing were also examined.

CONCLUSION

The Arabidopsis transcriptome contains a fairly large number of trans-NATs, whose possible functions include silencing of the corresponding sense transcripts or altering their splicing patterns. The interlaced relationships observed in some cis- and trans-NAT pairs suggest that antisense transcripts could be involved in complex regulatory networks in eukaryotes.

A large quantity of novel human antisense transcripts detected by LongSAGE

MOTIVATION

Taking advantage of the high sensitivity and specificity of LongSAGE tag for transcript detection and genome mapping, we analyzed the 632 813 unique human LongSAGE tags deposited in public databases to identify novel human antisense transcripts.

RESULTS

Our study identified 45 321 tags that match the antisense strand of 9804 known mRNA sequences, 6606 of which contain antisense ESTs and 3198 are mapped only by SAGE tags. Quantitative analysis showed that the detected antisense transcripts are present at levels lower than their counterpart sense transcripts. Experimental results confirmed the presence of antisense transcripts detected by the antisense tags. We also constructed an antisense tag database that can be used to identify the antisense SAGE tags originated from the antisense strand of known mRNA sequences included in the RefSeq database.

CONCLUSIONS

Our study highlights the benefits of exploring SAGE data for comprehensive identification of human antisense transcripts and demonstrates the prevalence of antisense transcripts in the human genome.

Genome-wide in silico identification and analysis of cis natural antisense transcripts (cis-NATs) in ten species

We developed a fast, integrative pipeline to identify cis natural antisense transcripts (cis-NATs) at genome scale. The pipeline mapped mRNAs and ESTs in UniGene to genome sequences in GoldenPath to find overlapping transcripts and combining information from coding sequence, poly(A) signal, poly(A) tail and splicing sites to deduce transcription orientation. We identified cis-NATs in 10 eukaryotic species, including 7830 candidate sense–antisense (SA) genes in 3915 SA pairs in human. The abundance of SA genes is remarkably low in worm and does not seem to be caused by the prevalence of operons. Hundreds of SA pairs are conserved across different species, even maintaining the same overlapping patterns. The convergent SA class is prevalent in fly, worm and sea squirt, but not in human or mouse as reported previously. The percentage of SA genes among imprinted genes in human and mouse is 24–47%, a range between the two previous reports. There is significant shortage of SA genes on Chromosome X in human and mouse but not in fly or worm, supporting X-inactivation in mammals as a possible cause. SA genes are over-represented in the catalytic activities and basic metabolism functions. All candidate cis-NATs can be downloaded from .

The complexity of antisense transcription revealed by the study of developing male germ cells

Computational analyses have identified the widespread occurrence of antisense transcripts in the human and the mouse genome. However, the structure and the origin of the majority of the antisense transcripts are unknown. The presence of antisense transcripts for 19 of 64 differentially expressed genes during mouse spermatogenesis was demonstrated with orientation-specific RT-PCR. These antisense transcripts were derived from a wide variety of origins, including processed sense transcripts, intronic and exonic sequences of a single gene or multiple genes, intergenic sequences, and pseudogenes. They underwent normal and alternative splicing, 5' capping, and 3' polyadenylation, similar to the sense transcripts. There were also antisense transcripts that were not capped and/or polyadenylated. The testicular levels of the sense transcripts were higher than those of the antisense transcripts in all cases, while the relative expression in nontesticular tissues was variable. Thus antisense transcripts have complex origins and structures and the sense and antisense transcripts can be regulated independently.

The silence RNA keeps: cis mechanisms of RNA mediated epigenetic silencing in mammals

One of the fundamental questions of modern biology is to unravel how genes are switched on and off at the right time and in the correct tissues. It is well recognized that gene regulation depends on a dynamic balance between activating and repressing forces, and multiple mechanisms are involved in both gene silencing and activation. Work over the last decade has revealed that in some cases transcriptional silencing of specific genes is mediated by RNAs that specifically recruit repressing complexes to homologous DNA sequences. Examples of both cis and trans RNA driven transcriptional silencing have been reported. This review focuses on those examples of transcriptional gene silencing in which the RNA component seems to act uniquely in cis. Speculative models of how such cis acting transcripts may trigger transcriptional silencing are proposed. Future experimental testing of these and other mechanisms is important to gain a fuller understanding of how genes are regulated and to identify instances in which such mechanisms are defective, leading to disease. Understanding the basic molecular basis of these phenomena will provide us with invaluable tools for the future development of targeted therapies and drugs for those diseases in which they are faulty.

Accumulation of sense-antisense transcripts of the rice catalase gene CatB under dark conditions requires signals from shoots

The amount of mRNA of the Oryza sativa L., cv. Nipponbare (rice) catalase gene, CatB, was decreased in the roots of intact seedlings kept in continuous darkness (DD). In contrast, sense and antisense unspliced CatB transcripts accumulated in the same tissue. Both strands cover the entire CatB-coding region, and form double-stranded RNA (dsRNA). The results of RNA dot-blot hybridization analysis using low molecular weight RNAs suggested that the sense and antisense CatB transcripts were more stable under DD conditions than under a light-dark regimen (LD). After shifting the lighting conditions from DD to LD, both the sense and antisense CatB transcripts were hardly detected, and the amount of CatB mRNA was restored. From these results, the antisense CatB transcripts might play a role in suppressing the normal processing of sense CatB transcript and also CatB protein synthesis by dsRNA formation, under conditions unsuitable for plant growth such as DD. This study indicates that signals transmitted from shoots are associated with the accumulation of sense and antisense CatB transcripts in roots under DD conditions, and that auxin is one of the possible signals.

A transient heterochromatic state in Xist preempts X inactivation choice without RNA stabilization

X chromosome inactivation (XCI) depends on a noncoding sense-antisense transcript pair, Xist and Tsix. At the onset of XCI, Xist RNA accumulates on one of two Xs, coating and silencing the chromosome in cis. The molecular basis for monoallelic Xist upregulation is not known, though evidence predominantly supports a posttranscriptional mechanism through RNA stabilization. Here, we test whether Tsix RNA destabilizes Xist RNA. Unexpectedly, we find that Xist upregulation is not based on transcript stabilization at all but is instead controlled by transcription in a sex-specific manner. Tsix directly regulates its transcription. On the future inactive X, Tsix downregulation induces a transient heterochromatic state in Xist, followed paradoxically by high-level Xist expression. A Tsix-deficient X chromosome adopts the heterochromatic state in pre-XCI cells. This state persists through XCI establishment and "reverts" to a euchromatic state during XCI maintenance. We have therefore identified chromatin marks that preempt and predict asymmetric Xist expression.

Complex Loci in human and mouse genomes

Mammalian genomes harbor a larger than expected number of complex loci, in which multiple genes are coupled by shared transcribed regions in antisense orientation and/or by bidirectional core promoters. To determine the incidence, functional significance, and evolutionary context of mammalian complex loci, we identified and characterized 5,248 cis-antisense pairs, 1,638 bidirectional promoters, and 1,153 chains of multiple cis-antisense and/or bidirectionally promoted pairs from 36,606 mouse transcriptional units (TUs), along with 6,141 cis-antisense pairs, 2,113 bidirectional promoters, and 1,480 chains from 42,887 human TUs. In both human and mouse, 25% of TUs resided in cis-antisense pairs, only 17% of which were conserved between the two organisms, indicating frequent species specificity of antisense gene arrangements. A sampling approach indicated that over 40% of all TUs might actually be in cis-antisense pairs, and that only a minority of these arrangements are likely to be conserved between human and mouse. Bidirectional promoters were characterized by variable transcriptional start sites and an identifiable midpoint at which overall sequence composition changed strand and the direction of transcriptional initiation switched. In microarray data covering a wide range of mouse tissues, genes in cis-antisense and bidirectionally promoted arrangement showed a higher probability of being coordinately expressed than random pairs of genes. In a case study on homeotic loci, we observed extensive transcription of nonconserved sequences on the noncoding strand, implying that the presence rather than the sequence of these transcripts is of functional importance. Complex loci are ubiquitous, host numerous nonconserved gene structures and lineage-specific exonification events, and may have a cis-regulatory impact on the member genes.

Complete Sequence Assembly and Characterization of the C57BL/6 Mouse Ig Heavy Chain V Region

The mechanisms that regulate variable (V) gene selection during the development of the mouse IgH repertoire are not fully understood, due in part to the absence of the complete locus sequence. To better understand these processes, we have assembled the entire 2.5-Mb mouse IgH (Igh) V region sequence of the C57BL/6 strain from public sequences and present the first complete annotated map of the region, including V genes, pseudogenes, repeats, and nonrepetitive intergenic sequences. In so doing, we have discovered a new V gene family, VH16. We have identified clusters of conserved region-specific intergenic sequences and have verified our assembly by genic and intergenic Southern blotting. We have observed that V pseudogenes are not evenly spread throughout the V region, but rather cluster together. The largest J558 family, which spans more than half of the locus, has two strikingly different domains, which suggest points of evolutionary divergence or duplication. The 5' end contains widely spaced J558 genes interspersed with 3609 genes and is pseudogene poor. The 3' end contains closely spaced J558 genes, no 3609 genes, and is pseudogene rich. Each occupies a different branch of the phylogenetic tree. Detailed analysis of 500-bp upstream of all functional genes has revealed several conserved binding sites, general and B cell-specific, as well as key differences between families. This complete and definitive assembly of the mouse Igh V region will facilitate detailed study of promoter function and large-scale mechanisms associated with V(D)J recombination including locus contraction and antisense intergenic transcription.

In silico discovery of human natural antisense transcripts

Background

Several high-throughput searches for ppotential natural antisense transcripts (NATs) have been performed recently, but most of the reports were focused on cis type. A thorough in silico analysis of human transcripts will help expand our knowledge of NATs.

Results

We have identified 568 NATs from human RefSeq RNA sequences. Among them, 403 NATs are reported for the first time, and at least 157 novel NATs are trans type. According to the pairing region of a sense and antisense RNA pair, hNATs are divided into 6 classes, of which about 87% involve 5' or 3' UTR sequences, supporting the regulatory role of UTRs. Among a total of 535 NAT pairs related with splice variants, 77.4% (414/535) have their pairing regions affected or completely eliminated by alternative splicing, suggesting significant relationship of alternative splicing and antisense-directed regulation. The extensive occurrence of splice variants in hNATs and other multiple pairing patterns results in a one-to-many relationship, allowing the formation of complex regulation networks. Based on microarray data from Stanford Microarray Database, two hNAT pairs were found to display significant inverse expression patterns before and after insulin injection.

Conclusion

NATs might carry out more extensive and complex functions than previously thought. Combined with endogenous micro RNAs, hNATs could be regarded as a special group of transcripts contributing to the complex regulation networks.

Regulatory RNAs in mammals

Recent years have brought a dramatic change in our understanding of the role of ribonucleic acids (RNAs) within the cell. In addition to the already well-known classes of RNAs that take part in the transmission of genetic information from DNA to proteins, a new highly heterogeneous group of RNA molecules has emerged. The regulatory nonprotein-coding RNAs (npcRNAs) have been shown to be involved in modulation of gene expression on both the transcriptional and post-transcriptional level. They participate in mechanisms of chromatin modification, regulation of transcription factor activity, and influencing mRNA stability, processing, and translation. npcRNAs are key factors in genetic imprinting, dosage compensation of X-chromosome-linked genes, and many processes of differentiation and development.

Natural expression of immature Ucn antisense RNA in the rat brain. Evidence favoring bidirectional transcription of the Ucn gene locus

Recently, it has been shown the endogenous expression of an antisense urocortin (Ucn) transcript in the rat brain and other tissues. In the present work, by means of two complementary techniques, specific-strand RT-PCR and in situ hybridization, we showed the natural expression of a second novel antisense Ucn RNA of higher size. Specific-strand RT-PCR of total RNA, cloning and sequence analysis together with the different subcellular localization observed for both antisense Ucn RNAs indicated that this novel antisense Ucn transcript corresponded to the immature form of the previously described antisense Ucn RNA. Sequence analysis indicated that this immature antisense Ucn transcript uses non-consensus CT-AC splice sites, exactly complementary to its sense counterpart. The mature antisense Ucn transcript was also amplified after specific-strand RT-PCR of poly(A)-RNA, suggesting that the mature antisense Ucn transcript is polyadenylated. We also proved that the region complementary to the promoter of sense Ucn RNA, including the TATA box, is part of the antisense Ucn RNA. Finally, we showed that the region complementary to the 3'-end of Ucn mRNA behaves as a functional promoter for the transcription of antisense Ucn RNA. Thus, the results indicate that the 3'-ends of both sense and antisense Ucn RNAs are the only non-complementary sequences between them. In conclusion, the present findings suggest that the Ucn gene locus naturally undergoes bidirectional transcription yielding a sense and an antisense RNA expanding the spectrum of antisense RNAs originated from the same genomic loci to antisense transcripts that are spliced using these non-consensus CT-AC splice sites.

Abundant novel transcriptional units and unconventional gene pairs on human chromosome 22

Novel transcriptional units (TUs) are EST-supported transcribed features not corresponding to known genes. Unconventional gene pairs (UGPs) are pairs of genes and/or TUs sharing exon-to-exon cis-antisense overlaps or putative bidirectional promoters. Computational TU and UGP discovery followed by manual curation was performed in the entire published 34.9-Mb human chromosome 22 euchromatic sequence. Novel TUs (n = 517) were as abundant as known genes (n = 492) and typically did not have nonprimate DNA and protein homologies. One hundred seventy-one (33%) of TUs, but only 13 (3%) of genes, both lacked nonprimate conservation and localized to gaps in the human-mouse BLASTZ alignment. Novel TUs were richer in exonic primate-specific interspersed repetitive elements (P = 0.001) and were more likely to rely on splice junctions provided by them, than were known genes: 19% of spliced TUs, versus 5% of spliced genes, had a splice site within a primate-specific repeat. Hence, novel TUs and known genes may represent different portions of the transcriptome. Two hundred nine (21%) of chromosome 22 transcripts participated in 77 cis-antisense and 42 promoter-sharing UGPs. Transcripts involved simultaneously in both UGP types were more common than was expected (P = 0.01). UGPs were nonrandomly distributed along the sequence: 89 (75%) clustered in distinct regions, the sum of which equaled 4.4 Mb (<13% of the chromosome). Eighty (67%) of the UGPs possessed significant locus structure differences between primates and rodents. Since some TUs may be functional noncoding transcripts and since the cis-regulatory potential of UGPs is well recognized, TUs and UGPs specific to the primate lineage may contribute to the genomic basis for primate-specific phenotypes.

Primate-specific endogenous cis-antisense transcription in the human 5q31 protocadherin gene cluster

Protocadherins (PCDH), localized to synaptic junctions, contribute to the formation of neuronal networks during brain development; thus, it is speculated that protocadherins may play a role in evolution of neuronal complexity. While protocadherin genes are highly conserved in vertebrates, EST evidence from the locus suggests apparently species-specific cis-antisense transcripts. Novel cis-antisense transcripts, which partially overlap the PCDHalpha12 variable exon, PCDHbeta3 single-exon gene, and PCDHpsi5 unprocessed pseudogene in the human 5q31 PCDHalpha/beta/gamma gene cluster and which are coexpressed with sense-strand transcripts in fetal and adult brain, were identified computationally and validated by gene-specific strand-specific reverse transcriptase PCR (SSRTPCR) and sequencing. Absence of antisense transcripts arising from equivalent genomic locations in mouse indicates that the antisense transcripts originated in the primates after the primate-rodent divergence. Furthermore, not all expected orthologues of human sense and antisense PCDH transcripts were detected in rhesus macaque brain, implying that protocadherin expression patterns differ between primate species. RT followed by quantitative real-time PCR (QPCR) analysis of the three genes in the brain of all three species, and of the PCDHbeta15 gene paralogous to PCDHpsi5 in human and rhesus, revealed that the presence of antisense transcripts was significantly associated with lower sense expression levels across all orthologues. This inverse relationship, along with the pattern of sense and antisense coexpression in the brain, is consistent with a regulatory role for the primate-specific PCDH cis-antisense transcripts, which may represent recent evolutionary inventions modulating the activity of this conserved gene cluster.

Rapid and efficient cDNA library screening by self-ligation of inverse PCR products (SLIP)

cDNA cloning is a central technology in molecular biology. cDNA sequences are used to determine mRNA transcript structures, including splice junctions, open reading frames (ORFs) and 5′- and 3′-untranslated regions (UTRs). cDNA clones are valuable reagents for functional studies of genes and proteins. Expressed Sequence Tag (EST) sequencing is the method of choice for recovering cDNAs representing many of the transcripts encoded in a eukaryotic genome. However, EST sequencing samples a cDNA library at random, and it recovers transcripts with low expression levels inefficiently. We describe a PCR-based method for directed screening of plasmid cDNA libraries. We demonstrate its utility in a screen of libraries used in our Drosophila EST projects for 153 transcription factor genes that were not represented by full-length cDNA clones in our Drosophila Gene Collection. We recovered high-quality, full-length cDNAs for 72 genes and variously compromised clones for an additional 32 genes. The method can be used at any scale, from the isolation of cDNA clones for a particular gene of interest, to the improvement of large gene collections in model organisms and the human. Finally, we discuss the relative merits of directed cDNA library screening and RT–PCR approaches.

Evidence for a preferential targeting of 3'-UTRs by cis-encoded natural antisense transcripts

Although both the 5'- and 3'-untranslated regions (5'- and 3'-UTRs) of eukaryotic mRNAs may play a crucial role in posttranscriptional gene regulation, we observe that cis-encoded natural antisense RNAs have a striking preferential complementarity to the 3'-UTRs of their target genes in mammalian (human and mouse) genomes. A null neutral model, evoking differences in the rate of 3'-UTR and 5'-UTR extension, could potentially explain high rates of 3'-to-3' overlap compared with 5'-to-5' overlap. However, employing a simulation model we show that this null model probably cannot explain the finding that 3'-to-3' overlapping pairs have a much higher probability (>5 times) of conservation in both mouse and human genomes with the same overlapping pattern than do 5'-to-5' overlaps. Furthermore, it certainly cannot explain the finding that overlapping pairs seen in both genomes have a significantly higher probability of having co-expression and inverse expression (i.e. characteristic of sense-antisense regulation) than do overlapping pairs seen in only one of the two species. We infer that the function of many 3'-to-3' overlaps is indeed antisense regulation. These findings underscore the preference for, and conservation of, 3'-UTR-targeted antisense regulation, and the importance of 3'-UTRs in gene regulation.

An abundance of RNA regulators

The importance of small, noncoding RNAs that act as regulators of transcription, of RNA modification or stability, and of mRNA translation is becoming increasingly apparent. Here we discuss current knowledge of regulatory RNA function and review how the RNAs have been identified in a variety of organisms. Many of the regulatory RNAs act through base-pairing interactions with target RNAs. The base-pairing RNAs can be grouped into two general classes: those that are encoded on the opposite strand of their target RNAs such that they contain perfect complementarity with their targets, and those that are encoded at separate locations on the chromosome and have imperfect base-pairing potential with their targets. Other regulatory RNAs act by modifying protein activity, in some cases by mimicking the structures of other RNA or DNA molecules.

Retention and repression: fates of hyperedited RNAs in the nucleus

Double-stranded RNA (dsRNA) is often formed in the nuclei of mammalian cells, but in this compartment it does not induce the effects characteristic of cytoplasmic dsRNA. Rather, recent work has suggested that nuclear dsRNA is a target for the ADAR class of enzymes, which deaminate adenosines to inosines. Further, there are a number of distinct fates of such edited RNA, including nuclear retention and perhaps also gene silencing.

Open reading frames provide a rich pool of potential natural antisense transcripts in fungal genomes

Natural antisense transcripts are reported from all kingdoms of life and several recent reports of genomewide screens indicate that they are widely distributed. These transcripts seem to be involved in various biological functions and may govern the expression of their respective sense partner. Very little, however, is known about the degree of evolutionary conservation of antisense transcripts. Furthermore, none of the earlier analyses has studied whether antisense relationships are solely dual or involved in more complex relationships. Here we present a systematic screen for cis- and trans-located antisense transcripts based on open reading frames (ORFs) from five fungal species. The relative number of ORFs involved in antisense relationships varies greatly between the five species. In addition, other significant differences are found between the species, such as the mean length of the antisense region. The majority of trans-located antisense transcripts is found to be involved in complex relationships, resulting in highly connected networks. The analysis of the degree of evolutionary conservation of antisense transcripts shows that most antisense transcripts have no ortholog in any other species. An annotation of antisense transcripts based on Gene Ontology directs to common terms and shows that proteins of genes involved in antisense relationships preferentially localize to the nucleus with common functions in the regulation or maintenance of nucleic acids.

Discovering functional transcription-factor combinations in the human cell cycle

With the completion of full genome sequences and advancement in high-throughput technologies, in silico methods have been successfully used to integrate diverse data sources toward unraveling the combinatorial nature of transcriptional regulation. So far, almost all of these studies are restricted to lower eukaryotes such as budding yeast. We describe here a computational search for functional transcription-factor (TF) combinations using phylogenetically conserved sequences and microarray-based expression data. Taking into account both orientational and positional constraints, we investigated the overrepresentation of binding sites in the vicinity of one another and whether these combinations result in more coherent expression profiles. Without any prior biological knowledge, the search led to the discovery of several experimentally established TF associations, as well as some novel ones. In particular, we identified a regulatory module controlling cell cycle-dependent transcription of G2-M genes and expanded its functional generality. We also detected many homotypic combinations, supporting the importance of binding-site density in transcriptional regulation of higher eukaryotes.

Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes

We have conducted a comprehensive search for conserved elements in vertebrate genomes, using genome-wide multiple alignments of five vertebrate species (human, mouse, rat, chicken, and Fugu rubripes). Parallel searches have been performed with multiple alignments of four insect species (three species of Drosophila and Anopheles gambiae), two species of Caenorhabditis, and seven species of Saccharomyces. Conserved elements were identified with a computer program called phastCons, which is based on a two-state phylogenetic hidden Markov model (phylo-HMM). PhastCons works by fitting a phylo-HMM to the data by maximum likelihood, subject to constraints designed to calibrate the model across species groups, and then predicting conserved elements based on this model. The predicted elements cover roughly 3%-8% of the human genome (depending on the details of the calibration procedure) and substantially higher fractions of the more compact Drosophila melanogaster (37%-53%), Caenorhabditis elegans (18%-37%), and Saccharaomyces cerevisiae (47%-68%) genomes. From yeasts to vertebrates, in order of increasing genome size and general biological complexity, increasing fractions of conserved bases are found to lie outside of the exons of known protein-coding genes. In all groups, the most highly conserved elements (HCEs), by log-odds score, are hundreds or thousands of bases long. These elements share certain properties with ultraconserved elements, but they tend to be longer and less perfectly conserved, and they overlap genes of somewhat different functional categories. In vertebrates, HCEs are associated with the 3' UTRs of regulatory genes, stable gene deserts, and megabase-sized regions rich in moderately conserved noncoding sequences. Noncoding HCEs also show strong statistical evidence of an enrichment for RNA secondary structure.

Detecting false expression signals in high-density oligonucleotide arrays by an in silico approach

High-density oligonucleotide arrays have become a popular assay for concurrent measurement of mRNA expression at the genome scale. Much effort has been devoted to the development of statistical analysis tools aimed at reducing experimental noise and normalizing experimental variation in gene expression analysis. However, these investigations do not detect or catalog systematic problems associated with specific oligonucleotide probes. Here, we present an investigation of problematic probes that yield consistent but inaccurate signals across multiple experiments. By evaluating data integrity among gene, probe sequence, and genomic structure we identified a total of 20,696 (10.5%) nonspecific probes that could cross-hybridize to multiple genes and a total of 18,363 (9.3%) probes that miss the target transcript sequences on the Affymetrix GeneChip U95A/Av2 array. The numbers of nonspecific and mistargeted probes on the U133A array are 29,405 (12.1%) and 19,717 (8.0%), respectively. The poor performance of the mistargeted probes was confirmed in two GeneChip experiments, in which these probes showed a 20-30% decrease in detecting present signals compared with normal probes. Comparison of qualitative expression signals obtained from SAGE and EST data with those from GeneChip arrays showed that the consistency of the two platforms is 30% lower in problematic probes than in normal probes. A Web application was developed to apply our results for improving the accuracy of expression analysis.

Diversity of antisense regulation in eukaryotes: multiple mechanisms, emerging patterns

High-throughput analysis of RNA molecules in multicellular eukaryotes has revealed an abundance of complementary antisense RNAs that are transcribed from separate or overlapping genes. In mammals these include many novel non-coding RNAs of unknown function. This unexpected complexity of the mammalian transcriptome suggests that expression of many genes is regulated post-transcriptionally by mechanisms mediated by RNA-RNA base pairing. The recent discovery of the widespread expression of microRNAs in animals and plants provides a prototypic example of such regulation in eukaryotes. However, there are likely to be numerous other types of antisense regulation in eukaryotes, many as yet uncharacterized.

Multiple sequence alignments of partially coding nucleic acid sequences

BACKGROUND

High quality sequence alignments of RNA and DNA sequences are an important prerequisite for the comparative analysis of genomic sequence data. Nucleic acid sequences, however, exhibit a much larger sequence heterogeneity compared to their encoded protein sequences due to the redundancy of the genetic code. It is desirable, therefore, to make use of the amino acid sequence when aligning coding nucleic acid sequences. In many cases, however, only a part of the sequence of interest is translated. On the other hand, overlapping reading frames may encode multiple alternative proteins, possibly with intermittent non-coding parts. Examples are, in particular, RNA virus genomes.

RESULTS

The standard scoring scheme for nucleic acid alignments can be extended to incorporate simultaneously information on translation products in one or more reading frames. Here we present a multiple alignment tool, codaln, that implements a combined nucleic acid plus amino acid scoring model for pairwise and progressive multiple alignments that allows arbitrary weighting for almost all scoring parameters. Resource requirements of codaln are comparable with those of standard tools such as ClustalW.

CONCLUSION

We demonstrate the applicability of codaln to various biologically relevant types of sequences (bacteriophage Levivirus and Vertebrate Hox clusters) and show that the combination of nucleic acid and amino acid sequence information leads to improved alignments. These, in turn, increase the performance of analysis tools that depend strictly on good input alignments such as methods for detecting conserved RNA secondary structure elements.

Naturally occurring antisense RNA of histone H2a in mouse cultured cell lines

Background

An antisense transcript of histone H2a that has no significant protein-coding region has been cloned from a mouse full-length cDNA library. In the present study, we evaluated this transcript by using RT-PCR and compared the expression patterns of the sense and antisense transcripts by using quantitative RT-PCR (qRT-PCR).

Results

This antisense RNA was expressed in three mouse cell lines. We call it ASH2a. ASH2a includes not only the complementary sequence of the transcript of Hist2h2aa2 (a replication-dependent histone H2a gene), but also that of the promoter of Hist2h2aa2. The upstream genomic sequence of the transcription start site of the ASH2a-coding gene (ASH2a) lacks both CCAAT and TATA boxes. This absence suggests that the regulation of ASH2a is different from that of the replication-dependent histone H2a genes. Findings from qRT-PCR indicated that the expression pattern of ASH2a was different from that of Hist2h2aa2. Expression of Hist2h2aa2 peaked at 2 to 4 h during S-phase, but that of ASH2a peaked at 1 h.

Conclusion

We showed the existence of ASH2a, a histone H2a antisense RNA, in mouse cultured cells. The expression pattern of ASH2a is different from that of the sense RNA.

Natural antisense mRNAs to hyaluronan synthase 2 inhibit hyaluronan biosynthesis and cell proliferation

We report the identification of a natural antisense mRNA of hyaluronan synthase 2 that we have chosen to designate as HASNT (for HA synthase 2 antisense) in human and mouse. HASNT is transcribed from the opposite strand of the HAS2 gene locus and is represented by several independent expressed sequence tags in human. Portions of the mouse Hasnt gene were identified through an exon-trapping approach. Sequence conservation is extremely low between human and mouse HASNT, and it is not clear whether these mRNAs contain functional open reading frames. HASNT has an alternate splice site in both human and mouse. This splice site is located at an identical position within the gene in both species and results in mRNAs of two different lengths. In each species, the antisense portion of the HASNT gene is complementary to the first exon of HAS2, which represents the 5'-untranslated region. To study the biological activity of HASNT, two human expressed sequence tag clones, representing long and short HASNT splice variants, were cloned into a tetracycline-inducible vector and were stably transfected into human osteosarcoma U2-OS Tet-on cells. The long and short HASNT-expressing cells had a reduction in HAS2 mRNA levels up to 94 and 86%, respectively, whereas hyaluronan biosynthesis was inhibited by 40 and 37%, respectively. Cell proliferation was reduced throughout the time frame of the experiment. Exogenous high molecular mass hyaluronan failed to rescue the suppressed cell proliferation, whereas adenoviral-mediated overexpression of hyaluronan synthase 3, which stimulated endogenous hyaluronan biosynthesis, was able to rescue. Collectively, our data suggest that natural antisense mRNAs of HAS2 are able to regulate HAS2 mRNA levels and hyaluronan biosynthesis in a cell culture model system and may have an important and novel regulatory role in the control of HAS2, HA biosynthesis, and HA-dependent cell functions in vivo.

RNA polymerase II synthesizes antisense RNA in Plasmodium falciparum

The recent identification of antisense RNA in the transcriptomes of many eukaryotes has generated enormous interest. The presence of antisense RNA in Plasmodium falciparum, the causative agent of severe malaria, remains controversial. Elucidation of the mechanism of antisense RNA in P. falciparum synthesis is critical in order to demonstrate the origin and function of these transcripts. Therefore, a systematic analysis of antisense and sense RNA synthesis was performed using direct labeling experiments. Nuclear run on experiments with single-stranded DNA probes demonstrated that antisense RNA is synthesized in the nucleus at several genomic loci. Antisense RNA synthesis is sensitive to the potent RNA polymerase II inhibitor alpha-amanitin. Antisense and sense transcription was also detected in nuclei isolated from synchronized parasites, suggesting concurrent synthesis. In summary, our experiments directly demonstrate that antisense RNA synthesis is a common transcriptional phenomenon in P. falciparum, and is catalyzed by RNA polymerase II.

Genome-wide prediction and identification of cis-natural antisense transcripts in Arabidopsis thaliana

A new computational method for predicting cis-encoded natural antisense transcripts (NATs) in Arabidopsis identified 1,340 potential NAT pairs. The expression of both sense and antisense transcripts of 957 NAT pairs was confirmed, and analysis of MPSS data suggested that for most pairs one of the two transcripts is predominantly expressed in a tissue-specific manner.

Background

Natural antisense transcripts (NAT) are a class of endogenous coding or non-protein-coding RNAs with sequence complementarity to other transcripts. Several lines of evidence have shown that cis- and trans-NATs may participate in a broad range of gene regulatory events. Genome-wide identification of cis-NATs in human, mouse and rice has revealed their widespread occurrence in eukaryotes. However, little is known about cis-NATs in the model plant Arabidopsis thaliana.

Results

We developed a new computational method to predict and identify cis-encoded NATs in Arabidopsis and found 1,340 potential NAT pairs. The expression of both sense and antisense transcripts of 957 NAT pairs was confirmed using Arabidopsis full-length cDNAs and public massively parallel signature sequencing (MPSS) data. Three known or putative Arabidopsis imprinted genes have cis-antisense transcripts. Sequences and the genomic arrangement of two Arabidopsis NAT pairs are conserved in rice.

Conclusion

We combined information from full-length cDNAs and Arabidopsis genome annotation in our NAT prediction work and reported cis-NAT pairs that could not otherwise be identified by using one of the two datasets only. Analysis of MPSS data suggested that for most Arabidopsis cis-NAT pairs, there is predominant expression of one of the two transcripts in a tissue-specific manner.

Naturally occurring antisense: transcriptional leakage or real overlap?

Naturally occurring antisense transcription is associated with the regulation of gene expression through a variety of biological mechanisms. Several recent genome-wide studies reported the identification of potential antisense transcripts for thousands of mammalian genes, many of them resulting from alternatively polyadenylated transcripts or heterogeneous transcription start sites. However, it is not clear whether this transcriptional plasticity is intentional, leading to regulated overlap between the transcripts, or, alternatively, represents a "leakage" of the RNA transcription machinery. To address this question through an evolutionary approach, we compared the genomic organization of genes, with or without antisense, between human, mouse, and the pufferfish Fugu rubripes. Our hypothesis was that if two neighboring genes overlap and have a sense-antisense relationship, we would expect negative selection acting on the evolutionary separation between them. We found that antisense gene pairs are twice as likely to preserve their genomic organization throughout vertebrates' evolution compared to nonantisense pairs, implying an overlap existence in the ancestral genome. In addition, we show that increasing the genomic distance between pairs of genes having a sense-antisense relationship is selected against. These findings indicate that, at least in part, the abundance of antisense transcripts observed in expressed data represents real overlap rather than transcriptional leakage. Moreover, our results imply that natural antisense transcription has considerably affected vertebrate genome evolution.

Tsix transcription- versus RNA-based mechanisms in Xist repression and epigenetic choice

Recent inquiries have revealed a surprisingly large number (>2500) of naturally occurring antisense transcripts, but their function remains largely undiscovered. A better understanding of antisense mechanisms is clearly needed because of their potentially diverse roles in gene regulation and disease. A well-documented case occurs in X inactivation, the mechanism by which X-linked gene expression is equalized between XX females and XY males. The antisense gene Tsix determines X chromosome choice and represses the noncoding silencer, Xist. In principle, Tsix action may involve RNA, the act of transcription, or local chromatin. Here, we create novel Tsix alleles to distinguish transcription- versus RNA-based mechanisms. When Tsix transcription is terminated before Xist (TsixTRAP), Tsix cannot block Xist upregulation, suggesting the importance of overlapping antisense transcription. To separate the act of transcription from RNA, we knocked in Tsix cDNA in the reverse orientation (Tsix(cDNA)) to restore RNA levels in cis without concurrent transcription across Xist. However, Tsix(cDNA) cannot complement TsixTRAP. Surprisingly, both mutations disrupt choice, indicating that this epigenetic step requires transcription. We conclude that the processed antisense RNA does not act alone and that Tsix function specifically requires antiparallel transcription through Xist. A mechanism of transcription-based feedback regulation is proposed.