Natural antisense transcripts with coding capacity in Arabidopsis may have a regulatory role that is not linked to double-stranded RNA degradation

Antisense Transcription in Plants: A Systematic Review and an Update on cis-NATs of Sugarcane

Initially, natural antisense transcripts (NATs, natRNAs, or asRNAs) were considered repressors; however, their functions in gene regulation are diverse. Positive, negative, or neutral correlations to the cognate gene expression have been noted. Although the first studies were published about 50 years ago, there is still much to be investigated regarding antisense transcripts in plants. A systematic review of scientific publications available in the Web of Science databases was conducted to contextualize how the studying of antisense transcripts has been addressed. Studies were classified considering three categories: “Natural antisense” (208), artificial antisense used in “Genetic Engineering” (797), or “Natural antisense and Genetic Engineering”-related publications (96). A similar string was used for a systematic search in the NCBI Gene database. Of the 1132 antisense sequences found for plants, only 0.8% were cited in PubMed and had antisense information confirmed. This value was the lowest when compared to fungi (2.9%), bacteria (2.3%), and mice (54.1%). Finally, we present an update for the cis-NATs identified in Saccharum spp. Of the 1413 antisense transcripts found in different experiments, 25 showed concordant expressions, 22 were discordant, 1264 did not correlate with the cognate genes, and 102 presented variable results depending on the experiment.

Long non-coding RNA IGFBP7-AS1 accelerates the odontogenic differentiation of stem cells from human exfoliated deciduous teeth by regulating IGFBP7 expression

Stem cells from human exfoliated deciduous teeth (SHED) are attractive seed cells for dental tissue engineering. We identified the effect of the long noncoding RNA insulin-like growth factor-binding protein 7 antisense RNA 1 (lncRNA IGFBP7-AS1) in vivo and its underlying mechanism during SHED odontogenic differentiation. IGFBP7-AS1 and insulin-like growth factor-binding protein 7 (IGFBP7) were overexpressed using lentiviruses. IGFBP7 expression was knocked down with small interfering RNA. The effect of IGFBP7-AS1 in vivo was confirmed by animal experiments. The effect of IGFBP7 on SHED odontogenic differentiation was assessed with alkaline phosphatase staining, alizarin red S staining, quantitative reverse transcription-PCR, and western blotting. The relationship between IGFBP7-AS1 and IGFBP7 was confirmed by quantitative reverse transcription–PCR and western blotting. IGFBP7-AS1 promoted SHED odontogenesis in vivo, and regulated the expression of the coding gene IGFBP7 positively. Inhibiting IGFBP7 led to suppress SHED odontogenic differentiation while IGFBP7 overexpression had the opposite effect. IGFBP7-AS1 enhanced the stability of IGFBP7. IGFBP7-AS1 promoted SHED odontogenic differentiation in vivo. The underlying mechanism may involve the enhancement of IGFBP7 stability. This may provide novel potential targets for dental tissue engineering.

Promoter switching in response to changing environment and elevated expression of protein-coding genes overlapping at their 5' ends

Despite the number of studies focused on sense-antisense transcription, the key question of whether such organization evolved as a regulator of gene expression or if this is only a byproduct of other regulatory processes has not been elucidated to date. In this study, protein-coding sense-antisense gene pairs were analyzed with a particular focus on pairs overlapping at their 5' ends. Analyses were performed in 73 human transcription start site libraries. The results of our studies showed that the overlap between genes is not a stable feature and depends on which TSSs are utilized in a given cell type. An analysis of gene expression did not confirm that overlap between genes causes downregulation of their expression. This observation contradicts earlier findings. In addition, we showed that the switch from one promoter to another, leading to genes overlap, may occur in response to changing environment of a cell or tissue. We also demonstrated that in transfected and cancerous cells genes overlap is observed more often in comparison with normal tissues. Moreover, utilization of overlapping promoters depends on particular state of a cell and, at least in some groups of genes, is not merely coincidental.

Inter-dependent Centrosomal Co-localization of the cen and ik2 cis-Natural Antisense mRNAs in Drosophila

Overlapping genes are prevalent in most genomes, but the extent to which this organization influences regulatory events operating at the post-transcriptional level remains unclear. Studying the cen and ik2 genes of Drosophila melanogaster, which are convergently transcribed as cis-natural antisense transcripts (cis-NATs) with overlapping 3' UTRs, we found that their encoded mRNAs strikingly co-localize to centrosomes. These transcripts physically interact in a 3' UTR-dependent manner, and the targeting of ik2 requires its 3' UTR sequence and the presence of cen mRNA, which serves as the main driver of centrosomal co-localization. The cen transcript undergoes localized translation in proximity to centrosomes, and its localization is perturbed by polysome-disrupting drugs. By interrogating global fractionation-sequencing datasets generated from Drosophila and human cellular models, we find that RNAs expressed as cis-NATs tend to co-localize to specific subcellular fractions. This work suggests that post-transcriptional interactions between RNAs with complementary sequences can dictate their localization fate in the cytoplasm.

Plant Long Noncoding RNAs: New Players in the Field of Post-Transcriptional Regulations

The first reference to the "C-value paradox" reported an apparent imbalance between organismal genome size and morphological complexity. Since then, next-generation sequencing has revolutionized genomic research and revealed that eukaryotic transcriptomes contain a large fraction of non-protein-coding components. Eukaryotic genomes are pervasively transcribed and noncoding regions give rise to a plethora of noncoding RNAs with undeniable biological functions. Among them, long noncoding RNAs (lncRNAs) seem to represent a new layer of gene expression regulation, participating in a wide range of molecular mechanisms at the transcriptional and post-transcriptional levels. In addition to their role in epigenetic regulation, plant lncRNAs have been associated with the degradation of complementary RNAs, the regulation of alternative splicing, protein sub-cellular localization, the promotion of translation and protein post-translational modifications. In this review, we report and integrate numerous and complex mechanisms through which long noncoding transcripts regulate post-transcriptional gene expression in plants.

Functional classification of plant long noncoding RNAs: a transcript is known by the company it keeps

The extraordinary maturation in high-throughput sequencing technologies has revealed the existence of a complex network of transcripts in eukaryotic organisms, including thousands of long noncoding (lnc) RNAs with little or no protein-coding capacity. Subsequent discoveries have shown that lncRNAs participate in a wide range of molecular processes, controlling gene expression and protein activity though direct interactions with proteins, DNA or other RNA molecules. Although significant advances have been achieved in the understanding of lncRNA biology in the animal kingdom, the functional characterization of plant lncRNAs is still in its infancy and remains a major challenge. In this review, we report emerging functional and mechanistic paradigms of plant lncRNAs and partner molecules, and discuss how cutting-edge technologies may help to identify and classify yet uncharacterized transcripts into functional groups.

