Antisense transcripts with rice full-length cDNAs

Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

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.

The Functional Meaning of 5'UTR in Protein-Coding Genes

As it is well known, messenger RNA has many regulatory regions along its sequence length. One of them is the 5' untranslated region (5'UTR), which itself contains many regulatory elements such as upstream ORFs (uORFs), internal ribosome entry sites (IRESs), microRNA binding sites, and structural components involved in the regulation of mRNA stability, pre-mRNA splicing, and translation initiation. Activation of the alternative, more upstream transcription start site leads to an extension of 5'UTR. One of the consequences of 5'UTRs extension may be head-to-head gene overlap. This review describes elements in 5'UTR of protein-coding transcripts and the functional significance of protein-coding genes 5' overlap with implications for transcription, translation, and disease.

LncPheDB: a genome-wide lncRNAs regulated phenotypes database in plants

LncPheDB (https://www.lncphedb.com/) is a systematic resource of genome-wide long non-coding RNAs (lncRNAs)-phenotypes associations for multiple species. It was established to display the genome-wide lncRNA annotations, target genes prediction, variant-trait associations, gene-phenotype correlations, lncRNA-phenotype correlations, and the similar non-coding regions of the queried sequence in multiple species. LncPheDB sorted out a total of 203,391 lncRNA sequences, 2000 phenotypes, and 120,271 variants of nine species (Zea mays L., Gossypium barbadense L., Triticum aestivum L., Lycopersicon esculentum Mille, Oryza sativa L., Hordeum vulgare L., Sorghum bicolor L., Glycine max L., and Cucumis sativus L.). By exploring the relationship between lncRNAs and the genomic position of variants in genome-wide association analysis, a total of 68,862 lncRNAs were found to be related to the diversity of agronomic traits. More importantly, to facilitate the study of the functions of lncRNAs, we analyzed the possible target genes of lncRNAs, constructed a blast tool for performing similar fragmentation studies in all species, linked the pages of phenotypic studies related to lncRNAs that possess similar fragments and constructed their regulatory networks. In addition, LncPheDB also provides a user-friendly interface, a genome visualization platform, and multi-level and multi-modal convenient data search engine. We believe that LncPheDB plays a crucial role in mining lncRNA-related plant data.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-022-00084-3.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara)

Background

Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown.

Results

Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding–noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding–noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding–noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare.

Conclusion

We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07754-4.

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.

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.

Establishing a System for Functional Characterization of Full-Length cDNAs of Camellia sinensis

Tea (Camellia sinensis) is enriched with bioactive secondary metabolites, and is one of the most popular nonalcoholic beverages globally. Two tea reference genomes have been reported; however, the functional analysis of tea genes has lagged, mainly due to tea’s recalcitrance to genetic transformation and the absence of alternative high throughput heterologous expression systems. A full-length cDNA collection with a streamlined cloning system is needed in this economically important woody crop species. RNAs were isolated from nine different vegetative tea tissues, pooled, then used to construct a normalized full-length cDNA library. The titer of unamplified and amplified cDNA library was 6.89 × 10⁶ and 1.8 × 10¹⁰ cfu/mL, respectively; the library recombinant rate was 87.2%. Preliminary characterization demonstrated that this collection can complement existing tea reference genomes and facilitate rare gene discovery. In addition, to streamline tea cDNA cloning and functional analysis, a binary vector (pBIG2113SF) was reengineered, seven tea cDNAs isolated from this library were successfully cloned into this vector, then transformed into Arabidopsis. One FL-cDNA, which encodes a putative P_1B-type ATPase 5 (CsHMA5), was characterized further as a proof of concept. We demonstrated that overexpression of CsHMA5 in Arabidopsis resulted in copper hyposensitivity. Thus, our data demonstrated that this represents an efficient system for rare gene discovery and functional characterization of tea genes. The integration of a tea FL-cDNA collection with efficient cloning and a heterologous expression system would facilitate functional annotation and characterization of tea genes.

OverGeneDB: a database of 5' end protein coding overlapping genes in human and mouse genomes

Gene overlap plays various regulatory functions on transcriptional and post-transcriptional levels. Most current studies focus on protein-coding genes overlapping with non-protein-coding counterparts, the so called natural antisense transcripts. Considerably less is known about the role of gene overlap in the case of two protein-coding genes. Here, we provide OverGeneDB, a database of human and mouse 5′ end protein-coding overlapping genes. The database contains 582 human and 113 mouse gene pairs that are transcribed using overlapping promoters in at least one analyzed library. Gene pairs were identified based on the analysis of the transcription start site (TSS) coordinates in 73 human and 10 mouse organs, tissues and cell lines. Beside TSS data, resources for 26 human lung adenocarcinoma cell lines also contain RNA-Seq and ChIP-Seq data for seven histone modifications and RNA Polymerase II activity. The collected data revealed that the overlap region is rarely conserved between the studied species and tissues. In ∼50% of the overlapping genes, transcription started explicitly in the overlap regions. In the remaining half of overlapping genes, transcription was initiated both from overlapping and non-overlapping TSSs. OverGeneDB is accessible at http://overgenedb.amu.edu.pl.

Long Non-Coding RNAs: Rising Regulators of Plant Reproductive Development

Long non-coding RNAs (lncRNAs) exert a multitude of functions in regulating numerous biological processes. Recent studies have uncovered a growing number of lncRNAs within the plant genome. These molecules show striking tissue-specific expression patterns, suggesting that they exert regulatory functions in the growth and development processes of plants. Plant reproductive development is tightly regulated by both environmental and endogenous factors. As plant reproductive development is a crucial aspect of crop breeding, lncRNAs that modulate reproductive development are now particularly worth regarding. Here, we summarize findings that implicate lncRNAs in the control of plant reproductive development, especially in flowering control. Additionally, we elaborate on the regulation mechanisms of lncRNAs, tools for research on their function and mechanism, and potential directions of future research.

LncRNA expression profile and ceRNA analysis in tomato during flowering

Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs that play essential regulatory roles in various developmental processes and stress responses. However, the functions of lncRNAs during the flowering period of tomato are largely unknown. To explore the lncRNA profiles and functions during flowering in tomato, we performed strand-specific paired-end RNA sequencing of tomato leaves, flowers and roots, with three biological replicates. We identified 10919 lncRNAs including 248 novel lncRNAs, of which 65 novel lncRNAs were significantly differentially expressed (DE) in the flowers, leaves, and roots. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out to identify the cis target gene of DE lncRNAs. The results showed that the lncRNAs might play an important role in the growth, development, and apoptosis of flowering tomato plant by regulating the formation of intima in flower tissues, binding to various molecules, influencing metabolic pathways, and inducing apoptosis. Moreover, we identified the interaction between 32, 78, and 397 kinds of miRNAs, lncRNAs, and mRNAs. The results suggest that the lncRNAs can regulate the expression of mRNA during flowering period in tomato by forming competitive endogenous RNA, and further regulate various biological metabolism pathways in tomato.

Transcriptome characterization of moso bamboo (Phyllostachys edulis) seedlings in response to exogenous gibberellin applications

BACKGROUND

Moso bamboo (Phyllostachys edulis) is a well-known bamboo species of high economic value in the textile industry due to its rapid growth. Phytohormones, which are master regulators of growth and development, serve as important endogenous signals. However, the mechanisms through which phytohormones regulate growth in moso bamboo remain unknown to date.

RESULTS

Here, we reported that exogenous gibberellins (GA) applications resulted in a significantly increased internode length and lignin condensation. Transcriptome sequencing revealed that photosynthesis-related genes were enriched in the GA-repressed gene class, which was consistent with the decrease in leaf chlorophyll concentrations and the lower rate of photosynthesis following GA treatment. Exogenous GA applications on seedlings are relatively easy to perform, thus we used 4-week-old whole seedlings of bamboo for GA- treatment followed by high throughput sequencing. In this study, we identified 932 cis-nature antisense transcripts (cis-NATs), and 22,196 alternative splicing (AS) events in total. Among them, 42 cis-nature antisense transcripts (cis-NATs) and 442 AS events were differentially expressed upon exposure to exogenous GA3, suggesting that post-transcriptional regulation might be also involved in the GA3 response. Targets of differential expression of cis-NATs included genes involved in hormone receptor, photosynthesis and cell wall biogenesis. For example, LAC4 and its corresponding cis-NATs were GA3-induced, and may be involved in the accumulation of lignin, thus affecting cell wall composition.

CONCLUSIONS

This study provides novel insights illustrating how GA alters post-transcriptional regulation and will shed light on the underlying mechanism of growth modulated by GA in moso bamboo.

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.

Cysteine peptidases of Eudiplozoon nipponicum: a broad repertoire of structurally assorted cathepsins L in contrast to the scarcity of cathepsins B in an invasive species of haematophagous monogenean of common carp

Background

Cysteine peptidases of clan CA, family C1 account for a major part of proteolytic activity in the haematophagous monogenean Eudiplozoon nipponicum. The full spectrum of cysteine cathepsins is, however, unknown and their particular biochemical properties, tissue localisation, and involvement in parasite-host relationships are yet to be explored.

Methods

Sequences of cathepsins L and B (EnCL and EnCB) were mined from E. nipponicum transcriptome and analysed bioinformatically. Genes encoding two EnCLs and one EnCB were cloned and recombinant proteins produced in vitro. The enzymes were purified by chromatography and their activity towards selected substrates was characterised. Antibodies and specific RNA probes were employed for localisation of the enzymes/transcripts in tissues of E. nipponicum adults.

Results

Transcriptomic analysis revealed a set of ten distinct transcripts that encode EnCLs. The enzymes are significantly variable in their active sites, specifically the S2 subsites responsible for interaction with substrates. Some of them display unusual structural features that resemble cathepsins B and S. Two recombinant EnCLs had different pH activity profiles against both synthetic and macromolecular substrates, and were able to hydrolyse blood proteins and collagen I. They were localised in the haematin cells of the worm’s digestive tract and in gut lumen. The EnCB showed similarity with cathepsin B2 of Schistosoma mansoni. It displays molecular features typical of cathepsins B, including an occluding loop responsible for its exopeptidase activity. Although the EnCB hydrolysed haemoglobin in vitro, it was localised in the vitelline cells of the parasite and not the digestive tract.

Conclusions

To our knowledge, this study represents the first complex bioinformatic and biochemical characterisation of cysteine peptidases in a monogenean. Eudiplozoon nipponicum adults express a variety of CLs, which are the most abundant peptidases in the worms. The properties and localisation of the two heterologously expressed EnCLs indicate a central role in the (partially extracellular?) digestion of host blood proteins. High variability of substrate-binding sites in the set of EnCLs suggests specific adaptation to a range of biological processes that require proteolysis. Surprisingly, a single cathepsin B is expressed by the parasite and it is not involved in digestion, but probably in vitellogenesis.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2666-2) contains supplementary material, which is available to authorized users.

