A transcriptome‐wide study on the microRNA‐ and the Argonaute 1‐enriched small RNA‐mediated regulatory networks involved in plant leaf senescence

Leaf senescence is an important physiological process during the plant life cycle. However, systemic studies on the impact of microRNAs (miRNAs) on the expression of senescence‐associated genes (SAGs) are lacking. Besides, whether other Argonaute 1 (AGO1)‐enriched small RNAs (sRNAs) play regulatory roles in leaf senescence remains unclear. In this study, a total of 5,123 and 1,399 AGO1‐enriched sRNAs, excluding miRNAs, were identified in Arabidopsis thaliana and rice (Oryza sativa), respectively. After retrieving SAGs from the Leaf Senescence Database, all of the AGO1‐enriched sRNAs and the miRBase‐registered miRNAs of these two plants were included for target identification. Supported by degradome signatures, 200 regulatory pairs involving 120 AGO1‐enriched sRNAs and 40 SAGs, and 266 regulatory pairs involving 64 miRNAs and 42 SAGs were discovered in Arabidopsis. Moreover, 13 genes predicted to interact with some of the above‐identified target genes at protein level were validated as regulated by 17 AGO1‐enriched sRNAs and ten miRNAs in Arabidopsis. In rice, only one SAG was targeted by three AGO1‐enriched sRNAs, and one SAG was targeted by miR395. However, five AGO1‐enriched sRNAs were conserved between Arabidopsis and rice. Target genes conserved between the two plants were identified for three of the above five sRNAs, pointing to the conserved roles of these regulatory pairs in leaf senescence or other developmental procedures. Novel targets were discovered for three of the five AGO1‐enriched sRNAs in rice, indicating species‐specific functions of these sRNA–target pairs. These results could advance our understanding of the sRNA‐involved molecular processes modulating leaf senescence.     

Identification of Novel Oryza sativa miRNAs in Deep Sequencing-Based Small RNA Libraries of Rice Infected with Rice Stripe Virus

MicroRNAs (miRNAs) play essential regulatory roles in the development of eukaryotes. Methods based on deep-sequencing have provided a powerful high-throughput strategy for identifying novel miRNAs and have previously been used to identify over 100 novel miRNAs from rice. Most of these reports are related to studies of rice development, tissue differentiation, or abiotic stress, but novel rice miRNAs related to viral infection have rarely been identified. In previous work, we constructed and pyrosequenced the small RNA (sRNA) libraries of rice infected with Rice stripe virus and described the character of the small interfering RNAs (siRNA) derived from the RSV RNA genome. We now report the identification of novel miRNAs from the abundant sRNAs (with a minimum of 100 sequencing reads) in the sRNA library of RSV-infected rice. 7 putative novel miRNAs (pn-miRNAs) whose precursor sequences have not previously been described were identified and could be detected by Northern blot or RT-PCR, and were recognized as novel miRNAs (n-miRNAs). Further analysis showed that 5 of the 7 n-miRNAs were up-expressed while the other 2 n-miRNAs were down-expressed in RSV-infected rice. In addition, 23 pn-miRNAs that were newly produced from 19 known miRNA precursors were also identified. This is first report of novel rice miRNAs produced from new precursors related to RSV infection.     

Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns

Plants feature a particularly diverse population of short (s)RNAs, the central component of all RNA silencing pathways. Next generation sequencing techniques enable deeper insights into this complex and highly conserved mechanism and allow identification and quantification of sRNAs. We employed deep sequencing to monitor the sRNAome of developing tomato fruits covering the period between closed flowers and ripened fruits by profiling sRNAs at 10 time-points. It is known that microRNAs (miRNAs) play an important role in development but very little information is available about the majority of sRNAs that are not miRNAs. Here we show distinctive patterns of sRNA expression that often coincide with stages of the developmental process such as flowering, early and late fruit maturation. Moreover, thousands of non-miRNA sRNAs are differentially expressed during fruit development and ripening. Some of these differentially expressed sRNAs derived from transposons but many derive from protein coding genes or regions that show homology to protein coding genes, several of which are known to play a role in flower and fruit development. These findings raise the possibility of a regulative role of these sRNAs during fruit onset and maturation in a crop species. We also identified six new miRNAs and experimentally validated two target mRNAs. These two mRNAs are targeted by the same miRNA but do not belong to the same gene family, which is rare for plant miRNAs. Expression pattern and putative function of these targets indicate a possible role in glutamate accumulation, which contributes to establishing the taste of the fruit.     