Additional Layer of Regulation via Convergent Gene Orientation in Yeasts

Convergent gene pairs can produce transcripts with complementary sequences. We had shown that mRNA duplexes form in vivo in Saccharomyces cerevisiae via interactions of mRNA overlapping 3′-ends and can lead to posttranscriptional regulatory events. Here we show that mRNA duplex formation is restricted to convergent genes separated by short intergenic distance, independently of their 3′-untranslated region (UTR) length. We disclose an enrichment in genes involved in biological processes related to stress among these convergent genes. They are markedly conserved in convergent orientation in budding yeasts, meaning that this mode of posttranscriptional regulation could be shared in these organisms, conferring an additional level for modulating stress response. We thus investigated the mechanistic advantages potentially conferred by 3′-UTR mRNA interactions. Analysis of genome-wide transcriptome data revealed that Pat1 and Lsm1 factors, having 3′-UTR binding preference and participating to the remodeling of messenger ribonucleoprotein particles, bind differently these messenger-interacting mRNAs forming duplexes in comparison to mRNAs that do not interact (solo mRNAs). Functionally, messenger-interacting mRNAs show limited translational repression upon stress. We thus propose that mRNA duplex formation modulates the regulation of mRNA expression by limiting their access to translational repressors. Our results thus show that posttranscriptional regulation is an additional factor that determines the order of coding genes.

Computational tools for plant small RNA detection and categorization

Small RNAs (sRNAs) are important short-length molecules with regulatory functions essential for plant development and plasticity. High-throughput sequencing of total sRNA populations has revealed that the largest share of sRNA remains uncategorized. To better understand the role of sRNA-mediated cellular regulation, it is necessary to create accurate and comprehensive catalogues of sRNA and their sequence features, a task that currently relies on nontrivial bioinformatic approaches. Although a large number of computational tools have been developed to predict features of sRNA sequences, these tools are mostly dedicated to microRNAs and none integrates the functionalities necessary to describe units from all sRNA pathways thus far discovered in plants. Here, we review the different classes of sRNA found in plants and describe available bioinformatics tools that can help in their detection and categorization.

Responsiveness and Adaptation to Salt Stress of the REDOX-RESPONSIVE TRANSCRIPTION FACTOR 1 (RRTF1) Gene are Controlled by its Promoter

The REDOX-RESPONSIVE TRANSCRIPTION FACTOR 1 (RRTF1) gene encodes a member of the ERF/AP2 transcription factor family involved in redox homeostasis. The RRTF1 gene shows tissue-specific responsiveness to various abiotic stress treatments including a response to salt stress in roots. An interesting feature of this response is an adaptation phase that follows its activation, when promoter levels revert to a base line level, even if salt stress is maintained. It is unclear if adaption is controlled by a switch in promoter activity or by changes in transcript levels. Here we show that the RRTF1 promoter is sufficient for the control of both activation and adaptation to salt stress. As constitutive expression of RRTF1 turned out to be detrimental to the plant, we propose that promoter-regulated adaptation evolved as a protection mechanism to balance the beneficial effects of short-term gene activation and the detrimental effects of long-term gene expression.

Arabidopsis PIRL6 Is Essential for Male and Female Gametogenesis and Is Regulated by Alternative Splicing

Plant intracellular Ras-group leucine-rich repeat (LRR) proteins (PIRLs) are related to Ras-interacting animal LRR proteins that participate in developmental cell signaling. Systematic knockout analysis has implicated some members of the Arabidopsis (Arabidopsis thaliana) PIRL family in pollen development. However, for PIRL6, no bona fide knockout alleles have been recovered, suggesting that it may have an essential function in both male and female gametophytes. To test this hypothesis, we investigated PIRL6 expression and induced knockdown by RNA interference. Knockdown triggered defects in gametogenesis, resulting in abnormal pollen and early developmental arrest in the embryo sac. Consistent with this, PIRL6 was expressed in gametophytes: functional transcripts were detected in wild-type flowers but not in sporocyteless (spl) mutant flowers, which do not produce gametophytes. A genomic PIRL6-GFP fusion construct confirmed expression in both pollen and the embryo sac. Interestingly, PIRL6 is part of a convergent overlapping gene pair, a scenario associated with an increased likelihood of alternative splicing. We detected multiple alternative PIRL6 mRNAs in vegetative organs and spl mutant flowers, tissues that lacked the functionally spliced transcript. cDNA sequencing revealed that all contained intron sequences and premature termination codons. These alternative mRNAs accumulated in the nonsense-mediated decay mutant upf3, indicating that they are normally subjected to degradation. Together, these results demonstrate that PIRL6 is required in both male and female gametogenesis and suggest that sporophytic expression is negatively regulated by unproductive alternative splicing. This posttranscriptional mechanism may function to minimize PIRL6 protein expression in sporophyte tissues while allowing the overlapping adjacent gene to remain widely transcribed.

A novel antisense long noncoding RNA, TWISTED LEAF, maintains leaf blade flattening by regulating its associated sense R2R3-MYB gene in rice

Natural antisense long noncoding RNAs (lncRNAs) are widespread in many organisms. However, their biological functions remain largely unknown, particularly in plants. We report the identification and characterization of an endogenous lncRNA, TWISTED LEAF (TL), which is transcribed from the opposite strand of the R2R3 MYB transcription factor gene locus, OsMYB60, in rice (Oryza sativa). TL and OsMYB60 were found to be coexpressed in many different tissues, and the expression level of TL was higher than that of OsMYB60. Downregulation of TL by RNA interference (RNAi) and overexpression of OsMYB60 resulted in twisted leaf blades in transgenic rice. The expression level of OsMYB60 was significantly increased in TL-RNAi transgenic plants. This suggests that TL may play a cis-regulatory role on OsMYB60 in leaf morphological development. We also determined that the antisense transcription suppressed the sense gene expression by mediating chromatin modifications. We further discovered that a C2H2 transcription factor, OsZFP7, is an OsMYB60 binding partner and involved in leaf development. Taken together, these findings reveal that the cis-natural antisense lncRNA plays a critical role in maintaining leaf blade flattening in rice. Our study uncovers a regulatory mechanism of lncRNA in plant leaf development.

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.

Plant organ evolution revealed by phylotranscriptomics in Arabidopsis thaliana

