Noncoding RNA for CR20, a cytokinin-repressed gene of cucumber

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

A stable tRNA-like molecule is generated from the long noncoding RNA GUT15 in Arabidopsis

The Arabidopsis GUT15 RNA belongs to a class of noncoding RNAs that are expressed from the intergenic regions of protein-coding genes. We show that the RNA polymerase II transcribed GUT15 transcript serves as a precursor for two stable RNA species, a tRNA-like molecule and GUT15-tRF-F5, which are both encoded by the final intron in the GUT15 gene. The GUT15-encoded tRNA-like molecule cannot be autonomously transcribed by RNA polymerase III. However, this molecule contains a CCA motif, suggesting that it may enter the tRNA maturation pathway. The GUT15-encoded tRNA-like sequence has an inhibiting effect on the splicing of its host intron. Moreover, we demonstrate that the canonical tRNA genes nested within introns do not affect the splicing patterns of their host protein-coding transcripts.

YeATS - a tool suite for analyzing RNA-seq derived transcriptome identifies a highly transcribed putative extensin in heartwood/sapwood transition zone in black walnut

The transcriptome provides a functional footprint of the genome by enumerating the molecular components of cells and tissues. The field of transcript discovery has been revolutionized through high-throughput mRNA sequencing (RNA-seq). Here, we present a methodology that replicates and improves existing methodologies, and implements a workflow for error estimation and correction followed by genome annotation and transcript abundance estimation for RNA-seq derived transcriptome sequences (YeATS - Yet Another Tool Suite for analyzing RNA-seq derived transcriptome). A unique feature of YeATS is the upfront determination of the errors in the sequencing or transcript assembly process by analyzing open reading frames of transcripts. YeATS identifies transcripts that have not been merged, result in broken open reading frames or contain long repeats as erroneous transcripts. We present the YeATS workflow using a representative sample of the transcriptome from the tissue at the heartwood/sapwood transition zone in black walnut. A novel feature of the transcriptome that emerged from our analysis was the identification of a highly abundant transcript that had no known homologous genes (GenBank accession: KT023102). The amino acid composition of the longest open reading frame of this gene classifies this as a putative extensin. Also, we corroborated the transcriptional abundance of proline-rich proteins, dehydrins, senescence-associated proteins, and the DNAJ family of chaperone proteins. Thus, YeATS presents a workflow for analyzing RNA-seq data with several innovative features that differentiate it from existing software.

YeATS - a tool suite for analyzing RNA-seq derived transcriptome identifies a highly transcribed putative extensin in heartwood/sapwood transition zone in black walnut

The transcriptome provides a functional footprint of the genome by enumerating the molecular components of cells and tissues. The field of transcript discovery has been revolutionized through high-throughput mRNA sequencing (RNA-seq). Here, we present a methodology that replicates and improves existing methodologies, and implements a workflow for error estimation and correction followed by genome annotation and transcript abundance estimation for RNA-seq derived transcriptome sequences (YeATS - Yet Another Tool Suite for analyzing RNA-seq derived transcriptome). A unique feature of YeATS is the upfront determination of the errors in the sequencing or transcript assembly process by analyzing open reading frames of transcripts. YeATS identifies transcripts that have not been merged, result in broken open reading frames or contain long repeats as erroneous transcripts. We present the YeATS workflow using a representative sample of the transcriptome from the tissue at the heartwood/sapwood transition zone in black walnut. A novel feature of the transcriptome that emerged from our analysis was the identification of a highly abundant transcript that had no known homologous genes (GenBank accession: KT023102). The amino acid composition of the longest open reading frame of this gene classifies this as a putative extensin. Also, we corroborated the transcriptional abundance of proline-rich proteins, dehydrins, senescence-associated proteins, and the DNAJ family of chaperone proteins. Thus, YeATS presents a workflow for analyzing RNA-seq data with several innovative features that differentiate it from existing software.

Identification and characterization of mRNA-like noncoding RNAs in Salvia miltiorrhiza

