Evidence for variation in abundance of antisense transcripts between multicellular animals but no relationship between antisense transcriptionand organismic complexity

Impacts of the feedback loop between sense-antisense RNAs in regulating circadian rhythms

Antisense transcripts are a unique group of non-coding RNAs that are transcribed from the opposite strand of a sense coding gene in an antisense orientation. Even though they do not encode a protein, these transcripts play a regulatory role in a variety of biological processes, including circadian rhythms. We and others found an antisense transcript, Per2AS , that is transcribed from the strand opposite the sense transcript Period2 ( Per2 ) and exhibits a rhythmic and antiphasic expression pattern compared to Per2 in mouse. By assuming that Per2AS and Per2 mutually repress each other, our previous mathematical model predicted that Per2AS regulates the robustness and the amplitude of circadian rhythms. In this study, we revised our previous model and developed a new mathematical model that mechanistically described the mutually repressive relationship between Per2 and Per2AS via transcriptional interference. We found that the simulation results are largely consistent with experimental observations including the counterintuitive ones that could not be fully explained by our previous model. These results indicate that our revised model serves as a foundation to build more detailed models in the future to better understand the impact of Per2AS-Per2 interaction in the mammalian circadian clock. Our mechanistic description of Per2AS-Per2 interaction can also be extended to other mathematical models that involve sense-antisense RNA pairs that mutually repress each other.

Natural antisense transcripts as versatile regulators of gene expression

Long non-coding RNAs (lncRNAs) are emerging as a major class of gene products that have central roles in cell and developmental biology. Natural antisense transcripts (NATs) are an important subset of lncRNAs that are expressed from the opposite strand of protein-coding and non-coding genes and are a genome-wide phenomenon in both eukaryotes and prokaryotes. In eukaryotes, a myriad of NATs participate in regulatory pathways that affect expression of their cognate sense genes. Recent developments in the study of NATs and lncRNAs and large-scale sequencing and bioinformatics projects suggest that whether NATs regulate expression, splicing, stability or translation of the sense transcript is influenced by the pattern and degrees of overlap between the sense-antisense pair. Moreover, epigenetic gene regulatory mechanisms prevail in somatic cells whereas mechanisms dependent on the formation of double-stranded RNA intermediates are prevalent in germ cells. The modulating effects of NATs on sense transcript expression make NATs rational targets for therapeutic interventions.

Natural antisense transcript of Period2, Per2AS, regulates the amplitude of the mouse circadian clock

In mammals, a set of core clock genes form transcription-translation feedback loops to generate circadian oscillations. We and others recently identified a novel transcript at the Period2 (Per2) locus that is transcribed from the antisense strand of Per2 This transcript, Per2AS, is expressed rhythmically and antiphasic to Per2 mRNA, leading to our hypothesis that Per2AS and Per2 mutually inhibit each other's expression and form a double negative feedback loop. By perturbing the expression of Per2AS, we found that Per2AS transcription, but not transcript, represses Per2 However, Per2 does not repress Per2AS, as Per2 knockdown led to a decrease in the Per2AS level, indicating that Per2AS forms a single negative feedback loop with Per2 and maintains the level of Per2 within the oscillatory range. Per2AS also regulates the amplitude of the circadian clock, and this function cannot be solely explained through its interaction with Per2, as Per2 knockdown does not recapitulate the phenotypes of Per2AS perturbation. Overall, our data indicate that Per2AS is an important regulatory molecule in the mammalian circadian clock machinery. Our work also supports the idea that antisense transcripts of core clock genes constitute a common feature of circadian clocks, as they are found in other organisms.

bsAS, an antisense long non-coding RNA, essential for correct wing development through regulation of blistered/DSRF isoform usage

Natural Antisense Transcripts (NATs) are long non-coding RNAs (lncRNAs) that overlap coding genes in the opposite strand. NATs roles have been related to gene regulation through different mechanisms, including post-transcriptional RNA processing. With the aim to identify NATs with potential regulatory function during fly development, we generated RNA-Seq data in Drosophila developing tissues and found bsAS, one of the most highly expressed lncRNAs in the fly wing. bsAS is antisense to bs/DSRF, a gene involved in wing development and neural processes. bsAS plays a crucial role in the tissue specific regulation of the expression of the bs/DSRF isoforms. This regulation is essential for the correct determination of cell fate during Drosophila development, as bsAS knockouts show highly aberrant phenotypes. Regulation of bs isoform usage by bsAS is mediated by specific physical interactions between the promoters of these two genes, which suggests a regulatory mechanism involving the collision of RNA polymerases transcribing in opposite directions. Evolutionary analysis suggests that bsAS NAT emerged simultaneously to the long-short isoform structure of bs, preceding the emergence of wings in insects.

Genomic determinants for initiation and length of natural antisense transcripts in Entamoeba histolytica