Bioinformatic analysis of regulation of natural antisense transcripts by transposable elements in human mRNA

Non-coding RNA is no longer considered to be "junk" DNA, based on evidence uncovered in recent decades. In particular, the important role played by natural antisense transcripts (NATs) in regulating the expression of genes is receiving increasing attention. However, the regulatory mechanisms of NATs remain incompletely understood. It is well-known that the insertion of transposable elements (TEs) can affect gene transcription. Using a bioinformatics approach, we identified NATs using human mRNA sequences from the UCSC Genome Browser Database. Our in silico analysis identified 1079 NATs and 700 sense-antisense gene pairs. We identified 179 NATs that showed evidence of having been affected by TEs during cellular gene expression. These findings may provide an understanding of the complex regulation mechanisms of NATs. If our understanding of NATs as modulators of gene expression is further enhanced, we can develop ways to control gene expression.

Identification of Gossypium hirsutum long non-coding RNAs (lncRNAs) under salt stress

BACKGROUND

Long non-coding RNAs (lncRNAs) represent a class of riboregulators that either directly act in long form or are processed into shorter microRNAs (miRNAs) and small interfering RNAs. Long noncoding RNAs (lncRNAs) are arbitrarily defined as RNA genes larger than 200 nt in length that have no apparent coding potential. lncRNAs have emerged as playing important roles in various biological regulatory processes and are expressed in a more tissue-specific manner than mRNA. Emerging evidence shows that lncRNAs participate in stress-responsive regulation.