Identification and characterization of 5,446 mlncRNAs from Salvia miltiorrhiza showed that the majority of identified mlncRNAs were stress responsive, providing a framework for elucidating mlncRNA functions in S. miltiorrhiza. mRNA-like noncoding RNAs (mlncRNAs) are transcribed by RNA polymerase II and are polyadenylated, capped and spliced. They play important roles in plant development and defense responses. However, there is no information available for mlncRNAs in Salvia miltiorrhiza Bunge, the first Chinese medicinal material entering the international market. To perform a transcriptome-wide identification of S. miltiorrhiza mlncRNAs, we assembled over 8 million RNA-seq reads from GenBank database and 5,624 ESTs from PlantGDB into 44422 unigenes. Using a computational identification pipeline, we identified 5446 S. miltiorrhiza mlncRNA candidates from the assembled unigenes. Of the 5446 mlncRNAs, 2 are primary transcripts of conserved miRNAs, and 2030 can be grouped into 470 families with at least two members in a family. Quantitative real-time PCR analysis of mlncRNAs with at least 900 nt showed that the majority were differentially expressed in roots, stems, leaves and flowers and responsive to methyl jasmonate (MeJA) treatment in S. miltiorrhiza. Analysis of published RNA-seq data showed that a total of 3,044 mlncRNAs were expressed in hairy roots of S. miltiorrhiza and the expression of 1,904 of the 3,044 mlncRNAs was altered by yeast extract and Ag(+) treatment. The results indicate that the majority of mlncRNAs are involved in plant response to stress. This study provides a framework for understanding the roles of mlncRNAs in S. miltiorrhiza.

Genome-wide identification, characterization and evolutionary analysis of long intergenic noncoding RNAs in cucumber

Long intergenic noncoding RNAs (lincRNAs) are intergenic transcripts with a length of at least 200 nt that lack coding potential. Emerging evidence suggests that lincRNAs from animals participate in many fundamental biological processes. However, the systemic identification of lincRNAs has been undertaken in only a few plants. We chose to use cucumber (Cucumis sativus) as a model to analyze lincRNAs due to its importance as a model plant for studying sex differentiation and fruit development and the rich genomic and transcriptome data available. The application of a bioinformatics pipeline to multiple types of gene expression data resulted in the identification and characterization of 3,274 lincRNAs. Next, 10 lincRNAs targeted by 17 miRNAs were also explored. Based on co-expression analysis between lincRNAs and mRNAs, 94 lincRNAs were annotated, which may be involved in response to stimuli, multi-organism processes, reproduction, reproductive processes, and growth. Finally, examination of the evolution of lincRNAs showed that most lincRNAs are under purifying selection, while 16 lincRNAs are under natural selection. Our results provide a rich resource for further validation of cucumber lincRNAs and their function. The identification of lincRNAs targeted by miRNAs offers new clues for investigations into the role of lincRNAs in regulating gene expression. Finally, evaluation of the lincRNAs suggested that some lincRNAs are under positive and balancing selection.

Genome-wide identification and characterization of long intergenic non-coding RNAs in Ganoderma lucidum

Ganoderma lucidum is a white-rot fungus best-known for its medicinal activities. We have previously sequenced its genome and annotated the protein coding genes. However, long non-coding RNAs in G. lucidum genome have not been analyzed. In this study, we have identified and characterized long intergenic non-coding RNAs (lincRNA) in G. lucidum systematically. We developed a computational pipeline, which was used to analyze RNA-Seq data derived from G. lucidum samples collected from three developmental stages. A total of 402 lincRNA candidates were identified, with an average length of 609 bp. Analysis of their adjacent protein-coding genes (apcGenes) revealed that 46 apcGenes belong to the pathways of triterpenoid biosynthesis and lignin degradation, or families of cytochrome P450, mating type B genes, and carbohydrate-active enzymes. To determine if lincRNAs and these apcGenes have any interactions, the corresponding pairs of lincRNAs and apcGenes were analyzed in detail. We developed a modified 3' RACE method to analyze the transcriptional direction of a transcript. Among the 46 lincRNAs, 37 were found unidirectionally transcribed, and 9 were found bidirectionally transcribed. The expression profiles of 16 of these 37 lincRNAs were found to be highly correlated with those of the apcGenes across the three developmental stages. Among them, 11 are positively correlated (r>0.8) and 5 are negatively correlated (r<-0.8). The co-localization and co-expression of lincRNAs and those apcGenes playing important functions is consistent with the notion that lincRNAs might be important regulators for cellular processes. In summary, this represents the very first study to identify and characterize lincRNAs in the genomes of basidiomycetes. The results obtained here have laid the foundation for study of potential lincRNA-mediated expression regulation of genes in G. lucidum.

BcMF11, a novel non-coding RNA gene from Brassica campestris, is required for pollen development and male fertility

