A global profiling of uncapped mRNAs under cold stress reveals specific decay patterns and endonucleolytic cleavages in Brachypodium distachyon

Degradome sequencing reveals an integrative miRNA-mediated gene interaction network regulating rice seed vigor

BACKGROUND

It is well known that seed vigor is essential for agricultural production and rice (Oryza sativa L.) is one of the most important crops in the world. Though we previously reported that miR164c regulates rice seed vigor, but whether and how other miRNAs cooperate with miR164c to regulate seed vigor is still unknown.

RESULTS

Based on degradome data of six RNA samples isolated from seeds of the wild-type (WT) indica rice cultivar 'Kasalath' as well as two modified lines in 'Kasalath' background (miR164c-silenced line [MIM164c] and miR164c overexpression line [OE164c]), which were subjected to either no aging treatment or an 8-day artificial aging treatment, 1247 different target transcripts potentially cleaved by 421 miRNAs were identified. The miRNA target genes were functionally annotated via GO and KEGG enrichment analyses. By STRING database assay, a miRNA-mediated gene interaction network regulating seed vigor in rice was revealed, which comprised at least four interconnected pathways: the miR5075-mediated oxidoreductase related pathway, the plant hormone related pathway, the miR164e related pathway, and the previously reported RPS27AA related pathway. Knockout and overexpression of the target gene Os02g0817500 of miR5075 decreased and enhanced seed vigor, respectively. By Y2H assay, the proteins encoded by five seed vigor-related genes, Os08g0295100, Os07g0633100, REFA1, OsPER1 and OsGAPC3, were identified to interact with Os02g0817500.

CONCLUSIONS

miRNAs cooperate to regulate seed vigor in rice via an integrative gene interaction network comprising miRNA target genes and other functional genes. The result provided a basis for fully understanding the molecular mechanisms of seed vigor regulation.

The complete mitochondrial genome of Choroterpes (Euthralus) yixingensis (Ephemeroptera: Leptophlebiidae) and its mitochondrial protein-coding gene expression under imidacloprid stress

As one of the most common benthic invertebrates in freshwater, mayflies are very sensitive to changes in water quality and have high requirements for the water environment to allow their nymphs to successfully live and grow. Neonicotinoids, such as imidacloprid, can enter fresh water and pollute the aquatic environment. The present study had two goals: (1) investigate imidacloprid effects on mayfly larvae Choroterpes (Euthralus) yixingensis, and (2) contribute to the phylogenetic status of Ephemeroptera that has always been controversial. Nymphs were collected from Jinhua, China and exposed to different concentrations imidacloprid (5, 10, 20, and 40 μg/L) in the laboratory. Survival of C. yixingensis nymphs decreased as a function of time and imidacloprid concentration with only ~ 55% survival after 72 h exposure to 40 μg/L imidacloprid. After culture under 40 μg/L imidacloprid for 24 h, the steady state transcript levels of mitochondrial COX3, ND4 and ND4L genes were reduced to just 0.07 ± 0.11, 0.30 ± 0.16, and 0.28 ± 0.13 as compared with respective control values (P < 0.01). Steady state transcript levels of ND4 and ND4L were also significantly reduced in a dose-dependent manner (P < 0.05), suggesting that the steady state transcript pattern of these genes in mayfly nymphs can change in response to different levels of environmental contamination. Hence, the mitochondrial protein-coding genes of mayflies could potentially be developed as biomarkers for water ecotoxicity monitoring in the future. In addition, we used the mitochondrial genome sequence of C. yixingensis for an assessment of the phylogenetic tree of Ephemeroptera. The monophyly of Leptophlebiidae was supported and showed that Leptophlebiidae was a sister group to the clade (Baetidae + Caenidae).

Global Analysis of the Human RNA Degradome Reveals Widespread Decapped and Endonucleolytic Cleaved Transcripts