RESULTS

In this study, in order to develop a comprehensive catalogue of lncRNAs in upland cotton under salt stress, we performed whole-transcriptome strand-specific RNA sequencing for three-leaf stage cotton seedlings treated with salt stress (S_NaCl) and controls (S_CK). In total we identified 1117 unique lncRNAs in this study and 44 differentially expressed RNAs were identified as potential non-coding RNAs. For the differentially expressed lncRNAs that were identified as intergenic lncRNAs (lincRNA), we analysed the gene ontology enrichment of cis targets and found that cis target protein-coding genes were mainly enriched in stress-related categories. Real-time quantitative PCR confirmed that all selected lincRNAs responsive to salt stress. We found lnc_388 was likely as regulator of Gh_A09G1182. And lnc_883 may participate in regulating tolerance to salt stress by modulating the expression of Gh_D03G0339 MS_channel. We then predicted the target mimics for miRNA in Gossypium. six miRNAs were identified, and the result of RT-qPCR with lncRNA and miRNA suggested that lnc_973 and lnc_253 may regulate the expression of ghr-miR399 and ghr-156e as a target mimic under salt stress.

CONCLUSIONS

We identified 44 lincRNAs that were differentially expressed under salt stress. These lincRNAs may target protein-coding genes via cis-acting regulation. We also discovered that specifically-expressed lincRNAs under salt stress may act as endogenous target mimics for conserved miRNAs. These findings extend the current view on lincRNAs as ubiquitous regulators under stress stress.

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.

Complex Regulation of Plant Phosphate Transporters and the Gap between Molecular Mechanisms and Practical Application: What Is Missing?

It has been almost 25 years since the first report of the gene encoding a high-affinity phosphate transporter (PT), PHO84, in yeast. Since then, an increasing number of yeast PHO84 homologs as well as other genes encoding proteins with phosphate (Pi) transport activities have been identified and functionally characterized in diverse plant species. Great progress has been made also in deciphering the molecular mechanism underlying the regulation of the abundance and/or activity of these genes and their products. The regulatory genes affect plant Pi homeostasis commonly through direct or indirect regulation of the abundance of PTs at different levels. However, little has been achieved in the use of PTs for developing genetically modified crops with high phosphorus use efficiency (PUE). This might be a consequence of overemphasizing Pi uptake from the rhizosphere and lack of knowledge about the roles of PTs in Pi transport and recycling within the plant that are required to optimize PUE. Here, we mainly focused on the genes encoding proteins with Pi transport activities and the emerging understanding of their regulation at the transcriptional, post-transcriptional, translational, and post-translational levels. In addition, we propose potential strategies for effective use of PTs in improving plant growth and development.

Functional and expression analyses of transcripts based on full-length cDNAs of Sorghum bicolor

Sorghum bicolor is one of the most important crops for food and bioethanol production. Its small diploid genome and resistance to environmental stress make sorghum an attractive model for studying the functional genomics of the Saccharinae and other C₄ grasses. We analyzed the domain-based functional annotation of the cDNAs using the gene ontology (GO) categories for molecular function to characterize all the genes cloned in the full-length cDNA library of sorghum. The sorghum cDNA library successfully captured a wide range of cDNA-encoded proteins with various functions. To characterize the protein function of newly identified cDNAs, a search of their deduced domains and comparative analyses in the Oryza sativa and Zea mays genomes were carried out. Furthermore, genes on the sense strand corresponding to antisense transcripts were classified based on the GO of molecular function. To add more information about these genes, we have analyzed the expression profiles using RNA-Seq of three tissues (spikelet, seed and stem) during the starch-filling phase. We performed functional analysis of tissue-specific genes and expression analysis of genes of starch biosynthesis enzymes. This functional analysis of sorghum full-length cDNAs and the transcriptome information will facilitate further analysis of the Saccharinae and grass families.

Functions of plants long non-coding RNAs

Long non-coding RNAs (lncRNAs) have been emerged as important players for various biological pathways, including dosage compensation, genomic imprinting, chromatin regulation, alternative splicing and nuclear organization. A large number of lncRNAs had already been identified by different approaches in plants, while the functions of only a few of them have been investigated. This review will summarize our current understanding of a wide range of plant lncRNAs functions, and highlight their roles in the regulation of diverse pathways in plants. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Background

Because of its size, allohexaploid nature, and high repeat content, the bread wheat genome is a good model to study the impact of the genome structure on gene organization, function, and regulation. However, because of the lack of a reference genome sequence, such studies have long been hampered and our knowledge of the wheat gene space is still limited. The access to the reference sequence of the wheat chromosome 3B provided us with an opportunity to study the wheat transcriptome and its relationships to genome and gene structure at a level that has never been reached before.