KEY MESSAGE : BcMF11 as a non-coding RNA gene has an essential role in pollen development, and might be useful for regulating the pollen fertility of crops by antisense RNA technology. We previously identified a 828-bp full-length cDNA of BcMF11, a novel pollen-specific non-coding mRNA-like gene from Chinese cabbage (Brassica campestris L. ssp. chinensis Makino). However, little information is known about the function of BcMF11 in pollen development. To investigate its exact biological roles in pollen development, the BcMF11 cDNA was antisense inhibited in transgenic Chinese cabbage under the control of a tapetum-specific promoter BcA9 and a constitutive promoter CaMV 35S. Antisense RNA transgenic plants displayed decreasing expression of BcMF11 and showed distinct morphological defects. Pollen germination test in vitro and in vivo of the transgenic plants suggested that inhibition of BcMF11 decreased pollen germination efficiency and delayed the pollen tubes' extension in the style. Under scanning electron microscopy, many shrunken and collapsed pollen grains were detected in the antisense BcMF11 transgenic Chinese cabbage. Further cytological observation revealed abnormal pollen development process in transgenic plants, including delayed degradation of tapetum, asynchronous separation of microspore, and aborted development of pollen grain. These results suggest that BcMF11, as a non-coding RNA, plays an essential role in pollen development and male fertility.

Analysis of antisense expression by whole genome tiling microarrays and siRNAs suggests mis-annotation of Arabidopsis orphan protein-coding genes

Background

MicroRNAs (miRNAs) and trans-acting small-interfering RNAs (tasi-RNAs) are small (20–22 nt long) RNAs (smRNAs) generated from hairpin secondary structures or antisense transcripts, respectively, that regulate gene expression by Watson-Crick pairing to a target mRNA and altering expression by mechanisms related to RNA interference. The high sequence homology of plant miRNAs to their targets has been the mainstay of miRNA prediction algorithms, which are limited in their predictive power for other kingdoms because miRNA complementarity is less conserved yet transitive processes (production of antisense smRNAs) are active in eukaryotes. We hypothesize that antisense transcription and associated smRNAs are biomarkers which can be computationally modeled for gene discovery.

Principal Findings

We explored rice (Oryza sativa) sense and antisense gene expression in publicly available whole genome tiling array transcriptome data and sequenced smRNA libraries (as well as C. elegans) and found evidence of transitivity of MIRNA genes similar to that found in Arabidopsis. Statistical analysis of antisense transcript abundances, presence of antisense ESTs, and association with smRNAs suggests several hundred Arabidopsis ‘orphan’ hypothetical genes are non-coding RNAs. Consistent with this hypothesis, we found novel Arabidopsis homologues of some MIRNA genes on the antisense strand of previously annotated protein-coding genes. A Support Vector Machine (SVM) was applied using thermodynamic energy of binding plus novel expression features of sense/antisense transcription topology and siRNA abundances to build a prediction model of miRNA targets. The SVM when trained on targets could predict the “ancient” (deeply conserved) class of validated Arabidopsis MIRNA genes with an accuracy of 84%, and 76% for “new” rapidly-evolving MIRNA genes.

Conclusions

Antisense and smRNA expression features and computational methods may identify novel MIRNA genes and other non-coding RNAs in plants and potentially other kingdoms, which can provide insight into antisense transcription, miRNA evolution, and post-transcriptional gene regulation.

Diamonds in the rough: mRNA-like non-coding RNAs

Non-coding RNAs are increasingly being identified as crucial regulators of gene expression and other cellular functions in plants. Experimental and computational methods have revealed the existence of mRNA-like non-coding RNAs (mlncRNAs), a class of non-coding RNAs that, in plants, are associated with tissue-specific expression, development and the phosphate-starvation response. Although their mechanisms of action are largely unknown, one can speculate that mlncRNAs act through secondary structures or specific sequences that bind to proteins or metabolites, or that have catalytic activity. This review summarizes the computational methods developed to identify candidate mlncRNAs, and the current experimental evidence regarding the function of several known mlncRNAs.

BcMF11, a putative pollen-specific non-coding RNA from Brassica campestris ssp. chinensis

A full-length cDNA, BcMF11, has been cloned from Chinese cabbage (Brassica campestris L. ssp. chinensis Makino) using rapid amplification of the cDNA ends (RACE) based on a pollen-specific cDNA fragment (DN237921). The BcMF11 cDNA has a total length of 828bp with poly (A) tail. Analysis of the sequence demonstrated that BcMF11 is a novel non-coding RNA which has no prominent open reading frame (ORF) or coding capacity. No significant similarities were observed between BcMF11 and previously published sequences in GenBank. Transcription analysis indicated that BcMF11 is a novel pollen-specific ncRNA and may be involved in the pollen development of Chinese cabbage.