Natural antisense transcripts (NAT) have been reported in prokaryotes and eukaryotes. While the functions of most reported NATs remain unknown, their potentials in regulating the transcription of their counterparts have been speculated. Entamoeba histolytica, which is a unicellular eukaryotic parasite, has a compact protein-coding genome with very short intronic and intergenic regions. The regulatory mechanisms of gene expression in this compact genome are under-described. In this study, by genome-wide mapping of RNA-Seq data in the genome of E. histolytica, we show that a substantial fraction of its protein-coding genes (28%) has significant transcription on their opposite strand (i.e. NAT). Intriguingly, we found the location of transcription start sites or polyadenylation sites of NAT are determined by the specific motifs encoded on the opposite strand of the gene coding sequences, thereby providing a compact regulatory system for gene transcription. Moreover, we demonstrated that NATs are globally up-regulated under various environmental conditions including temperature stress and pathogenicity. While NATs do not appear to be consequences of spurious transcription, they may play a role in regulating gene expression in E. histolytica, a hypothesis which needs to be tested.

First insights into the nature and evolution of antisense transcription in nematodes

BACKGROUND

The development of multicellular organisms is coordinated by various gene regulatory mechanisms that ensure correct spatio-temporal patterns of gene expression. Recently, the role of antisense transcription in gene regulation has moved into focus of research. To characterize genome-wide patterns of antisense transcription and to study their evolutionary conservation, we sequenced a strand-specific RNA-seq library of the nematode Pristionchus pacificus.

RESULTS

We identified 1112 antisense configurations of which the largest group represents 465 antisense transcripts (ASTs) that are fully embedded in introns of their host genes. We find that most ASTs show homology to protein-coding genes and are overrepresented in proteomic data. Together with the finding, that expression levels of ASTs and host genes are uncorrelated, this indicates that most ASTs in P. pacificus do not represent non-coding RNAs and do not exhibit regulatory functions on their host genes. We studied the evolution of antisense gene pairs across 20 nematode genomes, showing that the majority of pairs is lineage-specific and even the highly conserved vps-4, ddx-27, and sel-2 loci show abundant structural changes including duplications, deletions, intron gains and loss of antisense transcription. In contrast, host genes in general, are remarkably conserved and encode exceptionally long introns leading to unusually large blocks of conserved synteny.

CONCLUSIONS

Our study has shown that in P. pacificus antisense transcription as such does not define non-coding RNAs but is rather a feature of highly conserved genes with long introns. We hypothesize that the presence of regulatory elements imposes evolutionary constraint on the intron length, but simultaneously, their large size makes them a likely target for translocation of genomic elements including protein-coding genes that eventually end up as ASTs.

IAOseq: inferring abundance of overlapping genes using RNA-seq data

BACKGROUND

Overlapping transcription constitutes a common mechanism for regulating gene expression. A major limitation of the overlapping transcription assays is the lack of high throughput expression data.

RESULTS

We developed a new tool (IAOseq) that is based on reads distributions along the transcribed regions to identify the expression levels of overlapping genes from standard RNA-seq data. Compared with five commonly used quantification methods, IAOseq showed better performance in the estimation accuracy of overlapping transcription levels. For the same strand overlapping transcription, currently existing high-throughput methods are rarely available to distinguish which strand was present in the original mRNA template. The IAOseq results showed that the commonly used methods gave an average of 1.6 fold overestimation of the expression levels of same strand overlapping genes.

CONCLUSIONS

This work provides a useful tool for mining overlapping transcription levels from standard RNA-seq libraries. IAOseq could be used to help us understand the complex regulatory mechanism mediated by overlapping transcripts. IAOseq is freely available at http://lifecenter.sgst.cn/main/en/IAO_seq.jsp.

Development of genodynamic metrics for exploring the biophysics of DNA polymorphisms

Single nucleotide polymorphisms (SNPs) represent an important type of dynamic sites within the human genome. These common variants often locally correlate within more complex multi-SNP haploblocks that are maintained throughout generations in a stable population. Information encoded in the structure of SNPs and SNP haploblock variation can be characterized through a normalized information content metric. Genodynamics is being developed as the analogous "thermodynamics" characterizing the state variables for genomic populations that are stable under stochastic environmental stresses. Since living systems have not been found to develop in the absence of environmental influences, this paper describes the analogous genomic free energy metrics in a given environment. SNP haploblocks were constructed by Haploview v4.2 for five chromosomes from phase III HapMap data, and the genomic state variables for each chromosome were calculated. An in silico analysis was performed on SNP haploblocks with the lowest genomic energy measures. Highly favorable genomic energy measures were found to correlate with highly conserved SNP haploblocks. Moreover, the most conserved haploblocks were associated with an evolutionarily conserved regulatory element and domain.

Examining the condition-specific antisense transcription in S. cerevisiae and S. paradoxus