The evolution of phenotypes occurs through changes both in protein sequence and gene expression levels. Though much of plant morphological evolution can be explained by changes in gene expression, examining its evolution has challenges. To gain a new perspective on organ evolution in plants, we applied a phylotranscriptomics approach. We combined a phylostratigraphic approach with gene expression based on the strand-specific RNA-seq data from seedling, floral bud, and root of 19 Arabidopsis thaliana accessions to examine the age and sequence divergence of transcriptomes from these organs and how they adapted over time. Our results indicate that, among the sense and antisense transcriptomes of these organs, the sense transcriptomes of seedlings are the evolutionarily oldest across all accessions and are the most conserved in amino acid sequence for most accessions. In contrast, among the sense transcriptomes from these same organs, those from floral bud are evolutionarily youngest and least conserved in sequence for most accessions. Different organs have adaptive peaks at different stages in their evolutionary history; however, all three show a common adaptive signal from the Magnoliophyta to Brassicale stage. Our research highlights how phylotranscriptomic analyses can be used to trace organ evolution in the deep history of plant species.

From intronization to intron loss: How the interplay between mRNA-associated processes can shape the architecture and the expression of eukaryotic genes

Transcription-coupled processes such as capping, splicing, and cleavage/polyadenylation participate in the journey from genes to proteins. Although they are traditionally thought to serve only as steps in the generation of mature mRNAs, a synthesis of available data indicates that these processes could also act as a driving force for the evolution of eukaryotic genes. A theoretical framework for how mRNA-associated processes may shape gene structure and expression has recently been proposed. Factors that promote splicing and cleavage/polyadenylation in this framework compete for access to overlapping or neighboring signals throughout the transcription cycle. These antagonistic interactions allow mechanisms for intron gain and splice site recognition as well as common trends in eukaryotic gene structure and expression to be coherently integrated. Here, I extend this framework further. Observations that largely (but not exclusively) revolve around the formation of DNA-RNA hybrid structures, called R loops, and promoter directionality are integrated. Additionally, the interplay between splicing factors and cleavage/polyadenylation factors is theorized to also affect the formation of intragenic DNA double-stranded breaks thereby contributing to intron loss. The most notable prediction in this proposition is that RNA molecules can mediate intron loss by serving as a template to repair DNA double-stranded breaks. The framework presented here leverages a vast body of empirical observations, logically extending previous suggestions, and generating verifiable predictions to further substantiate the view that the intracellular environment plays an active role in shaping the structure and the expression of eukaryotic genes.

Small interfering RNAs from bidirectional transcripts of GhMML3_A12 regulate cotton fiber development

Natural antisense transcripts (NATs) are commonly observed in eukaryotic genomes, but only a limited number of such genes have been identified as being involved in gene regulation in plants. In this research, we investigated the function of small RNA derived from a NAT in fiber cell development. Using a map-based cloning strategy for the first time in tetraploid cotton, we cloned a naked seed mutant gene (N1 ) encoding a MYBMIXTA-like transcription factor 3 (MML3)/GhMYB25-like in chromosome A12, GhMML3_A12, that is associated with fuzz fiber development. The extremely low expression of GhMML3_A12 in N1 is associated with NAT production, driven by its 3' antisense promoter, as indicated by the promoter-driven histochemical staining assay. In addition, small RNA deep sequencing analysis suggested that the bidirectional transcriptions of GhMML3_A12 form double-stranded RNAs and generate 21-22 nt small RNAs. Therefore, in a fiber-specific manner, small RNA derived from the GhMML3_A12 locus can mediate GhMML3_A12 mRNA self-cleavage and result in the production of naked seeds followed by lint fiber inhibition in N1 plants. The present research reports the first observation of gene-mediated NATs and siRNA directly controlling fiber development in cotton.

Small RNAs in Plant Responses to Abiotic Stresses: Regulatory Roles and Study Methods

To survive under abiotic stresses in the environment, plants trigger a reprogramming of gene expression, by transcriptional regulation or translational regulation, to turn on protective mechanisms. The current focus of research on how plants cope with abiotic stresses has transitioned from transcriptomic analyses to small RNA investigations. In this review, we have summarized and evaluated the current methodologies used in the identification and validation of small RNAs and their targets, in the context of plant responses to abiotic stresses.

Cytoplasmic Control of Sense-Antisense mRNA Pairs

Transcriptome analyses have revealed that convergent gene transcription can produce many 3'-overlapping mRNAs in diverse organisms. Few studies have examined the fate of 3'-complementary mRNAs in double-stranded RNA-dependent nuclear phenomena, and nothing is known about the cytoplasmic destiny of 3'-overlapping messengers or their impact on gene expression. Here, we demonstrate that the complementary tails of 3'-overlapping mRNAs can interact in the cytoplasm and promote post-transcriptional regulatory events including no-go decay (NGD) in Saccharomyces cerevisiae. Genome-wide experiments confirm that these messenger-interacting mRNAs (mimRNAs) form RNA duplexes in wild-type cells and thus have potential roles in modulating the mRNA levels of their convergent gene pattern under different growth conditions. We show that the post-transcriptional fate of hundreds of mimRNAs is controlled by Xrn1, revealing the extent to which this conserved 5'-3' cytoplasmic exoribonuclease plays an unexpected but key role in the post-transcriptional control of convergent gene expression.

Molecular characterization of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) gene family from Citrus and the effect of fruit load on their expression

We recently identified a Citrus gene encoding SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor that contained a sequence complementary to miR156. Genes of the SPL family are known to play a role in flowering regulation and phase transition. In Citrus, the mRNA levels of the gene were significantly altered by fruit load in buds; under heavy fruit load (ON-Crop trees), known to suppress next year flowering, the mRNA levels were down-regulated, while fruit removal (de-fruiting), inducing next-year flowering, resulted in its up-regulation. In the current work, we set on to study the function of the gene. We showed that the Citrus SPL was able promote flowering independently of photoperiod in Arabidopsis, while miR156 repressed its flowering-promoting activity. In order to find out if fruit load affected the expression of additional genes of the SPL family, we identified and classified all SPL members in the Citrus genome, and studied their seasonal expression patterns in buds and leaves, and in response to de-fruiting. Results showed that two additional SPL-like genes and miR172, known to be induced by SPLs in Arabidopsis, were altered by fruit load. The relationships between these factors in relation to the fruit-load effect on Citrus flowering are discussed.

Rapid immunopurification of ribonucleoprotein complexes of plants

Hundreds of RNA binding proteins posttranscriptionally regulate gene expression, but relatively few have been characterized in plants. One successful approach to determine protein function has been to identify interacting molecules and the conditions of their association. The ribonucleoprotein immunopurification (RIP) assay facilitates the identification and quantitative comparison of RNA association to specific proteins under different experimental conditions. A variety of molecular techniques can be used to analyze the enriched RNAs, whether few as in the case of highly specific interactions, or many. Identification of associated RNAs can inform hypothesis generation about the processes or pathways regulated by the target protein. Downstream analysis of associated RNA sequences can lead to the identification of candidate motifs or features that mediate the protein-RNA interaction. We present a rapid method for RIP from tissues of plants that is suitable for experiments that require immediate tissue cryopreservation, such as monitoring a rapid response to an environmental stimulus.