RNA decay is an important regulatory mechanism for gene expression at the posttranscriptional level. Although the main pathways and major enzymes that facilitate this process are well defined, global analysis of RNA turnover remains under-investigated. Recent advances in the application of next-generation sequencing technology enable its use in order to examine various RNA decay patterns at the genome-wide scale. In this study, we investigated human RNA decay patterns using parallel analysis of RNA end-sequencing (PARE-seq) data from XRN1-knockdown HeLa cell lines, followed by a comparison of steady state and degraded mRNA levels from RNA-seq and PARE-seq data, respectively. The results revealed 1103 and 1347 transcripts classified as stable and unstable candidates, respectively. Of the unstable candidates, we found that a subset of the replication-dependent histone transcripts was polyadenylated and rapidly degraded. Additionally, we identified 380 endonucleolytically cleaved candidates by analyzing the most abundant PARE sequence on a transcript. Of these, 41.4% of genes were classified as unstable genes, which implied that their endonucleolytic cleavage might affect their mRNA stability. Furthermore, we identified 1877 decapped candidates, including HSP90B1 and SWI5, having the most abundant PARE sequences at the 5′-end positions of the transcripts. These results provide a useful resource for further analysis of RNA decay patterns in human cells.

Early dynamics of photosynthetic Lhcf2 and Lhcf15 transcription and mRNA stabilities in response to herbivory-related decadienal in Phaeodactylum tricornutum

Abiotic and biotic stresses widely reduce light harvesting complex (LHC) gene expression in higher plants and algae. However, control mechanisms and functions of these changes are not well understood. During herbivory, marine diatom species release oxylipins that impair grazer reproduction and serve as signaling molecules to nearby undamaged diatoms. To examine LHC mRNA regulation by oxylipin exposure, the diatom Phaeodactylum tricornutum was treated with a sublethal concentration of trans,trans-2,4-decadienal (DD) during the light cycle. Transcriptome analyses revealed extensive suppression of LHC mRNAs and a smaller set of up-regulated LHC mRNAs at 3 h. For two divergently regulated LHCF antennae family mRNAs, in vivo 4-thiouracil metabolic labeling was used to distinguish synthesis and degradation rates. Within 3 h of DD exposure, Lhcf2 mRNA levels and transcription were strongly suppressed and its mRNA half-life decreased. In contrast, Lhcf15 mRNA mainly accumulated between 3-9 h, its transcription increased and its mRNA was highly stabilized. Hence, DD-treated cells utilized transcriptional and mRNA stability control mechanisms which were likely major factors in the differing Lhcf2 and Lhcf15 expression patterns. Widespread LHC mRNA regulation and possible effects on photosynthesis may contribute to enhanced fitness in cells impacted by herbivory and other stresses.

RNA Regulation in Plant Cold Stress Response

In addition to plants, all organisms react to environmental stimuli via the perception of signals and subsequently respond through alterations of gene expression. However, genes/mRNAs are usually not the functional unit themselves, and instead, resultant protein products with individual functions result in various acquired phenotypes. In order to fully characterize the adaptive responses of plants to environmental stimuli, it is essential to determine the level of proteins, in addition to the regulation of mRNA expression. This regulatory step, which is referred to as "mRNA posttranscriptional regulation," occurs subsequent to mRNA transcription and prior to translation. Although these RNA regulatory mechanisms have been well-studied in many organisms, including plants, it is not fully understood how plants respond to environmental stimuli, such as cold stress, via these RNA regulations.A recent study described several RNA regulatory factors in relation to environmental stress responses, including plant cold stress tolerance. In this chapter, the functions of RNA regulatory factors and comprehensive analyses related to the RNA regulations involved in cold stress response are summarized, such as mRNA maturation, including capping, splicing, polyadenylation of mRNA, and the quality control system of mRNA; mRNA degradation, including the decapping step; and mRNA stabilization. In addition, the putative roles of messenger ribonucleoprotein (mRNP) granules, such as processing bodies (PBs) and stress granules (SGs), which are cytoplasmic particles, are described in relation to RNA regulations under stress conditions. These RNA regulatory systems are important for adjusting or fine-tuning and determining the final levels of mRNAs and proteins in order to adapt or respond to environmental stresses. Collectively, these new areas of study revealed that plants possess precise novel regulatory mechanisms which specifically function in the response to cold stress.