Background

Recent studies have demonstrated that antisense transcription is pervasive in budding yeasts and is conserved between Saccharomyces cerevisiae and S. paradoxus. While studies have examined antisense transcripts of S. cerevisiae for inverse expression in stationary phase and stress conditions, there is a lack of comprehensive analysis of the conditional specific evolutionary characteristics of antisense transcription between yeasts. Here we attempt to decipher the evolutionary relationship of antisense transcription of S. cerevisiae and S. paradoxus cultured in mid log, early stationary phase, and heat shock conditions.

Results

Massively parallel sequencing of sequence strand-specific cDNA library was performed from RNA isolated from S. cerevisiae and S. paradoxus cells at mid log, stationary phase and heat shock conditions. We performed this analysis using a stringent set of sense ORF transcripts and non-coding antisense transcripts that were expressed in all the three conditions, as well as in both species. We found the divergence of the condition-specific anti-sense transcription levels is higher than that in condition-specific sense transcription levels, suggesting that antisense transcription played a potential role in adapting to different conditions. Furthermore, 43% of sense-antisense pairs demonstrated inverse expression in either stationary phase or heat shock conditions relative to the mid log conditions. In addition, a large part of sense-antisense pairs (67%), which demonstrated inverse expression, were highly conserved between the two species. Our results were also concordant with known functional analyses from previous studies and with the evidence from mechanistic experiments of role of individual genes.

Conclusions

By performing a genome-scale computational analysis, we have tried to evaluate the role of antisense transcription in mediating sense transcription under different environmental conditions across and in two related yeast species. Our findings suggest that antisense regulation could control expression of the corresponding sense transcript via inverse expression under a range of different circumstances.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-521) contains supplementary material, which is available to authorized users.

Sequence signatures of genes with accompanying antisense transcripts in Saccharomyces cerevisiae

Recent studies have found many antisense non-coding transcripts at the opposite strand of some protein-coding genes. In yeast, it was reported that such antisense transcripts play regulatory roles for their partner genes by forming a feedback loop with the protein-coding genes. Since not all coding genes have accompanying antisense transcripts, it would be interesting to know whether there are sequence signatures in a coding gene that are decisive or associated with the existence of such antisense partners. We collected all the annotated antisense transcripts in the yeast Saccharomyces cerevisiae, analyzed sequence motifs around the genes with antisense partners, and classified genes with and without accompanying antisense transcripts by using machine learning methods. Some weak but statistically significant sequence features are detected, which indicates that there are sequence signatures around the protein-coding genes that may be decisive or indicative for the existence of accompanying antisense transcripts.

Emerging epigenetic mechanisms of long non-coding RNAs

Long non-coding RNAs (lncRNAs) have been increasingly appreciated as an integral component of gene regulatory networks. Genome-wide features of their origin and expression patterns ascribed a prominent role for lncRNAs to the regulation of protein-coding genes, and also suggest a potential link to many human diseases. Recent studies have begun to unravel the intricate regulatory mechanism of lncRNAs occurring at multiple levels. The brain is one of the richest sources of lncRNAs, many of which have already shown a close relationship with genes or genetic loci implicated in a wide range of neurological disorders. This review describes recently emerging mechanistic principles of lncRNA functions to provide neuroscientists with molecular insights that will help future research on lncRNAs in the brain.

Integrated detection of natural antisense transcripts using strand-specific RNA sequencing data

Pairs of RNA molecules transcribed from partially or entirely complementary loci are called cis-natural antisense transcripts (cis-NATs), and they play key roles in the regulation of gene expression in many organisms. A promising experimental tool for profiling sense and antisense transcription is strand-specific RNA sequencing (ssRNA-seq). To identify cis-NATs using ssRNA-seq, we developed a new computational method based on a model comparison framework that incorporates the inherent variable efficiency of generating perfectly strand-specific libraries. Applying the method to new ssRNA-seq data from whole-root and cell-type–specific Arabidopsis libraries confirmed most of the known cis-NAT pairs and identified 918 additional cis-NAT pairs. Newly identified cis-NAT pairs are supported by polyadenylation data, alternative splicing patterns, and RT-PCR validation. We found 209 cis-NAT pairs that have opposite expression levels in neighboring cell types, implying cell-type–specific roles for cis-NATs. By integrating a genome-wide epigenetic profile of Arabidopsis, we identified a unique chromatin signature of cis-NATs, suggesting a connection between cis-NAT transcription and chromatin modification in plants. An analysis of small-RNA sequencing data showed that ∼4% of cis-NAT pairs produce putative cis-NAT–induced siRNAs. Taken together, our data and analyses illustrate the potential for multifaceted regulatory roles of plant cis-NATs.

Mechanisms of small RNA generation from cis-NATs in response to environmental and developmental cues