Regulation of apoptosis by a prostate-specific and prostate cancer-associated noncoding gene, PCGEM1

PCGEM1 is a prostate tissue-specific, and prostate cancer-associated noncoding RNA (ncRNA) gene. Previous results revealed a significant association of elevated PCGEM1 expression levels in prostate cancer cells of African-American patients, whose mortality rate is the highest among prostate cancer patients. Functional study of PCGEM1 demonstrated a marked increase in colony formation in LNCaP prostate cancer cells and NIH3T3 mouse fibroblast cells. This study demonstrates that PCGEM1 overexpression in LNCaP cell culture model results in the inhibition of apoptosis induced by doxorubicin (DOX). Induction of p53 and p21(Waf1/Cip1) by DOX were delayed in LNCaP cells stably overexpressing PCGEM1 (LNCaP-PCGEM1 cells) compared to control LNCaP cells. The protein levels of cleaved caspase 7, and cleaved PARP were attenuated in DOXtreated LNCaP-PCGEM1 cells compared to control LNCaP cells. Similar results were observed in LNCaP cells transiently overexpressing PCGEM1. The inhibition of PARP cleavage by PCGEM1 overexpression was also observed in LNCaP-PCGEM1 cells incubated with etoposide and sodium selenite. Fluorescence-Activated Cell Sorter Annexin-V analysis revealed significantly lower percentage of apoptotic cells in DOX-treated LNCaP-PCGEM1 cells compared to control LNCaP cells. The attenuation of apoptic response appears to be androgen dependent in this experimental model, as androgen-independent variants of LNCaP cells did not exhibit this response. In summary, this study provides new insights into cell biologic functions and novel features of an ncRNA. Further, these data unravel biological mechanisms of cell growth/cell survival-associated functions of this ncRNA in a widely used prostate cancer cell culture model.

Genomic organization of the S core region and the S flanking regions of a class-II S haplotype in Brassica rapa

The nucleotide sequence of an 86.4-kb region that includes the SP11, SRK, and SLG genes of Brassica rapa S-60 (a class-II S haplotype) was determined. In the sequenced region, 13 putative genes were found besides SP11-60, SRK-60, and SLG-60. Five of these sequences were isolated as cDNAs, five were homologues of known genes, cDNAs, or ORFs, and three are hypothetical ORFs. Based on their nucleotide sequences, however, some of them are thought to be non-functional. Two regions of colinearity between the class-II S-60 and Brassica class-I S haplotypes were identified, i.e., S flanking region 1 which shows partial colinearity of non-genic sequences and S flanking region 2 which shows a high level of colinearity. The observed colinearity made it possible to compare the order of SP-11, SRK, and SLG genes in the S locus between the five sequenced S haplotypes. It emerged that the order of SRK and SLG in class-II S-60 is the reverse of that in the four class-I S haplotypes reported so far, and the order of SP11, SRK and SLG is the opposite of that in the class-I haplotype S-910. The possible gene designated as SAN1 (S locus Anther-expressed Non-coding RNA like-1), which is located in the region between SP11-60 and SRK-60, has features reminiscent of genes for non-coding RNAs (ncRNAs), but no homologous sequences were found in the databases. This sequence is transcribed in anthers but not in stigmas or leaves. These features of the genomic structure of S-60 are discussed with special reference to the characteristics of class-II S haplotypes.

Non-coding ribonucleic acids--a class of their own?

Non-coding ribonucleic acids (RNAs) do not contain a peptide-encoding open reading frame and are therefore not translated into proteins. They are expressed in all phyla, and in eukaryotic cells they are found in the nucleus, cytoplasm, and mitochondria. Non-coding RNAs either can exert structural functions, as do transfer and ribosomal RNAs, or they can regulate gene expression. Non-coding RNAs with regulatory functions differ in size ranging from a few nucleotides to over 100 kb and have diverse cell- or development-specific functions. Some of the non-coding RNAs associate with human diseases. This chapter summarizes the current knowledge about regulatory non-coding RNAs.

Beyond the proteome: non-coding regulatory RNAs

A variety of RNA molecules have been found over the last 20 years to have a remarkable range of functions beyond the well-known roles of messenger, ribosomal and transfer RNAs. Here, we present a general categorization of all non-coding RNAs and briefly discuss the ones that affect transcription, translation and protein function.

Identification and analysis of Arabidopsis expressed sequence tags characteristic of non-coding RNAs