NERF encodes a RING E3 ligase important for drought resistance and enhances the expression of its antisense gene NFYA5 in Arabidopsis

NFYA5 is an important drought-stress inducible transcription factor gene that is targeted by miR169 in Arabidopsis. We show here that the cis-natural antisense transcript gene of NFYA5, NFYA5 Enhancing RING FINGER (NERF), can produce siRNAs from their overlapping region (OR) and affect NFYA5 transcripts by functioning together with miR169. The NERF protein functions as an E3 ligase for ubiquitination. Overexpression of NERF or OR cDNA leads to siRNANERF accumulation, miR169 repression, and NFYA5 transcript enhancement; knock-down of NERF transcripts by an artificial miRNA enhances miR169 abundance and reduces NFYA5 transcripts. Overexpression of NFYA5 does not affect the NERF mRNA level. Deep sequencing of the small RNA library from 35S::OR plants identifies 960 sequences representing 323 unique siRNAs that originate from OR; the sequences of some siRNANERF are similar/complementary to those of miR169. Overexpression of the 195- to 280-bp OR cDNA-containing siRNAs similar/complementary to miR169 also leads to the accumulation of NFYA5 transcripts. Analysis of NERF knock-down plants and NERF overexpression lines showed that, like NFYA5, NERF is important for controlling stomatal aperture and drought resistance. This regulatory model might apply to other natural antisense transcripts with positively correlated expression patterns.

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.

Co-option of the piRNA pathway for germline-specific alternative splicing of C. elegans TOR

Many eukaryotic genes contain embedded antisense transcripts and repetitive sequences of unknown function. We report that male germline-specific expression of an antisense transcript contained in an intron of C. elegans Target of Rapamycin (TOR, let-363) is associated with (1) accumulation of endo-small interfering RNAs (siRNAs) against an embedded Helitron transposon and (2) activation of an alternative 3' splice site of TOR. The germline-specific Argonaute proteins PRG-1 and CSR-1, which participate in self/nonself RNA recognition, antagonistically regulate the generation of these endo-siRNAs, TOR mRNA levels, and 3' splice-site selection. Supply of exogenous double-stranded RNA against the region of sense/antisense overlap reverses changes in TOR expression and splicing and suppresses the progressive multigenerational sterility phenotype of prg-1 mutants. We propose that recognition of a "nonself" intronic transposon by endo-siRNAs/the piRNA system provides physiological regulation of expression and alternative splicing of a host gene that, in turn, contributes to the maintenance of germline function across generations.

EXPRSS: an Illumina based high-throughput expression-profiling method to reveal transcriptional dynamics

Background

Next Generation Sequencing technologies have facilitated differential gene expression analysis through RNA-seq and Tag-seq methods. RNA-seq has biases associated with transcript lengths, lacks uniform coverage of regions in mRNA and requires 10–20 times more reads than a typical Tag-seq. Most existing Tag-seq methods either have biases or not high throughput due to use of restriction enzymes or enzymatic manipulation of 5’ ends of mRNA or use of RNA ligations.

Results

We have developed EXpression Profiling through Randomly Sheared cDNA tag Sequencing (EXPRSS) that employs acoustic waves to randomly shear cDNA and generate sequence tags at a relatively defined position (~150-200 bp) from the 3′ end of each mRNA. Implementation of the method was verified through comparative analysis of expression data generated from EXPRSS, NlaIII-DGE and Affymetrix microarray and through qPCR quantification of selected genes. EXPRSS is a strand specific and restriction enzyme independent tag sequencing method that does not require cDNA length-based data transformations. EXPRSS is highly reproducible, is high-throughput and it also reveals alternative polyadenylation and polyadenylated antisense transcripts. It is cost-effective using barcoded multiplexing, avoids the biases of existing SAGE and derivative methods and can reveal polyadenylation position from paired-end sequencing.

Conclusions

EXPRSS Tag-seq provides sensitive and reliable gene expression data and enables high-throughput expression profiling with relatively simple downstream analysis.

Electronic supplementary material

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

Natural antisense transcripts and long non-coding RNA in Neurospora crassa

The prevalence of long non-coding RNAs (lncRNA) and natural antisense transcripts (NATs) has been reported in a variety of organisms. While a consensus has yet to be reached on their global importance, an increasing number of examples have been shown to be functional, regulating gene expression at the transcriptional and post-transcriptional level. Here, we use RNA sequencing data from the ABI SOLiD platform to identify lncRNA and NATs obtained from samples of the filamentous fungus Neurospora crassa grown under different light and temperature conditions. We identify 939 novel lncRNAs, of which 477 are antisense to annotated genes. Across the whole dataset, the extent of overlap between sense and antisense transcripts is large: 371 sense/antisense transcripts are complementary over 500 nts or more and 236 overlap by more than 1000 nts. Most prevalent are 3′ end overlaps between convergently transcribed sense/antisense pairs, but examples of divergently transcribed pairs and nested transcripts are also present. We confirm the expression of a subset of sense/antisense transcript pairs by qPCR. We examine the size, types of overlap and expression levels under the different environmental stimuli of light and temperature, and identify 11 lncRNAs that are up-regulated in response to light. We also find differences in transcript length and the position of introns between protein-coding transcripts that have antisense expression and transcripts with no antisense expression. These results demonstrate the ability of N. crassa lncRNAs and NATs to be regulated by different environmental stimuli and provide the scope for further investigation into the function of NATs.

Characterization and identification of cis-regulatory elements in Arabidopsis based on single-nucleotide polymorphism information

Identifying regulatory elements and revealing their role in gene expression regulation remains a central goal of plant genome research. We exploited the detailed genomic sequencing information of a large number of Arabidopsis (Arabidopsis thaliana) accessions to characterize known and to identify novel cis-regulatory elements in gene promoter regions of Arabidopsis by relying on conservation as the hallmark signal of functional relevance. Based on the genomic layout and the obtained density profiles of single-nucleotide polymorphisms (SNPs) in sequence regions upstream of transcription start sites, the average length of promoter regions in Arabidopsis could be established at 500 bp. Genes associated with high degrees of variability of their respective upstream regions are preferentially involved in environmental response and signaling processes, while low levels of promoter SNP density are common among housekeeping genes. Known cis-elements were found to exhibit a decreased SNP density than sequence regions not associated with known motifs. For 15 known cis-element motifs, strong positional preferences relative to the transcription start site were detected based on their promoter SNP density profiles. Five novel candidate cis-element motifs were identified as consensus motifs of 17 sequence hexamers exhibiting increased sequence conservation combined with evidence of positional preferences, annotation information, and functional relevance for inducing correlated gene expression. Our study demonstrates that the currently available resolution of SNP data offers novel ways for the identification of functional genomic elements and the characterization of gene promoter sequences.