A large proportion of eukaryotic genomes is transcribed from both positive and negative strands of DNA and thus may generate overlapping sense and antisense transcripts. Some of these so-called natural antisense transcripts (NATs) are possibly co-regulated. When the overlapping sense and antisense transcripts are expressed at the same time in the same cell in response to various developmental and environmental cues; they may form double-stranded RNAs, which could be recognized by the small RNA biogenesis machinery and processed into small interfering RNAs (siRNAs). cis-NAT-derived siRNAs (nat-siRNAs) are present in plants, animals, and fungi. In plants, the presence of nat-siRNAs is supported not only by Northern blot and genetic analyses, but also by the fact that there is an overall sixfold enrichment of siRNAs in the overlapping regions of cis-NATs and 19%-29% of the siRNA-generating cis-NATs in plants give rise to siRNAs only in their overlapping regions. Silencing mediated by nat-siRNAs is one of the mechanisms for regulating the expression of the cis-NATs. This review focuses on challenging issues related to the biogenesis mechanisms as well as regulation and detection of nat-siRNAs. The advantages and limitations of new technologies for detecting cis-NATs, including direct RNA sequencing and strand-specific RNA sequencing, are also discussed.

A genome-wide investigation of expression characteristics of natural antisense transcripts in liver and muscle samples of pigs

Natural antisense transcripts are endogenous transcripts that are complementary to the sense-strand of DNA. These transcripts have been identified in various eukaryotic species and are involved in a broad range of regulatory events and biological processes. However, their general biological functions, expression characteristics and regulatory mechanisms are still unclear. In this study, 497 liver and 586 muscle samples were harvested from a White Duroc×Erhualian F₂ resource population. The expression profiles of sense and antisense transcripts were determined by tag-based RNA sequencing. We identified 33.7% and 20.4% of transcripts having both sense and antisense expression, and 12.5% and 6.1% of transcripts only expressing antisense transcripts in liver and muscle, respectively. More than 32.2% of imprinting or predicted imprinting genes in the geneimprint database were detected with both sense and antisense expression. The correlations between sense and antisense expression in sense-antisense pairs were diverse in both liver and muscle, showing positive, negative or absent correlation. Antisense expression increases gene expression variability. More interestingly, compared to eQTL mapping of sense transcripts in which more than one eQTL was mapped for a transcript, only one eQTL was identified for each antisense transcript, and the percentage of cis-eQTL in antisense eQTL was higher than that in sense eQTL. This suggests that the expressions of antisense transcripts tend to be cis-regulated by a single genomic locus. To our knowledge, this study is the first systematical investigation of antisense transcription in pigs. The findings improve our understanding of the complexity of porcine transcriptome.

Transcriptional interference networks coordinate the expression of functionally related genes clustered in the same genomic loci

The regulation of gene expression is essential for normal functioning of biological systems in every form of life. Gene expression is primarily controlled at the level of transcription, especially at the phase of initiation. Non-coding RNAs are one of the major players at every level of genetic regulation, including the control of chromatin organization, transcription, various post-transcriptional processes, and translation. In this study, the Transcriptional Interference Network (TIN) hypothesis was put forward in an attempt to explain the global expression of antisense RNAs and the overall occurrence of tandem gene clusters in the genomes of various biological systems ranging from viruses to mammalian cells. The TIN hypothesis suggests the existence of a novel layer of genetic regulation, based on the interactions between the transcriptional machineries of neighboring genes at their overlapping regions, which are assumed to play a fundamental role in coordinating gene expression within a cluster of functionally linked genes. It is claimed that the transcriptional overlaps between adjacent genes are much more widespread in genomes than is thought today. The Waterfall model of the TIN hypothesis postulates a unidirectional effect of upstream genes on the transcription of downstream genes within a cluster of tandemly arrayed genes, while the Seesaw model proposes a mutual interdependence of gene expression between the oppositely oriented genes. The TIN represents an auto-regulatory system with an exquisitely timed and highly synchronized cascade of gene expression in functionally linked genes located in close physical proximity to each other. In this study, we focused on herpesviruses. The reason for this lies in the compressed nature of viral genes, which allows a tight regulation and an easier investigation of the transcriptional interactions between genes. However, I believe that the same or similar principles can be applied to cellular organisms too.

Multiple insert size paired-end sequencing for deconvolution of complex transcriptomes