Analysis of Brachypodium miRNA targets: evidence for diverse control during stress and conservation in bioenergy crops

BACKGROUND

Since the proposal of Brachypodium distachyon as a model for the grasses, over 500 Bdi-miRNAs have been annotated in miRBase making Brachypodium second in number only to rice. Other monocots, such as switchgrass, are completely absent from the miRBase database. While a significant number of miRNAs have been identified which are highly conserved across plants, little research has been done with respect to the conservation of miRNA targets. Plant responses to abiotic stresses are regulated by diverse pathways many of which involve miRNAs; however, it can be difficult to identify miRNA guided gene regulation when the miRNA is not the primary regulator of the target mRNA.

RESULTS

To investigate miRNA target conservation and stress response involvement, a set of PARE (Parallel Analysis of RNA Ends) libraries totaling over two billion reads was constructed and sequenced from Brachypodium, switchgrass, and sorghum representing the first report of RNA degradome data from the latter two species. Analysis of this data provided not only PARE evidence for miRNA guided cleavage of over 7000 predicted target mRNAs in Brachypodium, but also evidence for miRNA guided cleavage of over 1000 homologous transcripts in sorghum and switchgrass. A pipeline was constructed to compare RNA-seq and PARE data made from Brachypodium plants exposed to various abiotic stress conditions. This resulted in the identification of 44 miRNA targets which exhibit stress regulated cleavage. Time course experiments were performed to reveal the relationship between miR393ab, miR169a, miR394ab, and their respective targets throughout the first 36 h of the cold stress response in Brachypodium.

CONCLUSIONS

Knowledge gained from this study provides considerable insight into the RNA degradomes and the breadth of miRNA target conservation among these three species. Additionally, associations of a number of miRNAs and target mRNAs with the stress responses have been revealed which could aid in the development of stress tolerant transgenic crops.

Rapid Recovery Gene Downregulation during Excess-Light Stress and Recovery in Arabidopsis

Stress recovery may prove to be a promising approach to increase plant performance and, theoretically, mRNA instability may facilitate faster recovery. Transcriptome (RNA-seq, qPCR, sRNA-seq, and PARE) and methylome profiling during repeated excess-light stress and recovery was performed at intervals as short as 3 min. We demonstrate that 87% of the stress-upregulated mRNAs analyzed exhibit very rapid recovery. For instance, HSP101 abundance declined 2-fold every 5.1 min. We term this phenomenon rapid recovery gene downregulation (RRGD), whereby mRNA abundance rapidly decreases promoting transcriptome resetting. Decay constants (k) were modeled using two strategies, linear and nonlinear least squares regressions, with the latter accounting for both transcription and degradation. This revealed extremely short half-lives ranging from 2.7 to 60.0 min for 222 genes. Ribosome footprinting using degradome data demonstrated RRGD loci undergo cotranslational decay and identified changes in the ribosome stalling index during stress and recovery. However, small RNAs and 5'-3' RNA decay were not essential for recovery of the transcripts examined, nor were any of the six excess light-associated methylome changes. We observed recovery-specific gene expression networks upon return to favorable conditions and six transcriptional memory types. In summary, rapid transcriptome resetting is reported in the context of active recovery and cellular memory.

Co-ordinated Regulations of mRNA Synthesis and Decay during Cold Acclimation in Arabidopsis Cells

Plants possess a cold acclimation system to acquire freezing tolerance through pre-exposure to non-freezing low temperatures. The transcriptional cascade of C-repeat-binding factors (CBFs)/dehydration response element-binding factors (DREBs) is considered a major transcriptional regulatory pathway during cold acclimation. However, little is known regarding the functional significance of mRNA stability regulation in the response of gene expression to cold stress. The actual level of individual mRNAs is determined by a balance between mRNA synthesis and degradation. Therefore, it is important to assess the regulatory steps to increase our understanding of gene regulation. Here, we analyzed temporal changes in mRNA amounts and half-lives in response to cold stress in Arabidopsis cell cultures based on genome-wide analysis. In this mRNA decay array method, mRNA half-life measurements and microarray analyses were combined. In addition, temporal changes in the integrated value of transcription rates were estimated from the above two parameters using a mathematical approach. Our results showed that several cold-responsive genes, including Cold-regulated 15a, were relatively destabilized, whereas the mRNA amounts were increased during cold treatment by accelerating the transcription rate to overcome the destabilization. Considering the kinetics of mRNA synthesis and degradation, this apparently contradictory result supports that mRNA destabilization is advantageous for the swift increase in CBF-responsive genes in response to cold stress.