Sequencing of the Arabidopsis genome has led to the identification of thousands of new putative genes based on the predicted proteins they encode. Genes encoding tRNAs, ribosomal RNAs, and small nucleolar RNAs have also been annotated; however, a potentially important class of genes has largely escaped previous annotation efforts. These genes correspond to RNAs that lack significant open reading frames and encode RNA as their final product. Accumulating evidence indicates that such "non-coding RNAs" (ncRNAs) can play critical roles in a wide range of cellular processes, including chromosomal silencing, transcriptional regulation, developmental control, and responses to stress. Approximately 15 putative Arabidopsis ncRNAs have been reported in the literature or have been annotated. Although several have homologs in other plant species, all appear to be plant specific, with the exception of signal recognition particle RNA. Conversely, none of the ncRNAs reported from yeast or animal systems have homologs in Arabidopsis or other plants. To identify additional genes that are likely to encode ncRNAs, we used computational tools to filter protein-coding genes from genes corresponding to 20,000 expressed sequence tag clones. Using this strategy, we identified 19 clones with characteristics of ncRNAs, nine putative peptide-coding RNAs with open reading frames smaller than 100 amino acids, and 11 that could not be differentiated between the two categories. Again, none of these clones had homologs outside the plant kingdom, suggesting that most Arabidopsis ncRNAs are likely plant specific. These data indicate that ncRNAs represent a significant and underdeveloped aspect of Arabidopsis genomics that deserves further study.

Non-coding RNAs: the architects of eukaryotic complexity

Around 98% of all transcriptional output in humans is non-coding RNA. RNA-mediated gene regulation is widespread in higher eukaryotes and complex genetic phenomena like RNA interference, co-suppression, transgene silencing, imprinting, methylation, and possibly position-effect variegation and transvection, all involve intersecting pathways based on or connected to RNA signaling. I suggest that the central dogma is incomplete, and that intronic and other non-coding RNAs have evolved to comprise a second tier of gene expression in eukaryotes, which enables the integration and networking of complex suites of gene activity. Although proteins are the fundamental effectors of cellular function, the basis of eukaryotic complexity and phenotypic variation may lie primarily in a control architecture composed of a highly parallel system of trans-acting RNAs that relay state information required for the coordination and modulation of gene expression, via chromatin remodeling, RNA-DNA, RNA-RNA and RNA-protein interactions. This system has interesting and perhaps informative analogies with small world networks and dataflow computing.

Translational and structural requirements of the early nodulin gene enod40, a short-open reading frame-containing RNA, for elicitation of a cell-specific growth response in the alfalfa root cortex

A diversity of mRNAs containing only short open reading frames (sORF-RNAs; encoding less than 30 amino acids) have been shown to be induced in growth and differentiation processes. The early nodulin gene enod40, coding for a 0.7-kb sORF-RNA, is expressed in the nodule primordium developing in the root cortex of leguminous plants after infection by symbiotic bacteria. Ballistic microtargeting of this gene into Medicago roots induced division of cortical cells. Translation of two sORFs (I and II, 13 and 27 amino acids, respectively) present in the conserved 5' and 3' regions of enod40 was required for this biological activity. These sORFs may be translated in roots via a reinitiation mechanism. In vitro translation products starting from the ATG of sORF I were detectable by mutating enod40 to yield peptides larger than 38 amino acids. Deletion of a Medicago truncatula enod40 region between the sORFs, spanning a predicted RNA structure, did not affect their translation but resulted in significantly decreased biological activity. Our data reveal a complex regulation of enod40 action, pointing to a role of sORF-encoded peptides and structured RNA signals in developmental processes involving sORF-RNAs.

Collection of mRNA-like non-coding RNAs

In last few years much data has accumulated which shows that in different cells various RNA transcripts are synthesized. They lack protein coding capacity and do not produce mature protein. It seems that they work mainly or exclusively on the RNA level. Their function and mechanism of action is poorly understood. In this paper we have collected all known RNA transcript and prepared a database for further structural and functional studies. This is the first collection of the nucleotide sequences of RNAs of this kind. The data can be accessed via WWW at: http://www.man.poznan.pl/5SData/ncRNA/inde x.html

Plants embrace a stepchild: the discovery of peptide growth regulators

Over the past decade, peptides have been added to the collection of signalling molecules in plants. As the impact of peptide hormones in non-plants is enormous, a comparison of plant and non-plant peptide signal molecules at this stage deserves our attention-not only to reveal common and unique features, but also to point to new avenues of future research on plant hormones.