Deep sequencing of transcriptomes allows quantitative and qualitative analysis of many RNA species in a sample, with parallel comparison of expression levels, splicing variants, natural antisense transcripts, RNA editing and transcriptional start and stop sites the ideal goal. By computational modeling, we show how libraries of multiple insert sizes combined with strand-specific, paired-end (SS-PE) sequencing can increase the information gained on alternative splicing, especially in higher eukaryotes. Despite the benefits of gaining SS-PE data with paired ends of varying distance, the standard Illumina protocol allows only non-strand-specific, paired-end sequencing with a single insert size. Here, we modify the Illumina RNA ligation protocol to allow SS-PE sequencing by using a custom pre-adenylated 3' adaptor. We generate parallel libraries with differing insert sizes to aid deconvolution of alternative splicing events and to characterize the extent and distribution of natural antisense transcription in C. elegans. Despite stringent requirements for detection of alternative splicing, our data increases the number of intron retention and exon skipping events annotated in the Wormbase genome annotations by 127% and 121%, respectively. We show that parallel libraries with a range of insert sizes increase transcriptomic information gained by sequencing and that by current established benchmarks our protocol gives competitive results with respect to library quality.

Regulation of chromatin structure by long noncoding RNAs: focus on natural antisense transcripts

In the decade following the publication of the Human Genome, noncoding RNAs (ncRNAs) have reshaped our understanding of the broad landscape of genome regulation. During this period, natural antisense transcripts (NATs), which are transcribed from the opposite strand of either protein or non-protein coding genes, have vaulted to prominence. Recent findings have shown that NATs can exert their regulatory functions by acting as epigenetic regulators of gene expression and chromatin remodeling. Here, we review recent work on the mechanisms of epigenetic modifications by NATs and their emerging role as master regulators of chromatin states. Unlike other long ncRNAs, antisense RNAs usually regulate their counterpart sense mRNA in cis by bridging epigenetic effectors and regulatory complexes at specific genomic loci. Understanding the broad range of effects of NATs will shed light on the complex mechanisms that regulate chromatin remodeling and gene expression in development and disease.

The prevalence and regulation of antisense transcripts in Schizosaccharomyces pombe

A strand-specific transcriptome sequencing strategy, directional ligation sequencing or DeLi-seq, was employed to profile antisense transcriptome of Schizosaccharomyces pombe. Under both normal and heat shock conditions, we found that polyadenylated antisense transcripts are broadly expressed while distinct expression patterns were observed for protein-coding and non-coding loci. Dominant antisense expression is enriched in protein-coding genes involved in meiosis or stress response pathways. Detailed analyses further suggest that antisense transcripts are independently regulated with respect to their sense transcripts, and diverse mechanisms might be potentially involved in the biogenesis and degradation of antisense RNAs. Taken together, antisense transcription may have profound impacts on global gene regulation in S. pombe.

Intronic AT skew is a defendable proxy for germline transcription but does not predict crossing-over or protein evolution rates in Drosophila melanogaster

Recent evidence suggests that germline transcription may affect both protein evolutionary rates, possibly mediated by repair processes, and recombination rates, possibly mediated by chromatin and epigenetic modification. Here, we test these propositions in Drosophila melanogaster. The challenge for such analyses is to provide defendable measures of germline gene expression. Intronic AT skew is a good candidate measure as it is thought to be a consequence, at least in part, of transcription-coupled repair. Prior evidence suggests that intronic AT skew in D. melanogaster is not affected by proximity to intron extremities and differs between transcribed DNA and flanking sequence. We now also establish that intronic AT skew is a defendable proxy for germline expression as (a) it is more similar than expected by chance between introns of the same gene (which is not accounted for by physical proximity), (b) is correlated with male germline expression, and (c) is more pronounced in broadly expressed genes. Furthermore, (d) a trend for intronic skew to differ between 3' and 5' ends of genes is particular to broadly expressed genes. Finally, (e) controlling for physical distance, introns of proximate genes are most different in skew if they have different tissue specificity. We find that intronic AT skew, employed as a proxy for germline transcription, correlates neither with recombination rates nor with the rate of protein evolution. We conclude that there is no prima facie evidence that germline expression modulates recombination rates or monotonically affects protein evolution rates in D. melanogaster.

The elusive concept of the gene

In recent years geneticists have witnessed many significant observations which have seriously shaken the traditional concept of the gene. These specifically include the facts that (1) the boundaries of transcriptional units are far from clear; in fact, whole chromosomes if not the whole genome seem to be continuums of genetic transcription, (2) many examples of gene fusion are known, (3) likewise many examples of so-called encrypted genes are known in the organelle genomes of microbial eukaryotes and in prokaryotes, and (4) in addition to the structure of the gene, its functional status can also be inheritable, and, further, (5) epigenetic extra-genomic modes of inheritance, called genetic restoration, seem to be a rather common phenomenon, meaning that organisms can sometimes rewrite their DNA on the basis of RNA messages inherited from generations past. I will briefly review these observations and discuss the difficulties of defining the gene, and then formulate a new view, which is called the relational or systemic concept of the gene. It has to be noted that genes assume their information content characteristics in the Shannonian sense as nucleotide sequences of DNA (or RNA). However, on the basis of this we cannot say anything about their information content in the semantic sense. The semantic information content of genes is context-dependent. Genes namely assume their biochemical characteristics usually only within living cells, their developmental characteristics only within living organisms, and their evolutionary characteristics only within populations of living organisms.