Genome-Wide Identification of Reverse Complementary microRNA Genes in Plants

MicroRNAs (miRNAs) are ∼21-nucleotide small RNAs (sRNAs) with essential regulatory roles in plants. They are generated from stem-loop-structured precursors through two sequential Dicer-like 1 (DCL1)-mediated cleavages. To date, hundreds of plant miRNAs have been uncovered. However, the question, whether the sequences reverse complementary (RC) to the miRNA precursors could form hairpin-like structures and produce sRNA duplexes similar to the miRNA/miRNA* pairs has not been solved yet. Here, we interrogated this possibility in 16 plant species based on sRNA high-throughput sequencing data and secondary structure prediction. A total of 59 RC sequences with great potential to form stem-loop structures and generate miRNA/miRNA*-like duplexes were identified in ten plants, which were named as RC-miRNA precursors. Unlike the canonical miRNAs, only a few cleavage targets of the RC-miRNAs were identified in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), and none in Soybean (Glycine max) based on degradome data. Surprisingly, the genomic regions surrounding some of the RC-miRNA target recognition sites were observed to be specifically methylated in both Arabidopsis and rice. Taken together, we reported a new class of miRNAs, called RC-miRNAs, which were generated from the antisense strands of the miRNA precursors. Based on the results, we speculated that the mature RC-miRNAs might have subtle regulatory activity through target cleavages, but might possess short interfering RNA-like activity by guiding sequence-specific DNA methylation.     

Small RNA in rice genome

Rice has many characteristics of a model plant. The recent completion of the draft of the rice genome represents an important advance in our knowledge of plant biology and also has an important contribution to the understanding of general genomic evolution. Besides the rice genome finishing map, the next urgent step for rice researchers is to annotate the genes and noncoding functional sequences. The recent work shows that noncoding RNAs (ncRNAs) play significant roles in biological systems. We have explored all the known small RNAs (a kind of ncRNA) within rice genome and other six species sequences, including Arabidopsis, maize, yeast, worm, mouse and pig. As a result we find 160 out of 552 small RNAs (sRNAs) in database have homologs in 108 rice scaffolds, and almost all of them (99.41%) locate in intron regions of rice by gene predication. 19 sRNAs only appear in rice. More importantly, we find two special LJ14 sRNAs: one is located in a set of sRNA ZMU14SNR9(s) which only appears in three plants, 86% sequences of them can be compared as the same sequence in rice, Arabidopsis and maize; the other conserved sRNA XLHS7CU14 has a segment which appears in almost all these species from plants to animals. All these results indicate that sRNA do not have evident borderline between plants and animals.     

Small RNAs and Gene Network in a Durable Disease Resistance Gene—Mediated Defense Responses in Rice

Accumulating data have suggested that small RNAs (sRNAs) have important functions in plant responses to pathogen invasion. However, it is largely unknown whether and how sRNAs are involved in the regulation of rice responses to the invasion of Xanthomonas oryzae pv. oryzae (Xoo), which causes bacterial blight, the most devastating bacterial disease of rice worldwide. We performed simultaneous genome-wide analyses of the expression of sRNAs and genes during early defense responses of rice to Xoo mediated by a major disease resistance gene, Xa3/Xa26, which confers durable and race-specific qualitative resistance. A large number of sRNAs and genes showed differential expression in Xa3/Xa26-mediated resistance. These differentially expressed sRNAs include known microRNAs (miRNAs), unreported miRNAs, and small interfering RNAs. The candidate genes, with expression that was negatively correlated with the expression of sRNAs, were identified, indicating that these genes may be regulated by sRNAs in disease resistance in rice. These results provide a new perspective regarding the putative roles of sRNA candidates and their putative target genes in durable disease resistance in rice.     