Transcriptome and Degradome Sequencing Reveals Dormancy Mechanisms of Cunninghamia lanceolata Seeds

Seeds with physiological dormancy usually experience primary and secondary dormancy in the nature; however, little is known about the differential regulation of primary and secondary dormancy. We combined multiple approaches to investigate cytological changes, hormonal levels, and gene expression dynamics in Cunninghamia lanceolata seeds during primary dormancy release and secondary dormancy induction. Light microscopy and transmission electron microscopy revealed that protein bodies in the embryo cells coalesced during primary dormancy release and then separated during secondary dormancy induction. Transcriptomic profiling demonstrated that expression of genes negatively regulating gibberellic acid (GA) sensitivity reduced specifically during primary dormancy release, whereas the expression of genes positively regulating abscisic acid (ABA) biosynthesis increased during secondary dormancy induction. Parallel analysis of RNA ends revealed uncapped transcripts for ∼55% of all unigenes. A negative correlation between fold changes in expression levels of uncapped versus capped mRNAs was observed during primary dormancy release. However, this correlation was loose during secondary dormancy induction. Our analyses suggest that the reversible changes in cytology and gene expression during dormancy release and induction are related to ABA/GA balance. Moreover, mRNA degradation functions as a critical posttranscriptional regulator during primary dormancy release. These findings provide a mechanistic framework for understanding physiological dormancy in seeds.

Exploring drought stress-regulated genes in senna (Cassia angustifolia Vahl.): a transcriptomic approach

De novo assembly of reads produced by next-generation sequencing (NGS) technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. Senna (Cassia angustifolia Vahl.) is a drought-tolerant annual undershrub of Caesalpiniaceae, a subfamily of Fabaceae. There are insufficient transcriptomic and genomic data in public databases for understanding the molecular mechanism underlying the drought tolerance of senna. Therefore, the main purpose of this study was to know the transcriptome profile of senna, with special reference to drought stress. RNA from two different stages of leaf development was extracted and sequenced separately using the Illumina technology. A total of 200 million reads were generated, and a de novo assembly of processed reads in the pooled transcriptome using Trinity yielded 43,413 transcripts which were further annotated using NCBI BLAST with "green plant database (txid 33090)," Swiss Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). Out of the total transcripts, 42,280 (95.0 %) were annotated by BLASTX against the green plant database of NCBI. Senna transcriptome showed the highest similarity to Glycine max (41 %), followed by Phaseolus vulgaris (16 %), Cicer arietinum (15 %), and Medicago trancatula (5 %). The highest number of GO terms were enriched for the molecular functions category; of these "catalytic activity" (GO: 0003824) (25.10 %) and "binding activity" (GO: 0005488) (20.10 %) were most abundantly represented. We used InterProscan to see protein similarity at domain level; a total of 33,256 transcripts were annotated against the Pfam domains. The transcripts were assigned with various KEGG pathways. Coding DNA sequences (CDS) encoding various drought stress-regulated pathways such as signaling factors, protein-modifying/degrading enzymes, biosynthesis of phytohormone, phytohormone signaling, osmotically active compounds, free radical scavengers, chlorophyll metabolism, leaf cuticular wax, polyamines, and protective proteins were identified through BLASTX search. The lucine-rich repeat kinase family was the most abundantly found group of protein kinases. Orphan, bHLH, and bZIP family TFs were the most abundantly found in senna. Six genes encoding MYC2 transcription factor, 9-cis-epoxycarotenoid dioxygenase (NCED), l -ascorbate peroxidase (APX), aminocyclopropane carboxylate oxidase (ACO), abscisic acid 8'-hydroxylase (ABA), and WRKY transcription factor were confirmed through reverse transcriptase-PCR (RT-PCR) and Sanger sequencing for the first time in senna. The potential drought stress-related transcripts identified in this study provide a good start for further investigation into the drought adaptation in senna. Additionally, our transcriptome sequences are the valuable resource for accelerated genomics-assisted genetic improvement programs and facilitate manipulation of biochemical pathways for developing drought-tolerant genotypes of crop plants.