Consistent deregulation of gene expression between human and murine MLL rearrangement leukemias

Important biological and pathologic properties are often conserved across species. Although several mouse leukemia models have been well established, the genes deregulated in both human and murine leukemia cells have not been studied systematically. We performed a serial analysis of gene expression in both human and murine MLL-ELL or MLL-ENL leukemia cells and identified 88 genes that seemed to be significantly deregulated in both types of leukemia cells, including 57 genes not reported previously as being deregulated in MLL-associated leukemias. These changes were validated by quantitative PCR. The most up-regulated genes include several HOX genes (e.g., HOX A5, HOXA9, and HOXA10) and MEIS1, which are the typical hallmark of MLL rearrangement leukemia. The most down-regulated genes include LTF, LCN2, MMP9, S100A8, S100A9, PADI4, TGFBI, and CYBB. Notably, the up-regulated genes are enriched in gene ontology terms, such as gene expression and transcription, whereas the down-regulated genes are enriched in signal transduction and apoptosis. We showed that the CpG islands of the down-regulated genes are hypermethylated. We also showed that seven individual microRNAs (miRNA) from the mir-17-92 cluster, which are overexpressed in human MLL rearrangement leukemias, are also consistently overexpressed in mouse MLL rearrangement leukemia cells. Nineteen possible targets of these miRNAs were identified, and two of them (i.e., APP and RASSF2) were confirmed further by luciferase reporter and mutagenesis assays. The identification and validation of consistent changes of gene expression in human and murine MLL rearrangement leukemias provide important insights into the genetic base for MLL-associated leukemogenesis.

Bidirectional transcription directs both transcriptional gene activation and suppression in human cells

Small RNAs targeted to gene promoters in human cells have been shown to modulate both transcriptional gene suppression and activation. However, the mechanism involved in transcriptional activation has remained poorly defined, and an endogenous RNA trigger for transcriptional gene silencing has yet to be identified. Described here is an explanation for siRNA-directed transcriptional gene activation, as well as a role for non-coding antisense RNAs as effector molecules driving transcriptional gene silencing. Transcriptional activation of p21 gene expression was determined to be the result of Argonaute 2–dependent, post-transcriptional silencing of a p21-specific antisense transcript, which functions in Argonaute 1–mediated transcriptional control of p21 mRNA expression. The data presented here suggest that in human cells, bidirectional transcription is an endogenous gene regulatory mechanism whereby an antisense RNA directs epigenetic regulatory complexes to a sense promoter, resulting in RNA-directed epigenetic gene regulation. The observations presented here support the notion that epigenetic silencing of tumor suppressor genes, such as p21, may be the result of an imbalance in bidirectional transcription levels. This imbalance allows the unchecked antisense RNA to direct silent state epigenetic marks to the sense promoter, resulting in stable transcriptional gene silencing.

Non-coding RNAs have been shown to modulate transcriptional expression of genes in human cells. This form of gene regulation has been shown to be the result of RNA directing silent state epigenetic changes to the targeted gene promoter. Shortly after this seminal observation, small RNAs targeted to AT-rich regions of gene promoters were shown to modulate gene activation, termed RNA activation. While much is known regarding non-coding RNA-mediated transcriptional gene silencing, the mechanism of RNA activation has remained elusive. Here, we present evidence that RNA activation is the result of deregulation of endogenous bidirectional transcription. The antisense transcript in the bidirectionally transcribed gene is shown here to be operative in directing silent state epigenetic marks to the sense gene promoter. Suppression of the antisense transcript results in gene activation. Overall, these data support the notion that bidirectional transcription is an endogenous mechanism whereby RNA-directed gene regulation is operative and that RNA activation is the result of a disruption of this endogenous pathway. An understanding of this mode of RNA-based regulation will prove exceptionally useful in understanding gene expression as well as potential therapeutic approaches to controlling gene expression in human cells.

Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms

The hippocampal expression profiles of wild-type mice and mice transgenic for δC-doublecortin-like kinase were compared with Solexa/Illumina deep sequencing technology and five different microarray platforms. With Illumina's digital gene expression assay, we obtained ∼2.4 million sequence tags per sample, their abundance spanning four orders of magnitude. Results were highly reproducible, even across laboratories. With a dedicated Bayesian model, we found differential expression of 3179 transcripts with an estimated false-discovery rate of 8.5%. This is a much higher figure than found for microarrays. The overlap in differentially expressed transcripts found with deep sequencing and microarrays was most significant for Affymetrix. The changes in expression observed by deep sequencing were larger than observed by microarrays or quantitative PCR. Relevant processes such as calmodulin-dependent protein kinase activity and vesicle transport along microtubules were found affected by deep sequencing but not by microarrays. While undetectable by microarrays, antisense transcription was found for 51% of all genes and alternative polyadenylation for 47%. We conclude that deep sequencing provides a major advance in robustness, comparability and richness of expression profiling data and is expected to boost collaborative, comparative and integrative genomics studies.

