Sense and antisense transcripts of convergent gene pairs in Arabidopsis thaliana can share a common polyadenylation region

CCIVR2 facilitates comprehensive identification of both overlapping and non-overlapping antisense transcripts within specified regions

Pairs of sense and antisense transcriptions that are adjacent at their 5' and 3' regions are called divergent and convergent transcription, respectively. However, the structural properties of divergent/convergent transcription in different species or RNA biotypes are poorly characterized. Here, we developed CCIVR2, a program that facilitates identification of both overlapping and non-overlapping antisense transcripts produced from divergent/convergent transcription whose transcription start sites (TSS) or transcript end sites (TES) are located within a specified region. We used CCIVR2 to analyze antisense transcripts starting around the sense TSS (from divergent transcription) or ending around the sense TES (from convergent transcription) in 11 different species and found species- and RNA biotype-specific features of divergent/convergent transcription. Furthermore, we confirmed that CCIVR2 enables the identification of multiple sense/antisense transcript pairs from divergent transcription, including those with known functions in processes such as embryonic stem cell differentiation and TGFβ stimulation. CCIVR2 is therefore a valuable bioinformatics tool that facilitates the characterization of divergent/convergent transcription in different species and aids the identification of functional sense/antisense transcript pairs from divergent transcription in specified biological processes.

Novel Stress-Inducible Antisense RNAs of Protein-Coding Loci Are Synthesized by RNA-Dependent RNA Polymerase

Our previous study identified approximately 6,000 abiotic stress-responsive noncoding transcripts existing on the antisense strand of protein-coding genes and implied that a type of antisense RNA was synthesized from a sense RNA template by RNA-dependent RNA polymerase (RDR). Expression analyses revealed that the expression of novel abiotic stress-induced antisense RNA on 1,136 gene loci was reduced in the rdr1/2/6 mutants. RNase protection indicated that the RD29A antisense RNA and other RDR1/2/6-dependent antisense RNAs are involved in the formation of dsRNA. The accumulation of stress-inducible antisense RNA was decreased and increased in dcp5 and xrn4, respectively, but not changed in dcl2/3/4, nrpd1a and nrpd1b RNA-seq analyses revealed that the majority of the RDR1/2/6-dependent antisense RNA loci did not overlap with RDR1/2/6-dependent 20-30 nt RNA loci. Additionally, rdr1/2/6 mutants decreased the degradation rate of the sense RNA and exhibited arrested root growth during the recovery stage following a drought stress, whereas dcl2/3/4 mutants did not. Collectively, these results indicate that RDRs have stress-inducible antisense RNA synthesis activity and a novel biological function that is different from the known endogenous small RNA pathways from protein-coding genes. These data reveal a novel mechanism of RNA regulation during abiotic stress response that involves complex RNA degradation pathways.

Differential genomic arrangements in Caryophyllales through deep transcriptome sequencing of A. hypochondriacus

Genome duplication event in edible dicots under the orders Rosid and Asterid, common during the oligocene period, is missing for species under the order Caryophyllales. Despite this, grain amaranths not only survived this period but display many desirable traits missing in species under rosids and asterids. For example, grain amaranths display traits like C4 photosynthesis, high-lysine seeds, high-yield, drought resistance, tolerance to infection and resilience to stress. It is, therefore, of interest to look for minor genome rearrangements with potential functional implications that are unique to grain amaranths. Here, by deep sequencing and assembly of 16 transcriptomes (86.8 billion bases) we have interrogated differential genome rearrangement unique to Amaranthus hypochondriacus with potential links to these phenotypes. We have predicted 125,581 non-redundant transcripts including 44,529 protein coding transcripts identified based on homology to known proteins and 13,529 predicted as novel/amaranth specific coding transcripts. Of the protein coding de novo assembled transcripts, we have identified 1810 chimeric transcripts. More than 30% and 19% of the gene pairs within the chimeric transcripts are found within the same loci in the genomes of A. hypochondriacus and Beta vulgaris respectively and are considered real positives. Interestingly, one of the chimeric transcripts comprises two important genes, namely DHDPS1, a key enzyme implicated in the biosynthesis of lysine, and alpha-glucosidase, an enzyme involved in sucrose catabolism, in close proximity to each other separated by a distance of 612 bases in the genome of A. hypochondriacus in a convergent configuration. We have experimentally validated that transcripts of these two genes are also overlapping in the 3' UTR with their expression negatively correlated from bud to mature seed, suggesting a potential link between the high seed lysine trait and unique genome organization.

Post-transcriptional gene silencing triggered by sense transgenes involves uncapped antisense RNA and differs from silencing intentionally triggered by antisense transgenes