Ultra-fast alterations in mRNA levels uncover multiple players in light stress acclimation in plants

The acclimation of plants to changes in light intensity requires rapid responses at several different levels. These include biochemical and biophysical responses as well as alterations in the steady-state level of different transcripts and proteins. Recent studies utilizing promoter::reporter constructs suggested that transcriptional responses to changes in light intensity could occur within seconds, rates for which changes in mRNA expression are not routinely measured or functionally studied. To identify and characterize rapid changes in the steady-state level of different transcripts in response to light stress we performed RNA sequencing analysis of Arabidopsis thaliana plants subjected to light stress. Here we report that mRNA accumulation of 731 transcripts occurs as early as 20-60 sec following light stress application, and that at least five of these early response transcripts play an important biological role in the acclimation of plants to light stress. More than 20% of transcripts accumulating in plants within 20-60 sec of initiation of light stress are H2 O2 - and ABA-response transcripts, and the accumulation of several of these transcripts is inhibited by transcriptional inhibitors. In accordance with the association of rapid response transcripts with H2 O2 and ABA signaling, a mutant impaired in ABA sensing (abi-1) was found to be more tolerant to light stress, and the response of several of the rapid response transcripts was altered in mutants impaired in reactive oxygen metabolism. Our findings reveal that transcriptome reprogramming in plants could occur within seconds of initiation of abiotic stress and that this response could invoke known as well as unknown proteins and pathways.

Global analysis of uncapped mRNA changes under drought stress and microRNA-dependent endonucleolytic cleavages in foxtail millet

BACKGROUND

mRNA degradation plays an important role in the determination of mRNA abundance and can quickly regulate gene expression. The production of uncapped mRNAs, an important mechanism of mRNA degradation, can be initiated by decapping enzymes, endonucleases or small RNAs such as microRNAs (miRNAs). Little is known, however, about the role of uncapped mRNAs in plants under environmental stress.

RESULTS

Using a novel approach called parallel analysis of RNA ends (PARE), we performed a global study of uncapped mRNAs under drought stress in foxtail millet (Setaria italica [L.] P. Beauv.). When both gene degradation (PARE) and gene transcription (RNA-sequencing) data were considered, four types of mRNA decay patterns were identified under drought stress. In addition, 385 miRNA-target interactions were identified in the PARE data using PAREsnip. The PARE analysis also suggested that two miRNA hairpin processing mechanisms--loop-last and loop-first processing--operate in foxtail millet, with both miR319 and miR156 gene families undergoing precise processing via the unusual loop-first mechanism. Finally, we found 11 C4 photosynthesis-related enzymes encoded by drought-responsive genes.

CONCLUSIONS

We performed a global analysis of mRNA degradation under drought stress and uncovered diverse drought-response mechanisms in foxtail millet. This information will deepen our understanding of mRNA expression under stressful environmental conditions in gramineous plants. In addition, PARE analysis identified many miRNA targets and revealed miRNA-precursor processing modes in foxtail millet.

Mechanism of cytoplasmic mRNA translation

Protein synthesis is a fundamental process in gene expression that depends upon the abundance and accessibility of the mRNA transcript as well as the activity of many protein and RNA-protein complexes. Here we focus on the intricate mechanics of mRNA translation in the cytoplasm of higher plants. This chapter includes an inventory of the plant translational apparatus and a detailed review of the translational processes of initiation, elongation, and termination. The majority of mechanistic studies of cytoplasmic translation have been carried out in yeast and mammalian systems. The factors and mechanisms of translation are for the most part conserved across eukaryotes; however, some distinctions are known to exist in plants. A comprehensive understanding of the complex translational apparatus and its regulation in plants is warranted, as the modulation of protein production is critical to development, environmental plasticity and biomass yield in diverse ecosystems and agricultural settings.