Results

By combining this sequence with RNA-seq data, we construct a fine transcriptome map of the chromosome 3B. More than 8,800 transcription sites are identified, that are distributed throughout the entire chromosome. Expression level, expression breadth, alternative splicing as well as several structural features of genes, including transcript length, number of exons, and cumulative intron length are investigated. Our analysis reveals a non-monotonic relationship between gene expression and structure and leads to the hypothesis that gene structure is determined by its function, whereas gene expression is subject to energetic cost. Moreover, we observe a recombination-based partitioning at the gene structure and function level.

Conclusions

Our analysis provides new insights into the relationships between gene and genome structure and function. It reveals mechanisms conserved with other plant species as well as superimposed evolutionary forces that shaped the wheat gene space, likely participating in wheat adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0601-9) contains supplementary material, which is available to authorized users.

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.

The expression analysis of silk gland-enriched intermediate-size non-coding RNAs in silkworm Bombyx mori

Small non-protein coding RNAs (ncRNAs) play important roles in development, stress response and other cellular processes. Silkworm is an important model for studies on insect genetics and control of Lepidopterous pests. We have previously identified 189 novel intermediate-size ncRNAs in silkworm Bombyx mori, including 40 ncRNAs that showed altered expression in different developmental stages. Here we characterized the functions of these 40 ncRNAs by measuring their expressions in six tissues of the fifth instar larvae using Northern blot and real-time polymerase chain reaction assays. We identified nine ncRNAs (four small nucleolar RNAs and five unclassified ncRNAs) that were enriched in silk gland, including four ncRNAs that showed silk gland-specific expression. We further showed that three of nine silk gland-enriched ncRNAs were predominantly expressed in the anterior silk gland, whereas another three ncRNAs were highly accumulated in the posterior silk gland, suggesting that they may play different roles in fibroin synthesis. Furthermore, an unclassified ncRNA, Bm-152, exhibited converse expression pattern with its antisense host gene gartenzwerg in diverse tissues, and might regulate the expression of gartenzwerg through RNA-protein complex. In addition, two silk gland-enriched ncRNAs Bm-102 and Bm-159 can be found in histone modification complex, which indicated that they might play roles through epigenetic modifications. Taken together, we provided the first expression and preliminary functional analysis of silk gland-enriched ncRNAs, which will help understand the molecular mechanism of silk gland-development and fibroin synthesis.

Characterization of natural antisense transcript, sclerotia development and secondary metabolism by strand-specific RNA sequencing of Aspergillus flavus

Aspergillus flavus has received much attention owing to its severe impact on agriculture and fermented products induced by aflatoxin. Sclerotia morphogenesis is an important process related to A. flavus reproduction and aflatoxin biosynthesis. In order to obtain an extensive transcriptome profile of A. flavus and provide a comprehensive understanding of these physiological processes, the isolated mRNA of A. flavus CA43 cultures was subjected to high-throughput strand-specific RNA sequencing (ssRNA-seq). Our ssRNA-seq data profiled widespread transcription across the A. flavus genome, quantified vast transcripts (73% of total genes) and annotated precise transcript structures, including untranslated regions, upstream open reading frames (ORFs), alternative splicing variants and novel transcripts. We propose natural antisense transcripts in A. flavus might regulate gene expression mainly on the post-transcriptional level. This regulation might be relevant to tune biological processes such as aflatoxin biosynthesis and sclerotia development. Gene Ontology annotation of differentially expressed genes between the mycelia and sclerotia cultures indicated sclerotia development was related closely to A. flavus reproduction. Additionally, we have established the transcriptional profile of aflatoxin biosynthesis and its regulation model. We identified potential genes linking sclerotia development and aflatoxin biosynthesis. These genes could be used as targets for controlled regulation of aflatoxigenic strains of A. flavus.

Rice germline-specific Argonaute MEL1 protein binds to phasiRNAs generated from more than 700 lincRNAs

Small RNAs that interact with Argonaute (AGO) proteins play central roles in RNA-mediated silencing. MEIOSIS ARRESTED AT LEPTOTENE1 (MEL1), a rice AGO, has specific functions in the development of pre-meiotic germ cells and the progression of meiosis. Here, we show that MEL1, which is located mostly in the cytoplasm of germ cells, associates preferentially with 21-nucleotide phased small interfering RNAs (phasiRNAs) that bear a 5'-terminal cytosine. Most phasiRNAs are derived from 1171 intergenic clusters distributed on all rice chromosomes. From these clusters, over 700 large intergenic, non-coding RNAs (lincRNAs) that contain the consensus sequence complementary to miR2118 are transcribed specifically in inflorescences, and cleaved within the miR2118 site. Cleaved lincRNAs are processed via DICER-LIKE4 (DCL4) protein, resulting in production of phasiRNAs. This study provides the evidence that the miR2118-dependent and the DCL4-dependent pathways are both required for biogenesis of 21-nt phasiRNAs associated with germline-specific MEL1 AGO in rice, and over 700 lincRNAs are key factors for induction of this biogenesis during reproductive-specific stages.

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.

Global analysis of cis-natural antisense transcripts and their heat-responsive nat-siRNAs in Brassica rapa

BACKGROUND

Brassica rapa includes several important leaf vegetable crops whose production is often damaged by high temperature. Cis-natural antisense transcripts (cis-NATs) and cis-NATs-derived small interfering RNAs (nat-siRNAs) play important roles in plant development and stress responses. However, genome-wide cis-NATs in B. rapa are not known. The NATs and nat-siRNAs that respond to heat stress have never been well studied in B. rapa. Here, we took advantage of RNA-seq and small RNA (sRNA) deep sequencing technology to identify cis-NATs and heat responsive nat-siRNAs in B. rapa.