RNA-Seq of the Nucleolus Reveals Abundant SNORD44-Derived Small RNAs

Small non-coding RNAs represent RNA species that are not translated to proteins, but which have diverse and broad functional activities in physiological and pathophysiological states. The knowledge of these small RNAs is rapidly expanding in part through the use of massive parallel (deep) sequencing efforts. We present here the first deep sequencing of small RNomes in subcellular compartments with particular emphasis on small RNAs (sRNA) associated with the nucleolus. The vast majority of the cellular, cytoplasmic and nuclear sRNAs were identified as miRNAs. In contrast, the nucleolar sRNAs had a unique size distribution consisting of 19–20 and 25 nt RNAs, which were predominantly composed of small snoRNA-derived box C/D RNAs (termed as sdRNA). Sequences from 47 sdRNAs were identified, which mapped to both 5′ and 3′ ends of the snoRNAs, and retained conserved box C or D motifs. SdRNA reads mapping to SNORD44 comprised 74% of all nucleolar sdRNAs, and were confirmed by Northern blotting as comprising both 20 and 25 nt RNAs. A novel 120 nt SNORD44 form was also identified. The expression of the SNORD44 sdRNA and 120 nt form was independent of Dicer/Drosha–mediated processing pathways but was dependent on the box C/D snoRNP proteins/sno-ribonucleoproteins fibrillarin and NOP58. The 120 nt SNORD44-derived RNA bound to fibrillarin suggesting that C/D sno-ribonucleoproteins are involved in regulating the stability or processing of SNORD44. This study reveals sRNA cell-compartment specific expression and the distinctive unique composition of the nucleolar sRNAs.     

Comparative genomic analysis of CIPK gene family in <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Populus</Emphasis>

Calcium serves as a second messenger in various signal transduction pathways in plants. CBL-interacting protein kinases (CIPKs), which have a variety of functions, are involved in calcium signal transduction. Previous, the studies on CIPK family members focused on Arabidopsis and rice. Here, we present a comparative genomic analysis of the CIPK gene family in Arabidopsis and poplar, a model tree species. Twenty-seven potential CIPKs were identified from poplar using genome-wide analysis. Like the CIPK gene family from Arabidopsis, CIPK genes from poplar were also divided into intron-free and intron-harboring groups. In the intron-harboring group, the intron distribution of CIPKs is rather conserved during the genome evolutionary process. Many homologous gene pairs were found in the CIPK gene family, indicating duplication events might contribute to the amplification of this gene family. The phylogenetic comparison of CIPKs in combination with intron distribution analysis revealed that CIPK genes from both Arabidopsis and poplar might have an ancient origin, which formed earlier than the separation of these two eudicot species. Our genomic and bioinformatic analysis will provide an important foundation for further functional dissection of the CBL-CIPK signaling network in poplars. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Sex specific expression and distribution of small RNAs in papaya

Background

Regulatory function of small non-coding RNAs (sRNA) in response to environmental and developmental cues has been established. Additionally, sRNA, also plays an important role in maintaining the heterochromatin and centromere structures of the chromosome. Papaya, a trioecious species with recently evolved sex chromosomes, has emerged as an excellent model system to study sex determination and sex chromosome evolution in plants. However, role of small RNA in papaya sex determination is yet to be explored.

Results