Tudor staphylococcal nuclease links formation of stress granules and processing bodies with mRNA catabolism in Arabidopsis

Tudor Staphylococcal Nuclease (TSN or Tudor-SN; also known as SND1) is an evolutionarily conserved protein involved in the transcriptional and posttranscriptional regulation of gene expression in animals. Although TSN was found to be indispensable for normal plant development and stress tolerance, the molecular mechanisms underlying these functions remain elusive. Here, we show that Arabidopsis thaliana TSN is essential for the integrity and function of cytoplasmic messenger ribonucleoprotein (mRNP) complexes called stress granules (SGs) and processing bodies (PBs), sites of posttranscriptional gene regulation during stress. TSN associates with SGs following their microtubule-dependent assembly and plays a scaffolding role in both SGs and PBs. The enzymatically active tandem repeat of four SN domains is crucial for targeting TSN to the cytoplasmic mRNA complexes and is sufficient for the cytoprotective function of TSN during stress. Furthermore, our work connects the cytoprotective function of TSN with its positive role in stress-induced mRNA decapping. While stress led to a pronounced increase in the accumulation of uncapped mRNAs in wild-type plants, this increase was abrogated in TSN knockout plants. Taken together, our results establish TSN as a key enzymatic component of the catabolic machinery responsible for the processing of mRNAs in the cytoplasmic mRNP complexes during stress.

Identification of SMG6 cleavage sites and a preferred RNA cleavage motif by global analysis of endogenous NMD targets in human cells

In metazoans, cleavage by the endoribonuclease SMG6 is often the first degradative event in non-sense-mediated mRNA decay (NMD). However, the exact sites of SMG6 cleavage have yet to be determined for any endogenous targets, and most evidence as to the identity of SMG6 substrates is indirect. Here, we use Parallel Analysis of RNA Ends to specifically identify the 5′ termini of decay intermediates whose production is dependent on SMG6 and the universal NMD factor UPF1. In this manner, the SMG6 cleavage sites in hundreds of endogenous NMD targets in human cells have been mapped at high resolution. In addition, a preferred sequence motif spanning most SMG6 cleavage sites has been discovered and validated by mutational analysis. For many SMG6 substrates, depletion of SMG6 resulted in the accumulation of decapped transcripts, an effect indicative of competition between SMG6-dependent and SMG6-independent NMD pathways. These findings provide key insights into the mechanisms by which mRNAs targeted by NMD are degraded.

De novo transcriptome assembly and comparative analysis of differentially expressed genes in Prunus dulcis Mill. in response to freezing stress

Almond (Prunus dulcis Mill.), one of the most important nut crops, requires chilling during winter to develop fruiting buds. However, early spring chilling and late spring frost may damage the reproductive tissues leading to reduction in the rate of productivity. Despite the importance of transcriptional changes and regulation, little is known about the almond’s transcriptome under the cold stress conditions. In the current reserch, we used RNA-seq technique to study the response of the reporuductive tissues of almond (anther and ovary) to frost stress. RNA sequencing resulted in more than 20 million reads from anther and ovary tissues of almond, individually. About 40,000 contigs were assembled and annotated denovo in each tissue. Profile of gene expression in ovary showed significant alterations in 5,112 genes, whereas in anther 6,926 genes were affected by freezing stress. Around two thousands of these genes were common altered genes in both ovary and anther libraries. Gene ontology indicated the involvement of differentially expressed (DE) genes, responding to freezing stress, in metabolic and cellular processes. qRT-PCR analysis verified the expression pattern of eight genes randomley selected from the DE genes. In conclusion, the almond gene index assembled in this study and the reported DE genes can provide great insights on responses of almond and other Prunus species to abiotic stresses. The obtained results from current research would add to the limited available information on almond and Rosaceae. Besides, the findings would be very useful for comparative studies as the number of DE genes reported here is much higher than that of any previous reports in this plant.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon

BACKGROUND

The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon.

RESULTS

B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner.

CONCLUSIONS

B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.