RESULTS

Analyses of four RNA sequencing datasets revealed 1031 cis-NATs B. rapa ssp. chinensis cv Wut and B. rapa ssp. pekinensis cv. Bre. Based on sequence homology between Arabidopsis thaliana and B. rapa, 303 conserved cis-NATs in B. rapa were found to correspond to 280 cis-NATs in Arabidopsis; the remaining 728 novel cis-NATs were identified as Brassica-specific ones. Using six sRNA libraries, 4846 nat-siRNAs derived from 150 cis-NATs were detected. Differential expression analysis revealed that nat-siRNAs derived from 12 cis-NATs were responsive to heat stress, and most of them showed strand bias. Real-time PCR indicated that most of the transcripts generating heat-responsive nat-siRNAs were upregulated under heat stress, while the transcripts from the opposite strands of the same loci were downregulated.

CONCLUSIONS

Our results provide the first subsets of genome-wide cis-NATs and heat-responsive nat-siRNAs in B. rapa; these sRNAs are potentially useful for the genetic improvement of heat tolerance in B. rapa and other crops.

A rice cis-natural antisense RNA acts as a translational enhancer for its cognate mRNA and contributes to phosphate homeostasis and plant fitness

cis-natural antisense transcripts (cis-NATs) are widespread in plants and are often associated with downregulation of their associated sense genes. We found that a cis-NAT positively regulates the level of a protein critical for phosphate homeostasis in rice (Oryza sativa). PHOSPHATE1;2 (PHO1;2), a gene involved in phosphate loading into the xylem in rice, and its associated cis-NATPHO1;2 are both controlled by promoters active in the vascular cylinder of roots and leaves. While the PHO1;2 promoter is unresponsive to the plant phosphate status, the cis-NATPHO1;2 promoter is strongly upregulated under phosphate deficiency. Expression of both cis-NATPHO1;2 and the PHO1;2 protein increased in phosphate-deficient plants, while the PHO1;2 mRNA level remained stable. Downregulation of cis-NATPHO1;2 expression by RNA interference resulted in a decrease in PHO1;2 protein, impaired the transfer of phosphate from root to shoot, and decreased seed yield. Constitutive overexpression of NATPHO1;2 in trans led to a strong increase of PHO1;2, even under phosphate-sufficient conditions. Under all conditions, no changes occurred in the level of expression, sequence, or nuclear export of PHO1;2 mRNA. However, expression of cis-NATPHO1;2 was associated with a shift of both PHO1;2 and cis-NATPHO1;2 toward the polysomes. These findings reveal an unexpected role for cis-NATPHO1;2 in promoting PHO1;2 translation and affecting phosphate homeostasis and plant fitness.

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.

Safety assessment of food and feed from biotechnology-derived crops employing RNA-mediated gene regulation to achieve desired traits: a scientific review

Gene expression can be modulated in plants to produce desired traits through agricultural biotechnology. Currently, biotechnology-derived crops are compared to their conventional counterparts, with safety assessments conducted on the genetic modification and the intended and unintended differences. This review proposes that this comparative safety assessment paradigm is appropriate for plants modified to express mediators of RNA-mediated gene regulation, including RNA interference (RNAi), a gene suppression mechanism that naturally occurs in plants and animals. The molecular mediators of RNAi, including long double-stranded RNAs (dsRNA), small interfering RNAs (siRNA), and microRNAs (miRNA), occur naturally in foods; therefore, there is an extensive history of safe consumption. Systemic exposure following consumption of plants containing dsRNAs that mediate RNAi is limited in higher organisms by extensive degradation of ingested nucleic acids and by biological barriers to uptake and efficacy of exogenous nucleic acids. A number of mammalian RNAi studies support the concept that a large margin of safety will exist for any small fraction of RNAs that might be absorbed following consumption of foods from biotechnology-derived plants that employ RNA-mediated gene regulation. Food and feed derived from these crops utilizing RNA-based mechanisms is therefore expected to be as safe as food and feed derived through conventional plant breeding.

Strand-specific RNA-seq reveals widespread occurrence of novel cis-natural antisense transcripts in rice

Background

Cis-natural antisense transcripts (cis-NATs) are RNAs transcribed from the antisense strand of a gene locus, and are complementary to the RNA transcribed from the sense strand. Common techniques including microarray approach and analysis of transcriptome databases are the major ways to globally identify cis-NATs in various eukaryotic organisms. Genome-wide in silico analysis has identified a large number of cis-NATs that may generate endogenous short interfering RNAs (nat-siRNAs), which participate in important biogenesis mechanisms for transcriptional and post-transcriptional regulation in rice. However, the transcriptomes are yet to be deeply sequenced to comprehensively investigate cis-NATs.

Results

We applied high-throughput strand-specific complementary DNA sequencing technology (ssRNA-seq) to deeply sequence mRNA for assessing sense and antisense transcripts that were derived under salt, drought and cold stresses, and normal conditions, in the model plant rice (Oryza sativa). Combined with RAP-DB genome annotation (the Rice Annotation Project Database build-5 data set), 76,013 transcripts corresponding to 45,844 unique gene loci were assembled, in which 4873 gene loci were newly identified. Of 3819 putative rice cis-NATs, 2292 were detected as expressed and giving rise to small RNAs from their overlapping regions through integrated analysis of ssRNA-seq data and small RNA data. Among them, 503 cis-NATs seemed to be associated with specific conditions. The deep sequence data from isolated epidermal cells of rice seedlings further showed that 54.0% of cis-NATs were expressed simultaneously in a population of homogenous cells. Nearly 9.7% of rice transcripts were involved in one-to-one or many-to-many cis-NATs formation. Furthermore, only 17.4-34.7% of 223 many-to-many cis-NAT groups were all expressed and generated nat-siRNAs, indicating that only some cis-NAT groups may be involved in complex regulatory networks.