Although post-transcriptional gene silencing (PTGS) has been studied for more than a decade, there is still a gap in our understanding of how de novo silencing is initiated against genetic elements that are not supposed to produce double-stranded (ds)RNA. Given the pervasive transcription occurring throughout eukaryote genomes, we tested the hypothesis that unintended transcription could produce antisense (as)RNA molecules that participate to the initiation of PTGS triggered by sense transgenes (S-PTGS). Our results reveal a higher level of asRNA in Arabidopsis thaliana lines that spontaneously trigger S-PTGS than in lines that do not. However, PTGS triggered by antisense transgenes (AS-PTGS) differs from S-PTGS. In particular, a hypomorphic ago1 mutation that suppresses S-PTGS prevents the degradation of asRNA but not sense RNA during AS-PTGS, suggesting a different treatment of coding and non-coding RNA by AGO1, likely because of AGO1 association to polysomes. Moreover, the intended asRNA produced during AS-PTGS is capped whereas the asRNA produced during S-PTGS derives from 3′ maturation of a read-through transcript and is uncapped. Thus, we propose that uncapped asRNA corresponds to the aberrant RNA molecule that is converted to dsRNA by RNA-DEPENDENT RNA POLYMERASE 6 in siRNA-bodies to initiate S-PTGS, whereas capped asRNA must anneal with sense RNA to produce dsRNA that initiate AS-PTGS.

RNA around the clock - regulation at the RNA level in biological timing

The circadian timing system in plants synchronizes their physiological functions with the environment. This is achieved by a global control of gene expression programs with a considerable part of the transcriptome undergoing 24-h oscillations in steady-state abundance. These circadian oscillations are driven by a set of core clock proteins that generate their own 24-h rhythm through periodic feedback on their own transcription. Additionally, post-transcriptional events are instrumental for oscillations of core clock genes and genes in clock output. Here we provide an update on molecular events at the RNA level that contribute to the 24-h rhythm of the core clock proteins and shape the circadian transcriptome. We focus on the circadian system of the model plant Arabidopsis thaliana but also discuss selected regulatory principles in other organisms.

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.

Association of large noncoding RNA HOTAIR expression and its downstream intergenic CpG island methylation with survival in breast cancer

Large noncoding RNA HOTAIR, transcribed from the antisense strand of HOXC12, interacts with Polycomb Repressive Complex 2 (PRC2) in the regulation of gene activities. Recent work suggests that it may have effects on breast cancer progression and survival. We evaluated HOTAIR expression and the methylation status of its downstream intergenic CpG island in primary breast cancers, and examined associations of these factors with clinical and pathologic features and patient survival. HOTAIR expression and DNA methylation were analyzed in tissue from 348 primary breast cancers with quantitative RT-PCR and methylation-specific PCR, respectively. HOTAIR expression and methylation varied widely in the tissues. A positive correlation was found between DNA methylation and HOTAIR expression. Methylation was associated with unfavorable disease characteristics, whereas no significant associations were found between HOTAIR expression and clinical or pathologic features. In multivariate, but not in univariate, Cox proportional hazard regression models, patients with high HOTAIR expression had lower risks of relapse and mortality than those with low HOTAIR expression. These findings suggest that the intergenic DNA methylation may have important biologic relevance in regulating HOTAIR expression, and that HOTAIR expression may not be an independent prognostic marker in breast cancer, but needs further validation in independent studies.

De novo DNA methylation activity of methyltransferase 1 (MET1) partially restores body methylation in Arabidopsis thaliana

Arabidopsis METHYLTRANSFERASE 1 (MET1) controls faithful maintenance of cytosine methylation at CG sites in repetitive regions and central body regions of active genes. If MET1 is removed in a mutant background, CG methylation is lost and is only restored in specific heterochromatic regions that have maintained competence for re-methylation due to the presence of small RNAs and the RNA-directed DNA methylation pathway that controls de novo DNA methylation functions. We analysed re-methylation at a locus that loses body methylation in an met1 mutant. We found that body methylation at this locus is at least partially restored when MET1 is re-introduced into the met1 mutant background, either via genetic cross or DNA transfer. Re-methylation is region-specific but random with respect to individual CG targets, does not require passage through the germline, and its efficiency appears to be influenced by transcription. This suggests that, at least at some loci, MET1 has de novo methylation activity that can restore lost body methylation patterns. We propose that this activity helps to stabilize body methylation patterns, and the random target site selection probably also enhances the variability of body methylation patterns.

A pair of partially overlapping Arabidopsis genes with antagonistic circadian expression