Long noncoding RNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana

Short noncoding RNAs have been demonstrated to play important roles in regulation of gene expression and stress responses, but the repertoire and functions of long noncoding RNAs (lncRNAs) remain largely unexplored, particularly in plants. To explore the role of lncRNAs in disease resistance, we used a strand-specific RNA-sequencing approach to identify lncRNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana. Antisense transcription was found in c. 20% of the annotated A. thaliana genes. Several noncoding natural antisense transcripts responsive to F. oxysporum infection were found in genes implicated in disease defense. While the majority of the novel transcriptionally active regions (TARs) were adjacent to annotated genes and could be an extension of the annotated transcripts, 159 novel intergenic TARs, including 20 F. oxysporum-responsive lncTARs, were identified. Ten F. oxysporum-induced lncTARs were functionally characterized using T-DNA insertion or RNA-interference knockdown lines, and five were demonstrated to be related to disease development. Promoter analysis suggests that some of the F. oxysporum-induced lncTARs are direct targets of transcription factor(s) responsive to pathogen attack. Our results demonstrated that strand-specific RNA sequencing is a powerful tool for uncovering hidden levels of transcriptome and that IncRNAs are important components of the antifungal networks in A. thaliana.

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.

Gene silencing in plants: a diversity of pathways

Eukaryotic organisms have evolved a variety of gene silencing pathways in which small RNAs, 20- to 30-nucleotides in length, repress the expression of sequence homologous genes at the transcriptional or post-transcriptional levels. In plants, RNA silencing pathways play important roles in regulating development and response to both biotic and abiotic stresses. The molecular basis of these complex and interconnected pathways has emerged only in recent years with the identification of many of the genes necessary for the biogenesis and action of small RNAs. This review covers the diversity of RNA silencing pathways identified in plants.

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.

Natural antisense transcripts in plants: a review and identification in soybean infected with Phakopsora pachyrhizi SuperSAGE library

Natural antisense ranscripts (NAT) are RNA molecules complementary to other endogenous RNAs. They are capable of regulating the expression of target genes at different levels (transcription, mRNA stability, translation, etc.). Such a property makes them ideal for interventions in organisms' metabolism. The present study reviewed plant NAT aspects, including features, availability and genesis, conservation and distribution, coding capacity, NAT pair expression, and functions. Besides, an in silico identification of NATs pairs was presented, using deepSuperSAGE libraries of soybean infected or not with Phakopsora pachyrhizi. Results showed that around 1/3 of the 77,903 predicted trans-NATs (by PlantsNATsDB database) detected had unitags mapped in both sequences of each pair. The same 1/3 of the 436 foreseen cis-NATs showed unitags anchored in both sequences of the related pairs. For those unitags mapped in NAT pairs, a modulation expression was assigned as upregulated, downregulated, or constitutive, based on the statistical analysis (P < 0.05). As a result, the infected treatment promoted the expression of 2,313 trans-NATs pairs comprising unitags exclusively from that library (1,326 pairs had unitags only found in the mock library). To understand the regulation of these NAT pairs could be a key aspect in the ASR plant response.

Protein-coding cis-natural antisense transcripts have high and broad expression in Arabidopsis

Pairs of genes within eukaryotic genomes are often located on opposite DNA strands such that transcription generates cis-natural sense antisense transcripts (cis-NATs). This orientation of genes has been associated with the biogenesis of splice variants and natural antisense small RNAs. Here, in an analysis of currently available data, we report that within Arabidopsis (Arabidopsis thaliana), protein-coding cis-NATs are also characterized by high abundance, high coexpression, and broad expression. Our results suggest that a permissive chromatin environment may have led to the proximity of these genes. Compared with other genes, cis-NAT-encoding genes have enriched low-nucleosome-density regions, high levels of histone H3 lysine-9 acetylation, and low levels of H3 lysine-27 trimethylation. Promoters associated with broadly expressed genes are preferentially found in the 5' regulatory sequences of cis-NAT-encoding genes. Our results further suggest that natural antisense small RNA production from cis-NATs is limited. Small RNAs sequenced from natural antisense small RNA biogenesis mutants including dcl1, dcl2, dcl3, and rdr6 map to cis-NATs as frequently as small RNAs sequenced from wild-type plants. Future work will investigate if the positive transcriptional regulation of overlapping protein-coding genes contributes to the prevalence of these genes within other eukaryotic genomes.

Natural antisense transcripts of Trifolium repens dehydrins

The recently described complex nature of some dehydrin-coding sequences in Trifolium repens could explain the considerable variability among transcripts originating from a single gene.1 For some of the sequences the existence of natural antisense transcripts (NAT s), which could form sense-antisense (SAS) pairs, was predicted. The present study demonstrates that cis-natural antisense transcripts of 2 dehydrin types (YnKn and YnSKn) accumulate in white clover plants subjected to treatments with polyethylene glycol (PEG), abscisic acid (ABA), and high salt concentration. The isolated YnKn cis-NAT s mapped to sequence site enriched in alternative start codons. Some of the sense-antisense pairs exhibited inverse expression with differing profiles which depended on the applied stress. A natural antisense transcript coding for an ABC F family protein (a trans-NAT ) which shares short sequence homology with YnSKn dehydrin was identified in plants subjected to salt stress. Forthcoming experiments will evaluate the impact of NAT s on transcript abundances, elucidating the role of transcriptional and post-transcriptional interferences in the regulation of dehydrin levels under various abiotic stresses.

Transcriptome-wide characterization of miRNA-directed and non-miRNA-directed endonucleolytic cleavage using Degradome analysis under low ambient temperature in Phalaenopsis aphrodite subsp. formosana