The effect of temperature on Natural Antisense Transcript (NAT) expression in Aspergillus flavus

Naturally occurring Antisense Transcripts (NATs) compose an emerging group of regulatory RNAs. These regulatory elements appear in all organisms examined, but little is known about global expression of NATs in fungi. Analysis of currently available EST sequences suggests that 352 cis NATs are present in Aspergillus flavus. An Affymetrix GeneChip microarray containing probes for these cis NATs, as well as all predicted genes in A. flavus, allowed a whole genome expression analysis of these elements in response to two ecologically important temperatures for the fungus. RNA expression analysis showed that 32 NATs and 2,709 genes were differentially expressed between 37 degrees C, the optimum temperature for growth, and 28 degrees C, the conducive temperature for the biosynthesis of aflatoxin (AF) and many other secondary metabolites. These NATs correspond to sense genes with diverse functions including transcription initiation, carbohydrate processing and binding, temperature sensitive morphogenesis, and secondary metabolism. This is the first report of a whole genome transcriptional analysis of NAT expression in a fungus.

Arabidopsis transcriptome analysis under drought, cold, high-salinity and ABA treatment conditions using a tiling array

Plants respond and adapt to drought, cold and high-salinity stresses in order to survive. In this study, we applied Arabidopsis Affymetrix tiling arrays to study the whole genome transcriptome under drought, cold, high-salinity and ABA treatment conditions. The bioinformatic analysis using the tiling array data showed that 7,719 non-AGI transcriptional units (TUs) exist in the unannotated "intergenic" regions of Arabidopsis genome. These include 1,275 and 181 TUs that are induced and downregulated, respectively, by the stress or ABA treatments. Most of the non-AGI TUs are hypothetical non-protein-coding RNAs. About 80% of the non-AGI TUs belong to pairs of the fully overlapping sense-antisense transcripts (fSATs). Significant linear correlation between the expression ratios (treated/untreated) of the sense TUs and the ratios of the antisense TUs was observed in the SATs of AGI code/non-AGI TU. We studied the biogenesis mechanisms of the stress- or ABA-inducible antisense RNAs and found that the expression of sense TUs is necessary for the stress- or ABA-inducible expression of the antisense TUs in the fSATs (AGI code/non-AGI TU).

Synergistic upregulation of erythropoietin receptor (EPO-R) expression by sense and antisense EPO-R transcripts in the canine lung

We previously found increased erythropoietin receptor (EPO-R) protein levels in vigorously growing canine lungs after pneumonectomy (PNX), suggesting a role for paracrine EPO signaling in lung growth and remodeling. Now we find that sense and antisense EPO-R transcripts (sEPO-R and asEPO-R, respectively) are concordantly up-regulated in the post-PNX remaining lung, leading to the hypothesis that sEPO-R and asEPO-R interactions enhance EPO signaling during lung growth. We cloned a canine asEPO-R cDNA, which is fully complementary to the sense strand of the EPO-R gene from 2.5kb 3' to the sense stop codon, and extends into the 5' UTR of the sEPO-R transcript. Both asEPO-R and sEPO-R transcripts colocalize with EPO-R protein in the same lung cells. In cultured human embryonic kidney (HEK293) cells, transfection with sEPO-R (+FLAG tag) cDNA alone increased EPO-R protein expression (anti-EPO-R and anti-FLAG). At constant sEPO-R cDNA levels, cotransfection with escalating asEPO-R cDNA further increased recombinant EPO-R protein expression. The asEPO-R transcript harbors two putative opening reading frames (ORFs). Separate transfection of each asEPO-R ORF cDNA resulted in differential stimulatory effects on EPO-R protein expression. We conclude that both sEPO-R and asEPO-R transcripts contribute to in vivo up-regulation of EPO-R protein expression in the post-PNX remaining lung. This demonstrates synergism between sense-antisense EPO-R transcripts in response to physiological stimulation in a robust model of induced lung growth.

Non-random retention of protein-coding overlapping genes in Metazoa

Background

Although the overlap of transcriptional units occurs frequently in eukaryotic genomes, its evolutionary and biological significance remains largely unclear. Here we report a comparative analysis of overlaps between genes coding for well-annotated proteins in five metazoan genomes (human, mouse, zebrafish, fruit fly and worm).

Results