We analyzed the high throughput sRNAs reads in the Illumina libraries prepared from male, female, and hermaphrodite flowers of papaya. Using the sRNA reads, we identified 29 miRNAs that were not previously reported from papaya. Including this and two previous studies, a total of 90 miRNAs has been identified in papaya. We analyzed the expression of these miRNAs in each sex types. A total of 65 miRNAs, including 31 conserved and 34 novel mirNA, were detected in at least one library. Fourteen of the 65 miRNAs were differentially expressed among different sex types. Most of the miRNA expressed higher in male flowers were related to the auxin signaling pathways, whereas the miRNAs expressed higher in female flowers were the potential regulators of the apical meristem identity genes. Aligning the sRNA reads identified the sRNA hotspots adjacent to the gaps of the X and Y chromosomes. The X and Y chromosomes sRNA hotspots has a 7.8 and 4.4 folds higher expression of sRNA, respectively, relative to the chromosome wide average. Approximately 75% of the reads aligned to the X chromosome hotspot was identical to that of the Y chromosome hotspot.

Conclusion

By analyzing the large-scale sRNA sequences from three sex types, we identified the sRNA hotspots flanking the gaps of papaya X, Y, and Y^h chromosome. The sRNAs expression patterns in these regions were reminiscent of the pericentromeric region indicating that the only remaining gap in each of these chromosomes is likely the centromere. We also identified 14 differentially expressed miRNAs in male, female and hermaphrodite flowers of papaya. Our results provide valuable information toward understanding the papaya sex determination.

Electronic supplementary material

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

Analysis of plant-derived miRNAs in animal small RNA datasets

ABSTRACT: BACKGROUND: Plants contain significant quantities of small RNAs (sRNAs) derived from various sRNA biogenesis pathways. Many of these sRNAs play regulatory roles in plants. Previous analysis revealed that numerous sRNAs in corn, rice and soybean seeds have high sequence similarity to animal genes. However, exogenous RNA is considered to be unstable within the gastrointestinal tract of many animals, thus limiting potential for any adverse effects from consumption of dietary RNA. A recent paper reported that putative plant miRNAs were detected in animal plasma and serum, presumably acquired through ingestion, and may have a functional impact in the consuming organisms. RESULTS: To address the question of how common this phenomenon could be, we searched for plant miRNAs sequences in public sRNA datasets from various tissues of mammals, chicken and insects. Our analyses revealed that plant miRNAs were present in the animal sRNA datasets, and significantly miR168 was extremely over-represented. Furthermore, all or nearly all (>96%) miR168 sequences were monocot derived for most datasets, including datasets for two insects reared on dicot plants in their respective experiments. To investigate if plant-derived miRNAs, including miR168, could accumulate and move systemically in insects, we conducted insect feeding studies for three insects including corn rootworm, which has been shown to be responsive to plant-produced long double-stranded RNAs. CONCLUSIONS: Our analyses suggest that the observed plant miRNAs in animal sRNA datasets can originate in the process of sequencing, and that accumulation of plant miRNAs via dietary exposure is not universal in animals.     

Photooxidative stress‐induced and abundant small RNAs in Rhodobacter sphaeroides

Exposure to oxygen and light generates photooxidative stress by the bacteriochlorophyll a mediated formation of singlet oxygen (¹O₂) in Rhodobacter sphaeroides. Our study reports the genome‐wide search for small RNAs (sRNAs) involved in the regulatory response to ¹O₂. By using 454 pyrosequencing and Northern blot analysis, we identified 20 sRNAs from R. sphaeroides aerobic cultures or following treatment with ¹O₂ or superoxide (O^–₂). One sRNA was specifically induced by ¹O₂ and its expression depends on the extracytoplasmic function sigma factor RpoE. Two sRNAs induced by ¹O₂ and O^–₂ were cotranscribed with upstream genes preceded by promoters with target sequences for the alternative sigma factors RpoH_I and RpoH_II. The most abundant sRNA was processed in the presence of ¹O₂ but not by O^–₂. From this and a second sRNA a conserved 3′‐segment accumulated from a larger precursor. Absence of the RNA chaperone Hfq changed the half‐lives, abundance and processing of ¹O₂‐affected sRNAs. Orthologues of three sRNA genes are present in different alpha‐proteobacteria, but the majority was unique to R. sphaeroides or Rhodobacterales species. Our discovery that abundant sRNAs are affected by ¹O₂ exposure extends the knowledge on the role of sRNAs and Hfq in the regulatory response to oxidative stress.