Plant microRNAs (miRNAs) regulate gene expression through post-transcriptional gene silencing. Phalaenopsis aphrodite subsp. formosana is an orchid species native to Taiwan, which has high economic value and a high frequency of floral polymorphism. To date, few studies have focused on the regulatory roles of miRNAs and functional small RNAs (sRNAs) in orchids although understanding the regulation of flower development and flowering time is potentially important. Here, we combined analyses of the transcriptome, sRNAs and the degradome to identify sRNA-directed transcript cleavages in Phalaenopsis. Degradome analysis provided large-scale evidence of conserved and novel miRNA-directed cleavage of target transcripts, and 46 abundant sRNA groups and their target transcripts were identified. Low temperature-responsive sRNAs were validated with normalized reads from an sRNA library and quantitative stem-loop reverse transcription-PCR (RT-PCR) analysis. According to gene ontology (GO) categorization, target transcripts of the novel miRNAs and sRNAs are functionally involved in metabolic processes or responses to stress. One particular homologous gene, Allcontig28452, which encodes digalactosyldiacylglycerol synthase 2 (DGD2), was found to be targeted by natural antisense transcripts (NATs) unique to Phalaenopsis. In summary, comprehensive analyses of the transcriptome, sRNAs and degradome using deep sequencing technology provided a useful platform for investigating miRNA-directed and non-miRNA-directed endonucleolytic cleavage in a non-model plant, the orchid Phalaenopsis.

Development of a rapid and inexpensive method to reveal natural antisense transcripts

BACKGROUND

Natural antisense transcripts (NATs) are a group of RNAs encoded within a cell that have transcript complementarity to other RNA transcripts. NATs have been identified in multiple eukaryotes, including humans, mice, yeast and several plants, and are known to play crucial roles in gene regulation and modification via RNA interference, alternative splicing and genomic imprinting. NATs are also involved in several human diseases.

RESULTS

We describe a novel method to detect the occurrence of target NATs in specific plant tissues. This method differs from the others currently used in molecular biology laboratories for a number of reasons, particularly the simplicity and versatility of application, low cost and lower material requirement. We demonstrate that NATs can be detected by using diluted cDNA, avoiding the need for a large amount of RNA, thus differing from basic techniques, such as Northern blot hybridisation and reverse-transcription PCR amplification. Furthermore, our method also allows the precise detection of long NATs and their cloning into plasmid vectors for downstream applications. We also reported the first case of a tissue-specific NAT occurring in Oleaceae family and, the antisense orientation of this transcript, allows the splicing of two introns otherwise impossible in the sense orientation.

CONCLUSIONS

This method is the first that combines the polymerisation and cleavage activity of DNA polymerase and exonuclease enzymes, respectively, to discover NATs in living organisms. It may simplify the discovery of NATs in plants providing a new strategy for an easy identification and characterization of this group of RNA molecules. Furthermore, since NATs are found in multiple eukaryotes, our method can be easily applied to a wide range of organisms, including human, mice and yeast.

Direct sequencing of Arabidopsis thaliana RNA reveals patterns of cleavage and polyadenylation

It has recently been shown that RNA 3′ end formation plays a more widespread role in controlling gene expression than previously thought. In order to examine the impact of regulated 3′ end formation genome-wide we applied direct RNA sequencing to A. thaliana. Here we show the authentic transcriptome in unprecedented detail and how 3′ end formation impacts genome organization. We reveal extreme heterogeneity in RNA 3′ ends, discover previously unrecognized non-coding RNAs and propose widespread re-annotation of the genome. We explain the origin of most poly(A)+ antisense RNAs and identify cis-elements that control 3′ end formation in different registers. These findings are essential to understand what the genome actually encodes, how it is organized and the impact of regulated 3′ end formation on these processes.

Flowering time control: another window to the connection between antisense RNA and chromatin

A high proportion of all eukaryotic genes express antisense RNA (asRNA), which accumulates to varying degrees at different loci. Whether there is a general function for asRNA is unknown, but its widespread occurrence and frequent regulation by stress suggest an important role. The best-characterized plant gene exhibiting a complex antisense transcript pattern is the Arabidopsis floral regulator FLOWERING LOCUS C (FLC). Changes occur in the accumulation, splicing, and polyadenylation of this antisense transcript, termed COOLAIR, in different environments and genotypes. These changes are associated with altered chromatin regulation and differential FLC expression, provoking mechanistic comparisons with many well-studied loci in yeast and mammals. Detailed analysis of these specific examples may shed light on the complex interplay between asRNA and chromatin modifications in different genomes.

Identification of sense and antisense transcripts regulated by drought in sugarcane

Sugarcane is an important sugar and energy crop that can be used efficiently for biofuels production. The development of sugarcane cultivars tolerant to drought could allow for the expansion of plantations to sub-prime regions. Knowledge on the mechanisms underlying drought responses and its relationship with carbon partition would greatly help to define routes to increase yield. In this work we studied sugarcane responses to drought using a custom designed oligonucleotide array with 21,901 different probes. The oligoarrays were designed to contain probes that detect transcription in both sense and antisense orientation. We validated the results obtained using quantitative real-time PCR (qPCR). A total of 987 genes were differentially expressed in at least one sample of sugarcane plants submitted to drought for 24, 72 and 120 h. Among them, 928 were sense transcripts and 59 were antisense transcripts. Genes related to Carbohydrate Metabolism, RNA Metabolism and Signal Transduction were selected for gene expression validation by qPCR that indicated a validation percentage of 90%. From the probes presented on the array, 75% of the sense probes and 11.9% of the antisense probes have signal above background and can be classified as expressed sequences. Our custom sugarcane oligonucleotide array provides sensitivity and good coverage of sugarcane transcripts for the identification of a representative proportion of natural antisense transcripts (NATs) and sense-antisense transcript pairs (SATs). The antisense transcriptome showed, in most cases, co-expression with respective sense transcripts.

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.

Unexpected Diversity of Chloroplast Noncoding RNAs as Revealed by Deep Sequencing of the Arabidopsis Transcriptome

Noncoding RNAs (ncRNA) are widely expressed in both prokaryotes and eukaryotes. Eukaryotic ncRNAs are commonly micro- and small-interfering RNAs (18-25 nt) involved in posttranscriptional gene silencing, whereas prokaryotic ncRNAs vary in size and are involved in various aspects of gene regulation. Given the prokaryotic origin of organelles, the presence of ncRNAs might be expected; however, the full spectrum of organellar ncRNAs has not been determined systematically. Here, strand-specific RNA-Seq analysis was used to identify 107 candidate ncRNAs from Arabidopsis thaliana chloroplasts, primarily encoded opposite protein-coding and tRNA genes. Forty-eight ncRNAs were shown to accumulate by RNA gel blot as discrete transcripts in wild-type (WT) plants and/or the pnp1-1 mutant, which lacks the chloroplast ribonuclease polynucleotide phosphorylase (cpPNPase). Ninety-eight percent of the ncRNAs detected by RNA gel blot had different transcript patterns between WT and pnp1-1, suggesting cpPNPase has a significant role in chloroplast ncRNA biogenesis and accumulation. Analysis of materials deficient for other major chloroplast ribonucleases, RNase R, RNase E, and RNase J, showed differential effects on ncRNA accumulation and/or form, suggesting specificity in RNase-ncRNA interactions. 5' end mapping demonstrates that some ncRNAs are transcribed from dedicated promoters, whereas others result from transcriptional read-through. Finally, correlations between accumulation of some ncRNAs and the symmetrically transcribed sense RNA are consistent with a role in RNA stability. Overall, our data suggest that this extensive population of ncRNAs has the potential to underpin a previously underappreciated regulatory mode in the chloroplast.