Conclusions

Our study profiles an abundance of cis-NATs and nat-siRNAs in rice. These data are valuable for gaining insight into the complex function of the rice transcriptome.

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

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

Identification of potential antisense transcripts in rice using conventional microarray

Natural antisense transcripts (NATs) are endogenous transcripts that contain reverse complementary sequences to other RNAs (usually called sense transcripts). NATs regulate the expression of sense transcripts in a wide range of species. The identification and analysis of NATs are the prerequisite to elucidate their functions. Microarray is a genome-wide method to detect gene expression. However, conventional microarrays do not contain the specific probes of NATs; thus, they cannot be utilized to detect NATs. In this article, we developed a novel method to identify potential NATs with the conventional microarrays. In this method of our study, we labeled the first strand cDNA from one sample with Cy5 and labeled the second strand cDNA from another sample with Cy3, and then hybridized these labeled samples with oligonucleotide microarray. Using this method, we identified 920 potential NATs in rice variety Nipponbare. Among these potential NATs, 88 of them were confirmed by either full-length cDNA or orientated ESTs (expressed sequence tags). This is the first time that a conventional oligonucleotide microarray was employed to identify NATs in rice.

Antisense RNAs and epigenetic regulation

Antisense RNA is the first noncoding RNA found to have a regulatory function. With the advances of biological science, it has been recognized that the function of antisense RNAs is not only limited to post-transcriptional regulation, but extends to transcriptional regulation of various important genes leading to epigenetic changes in DNA methylation and histone modifications. Gene regulation by antisense RNA is a general phenomenon observed in eukaryotic cells while genome-wide natural antisense transcripts have been identified in many animals and plants. Antisense RNAs are not only involved in X-chromosome inactivation and imprinted silencing in normal cells, but aberrantly expressed antisense RNAs can also induce epigenetic silencing of tumor suppressor genes in cancer cells and deletion-induced aberrant antisense RNAs lead to epigenetic silencing and diseases. While a general picture of the pathways involved in antisense RNA-mediated gene regulation has emerged, many questions remain. The mechanisms by which genes are regulated by antisense RNAs, antisense transcript itself is produced and aberrant antisense RNAs induce human diseases are all research focuses of the future.

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

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

C6orf176: a novel possible regulator of cAMP-mediated gene expression

cAMP mediates diverse cellular signals including prostaglandin (PG) E(2)-mediated intraocular pressure (IOP)-lowering activity in human ocular ciliary smooth muscle cells (hCSM). We have identified gene regulatory networks and key genes upon activation of the cAMP pathway in hCSM, using novel agonists highly selective for PGE(2) receptor subtypes EP2 or EP4, which are G protein-coupled receptors well known to activate cAMP signaling. Here we describe a novel, EP2/EP4-induced, primate-specific gene of hitherto unknown function, also known as C6orf176 (chromosome 6 open reading frame 176) and recently reclassified as noncoding RNA in NCBI's database. Its expression, as determined by quantitative real-time RT-PCR (qRT-PCR), is dramatically upregulated (>2,000-fold) subsequent to transduction of EP2/EP4/Gs/cAMP signaling not only in hCSM, but also in HEK cells overexpressing the recombinant receptors. Moreover, activation of other IOP lowering, Gs-coupled prostanoid receptors, such as DP1 and IP, as well as a direct activator of adenylyl cyclase, forskolin, also substantially upregulated C6orf176 in hCSM, while FP and TP, which are Gq-coupled prostanoid receptor subtypes, did not. Novel transcript variants carrying open reading frames, derived from an at least 67 kb genomic locus on chromosome 6q27 with putative alternative transcription start sites, were identified. Transcriptional upregulation of transcript variants as well as of two genes expressed in antisense orientation that partially overlap the transcribed C6orf176 region was observed, to varying degrees, subsequent to induction of cAMP signaling using various agonists. Small interfering RNA-mediated C6orf176 gene silencing experiments showed modulation of several cAMP-responsive genes. These transcriptional activities identify C6orf176 as a potential biomarker and/or therapeutic target in context with diseases linked to deregulated cAMP signaling. Also, the cAMP-inducible C6orf176 gene locus could be useful as a model system for studying transcriptional regulation by chromatin and RNA polymerase II.

PlantNATsDB: a comprehensive database of plant natural antisense transcripts

Natural antisense transcripts (NATs), as one type of regulatory RNAs, occur prevalently in plant genomes and play significant roles in physiological and pathological processes. Although their important biological functions have been reported widely, a comprehensive database is lacking up to now. Consequently, we constructed a plant NAT database (PlantNATsDB) involving approximately 2 million NAT pairs in 69 plant species. GO annotation and high-throughput small RNA sequencing data currently available were integrated to investigate the biological function of NATs. PlantNATsDB provides various user-friendly web interfaces to facilitate the presentation of NATs and an integrated, graphical network browser to display the complex networks formed by different NATs. Moreover, a ‘Gene Set Analysis’ module based on GO annotation was designed to dig out the statistical significantly overrepresented GO categories from the specific NAT network. PlantNATsDB is currently the most comprehensive resource of NATs in the plant kingdom, which can serve as a reference database to investigate the regulatory function of NATs. The PlantNATsDB is freely available at http://bis.zju.edu.cn/pnatdb/.

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.

Comprehensive sequence analysis of 24,783 barley full-length cDNAs derived from 12 clone libraries