A large number of plant genes are aligned with partially overlapping genes in antisense orientation. Transcription of both genes would therefore favour the formation of double-stranded RNA, providing a substrate for the RNAi machinery, and enhanced antisense transcription should therefore reduce sense transcript levels. We have identified a gene pair that resembles a model for antisense-based gene regulation as a T-DNA insertion into the antisense gene causes a reduction in antisense transcript levels and an increase in sense transcript levels. The same effect was, however, also observed when the two genes were inserted as transgenes into different chromosomal locations, independent of the sense and antisense gene being expressed individually or jointly. Our results therefore indicate that antagonistic changes in sense/antisense transcript levels do not necessarily reflect antisense-mediated regulation. More likely, the partial overlap of the two genes may have favoured the evolution of antagonistic expression patterns preventing RNAi effects.

The half-size ABC transporters STR1 and STR2 are indispensable for mycorrhizal arbuscule formation in rice

The central structure of the symbiotic association between plants and arbuscular mycorrhizal (AM) fungi is the fungal arbuscule that delivers minerals to the plant. Our earlier transcriptome analyses identified two half-size ABCG transporters that displayed enhanced mRNA levels in mycorrhizal roots. We now show specific transcript accumulation in arbusculated cells of both genes during symbiosis. Presently, arbuscule-relevant factors from monocotyledons have not been reported. Mutation of either of the Oryza sativa (rice) ABCG transporters blocked arbuscule growth of different AM fungi at a small and stunted stage, recapitulating the phenotype of Medicago truncatula stunted arbuscule 1 and 2 (str1 and str2) mutants that are deficient in homologous ABCG genes. This phenotypic resemblance and phylogenetic analysis suggest functional conservation of STR1 and STR2 across the angiosperms. Malnutrition of the fungus underlying limited arbuscular growth was excluded by the absence of complementation of the str1 phenotype by wild-type nurse plants. Furthermore, plant AM signaling was found to be intact, as arbuscule-induced marker transcript accumulation was not affected in str1 mutants. Strigolactones have previously been hypothesized to operate as intracellular hyphal branching signals and possible substrates of STR1 and STR2. However, full arbuscule development in the strigolactone biosynthesis mutants d10 and d17 suggested strigolactones to be unlikely substrates of STR1/STR2. Interestingly, rice STR1 is associated with a cis-natural antisense transcript (antiSTR1). Analogous to STR1 and STR2, at the root cortex level, the antiSTR1 transcript is specifically detected in arbusculated cells, suggesting unexpected modes of STR1 regulation in rice.

RNA-based regulation in the plant circadian clock

The circadian clock is an endogenous, approximately 24-h timer that enables plants to anticipate daily changes in their environment and regulates a considerable fraction of the transcriptome. At the core of the circadian system is the oscillator, made up of interconnected feedback loops, involving transcriptional regulation of clock genes and post-translational modification of clock proteins. Recently, it has become clear that post-transcriptional events are also critical for shaping rhythmic mRNA and protein profiles. This review covers regulation at the RNA level of both the core clock and output genes in Arabidopsis (Arabidopsis thaliana), with comparisons with other model organisms. We discuss the role of splicing, mRNA decay and translational regulation as well as recent insights into rhythms of noncoding regulatory RNAs.

Dual RNAs in plants

Plants have remarkable developmental plasticity, and the same genotype can result in different phenotypes depending on environmental variation. Indeed, abiotic stresses or biotic interactions affect organogenesis and post-embryonic growth and significantly affect gene regulation. The large diversity of non-protein-coding RNAs (npcRNAs) and genes containing only short open reading frames that are expressed during plant growth and development, contribute to the regulation of gene expression. Certain npcRNAs code for oligopeptides and may possess additional biological activity linked to the RNA moiety. The ENOD40 gene is a dual RNA that is activated during a symbiotic interaction leading to root nodule organogenesis. Both the oligopeptides encoded by ENOD40 and the structured regions of the ENOD40 RNA have been shown to interact with different proteins in the cell to control enzymatic activities or induce the relocalisation of ribonucleoproteins, respectively. Other npcRNAs encode for small signalling peptides or are the precursors of small RNAs involved in post-transcriptional or transcriptional gene silencing. They may have RNA-related activities or encode peptides (or even larger proteins), and therefore act as dual RNAs. In addition, long natural antisense RNAs with a coding function and a regulatory RNA-mediated action that are expressed in response to abiotic stress in plants have been identified. In certain cases, these RNAs lead to the synthesis of nat-siRNAs, that are small RNAs derived from the overlapping double-stranded RNA region of natural antisense RNAs, which facilitates the silencing of complementary mRNAs. Finally, the advent of deep sequencing technologies has identified a large number of non-protein-coding RNAs in plants, which could be a large reservoir for dual RNAs.