Double-stranded RNA binding proteins DRB2 and DRB4 have an antagonistic impact on polymerase IV-dependent siRNA levels in Arabidopsis

Biogenesis of the vast majority of plant siRNAs depends on the activity of the plant-specific RNA polymerase IV (PolIV) enzyme. As part of the RNA-dependent DNA methylation (RdDM) process, PolIV-dependent siRNAs (p4-siRNAs) are loaded onto an ARGONAUTE4-containing complex and guide de novo DNA methyltransferases to target loci. Here we show that the double-stranded RNA binding proteins DRB2 and DRB4 are required for proper accumulation of p4-siRNAs. In flowers, loss of DRB2 results in increased accumulation of p4-siRNAs but not ta-siRNAs, inverted repeat (IR)-derived siRNAs, or miRNA. Loss of DRB2 does not impair uniparental expression of p4-dependent siRNAs in developing endosperm, indicating that p4-siRNA increased accumulation is not the result of the activation of the polIV pathway in the male gametophyte. In contrast to drb2, drb4 mutants exhibit reduced p4-siRNA levels, but the extent of this reduction is variable, according to the nature and size of the p4-siRNAs. Loss of DRB4 also leads to a spectacular increase of p4-independent IR-derived 24-nt siRNAs, suggesting a reallocation of factors from p4-dependent to p4-independent siRNA pathways in drb4. Opposite effects of drb2 and drb4 mutations on the accumulation of p4-siRNAs were also observed in vegetative tissues. Moreover, transgenic plants overexpressing DRB2 mimicked drb4 mutants at the morphological and molecular levels, confirming the antagonistic roles of DRB2 and DRB4.

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.

Genome-wide landscape of polyadenylation in Arabidopsis provides evidence for extensive alternative polyadenylation

Alternative polyadenylation (APA) has been shown to play an important role in gene expression regulation in animals and plants. However, the extent of sense and antisense APA at the genome level is not known. We developed a deep-sequencing protocol that queries the junctions of 3'UTR and poly(A) tails and confidently maps the poly(A) tags to the annotated genome. The results of this mapping show that 70% of Arabidopsis genes use more than one poly(A) site, excluding microheterogeneity. Analysis of the poly(A) tags reveal extensive APA in introns and coding sequences, results of which can significantly alter transcript sequences and their encoding proteins. Although the interplay of intron splicing and polyadenylation potentially defines poly(A) site uses in introns, the polyadenylation signals leading to the use of CDS protein-coding region poly(A) sites are distinct from the rest of the genome. Interestingly, a large number of poly(A) sites correspond to putative antisense transcripts that overlap with the promoter of the associated sense transcript, a mode previously demonstrated to regulate sense gene expression. Our results suggest that APA plays a far greater role in gene expression in plants than previously expected.

Evolutionary origins of Brassicaceae specific genes in Arabidopsis thaliana

Background

All sequenced genomes contain a proportion of lineage-specific genes, which exhibit no sequence similarity to any genes outside the lineage. Despite their prevalence, the origins and functions of most lineage-specific genes remain largely unknown. As more genomes are sequenced opportunities for understanding evolutionary origins and functions of lineage-specific genes are increasing.

Results

This study provides a comprehensive analysis of the origins of lineage-specific genes (LSGs) in Arabidopsis thaliana that are restricted to the Brassicaceae family. In this study, lineage-specific genes within the nuclear (1761 genes) and mitochondrial (28 genes) genomes are identified. The evolutionary origins of two thirds of the lineage-specific genes within the Arabidopsis thaliana genome are also identified. Almost a quarter of lineage-specific genes originate from non-lineage-specific paralogs, while the origins of ~10% of lineage-specific genes are partly derived from DNA exapted from transposable elements (twice the proportion observed for non-lineage-specific genes). Lineage-specific genes are also enriched in genes that have overlapping CDS, which is consistent with such novel genes arising from overprinting. Over half of the subset of the 958 lineage-specific genes found only in Arabidopsis thaliana have alignments to intergenic regions in Arabidopsis lyrata, consistent with either de novo origination or differential gene loss and retention, with both evolutionary scenarios explaining the lineage-specific status of these genes. A smaller number of lineage-specific genes with an incomplete open reading frame across different Arabidopsis thaliana accessions are further identified as accession-specific genes, most likely of recent origin in Arabidopsis thaliana. Putative de novo origination for two of the Arabidopsis thaliana-only genes is identified via additional sequencing across accessions of Arabidopsis thaliana and closely related sister species lineages. We demonstrate that lineage-specific genes have high tissue specificity and low expression levels across multiple tissues and developmental stages. Finally, stress responsiveness is identified as a distinct feature of Brassicaceae-specific genes; where these LSGs are enriched for genes responsive to a wide range of abiotic stresses.

Conclusion

Improving our understanding of the origins of lineage-specific genes is key to gaining insights regarding how novel genes can arise and acquire functionality in different lineages. This study comprehensively identifies all of the Brassicaceae-specific genes in Arabidopsis thaliana and identifies how the majority of such lineage-specific genes have arisen. The analysis allows the relative importance (and prevalence) of different evolutionary routes to the genesis of novel ORFs within lineages to be assessed. Insights regarding the functional roles of lineage-specific genes are further advanced through identification of enrichment for stress responsiveness in lineage-specific genes, highlighting their likely importance for environmental adaptation strategies.

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

The Arabidopsis genome contains a large number of gene pairs that encode sense and antisense transcripts with overlapping 3′ regions, indicative for a potential role of natural antisense transcription in regulating sense gene expression or transcript processing. When we mapped poly(A) transcripts of three plant gene pairs with long overlapping antisense transcripts, we identified an unusual transcript composition for two of the three gene pairs. Both genes pairs encoded a class of long sense transcripts and a class of short sense transcripts that terminate within the same polyadenylation region as the antisense transcripts encoded by the opposite strand. We find that the presence of the short sense transcript was not dependent on the expression of an antisense transcript. This argues against the assumption that the common termination region for sense and antisense poly(A) transcripts is the result of antisense-specific regulation. We speculate that for some genes evolution may have especially favoured alternative polyadenylation events that shorten transcript length for gene pairs with overlapping sense/antisense transcription, if this reduces the likelihood for dsRNA formation and transcript degradation.