Full-length cDNA (FLcDNA) libraries consisting of 172,000 clones were constructed from a two-row malting barley cultivar (Hordeum vulgare 'Haruna Nijo') under normal and stressed conditions. After sequencing the clones from both ends and clustering the sequences, a total of 24,783 complete sequences were produced. By removing duplicates between these and publicly available sequences, 22,651 representative sequences were obtained: 17,773 were novel barley FLcDNAs, and 1,699 were barley specific. Highly conserved genes were found in the barley FLcDNA sequences for 721 of 881 rice (Oryza sativa) trait genes with 50% or greater identity. These FLcDNA resources from our Haruna Nijo cDNA libraries and the full-length sequences of representative clones will improve our understanding of the biological functions of genes in barley, which is the cereal crop with the fourth highest production in the world, and will provide a powerful tool for annotating the barley genome sequences that will become available in the near future.

Regulated antisense transcription controls expression of cell-type-specific genes in yeast

Transcriptome profiling studies have recently uncovered a large number of noncoding RNA transcripts (ncRNAs) in eukaryotic organisms, and there is growing interest in their role in the cell. For example, in haploid Saccharomyces cerevisiae cells, the expression of an overlapping antisense ncRNA, referred to here as RME2 (Regulator of Meiosis 2), prevents IME4 expression. In diploid cells, the a1-α2 complex represses the transcription of RME2, allowing IME4 to be induced during meiosis. In this study we show that antisense transcription across the IME4 promoter region does not block transcription factors from binding and is not required for repression. Mutational analyses found that sequences within the IME4 open reading frame (ORF) are required for the repression mediated by RME2 transcription. These results support a model where transcription of RME2 blocks the elongation of the full-length IME4 transcript but not its initiation. We have found that another antisense transcript, called RME3, represses ZIP2 in a cell-type-specific manner. These results suggest that regulated antisense transcription may be a widespread mechanism for the control of gene expression and may account for the roles of some of the previously uncharacterized ncRNAs in yeast.

Experimental RNomics and genomic comparative analysis reveal a large group of species-specific small non-message RNAs in the silkworm Bombyx mori

Accumulating evidences show that small non-protein coding RNAs (ncRNAs) play important roles in development, stress response and other cellular processes. The silkworm is an important model for studies on insect genetics and control of lepidopterous pests. Here, we have performed the first systematic identification and analysis of intermediate size ncRNAs (50–500 nt) in the silkworm. We identified 189 novel ncRNAs, including 141 snoRNAs, six snRNAs, three tRNAs, one SRP and 38 unclassified ncRNAs. Forty ncRNAs showed significantly altered expression during silkworm development or across specific stage transitions. Genomic comparisons revealed that 123 of these ncRNAs are potentially silkworm-specific. Analysis of the genomic organization of the ncRNA loci showed that 32.62% of the novel snoRNA loci are intergenic, and that all the intronic snoRNAs follow the pattern of one-snoRNA-per-intron. Target site analysis predicted a total of 95 2′-O-methylation and pseudouridylation modification sites of rRNAs, snRNAs and tRNAs. Together, these findings provide new clues for future functional study of ncRNA during insect development and evolution.

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 prevalence and regulation of antisense transcripts in Schizosaccharomyces pombe

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

Detection of antisense RNA transcripts by strand-specific RT-PCR

Comprehensive genome annotation requires extensive cDNA analysis. This analysis has identified natural antisense transcripts (NATs), which are distinct from the microRNAs, siRNAs, and piRNAs, in a number of diverse eukaryotes. This wide conservation supports the possibility of an important role for NATs in regulating cellular processes. Investigating their roles requires the confirmation of expressed sequence tag (EST) data and the detection of antisense transcripts in distinct cellular backgrounds. This chapter describes the use of a reverse transcription polymerase chain reaction (RT-PCR) method for the detection of antisense transcripts. The protocol was designed to reduce the number of first strand synthesis reactions during screening for antisense transcripts through the utilization of antisense directed primers and oligo dT to prime first strand synthesis. These results are further confirmed using sense and antisense directed primers in first strand synthesis. Results indicate that optimization of the screens requires proper controls to confirm removal of gDNA contamination and to rule out self-priming as a source of first strand products.

Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons

Phosphate homeostasis was studied in a monocotyledonous model plant through the characterization of the PHO1 gene family in rice (Oryza sativa). Bioinformatics and phylogenetic analysis showed that the rice genome has three PHO1 homologs, which cluster with the Arabidopsis (Arabidopsis thaliana) AtPHO1 and AtPHO1;H1, the only two genes known to be involved in root-to-shoot transfer of phosphate. In contrast to the Arabidopsis PHO1 gene family, all three rice PHO1 genes have a cis-natural antisense transcript located at the 5 ' end of the genes. Strand-specific quantitative reverse transcription-PCR analyses revealed distinct patterns of expression for sense and antisense transcripts for all three genes, both at the level of tissue expression and in response to nutrient stress. The most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. However, while the OsPHO1;2 sense transcript was relatively stable under various nutrient deficiencies, the antisense transcript was highly induced by inorganic phosphate (Pi) deficiency. Characterization of Ospho1;1 and Ospho1;2 insertion mutants revealed that only Ospho1;2 mutants had defects in Pi homeostasis, namely strong reduction in Pi transfer from root to shoot, which was accompanied by low-shoot and high-root Pi. Our data identify OsPHO1;2 as playing a key role in the transfer of Pi from roots to shoots in rice, and indicate that this gene could be regulated by its cis-natural antisense transcripts. Furthermore, phylogenetic analysis of PHO1 homologs in monocotyledons and dicotyledons revealed the emergence of a distinct clade of PHO1 genes in dicotyledons, which include members having roles other than long-distance Pi transport.

Rice MicroRNA effector complexes and targets