For all analyzed species the observed number of overlapping genes is always lower than expected assuming functional neutrality, suggesting that gene overlap is negatively selected. The comparison to the random distribution also shows that retained overlaps do not exhibit random features: antiparallel overlaps are significantly enriched, while overlaps lying on the same strand and those involving coding sequences are highly underrepresented. We confirm that overlap is mostly species-specific and provide evidence that it frequently originates through the acquisition of terminal, non-coding exons. Finally, we show that overlapping genes tend to be significantly co-expressed in a breast cancer cDNA library obtained by 454 deep sequencing, and that different overlap types display different patterns of reciprocal expression.

Conclusion

Our data suggest that overlap between protein-coding genes is selected against in Metazoa. However, when retained it may be used as a species-specific mechanism for the reciprocal regulation of neighboring genes. The tendency of overlaps to involve non-coding regions of the genes leads to the speculation that the advantages achieved by an overlapping arrangement may be optimized by evolving regulatory non-coding transcripts.

Mapping of small RNAs in the human ENCODE regions

The elucidation of the largely unknown transcriptome of small RNAs is crucial for the understanding of genome and cellular function. We report here the results of the analysis of small RNAs (< 50 nt) in the ENCODE regions of the human genome. Size-fractionated RNAs from four different cell lines (HepG2, HelaS3, GM06990, SK-N-SH) were mapped with the forward and reverse ENCODE high-density resolution tiling arrays. The top 1% of hybridization signals are termed SmRfrags (Small RNA fragments). Eight percent of SmRfrags overlap the GENCODE genes (CDS), given that the majority map to intergenic regions (34%), intronic regions (53%), and untranslated regions (UTRs) (5%). In addition, 9.6% and 16.8% of SmRfrags in the 5' UTR regions overlap significantly with His/Pol II/TAF250 binding sites and DNase I Hypersensitive sites, respectively (compared to the 5.3% and 9% expected). Interestingly, 17%-24% (depending on the cell line) of SmRfrags are sense-antisense strand pairs that show evidence of overlapping transcription. Only 3.4% and 7.2% of SmRfrags in intergenic regions overlap transcribed fragments (Txfrags) in HeLa and GM06990 cell lines, respectively. We hypothesized that a fraction of the identified SmRfrags corresponded to microRNAs. We tested by Northern blot a set of 15 high-likelihood predictions of microRNA candidates that overlap with smRfrags and validated three potential microRNAs ( approximately 20 nt length). Notably, most of the remaining candidates showed a larger hybridizing band ( approximately 100 nt) that could be a microRNA precursor. The small RNA transcriptome is emerging as an important and abundant component of the genome function.

antiCODE: a natural sense-antisense transcripts database

Background

Natural antisense transcripts (NATs) are endogenous RNA molecules that exhibit partial or complete complementarity to other RNAs, and that may contribute to the regulation of molecular functions at various levels. In recent years, large-scale NAT screens in several model organisms have produced much data, but there is no database to assemble all these data. AntiCODE intends to function as an integrated NAT database for this purpose.

Results

This release of antiCODE contains more than 30,000 non-redundant natural sense-antisense transcript pairs from 12 eukaryotic model organisms. In order to provide an integrated NAT research platform, efficient browser, search and Blast functions have been included to enable users to easily access information through parameters such as species, accession number, overlapping patterns, coding potential etc. In addition to the collected information, antiCODE also introduces a simple classification system to facilitate the study of natural antisense transcripts.

Conclusion

Though a few similar databases also dealing with NATs have appeared lately, antiCODE is the most comprehensive among these, comprising almost all currently detected NAT pairs.

How do synonymous mutations affect fitness?

While it has often been assumed that, in humans, synonymous mutations would have no effect on fitness, let alone cause disease, this position has been questioned over the last decade. There is now considerable evidence that such mutations can, for example, disrupt splicing and interfere with miRNA binding. Two recent publications suggest involvement of additional mechanisms: modification of protein abundance most probably mediated by alteration in mRNA stability and modification of protein structure and activity, probably mediated by induction of translational pausing. These case histories put a further nail into the coffin of the assumption that synonymous mutations must be neutral.

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.

Genome-wide natural antisense transcription: coupling its regulation to its different regulatory mechanisms

Many genomic loci contain transcription units on both strands, therefore two oppositely oriented transcripts can overlap. Often, one strand codes for a protein, whereas the transcript from the other strand is non-encoding. Such natural antisense transcripts (NATs) can negatively regulate the conjugated sense transcript. NATs are highly prevalent in a wide range of species--for example, around 15% of human protein-encoding genes have an associated NAT. The regulatory mechanisms by which NATs act are diverse, as are the means to control their expression. Here, we review the current understanding of NAT function and its mechanistic basis, which has been gathered from both individual gene cases and genome-wide studies. In parallel, we survey findings about the regulation of NAT transcription. Finally, we hypothesize that the regulation of antisense transcription might be tailored to its mode of action. According to this model, the observed relationship between the expression patterns of NATs and their targets might indicate the regulatory mechanism that is in action.

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.