RNA-Seq analysis of splicing in Plasmodium falciparum uncovers new splice junctions, alternative splicing and splicing of antisense transcripts

Over 50% of genes in Plasmodium falciparum, the deadliest human malaria parasite, contain predicted introns, yet experimental characterization of splicing in this organism remains incomplete. We present here a transcriptome-wide characterization of intraerythrocytic splicing events, as captured by RNA-Seq data from four timepoints of a single highly synchronous culture. Gene model-independent analysis of these data in conjunction with publically available RNA-Seq data with HMMSplicer, an in-house developed splice site detection algorithm, revealed a total of 977 new 5' GU-AG 3' and 5 new 5' GC-AG 3' junctions absent from gene models and ESTs (11% increase to the current annotation). In addition, 310 alternative splicing events were detected in 254 (4.5%) genes, most of which truncate open reading frames. Splicing events antisense to gene models were also detected, revealing complex transcriptional arrangements within the parasite's transcriptome. Interestingly, antisense introns overlap sense introns more than would be expected by chance, perhaps indicating a functional relationship between overlapping transcripts or an inherent organizational property of the transcriptome. Independent experimental validation confirmed over 30 new antisense and alternative junctions. Thus, this largest assemblage of new and alternative splicing events to date in Plasmodium falciparum provides a more precise, dynamic view of the parasite's transcriptome.

Long Nonprotein-Coding RNAs in Plants

In recent years, nonprotein-coding RNAs (or npcRNAs) have emerged as a major part of the eukaryotic transcriptome. Many new regulatory npcRNAs or riboregulators riboregulators have been discovered and characterized due to the advent of new genomic approaches. This growing number suggests that npcRNAs could play a more important role than previously believed and significantly contribute to the generation of evolutionary complexity in multicellular organisms. Regulatory npcRNAs range from small RNAs (si/miRNAs) to very large transcripts (or long npcRNAs) and play diverse functions in development and/or environmental stress responses. Small RNAs include an expanding number of 20-40 nt RNAs that function in the regulation of gene expression by affecting mRNA decay and translational inhibition or lead to DNA methylation and gene silencing. They generally involve double-stranded RNA or stem loops and imply transcriptional or posttranscriptional gene silencing (PTGS). RNA silencing besides small interfering RNA and microRNA, gene silencing in plants is also mediated by tasiRNAs (trans-acting siRNAs) and nat-siRNAs (natural antisense mediated siRNAs). In contrast to small RNAs, much less is known about the large and diverse population of long npcRNAs, and only a few have been implicated in diverse functions such as abiotic stress responses, nodulation and flower development, and sex chromosome-specific expression. Moreover, many long npcRNAs act as antisense transcripts or are substrates of the small RNA pathways, thus interfering with a variety of RNA-related metabolisms. An emerging hypothesis is that long npcRNAs, as shown for small si/miRNAs, integrate into ribonucleoprotein particles (RNPs) to modulate their function, localization, or stability to act on target mRNAs. As plants show a remarkable developmental plasticity to adapt their growth to changing environmental conditions, understanding how npcRNAs work may reveal novel mechanisms involved in growth control and differentiation and help to design new tools for biotechnological applications.

The Functions of RNA-Dependent RNA Polymerases in Arabidopsis

One recently identified mechanism that regulates mRNA abundance is RNA silencing, and pioneering work in Arabidopsis thaliana and other genetic model organisms helped define this process. RNA silencing pathways are triggered by either self-complementary fold-back structures or the production of double-stranded RNA (dsRNA) that gives rise to small RNAs (smRNAs) known as microRNAs (miRNAs) or small-interfering RNAs (siRNAs). These smRNAs direct sequence-specific regulation of various gene transcripts, repetitive sequences, viruses, and mobile elements via RNA cleavage, translational inhibition, or transcriptional silencing through DNA methylation and heterochromatin formation. Early genetic screens in Arabidopsis were instrumental in uncovering numerous proteins required for these important regulatory pathways. Among the factors identified by these studies were RNA-dependent RNA polymerases (RDRs), which are proteins that synthesize siRNA-producing dsRNA molecules using a single-stranded RNA (ssRNA) molecule as a template. Recently, a growing body of evidence has implicated RDR-dependent RNA silencing in many different aspects of plant biology ranging from reproductive development to pathogen resistance. Here, we focus on the specific functions of the six Arabidopsis RDRs in RNA silencing, their ssRNA substrates and resulting RDR-dependent smRNAs, and the numerous biological functions of these proteins in plant development and stress responses.

Abscisic acid response element binding factor 1 is required for establishment of Arabidopsis seedlings during winter

Abscisic acid (ABA) plays a crucial role in abiotic stress response apart from its influence on growth and development of a plant. Our studies on abscisic acid response element binding factor 1 (ABF1) gene in Arabidopsis demonstrate that it is required for seedling establishment during winter. ABF1 is also involved in regulating seed dormancy and seed germination to some extent. Analysis of transcriptional activity of ABF1 promoter reveals that ABF1 expresses specifically in trichomes of young leaves and constitutively in cotyledons, roots, older leaves and flowers. The expression is induced upon exposure to ABA, cold and heat. The alignment of cDNAs of ABF1 (At1g49720) and At1g49730 (encodes a protein kinase of unknown function), reveals an overlap of 88 bp at their 3' UTR region suggesting that they can potentially form natural cis-antisense mRNAs pair in a tail-to-tail manner. Analysis by Genevestigator microarray stress response viewer further supports the regulatory role of these genes. An inverse proportion is observed in the transcription the two loci in number of stress responses. The abf1 mutants do not show any seedling establishment defects when grown under standard growth conditions. The mutant seedlings exhibit growth defects during winter in the western Himalayan region. Our study also signifies the importance of functional analysis for mutant phenotypes in natural habitats by reverse genetic approaches, in order to identify specific function of particular gene/s whose expression level is altered upon exposure to changes in environmental cues such as temperature and light.

Alternative polyadenylation of antisense RNAs and flowering time control

Flowering time is controlled by precision in gene regulation mediated by different pathways. Two Arabidopsis thaliana components of the autonomous flowering pathway, FCA and FPA, function as genetically independent trans-acting regulators of alternative cleavage and polyadenylation. FCA and FPA directly associate with chromatin at the locus encoding the floral repressor FLC, but appear to control FLC transcription by mediating alternative polyadenylation of embedded non-coding antisense RNAs. These findings prompt the re-examination of how other factors control FLC expression, as it is formally possible that they function primarily to control alternative processing of antisense RNAs. As co-expressed sense and antisense gene pairs are widespread in eukaryotes, alternative processing of antisense RNAs may represent a significant form of gene regulation.