MicroRNAs (miRNAs) are small silencing RNAs with regulatory roles in gene expression. miRNAs interact with Argonaute (AGO) proteins to form effector complexes that cleave target mRNAs or repress translation. Rice (Oryza sativa) encodes four AGO1 homologs (AGO1a, AGO1b, AGO1c, and AGO1d). We used RNA interference (RNAi) to knock down the four AGO1s. The RNAi lines displayed pleiotropic developmental phenotypes and had increased accumulation of miRNA targets. AGO1a, AGO1b, and AGO1c complexes were purified and further characterized. The three AGO1s all have a strong preference for binding small RNAs (sRNAs) with 5' U and have Slicer activity. We cataloged the sRNAs in each AGO1 complex by deep sequencing and found that all three AGO1s predominantly bound known miRNAs. Most of the miRNAs were evenly distributed in the three AGO1 complexes, suggesting a redundant role for the AGO1s. Intriguingly, a subset of miRNAs were specifically incorporated into or excluded from one of the AGO1s, suggesting functional specialization among the AGO1s. Furthermore, we identified rice miRNA targets at a global level. The validated targets include transcription factors that control major stages of development and also genes involved in a variety of physiological processes, indicating a broad regulatory role for miRNAs in rice.

Large-scale analysis of antisense transcription in wheat using the Affymetrix GeneChip Wheat Genome Array

Background

Natural antisense transcripts (NATs) are transcripts of the opposite DNA strand to the sense-strand either at the same locus (cis-encoded) or a different locus (trans-encoded). They can affect gene expression at multiple stages including transcription, RNA processing and transport, and translation. NATs give rise to sense-antisense transcript pairs and the number of these identified has escalated greatly with the availability of DNA sequencing resources and public databases. Traditionally, NATs were identified by the alignment of full-length cDNAs or expressed sequence tags to genome sequences, but an alternative method for large-scale detection of sense-antisense transcript pairs involves the use of microarrays. In this study we developed a novel protocol to assay sense- and antisense-strand transcription on the 55 K Affymetrix GeneChip Wheat Genome Array, which is a 3' in vitro transcription (3'IVT) expression array. We selected five different tissue types for assay to enable maximum discovery, and used the 'Chinese Spring' wheat genotype because most of the wheat GeneChip probe sequences were based on its genomic sequence. This study is the first report of using a 3'IVT expression array to discover the expression of natural sense-antisense transcript pairs, and may be considered as proof-of-concept.

Results

By using alternative target preparation schemes, both the sense- and antisense-strand derived transcripts were labeled and hybridized to the Wheat GeneChip. Quality assurance verified that successful hybridization did occur in the antisense-strand assay. A stringent threshold for positive hybridization was applied, which resulted in the identification of 110 sense-antisense transcript pairs, as well as 80 potentially antisense-specific transcripts. Strand-specific RT-PCR validated the microarray observations, and showed that antisense transcription is likely to be tissue specific. For the annotated sense-antisense transcript pairs, analysis of the gene ontology terms showed a significant over-representation of transcripts involved in energy production. These included several representations of ATP synthase, photosystem proteins and RUBISCO, which indicated that photosynthesis is likely to be regulated by antisense transcripts.

Conclusion

This study demonstrated the novel use of an adapted labeling protocol and a 3'IVT GeneChip array for large-scale identification of antisense transcription in wheat. The results show that antisense transcription is relatively abundant in wheat, and may affect the expression of valuable agronomic phenotypes. Future work should select potentially interesting transcript pairs for further functional characterization to determine biological activity.

Characterization and expression profiles of miRNAs in rice seeds

Small RNAs (sRNAs) are common and effective modulators of gene expression in eukaryotic organisms. To characterize the sRNAs expressed during rice seed development, massively parallel signature sequencing (MPSS) was performed, resulting in the obtainment of 797 399 22-nt sequence signatures, of which 111 161 are distinct ones. Analysis on the distributions of sRNAs on chromosomes showed that most sRNAs originate from interspersed repeats that mainly consist of transposable elements, suggesting the major function of sRNAs in rice seeds is transposon silencing. Through integrative analysis, 26 novel miRNAs and 12 miRNA candidates were identified. Further analysis on the expression profiles of the known and novel miRNAs through hybridizing the generated chips revealed that most miRNAs were expressed preferentially in one or two rice tissues. Detailed comparison of the expression patterns of miRNAs and corresponding target genes revealed the negative correlation between them, while few of them are positively correlated. In addition, differential accumulations of miRNAs and corresponding miRNA*s suggest the functions of miRNA*s other than being passenger strands of mature miRNAs, and in regulating the miRNA functions.

rasiRNA pathway controls antisense expression of Drosophila telomeric retrotransposons in the nucleus

Telomeres in Drosophila are maintained by the specialized telomeric retrotransposons HeT-A, TART and TAHRE. Sense transcripts of telomeric retroelements were shown to be the targets of a specialized RNA-interference mechanism, a repeat-associated short interfering (rasi)RNA-mediated system. Antisense rasiRNAs play a key role in this mechanism, highlighting the importance of antisense expression in retrotransposon silencing. Previously, bidirectional transcription was reported for the telomeric element TART. Here, we show that HeT-A is also bidirectionally transcribed, and HeT-A antisense transcription in ovaries is regulated by a promoter localized within its 3' untranslated region. A remarkable feature of noncoding HeT-A antisense transcripts is the presence of multiple introns. We demonstrate that sense and antisense HeT-A-specific rasiRNAs are present in the same tissue, indicating that transcripts of both directions may be considered as natural targets of the rasiRNA pathway. We found that the expression of antisense transcripts of telomeric elements is regulated by the RNA silencing machinery, suggesting rasiRNA-mediated interplay between sense and antisense transcripts in the cell. Finally, this regulation occurs in the nucleus since disruption of the rasiRNA pathway leads to an accumulation of TART and HeT-A transcripts in germ cell nuclei.