Characterization of unique small RNA populations from rice grain

Dicer-like2b suppresses the wiry leaf phenotype in tomato induced by tobacco mosaic virus

Dicer-like (DCL) proteins are principal components of RNA silencing, a major defense mechanism against plant virus infections. However, their functions in suppressing virus-induced disease phenotypes remain largely unknown. Here, we identified a role for tomato (Solanum lycopersicum) DCL2b in regulating the wiry leaf phenotype during defense against tobacco mosaic virus (TMV). Knocking out SlyDCL2b promoted TMV accumulation in the leaf primordium, resulting in a wiry phenotype in distal leaves. Biochemical and bioinformatics analyses showed that 22-nt virus-derived small interfering RNAs (vsiRNAs) accumulated less abundantly in slydcl2b mutants than in wild-type plants, suggesting that SlyDCL2b-dependent 22-nt vsiRNAs are required to exclude virus from leaf primordia. Moreover, the wiry leaf phenotype was accompanied by upregulation of Auxin Response Factors (ARFs), resulting from a reduction in trans-acting siRNAs targeting ARFs (tasiARFs) in TMV-infected slydcl2b mutants. Loss of tasiARF production in the slydcl2b mutant was in turn caused by inhibition of miRNA390b function. Importantly, silencing SlyARF3 and SlyARF4 largely restored the wiry phenotype in TMV-infected slydcl2b mutants. Our work exemplifies the complex relationship between RNA viruses and the endogenous RNA silencing machinery, whereby SlyDCL2b protects the normal development of newly emerging organs by excluding virus from these regions and thus maintaining developmental silencing.

Targeted suppression of gibberellin biosynthetic genes ZmGA20ox3 and ZmGA20ox5 produces a short stature maize ideotype

Maize is one of the world's most widely cultivated crops. As future demands for maize will continue to rise, fields will face ever more frequent and extreme weather patterns that directly affect crop productivity. Development of environmentally resilient crops with improved standability in the field, like wheat and rice, was enabled by shifting the architecture of plants to a short stature ideotype. However, such architectural change has not been implemented in maize due to the unique interactions between gibberellin (GA) and floral morphology which limited the use of the same type of mutations as in rice and wheat. Here, we report the development of a short stature maize ideotype in commercial hybrid germplasm, which was generated by targeted suppression of the biosynthetic pathway for GA. To accomplish this, we utilized a dominant, miRNA-based construct expressed in a hemizygous state to selectively reduce expression of the ZmGA20ox3 and ZmGA20ox5 genes that control GA biosynthesis primarily in vegetative tissues. Suppression of both genes resulted in the reduction of GA levels leading to inhibition of cell elongation in internodal tissues, which reduced plant height. Expression of the miRNA did not alter GA levels in reproductive tissues, and thus, the reproductive potential of the plants remained unchanged. As a result, we developed a dominant, short-stature maize ideotype that is conducive for the commercial production of hybrid maize. We expect that the new maize ideotype would enable more efficient and more sustainable maize farming for a growing world population.

TRANS-ACTING SIRNA3-derived short interfering RNAs confer cleavage of mRNAs in rice

Plant TRANS-ACTING SIRNA3 (TAS3)-derived short interfering RNAs (siRNAs) include tasiR-AUXIN RESPONSE FACTORs (ARFs), which are functionally conserved in targeting ARF genes, and a set of non-tasiR-ARF siRNAs, which have rarely been studied. In this study, TAS3 siRNAs were systematically characterized in rice (Oryza sativa). Small RNA sequencing results showed that an overwhelming majority of TAS3 siRNAs belong to the non-tasiR-ARF group, while tasiR-ARFs occupy a diminutive fraction. Phylogenetic analysis of TAS3 genes across dicot and monocot plants revealed that the siRNA-generating regions were highly conserved in grass species, especially in the Oryzoideae. Target genes were identified for not only tasiR-ARFs but also non-tasiR-ARF siRNAs by analyzing rice Parallel Analysis of RNA Ends datasets, and some of these siRNA-target interactions were experimentally confirmed using tas3 mutants generated by genome editing. Consistent with the de-repression of target genes, phenotypic alterations were observed for mutants in three TAS3 loci in comparison to wild-type rice. The regulatory role of ribosomes in the TAS3 siRNA-target interactions was further revealed by the fact that TAS3 siRNA-mediated target cleavage, in particular tasiR-ARFs targeting ARF2/3/14/15, occurred extensively in rice polysome samples. Altogether, our study sheds light into TAS3 genes in plants and expands our knowledge about rice TAS3 siRNA-target interactions.

RNAi for Western Corn Rootworm Management: Lessons Learned, Challenges, and Future Directions

The western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte, is considered one of the most economically important pests of maize (Zea mays L.) in the United States (U.S.) Corn Belt with costs of management and yield losses exceeding USD ~1-2 billion annually. WCR management has proven challenging given the ability of this insect to evolve resistance to multiple management strategies including synthetic insecticides, cultural practices, and plant-incorporated protectants, generating a constant need to develop new management tools. One of the most recent developments is maize expressing double-stranded hairpin RNA structures targeting housekeeping genes, which triggers an RNA interference (RNAi) response and eventually leads to insect death. Following the first description of in planta RNAi in 2007, traits targeting multiple genes have been explored. In June 2017, the U.S. Environmental Protection Agency approved the first in planta RNAi product against insects for commercial use. This product expresses a dsRNA targeting the WCR snf7 gene in combination with Bt proteins (Cry3Bb1 and Cry34Ab1/Cry35Ab1) to improve trait durability and will be introduced for commercial use in 2022.

Comparative Analysis of miRNA Expression Profiles Between Heat-Tolerant and Heat-Sensitive Genotypes of Flowering Chinese Cabbage Under Heat Stress Using High-Throughput Sequencing

Heat stress disturbs cellular homeostasis, thus usually impairs yield of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). MicroRNAs (miRNAs) play a significant role in plant responses to different stresses by modulating gene expression at the post-transcriptional level. However, the roles that miRNAs and their target genes may play in heat tolerance of flowering Chinese cabbage remain poorly characterized. The current study sequenced six small RNA libraries generated from leaf tissues of flowering Chinese cabbage collected at 0, 6, and 12 h after 38 °C heat treatment, and identified 49 putative novel miRNAs and 43 known miRNAs that differentially expressed between heat-tolerant and heat-sensitive flowering Chinese cabbage. Among them, 14 novel and nine known miRNAs differentially expressed only in the heat-tolerant genotype under heat-stress, therefore, their target genes including disease resistance protein TAO1-like, RPS6, reticuline oxidase-like protein, etc. might play important roles in enhancing heat-tolerance. Gene Ontology (GO) analysis revealed that targets of these differentially expressed miRNAs may play key roles in responses to temperature stimulus, cell part, cellular process, cell, membrane, biological regulation, binding, and catalytic activities. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified their important functions in signal transduction, environmental adaptation, global and overview maps, as well as in stress adaptation and in MAPK signaling pathways such as cell death. These findings provide insight into the functions of the miRNAs in heat stress tolerance of flowering Chinese cabbage.

Safety Considerations for Humans and Other Vertebrates Regarding Agricultural Uses of Externally Applied RNA Molecules

The potential of double-stranded RNAs (dsRNAs) for use as topical biopesticides in agriculture was recently discussed during an OECD (Organisation for Economic Co-operation and Development) Conference on RNA interference (RNAi)-based pesticides. Several topics were presented and these covered different aspects of RNAi technology, its application, and its potential effects on target and non-target organisms (including both mammals and non-mammals). This review presents information relating to RNAi mechanisms in vertebrates, the history of safe RNA consumption, the biological barriers that contribute to the safety of its consumption, and effects related to humans and other vertebrates as discussed during the conference. We also review literature related to vertebrates exposed to RNA molecules and further consider human health safety assessments of RNAi-based biopesticides. This includes possible routes of exposure other than the ingestion of potential residual material in food and water (such as dermal and inhalation exposures during application in the field), the implications of different types of formulations and RNA structures, and the possibility of non-specific effects such as the activation of the innate immune system or saturation of the RNAi machinery.

Identification and Characterization of Stress-Responsive TAS3-Derived TasiRNAs in Melon

Small interfering RNAs (siRNA) are key regulators of gene expression that play essential roles in diverse biological processes. Trans-acting siRNAs (tasiRNAs) are a class of plant-endogenous siRNAs that lead the cleavage of nonidentical transcripts. TasiRNAs are usually involved in fine-tuning development. However, increasing evidence supports that tasiRNAs may be involved in stress response. Melon is a crop of great economic importance extensively cultivated in semiarid regions frequently exposed to changing environmental conditions that limit its productivity. However, knowledge of the precise role of siRNAs in general, and of tasiRNAs in particular, in regulating the response to adverse environmental conditions is limited. Here, we provide the first comprehensive analysis of computationally inferred melon-tasiRNAs responsive to two biotic (viroid-infection) and abiotic (cold treatment) stress conditions. We identify two TAS3-loci encoding to length (TAS3-L) and short (TAS3-S) transcripts. The TAS candidates predicted from small RNA-sequencing data were characterized according to their chromosome localization and expression pattern in response to stress. The functional activity of cmTAS genes was validated by transcript quantification and degradome assays of the tasiRNA precursors and their predicted targets. Finally, the functionality of a representative cmTAS3-derived tasiRNA (TAS3-S) was confirmed by transient assays showing the cleavage of ARF target transcripts.

Molecular insights into the regulation of rice kernel elongation

A large number of rice agronomic traits are complex, multi factorial and polygenic. As the mechanisms and genes determining grain size and yield are largely unknown, the identification of regulatory genes related to grain development remains a preeminent approach in rice genetic studies and breeding programs. Genes regulating cell proliferation and expansion in spikelet hulls and participating in endosperm development are the main controllers of rice kernel elongation and grain size. We review here and discuss recent findings on genes controlling rice grain size and the mechanisms, epialleles, epigenomic variation, and assessment of controlling genes using genome-editing tools relating to kernel elongation.

Deep sequencing identified potential miRNAs involved in defence response, stress and plant growth characteristics of wild genotypes of cardamom

Cardamom has long been used as a food flavouring agent and in ayurvedic medicines for mouth ulcers, digestive problems and even depression. Extensive occurrence of pests and diseases adversely affect its cultivation and result in substantial reductions in total production and productivity. Numerous studies revealed the significant role of miRNAs in plant biotic stress responses. In the current study, miRNA profiling of cultivar and wild cardamom genotypes was performed using an Ion Proton sequencer. We identified 161 potential miRNAs representing 42 families, including monocot/tissue-specific and 14 novel miRNAs in both genotypes. Significant differences in miRNA family abundance between the libraries were observed in read frequencies. A total of 19 miRNAs (from known miRNAs) displayed a twofold difference in expression between wild and cultivar genotypes. We found 1168 unique potential targets for 40 known miRNA families in wild and 1025 potential targets for 42 known miRNA families in cultivar genotypes. The differential expression analysis revealed that most miRNAs identified were highly expressed in cultivars and, furthermore, lower expression of miR169 and higher expression of miR529 in wild cardamom proved evidence that wild genotypes have stronger drought stress tolerance and floral development than cultivars. Potential targets predicted for the newly identified miRNAs from the miRNA libraries of wild and cultivar cardamom genotypes involved in metabolic and developmental processes and in response to various stimuli. qRT-PCR confirmed miRNAs were differentially expressed between wild and cultivar genotypes. Furthermore, four target genes were validated experimentally to confirm miRNA-mRNA target pairing using RNA ligase-mediated 5' Rapid Amplification of cDNA Ends (5'RLM-RACE) PCR.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Identification and validation of a virus-inducible ta-siRNA-generating TAS4 locus in tomato

Trans-acting small interfering RNAs (ta-siRNAs) are a class of endogenous small RNA, associated with posttranscriptional gene silencing. Their biogenesis requires an initial microRNA (miRNA)-mediated cleavage of precursor RNA. Around 20 different ta-siRNA-producing loci (TASs), whose sequences are conserved, are reported in plants. In tomato, two TAS gene families have been identified, which are found to target auxin response factor gene and bacterial spot disease resistance protein Bs4 gene. Using high-throughput computational and experimental approach, we identified a new locus-producing ta-siRNA in tomato. We have also identified the putative miRNA regulating the production of ta-siRNA from this locus. The ta-siRNAs generated from TAS4 were up-regulated upon infection with a DNA virus. The potential targets of ta-siRNAs were predicted to be variety of proteins including MYB transcription factors and cell cycle regulators for some of the ta-siRNAs produced.

Generation of different sizes and classes of small RNAs in barley is locus, chromosome and/or cultivar-dependent

Background

Various small RNA (sRNA) sizes and varieties have been identified, but their relationship as well as relationship with their origins and allocations have not been well understood or investigated.

Results

By comparing sRNAs generated from two barley cultivars, Golden Promise (GP) and Pallas, we identified that the generation of different sizes and types of sRNAs in barley was locus-, chromosome- and/or cultivar-dependent. 20-nt sRNAs mainly comprising miRNAs and chloroplast-derived sRNAs were significantly over-expressed in Pallas vs. GP on chromosomes 3H and 6H. MiRNAs-enriched 21-nt sRNAs were significantly over-expressed in Pallas vs. GP only on chromosome 4H. On chromosome 5H this size of sRNAs was significantly under-expressed in Pallas, so were 22-nt sRNAs mainly comprising miRNAs and repeat-derived sRNAs. 24-nt sRNAs mostly derived from repeats were evenly distributed in all chromosomes and expressed similarly between GP and Pallas. Unlike other sizes of sRNAs, 24-nt sRNAs were little conserved in other plant species. Abundant sRNAs were mostly generated from 3’ terminal regions of chromosome 1H and 5’ terminal regions of chromosome 5H. Over-expressed miRNAs in GP vs. Pallas primarily function in stress responses and iron-binding.

Conclusions

Our study indicates that 23−24-nt sRNAs may be linked to repressive chromatin modifications and function in genome stability while 20−21-nt sRNAs may be important for the cultivar specificity. This study provides a novel insight into the mechanism of sRNA expression and function in barley.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3023-5) contains supplementary material, which is available to authorized users.

Ecological risk assessment for DvSnf7 RNA: A plant-incorporated protectant with targeted activity against western corn rootworm

MON 87411 maize, which expresses DvSnf7 RNA, was developed to provide an additional mode of action to confer protection against corn rootworm (Diabrotica spp.). A critical step in the registration of a genetically engineered crop with an insecticidal trait is performing an ecological risk assessment to evaluate the potential for adverse ecological effects. For MON 87411, an assessment plan was developed that met specific protection goals by characterizing the routes and levels of exposure, and testing representative functional taxa that would be directly or indirectly exposed in the environment. The potential for toxicity of DvSnf7 RNA was evaluated with a harmonized battery of non-target organisms (NTOs) that included invertebrate predators, parasitoids, pollinators, soil biota as well as aquatic and terrestrial vertebrate species. Laboratory tests evaluated ecologically relevant endpoints such as survival, growth, development, and reproduction and were of sufficient duration to assess the potential for adverse effects. No adverse effects were observed with any species tested at, or above, the maximum expected environmental concentration (MEEC). All margins of exposure for NTOs were >10-fold the MEEC. Therefore, it is reasonable to conclude that exposure to DvSnf7 RNA, both directly and indirectly, is safe for NTOs at the expected field exposure levels.

Identification of microRNAs in rice root in response to nitrate and ammonium

Nitrate and ammonium are two major nitrogen (N) sources for higher plants, but they differ in utilization and signaling. MicroRNAs (miRNAs) play an essential role in N signal transduction; however, knowledge remains limited about the regulatory role of miRNAs responsive to different N sources, especially in crop plants. To get global overview on miRNAs involved in N response in rice, we performed high-throughput small RNA-sequencing under different nitrate and ammonium treatments. The results demonstrated that only 16 and 11 miRNAs were significantly induced by nitrate and ammonium under short-term treatment, respectively. However, 60 differentially expressed miRNAs were found between nitrate and ammonium under long-term cultivation. These results suggested that miRNA response greatly differentiates between nitrate and ammonium treatments. Furthermore, 44 miRNAs were found to be differentially expressed between high- and low-N conditions. Our study reveals comprehensive expression profiling of miRNAs responsive to different N sources and different N treatments, which advances our understanding on the regulation of different N signaling and homeostasis mediated by miRNAs.

High-resolution identification and abundance profiling of cassava (Manihot esculenta Crantz) microRNAs

BACKGROUND

Small RNAs (sRNAs) are endogenous sRNAs that play regulatory roles in plant growth, development, and biotic and abiotic stress responses. In plants, one subset of sRNAs, microRNAs (miRNAs) exhibit tissue-differential expression and regulate gene expression mainly through direct cleavage of mRNA or indirectly via production of secondary phased siRNAs (phasiRNAs) that silence cognate target transcripts in trans.

RESULTS

Here, we have identified cassava (Manihot esculenta Crantz) miRNAs using high resolution sequencing of sRNA libraries from leaf, stem, callus, male and female flower tissues. To analyze the data, we built a cassava genome database and, via sequence analysis and secondary structure prediction, 38 miRNAs not previously reported in cassava were identified. These new cassava miRNAs included two miRNAs not previously been reported in any plant species. The miRNAs exhibited tissue-differential accumulation as confirmed by quantitative RT-PCR and Northern blot analysis, largely reflecting levels observed in sequencing data. Some of the miRNAs identified were predicted to trigger production of secondary phased siRNAs (phasiRNAs) from 80 PHAS loci.

CONCLUSIONS

Cassava is a woody perennial shrub, grown principally for its starch-rich storage roots, which are rich in calories. In this study, new miRNAs were identified and their expression was validated using qRT-PCR of RNA from five different tissues. The data obtained expand the list of annotated miRNAs and provide additional new resources for cassava improvement research.

Genome-wide discovery and validation of Eucalyptus small RNAs reveals variable patterns of conservation and diversity across species of Myrtaceae

Background

Micro RNAs are a class of small non coding RNAs of 20–24 nucleotides transcribed as single stranded precursors from MIR gene loci. Initially described as post-transcriptional regulators involved in development, two decades ago, miRNAs have been proven to regulate a wide range of processes in plants such as germination, morphology and responses to biotic and abiotic stress. Despite wide conservation in plants, a number of miRNAs are lineage specific. We describe the first genome wide survey of Eucalyptus miRNAs based on high throughput sequencing.

Results

In addition to discovering small RNA sequences, MIR loci were mapped onto the reference genome and interspecific variability investigated. Sequencing was carried out for the two most world widely planted species, E. grandis and E. globulus. To maximize discovery, E. grandis samples were from BRASUZ1, the same tree whose genome provided the reference sequence. Interspecific analysis reinforces the variability in small RNA repertoire even between closely related species. Characterization of Eucalyptus small RNA sequences showed 95 orthologous to conserved miRNAs and 193 novel miRNAs. In silico target prediction confirmed 163 novel miRNAs and degradome sequencing experimentally confirmed several hundred targets. Experimental evidence based on the exclusive expression of a set of small RNAs across 16 species within Myrtaceae further highlighted variable patterns of conservation and diversity of these regulatory elements.

Conclusions

The description of miRNAs in Eucalyptus contributes to scientific knowledge of this vast genre, which is the most widely planted hardwood crop in the tropical and subtropical world, adding another important element to the annotation of Eucalyptus grandis reference genome.

Electronic supplementary material

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

Highly specific gene silencing in a monocot species by artificial microRNAs derived from chimeric miRNA precursors

Artificial microRNAs (amiRNAs) are used for selective gene silencing in plants. However, current methods to produce amiRNA constructs for silencing transcripts in monocot species are not suitable for simple, cost-effective and large-scale synthesis. Here, a series of expression vectors based on Oryza sativa MIR390 (OsMIR390) precursor was developed for high-throughput cloning and high expression of amiRNAs in monocots. Four different amiRNA sequences designed to target specifically endogenous genes and expressed from OsMIR390-based vectors were validated in transgenic Brachypodium distachyon plants. Surprisingly, amiRNAs accumulated to higher levels and were processed more accurately when expressed from chimeric OsMIR390-based precursors that include distal stem-loop sequences from Arabidopsis thaliana MIR390a (AtMIR390a). In all cases, transgenic plants displayed the predicted phenotypes induced by target gene repression, and accumulated high levels of amiRNAs and low levels of the corresponding target transcripts. Genome-wide transcriptome profiling combined with 5'-RLM-RACE analysis in transgenic plants confirmed that amiRNAs were highly specific.

Predicted trans-acting siRNAs in the human brain

Endogenous small non-coding RNAs play pivotal roles in regulating gene expression in eukaryotes. Many studies have investigated the function and molecular mechanism of microRNAs in the development and disease of various organisms via mRNA repression of protein-coding genes. Recent findings indicate microRNAs might trigger the generation of trans-acting small interfering RNAs (ta-siRNAs). The interaction among different types of small RNA molecules reveals an even more complicated and elaborate pattern of RNA regulation during gene expression than previously thought. We developed a method for mining ta-siRNA sequences and evaluated the performance of our novel method using data from Arabidopsis thaliana. Additionally, using small RNA and degradome data for the human brain, we identified 155 small RNAs that satisfied ta-siRNA characteristics. The DRAXIN and ATCAY genes, which are preferentially expressed in the human brain, were predicted to be the targets of 12 potential ta-siRNAs.

MicroRNA-based biotechnology for plant improvement

MicroRNAs (miRNAs) are an extensive class of newly discovered endogenous small RNAs, which negatively regulate gene expression at the post-transcription levels. As the application of next-generation deep sequencing and advanced bioinformatics, the miRNA-related study has been expended to non-model plant species and the number of identified miRNAs has dramatically increased in the past years. miRNAs play a critical role in almost all biological and metabolic processes, and provide a unique strategy for plant improvement. Here, we first briefly review the discovery, history, and biogenesis of miRNAs, then focus more on the application of miRNAs on plant breeding and the future directions. Increased plant biomass through controlling plant development and phase change has been one achievement for miRNA-based biotechnology; plant tolerance to abiotic and biotic stress was also significantly enhanced by regulating the expression of an individual miRNA. Both endogenous and artificial miRNAs may serve as important tools for plant improvement.

miR168 influences phase transition, leaf epinasty, and fruit development via SlAGO1s in tomato

In Arabidopsis thaliana, Argonaute1 (AGO1) interacts with miR168 to modulate the small RNA regulatory pathway. However, the underlying mechanism of regulation and relationship between AGO1 and miR168 is poorly understood in the cash crop Solanum lycopersicum (tomato). We previously found that SlAGO1A and SlAGO1B were cleaved by miR168 in tomato. In this study, we show that SlAGO1A and SlAGO1B accumulate in miR168-sponge transgenic plants, and that expression of miR168-resistant SlAGO1A (4m-SlAGO1A) and SlAGO1B (4m-SlAGO1B) in tomato results in a series of defects affecting growth rate, floral timing, leaves, and fruit. Accumulation of miR156 was found when 4m-SlAGO1A was at an early developmental stage compared to the wild type and original SlAGO1A transgenic plants, and miR172 was highly expressed in adult 4m-SlAGO1A compared to the controls. In addition, the expression of multiple small RNAs was altered in 4m-SlAGO1A. Taken together, our data provide novel insights into the interaction between SlAGO1s and miR168 in determining growth rate, phase change, leaf epinasty, fruit initiation and expansion, and other developmental processes in tomato.

A 28-day oral toxicity evaluation of small interfering RNAs and a long double-stranded RNA targeting vacuolar ATPase in mice

New biotechnology-derived crop traits have been developed utilizing the natural process of RNA interference (RNAi). However, plant-produced double stranded RNAs (dsRNAs) are not known to present a hazard to mammals because numerous biological barriers limit uptake and potential for activity. To evaluate this experimentally, dsRNA sequences matching the mouse vATPase gene (an established target for control of corn rootworms) were evaluated in a 28-day toxicity study with mice. Test groups were orally gavaged with escalating doses of either a pool of four 21-mer vATPase small interfering RNAs (siRNAs) or a 218-base pair vATPase dsRNA. There were no treatment-related effects on body weight, food consumption, clinical observations, clinical chemistry, hematology, gross pathology, or histopathology endpoints. The highest dose levels tested were considered to be the no observed adverse effect levels (NOAELs) for the 21-mer siRNAs (48 mg/kg/day) and the 218 bp dsRNA (64 mg/kg/day). As an additional exploratory endpoint, vATPase gene expression, was evaluated in selected gastrointestinal tract and systemic tissues. The results of this assay did not indicate treatment-related suppression of vATPase. The results of this study indicate that orally ingested dsRNAs, even those targeting a gene in the test species, do not produce adverse health effects in mammals.

Genomic methylation patterns in archaeological barley show de-methylation as a time-dependent diagenetic process

Genomic methylation is variable under biotic and abiotic stresses in plants. In particular, viral infection is thought to significantly increase genomic methylation with particularly high activity around transposable elements. Here we present the genomic methylation profiles of grains of archaeological barley (Hordeum vulgare) from several strata from a site in southern Egypt, from the Napatan to the Islamic periods (800 BCE - 1812 CE). One sample tested positive for viral infection and exhibits an unusually high degree of genomic methylation compared to the rest. A decreasing trend in global methylation levels according to deposition date shows in-situ de-methylation of 5-methylcytosine, which can be described as a diagenetic process. This is most likely a deamination mediated de-methylation process and is expected to lead to 5 mC > T base modifications in addition to the C > U modifications due to cytosine deamination, so represents a time-dependent process of DNA diagenesis in ancient DNA.

Whole-genome discovery of miRNAs and their targets in wheat (Triticum aestivum L.)

BACKGROUND

MicroRNAs (miRNAs) are small, non-coding RNAs playing essential roles in plant growth, development, and stress responses. Sequencing of small RNAs is a starting point for understanding their number, diversity, expression and possible roles in plants.

RESULTS

In this study, we conducted a genome-wide survey of wheat miRNAs from 11 tissues, characterizing a total of 323 novel miRNAs belonging to 276 families in wheat. A miRNA conservation analysis identified 191 wheat-specific miRNAs, 2 monocot-specific miRNAs, and 30 wheat-specific variants from 9 highly conserved miRNA families. To understand possible roles of wheat miRNAs, we determined 524 potential targets for 124 miRNA families through degradome sequencing, and cleavage of a subset of them was validated via 5' RACE. Based on the genome-wide identification and characterization of miRNAs and their associated target genes, we further identified 64 miRNAs preferentially expressing in developing or germinating grains, which could play important roles in grain development.

CONCLUSION

We discovered 323 wheat novel miRNAs and 524 target genes for 124 miRNA families in a genome-wide level, and our data will serve as a foundation for future research into the functional roles of miRNAs in wheat.

A diverse set of miRNAs responsive to begomovirus-associated betasatellite in Nicotiana benthamiana

BACKGROUND

Roles of microRNAs (miRNAs) and short interfering RNAs (siRNAs) in biotic stress responses, e.g., viral infection, have been demonstrated in plants by many studies. Tomato yellow leaf curl China virus (TYLCCNV) is a monopartite begomovirus that can systemically infect Solanaceae plants, and induces leaf curling, yellowing and enation symptoms when co-inoculated with a betasatellite (TYLCCNB). The released genome sequence of Nicotiana benthamiana provides an opportunity to identify miRNAs and siRNAs responsive to begomovirus-associated betasatellite in N. benthamiana.

RESULTS

miRNAs were identified in three small RNA libraries generated using RNA isolated from N. benthamiana plants systemically infected with TYLCCNV (Y10A) alone, co-infected with Y10A and its betasatellite TYLCCNB (Y10β) or a TYLCCNB mutant (Y10mβ) that contains a mutated βC1, the sole betasatellite-encoded protein. A total of 196 conserved miRNAs from 38 families and 197 novel miRNAs from 160 families were identified. Northern blot analysis confirmed that expression of species-specific miRNAs was much lower than that of conserved miRNAs. Several conserved and novel miRNAs were found to be responsive to co-infection of Y10A and Y10β but not to co-infection of Y10A and Y10mβ, suggesting that these miRNAs might play a role unique to interaction between Y10β and N. benthamiana. Additionally, we identified miRNAs that can trigger the production of phased secondary siRNAs (phasiRNAs).

CONCLUSIONS

Identification of miRNAs with differential expression profiles in N. benthamiana co-infected with Y10A and Y10β and co-infected with Y10A and Y10mβ indicates that these miRNAs are betasatellite-responsive. Our result also suggested a potential role of miRNA-mediated production of phasiRNAs in interaction between begomovirus and N. benthamiana.

Identification of novel and conserved miRNAs involved in pollen development in Brassica campestris ssp. chinensis by high-throughput sequencing and degradome analysis

Background

microRNAs (miRNAs) are endogenous, noncoding, small RNAs that have essential regulatory functions in plant growth, development, and stress response processes. However, limited information is available about their functions in sexual reproduction of flowering plants. Pollen development is an important process in the life cycle of a flowering plant and is a major factor that affects the yield and quality of crop seeds.

Results

This study aims to identify miRNAs involved in pollen development. Two independent small RNA libraries were constructed from the flower buds of the male sterile line (Bcajh97-01A) and male fertile line (Bcajh97-01B) of Brassica campestris ssp. chinensis. The libraries were subjected to high-throughput sequencing by using the Illumina Solexa system. Eight novel miRNAs on the other arm of known pre-miRNAs, 54 new conserved miRNAs, and 8 novel miRNA members were identified. Twenty-five pairs of novel miRNA/miRNA* were found. Among all the identified miRNAs, 18 differentially expressed miRNAs with over two-fold change between flower buds of male sterile line (Bcajh97-01A) and male fertile line (Bcajh97-01B) were identified. qRT-PCR analysis revealed that most of the differentially expressed miRNAs were preferentially expressed in flower buds of the male fertile line (Bcajh97-01B). Degradome analysis showed that a total of 15 genes were predicted to be the targets of seven miRNAs.

Conclusions

Our findings provide an overview of potential miRNAs involved in pollen development and interactions between miRNAs and their corresponding targets, which may provide important clues on the function of miRNAs in pollen development.

Identification of phasiRNAs in wild rice (Oryza rufipogon)

Plant miRNAs can trigger the production of phased, secondary siRNAs from either non-coding or protein-coding genes. In this study, at least 864 and 3,961 loci generating 21-nt and 24-nt phased siRNAs (phasiRNAs),respectively, were identified in three tissues from wild rice. Of these phasiRNA-producing loci, or PHAS genes, biogenesis of phasiRNAs in at least 160 of 21-nt and 254 of 24-nt loci could be triggered by interaction with miRNA(s). Developing seeds had more PHAS genes than leaves and roots. Genetic constrain on miRNA-triggered PHAS genes suggests that phasiRNAs might be one of the driving forces contributed to rice domestication.

Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks

Plant genomes are the source of large numbers of small RNAs, generated via a variety of genetically separable pathways. Several of these pathways converge in the production of phased, secondary, small interfering RNAs (phasiRNAs), originally designated as trans-acting small interfering RNAs or tasiRNAs. PhasiRNA biogenesis requires the involvement of microRNAs as well as the cellular machinery for the production of siRNAs. PhasiRNAs in Arabidopsis thaliana have been well described for their ability to function in trans to suppress target transcript levels. Plant genomic data from an expanding set of species have demonstrated that Arabidopsis is relatively sparing in its use of phasiRNAs, while other genomes contain hundreds or even thousands of phasiRNA-generating loci. In the dicots, targets of those phasiRNAs include several large or conserved families of genes, such as those encoding disease resistance proteins or transcription factors. Suppression of nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes by small RNAs is particularly unusual because of a high level of redundancy. In this review, we discuss plant phasiRNAs and the possible mechanistic significance of phasiRNA-based regulation of the NB-LRRs.

Polymerase IV occupancy at RNA-directed DNA methylation sites requires SHH1

DNA methylation is an epigenetic modification that has critical roles in gene silencing, development and genome integrity. In Arabidopsis, DNA methylation is established by DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2) and targeted by 24-nucleotide small interfering RNAs (siRNAs) through a pathway termed RNA-directed DNA methylation (RdDM). This pathway requires two plant-specific RNA polymerases: Pol-IV, which functions to initiate siRNA biogenesis, and Pol-V, which functions to generate scaffold transcripts that recruit downstream RdDM factors. To understand the mechanisms controlling Pol-IV targeting we investigated the function of SAWADEE HOMEODOMAIN HOMOLOG 1 (SHH1), a Pol-IV-interacting protein. Here we show that SHH1 acts upstream in the RdDM pathway to enable siRNA production from a large subset of the most active RdDM targets, and that SHH1 is required for Pol-IV occupancy at these same loci. We also show that the SHH1 SAWADEE domain is a novel chromatin-binding module that adopts a unique tandem Tudor-like fold and functions as a dual lysine reader, probing for both unmethylated K4 and methylated K9 modifications on the histone 3 (H3) tail. Finally, we show that key residues within both lysine-binding pockets of SHH1 are required in vivo to maintain siRNA and DNA methylation levels as well as Pol-IV occupancy at RdDM targets, demonstrating a central role for methylated H3K9 binding in SHH1 function and providing the first insights into the mechanism of Pol-IV targeting. Given the parallels between methylation systems in plants and mammals, a further understanding of this early targeting step may aid our ability to control the expression of endogenous and newly introduced genes, which has broad implications for agriculture and gene therapy.

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.

Identification and expression analysis of microRNAs at the grain filling stage in rice(Oryza sativa L.)via deep sequencing

MicroRNAs (miRNAs) have been shown to play crucial roles in the regulation of plant development. In this study, high-throughput RNA-sequencing technology was used to identify novel miRNAs, and to reveal miRNAs expression patterns at different developmental stages during rice (Oryza sativa L.) grain filling. A total of 434 known miRNAs (380, 402, 390 and 392 at 5, 7, 12 and 17 days after fertilization, respectively.) were obtained from rice grain. The expression profiles of these identified miRNAs were analyzed and the results showed that 161 known miRNAs were differentially expressed during grain development, a high proportion of which were up-regulated from 5 to 7 days after fertilization. In addition, sixty novel miRNAs were identified, and five of these were further validated experimentally. Additional analysis showed that the predicted targets of the differentially expressed miRNAs may participate in signal transduction, carbohydrate and nitrogen metabolism, the response to stimuli and epigenetic regulation. In this study, differences were revealed in the composition and expression profiles of miRNAs among individual developmental stages during the rice grain filling process, and miRNA editing events were also observed, analyzed and validated during this process. The results provide novel insight into the dynamic profiles of miRNAs in developing rice grain and contribute to the understanding of the regulatory roles of miRNAs in grain filling.

MicroRNAs in trees

MicroRNAs (miRNAs) are 20-24 nucleotide long molecules processed from a specific class of RNA polymerase II transcripts that mainly regulate the stability of mRNAs containing a complementary sequence by targeted degradation in plants. Many features of tree biology are regulated by miRNAs affecting development, metabolism, adaptation and evolution. MiRNAs may be modified and harnessed for controlled suppression of specific genes to learn about gene function, or for practical applications through genetic engineering. Modified (artificial) miRNAs act as dominant suppressors and are particularly useful in tree genetics because they bypass the generations of inbreeding needed for fixation of recessive mutations. The purpose of this review is to summarize the current status of information on miRNAs in trees and to guide future studies on the role of miRNAs in the biology of woody perennials and to illustrate their utility in directed genetic modification of trees.

Global analysis of RNA secondary structure in two metazoans

The secondary structure of RNA is necessary for its maturation, regulation, processing, and function. However, the global influence of RNA folding in eukaryotes is still unclear. Here, we use a high-throughput, sequencing-based, structure-mapping approach to identify the paired (double-stranded RNA [dsRNA]) and unpaired (single-stranded RNA [ssRNA]) components of the Drosophila melanogaster and Caenorhabditis elegans transcriptomes, which allows us to identify conserved features of RNA secondary structure in metazoans. From this analysis, we find that ssRNAs and dsRNAs are significantly correlated with specific epigenetic modifications. Additionally, we find key structural patterns across protein-coding transcripts that indicate that RNA folding demarcates regions of protein translation and likely affects microRNA-mediated regulation of mRNAs in animals. Finally, we identify and characterize 546 mRNAs whose folding pattern is significantly correlated between these metazoans, suggesting that their structure has some function. Overall, our findings provide a global assessment of RNA folding in animals.

Identification of novel maize miRNAs by measuring the precision of precursor processing

BACKGROUND

miRNAs are known to play important regulatory roles throughout plant development. Until recently, nearly all the miRNAs in maize were identified by comparative analysis to miRNAs sequences of other plant species, such as rice and Arabidopsis.

RESULTS

To find new miRNA in this important crop, small RNAs from mixed tissues were sequenced, resulting in over 15 million unique sequences. Our sequencing effort validated 23 of the 28 known maize miRNA families, including 49 unique miRNAs. Using a newly established criterion, based on the precision of miRNA processing from precursors, we identified 66 novel miRNAs in maize. These miRNAs can be grouped into 58 families, 54 of which have not been identified in any other species. Five new miRNAs were validated by northern blot. Moreover, we found targets for 23 of the 66 new miRNAs. The targets of two of these newly identified miRNAs were confirmed by 5'RACE.

CONCLUSION

We have implemented a novel method of identifying miRNA by measuring the precision of miRNA processing from precursors. Using this method, 66 novel miRNAs and 50 potential miRNAs have been identified in maize.

Small RNA diversity in plants and its impact in development

MicroRNAs are a class of non-coding RNAs involved in post-transcriptional control of gene expression, either via degradation or translational inhibition of target mRNAs. Both experimental and computational approaches have been used to identify miRNAs and their target genes. In plants, deep sequencing methods have recently allowed the analysis of small RNA diversity in different species and/or mutants. Most sequencing efforts have been concentrated on the identification of miRNAs and their mRNA targets have been predicted based on complementarity criteria. The recent demonstration that certain plant miRNAs could act partly via inhibition of protein translation certainly opens new fields of analysis for plant miRNA function on a broader group of targets. The roles of conserved miRNAs on target mRNA stability have been analysed in different species and defined common mechanisms in development and stress responses. In contrast, much less is known about expression patterns or functions of non-conserved miRNAs. In this review, we focus on the comparative analyses of plant small RNA diversity and the action of si/miRNAs in post-transcriptional regulation of some key genes involved in root development.

Characterization of the small RNA component of the transcriptome from grain and sweet sorghum stems

Background

Sorghum belongs to the tribe of the Andropogoneae that includes potential biofuel crops like switchgrass, Miscanthus and successful biofuel crops like corn and sugarcane. However, from a genomics point of view sorghum has compared to these other species a simpler genome because it lacks the additional rounds of whole genome duplication events. Therefore, it has become possible to generate a high-quality genome sequence. Furthermore, cultivars exists that rival sugarcane in levels of stem sugar so that a genetic approach can be used to investigate which genes are differentially expressed to achieve high levels of stem sugar.

Results

Here, we characterized the small RNA component of the transcriptome from grain and sweet sorghum stems, and from F2 plants derived from their cross that segregated for sugar content and flowering time. We found that variation in miR172 and miR395 expression correlated with flowering time whereas variation in miR169 expression correlated with sugar content in stems. Interestingly, genotypic differences in the ratio of miR395 to miR395* were identified, with miR395* species expressed as abundantly as miR395 in sweet sorghum but not in grain sorghum. Finally, we provided experimental evidence for previously annotated miRNAs detecting the expression of 25 miRNA families from the 27 known and discovered 9 new miRNAs candidates in the sorghum genome.

Conclusions

Sequencing the small RNA component of sorghum stem tissue provides us with experimental evidence for previously predicted microRNAs in the sorghum genome and microRNAs with a potential role in stem sugar accumulation and flowering time.

Evolution of MIR159/319 microRNA genes and their post-transcriptional regulatory link to siRNA pathways

Background

MicroRNAs (miRNAs) are prevalent and important endogenous gene regulators in eukaryotes. MiR159 and miR319 are highly conserved miRNAs essential for plant development and fertility. Despite high similarity in conservation pattern and mature miRNA sequences, miR159 and miR319 have distinct expression patterns, targets and functions. In addition, both MIR319 and MIR159 precursors produce multiple miRNAs in a phased loop-to-base manner. Thus, MIR159 and MIR319 appear to be related in origin and considerably diverged. However the phylogeny of MIR159 and MIR319 genes and why such unusual style of miRNA production has been conserved during evolution is not well understood.

Results

We reconstructed the phylogeny of MIR159/319 genes and analyzed their mature miRNA expression. The inferred phylogeny suggests that the MIR159/319 genes may have formed at least ten extant early-branching clades through gene duplication and loss. A series of duplications occurred in the common ancestor of seed plants leading to the original split of flowering plant MIR159 and MIR319. The results also indicate that the expression of MIR159/319 is regulated at post-transcriptional level to switch on the expression of alternative miRNAs during development in a highly spatio-temporal specific manner, and to selectively respond to the disruption of defensive siRNA pathways. Such intra-stem-loop regulation appears diverged across the early-branching clades of MIR159/319 genes.

Conclusions

Our results support that the MIR159 and MIR319 genes evolve from a common ancestor, which is likely to be a phased stem-loop small RNA. Through duplication and loss of genes this miRNA gene family formed clades specific to moss, lycopods, gymnosperms and angiosperms including the two major clades of flowering plants containing the founding members of MIR319 and MIR159 genes in A.thaliana. Our analyses also suggest that some MIR159/319 have evolved into unusual miRNA genes that are regulated at post-transcriptional level to express multiple mature products with variable proportions under different circumstances. Moreover, our analyses reveal conserved regulatory link of MIR159/319 genes to siRNA pathway through post-transcriptional regulation.

Domestication of transposable elements into MicroRNA genes in plants

Transposable elements (TE) usually take up a substantial portion of eukaryotic genome. Activities of TEs can cause genome instability or gene mutations that are harmful or even disastrous to the host. TEs also contribute to gene and genome evolution at many aspects. Part of miRNA genes in mammals have been found to derive from transposons while convincing evidences are absent for plants. We found that a considerable number of previously annotated plant miRNAs are identical or homologous to transposons (TE-MIR), which include a small number of bona fide miRNA genes that conform to generally accepted plant miRNA annotation rules, and hairpin derived siRNAs likely to be pre-evolved miRNAs. Analysis of these TE-MIRs indicate that transitions from the medium to high copy TEs into miRNA genes may undergo steps such as inverted repeat formation, sequence speciation and adaptation to miRNA biogenesis. We also identified initial target genes of the TE-MIRs, which contain homologous sequences in their CDS as consequence of cognate TE insertions. About one-third of the initial target mRNAs are supported by publicly available degradome sequencing data for TE-MIR sRNA induced cleavages. Targets of the TE-MIRs are biased to non-TE related genes indicating their penchant to acquire cellular functions during evolution. Interestingly, most of these TE insertions span boundaries between coding and non-coding sequences indicating their incorporation into CDS through alteration of splicing or translation start or stop signals. Taken together, our findings suggest that TEs in gene rich regions can form foldbacks in non-coding part of transcripts that may eventually evolve into miRNA genes or be integrated into protein coding sequences to form potential targets in a “temperate” manner. Thus, transposons may supply as resources for the evolution of miRNA-target interactions in plants.

Experimental design, preprocessing, normalization and differential expression analysis of small RNA sequencing experiments

Prior to the advent of new, deep sequencing methods, small RNA (sRNA) discovery was dependent on Sanger sequencing, which was time-consuming and limited knowledge to only the most abundant sRNA. The innovation of large-scale, next-generation sequencing has exponentially increased knowledge of the biology, diversity and abundance of sRNA populations. In this review, we discuss issues involved in the design of sRNA sequencing experiments, including choosing a sequencing platform, inherent biases that affect sRNA measurements and replication. We outline the steps involved in preprocessing sRNA sequencing data and review both the principles behind and the current options for normalization. Finally, we discuss differential expression analysis in the absence and presence of biological replicates. While our focus is on sRNA sequencing experiments, many of the principles discussed are applicable to the sequencing of other RNA populations.

Evolution and functional diversification of MIRNA genes

MicroRNAs (miRNAs) are small regulatory RNAs found in diverse eukaryotic lineages. In plants, a minority of annotated MIRNA gene families are conserved between plant families, while the majority are family- or species-specific, suggesting that most known MIRNA genes arose relatively recently in evolutionary time. Given the high proportion of young MIRNA genes in plant species, new MIRNA families are likely spawned and then lost frequently. Unlike highly conserved, ancient miRNAs, young miRNAs are often weakly expressed, processed imprecisely, lack targets, and display patterns of neutral variation, suggesting that young MIRNA loci tend to evolve neutrally. Genome-wide analyses from several plant species have revealed that variation in miRNA foldback expression, structure, processing efficiency, and miRNA size have resulted in the unique functionality of MIRNA loci and resulting miRNAs. Additionally, some miRNAs have evolved specific properties and functions that regulate other transcriptional or posttranscriptional silencing pathways. The evolution of miRNA processing and functional diversity underscores the dynamic nature of miRNA-based regulation in complex regulatory networks.

Complex regulation of two target genes encoding SPX-MFS proteins by rice miR827 in response to phosphate starvation

Here we report on the characterization of rice osa-miR827 and its two target genes, OsSPX-MFS1 and OsSPX-MFS2, which encode SPX-MFS proteins predicted to be implicated in phosphate (Pi) sensing or transport. We first show by Northern blot analysis that osa-miR827 is strongly induced by Pi starvation in both shoots and roots. Hybridization of osa-miR827 in situ confirms its strong induction by Pi starvation, with signals concentrated in mesophyll, epidermis and ground tissues of roots. In parallel, we analyzed the responses of the two OsSPX-MFS1 and OsSPX-MFS2 gene targets to Pi starvation. OsSPX-MFS1 mRNA is mainly expressed in shoots under sufficient Pi supply while its expression is reduced on Pi starvation, revealing a direct relationship between induction of osa-miR827 and down-regulation of OsSPX-MFS1. In contrast, OsSPX-MFS2 responds in a diametrically opposed manner to Pi starvation. The accumulation of OsSPX-MFS2 mRNA is dramatically enhanced under Pi starvation, suggesting the involvement of complex regulation of osa-miR827 and its two target genes. We further produced transgenic rice lines overexpressing osa-miR827 and T-DNA knockout mutant lines in which the expression of osa-miR827 is abolished. Compared with wild-type controls, both target mRNAs exhibit similar changes, their expression being reduced and increased in overexpressing and knockout lines, respectively. This suggests that OsSPX-MFS1 and OsSPX-MFS2 are both negatively regulated by osa-miR827 abundance although they respond differently to external Pi conditions. We propose that this is a complex mechanism comprising fine tuning of spatial or temporal regulation of both targets by osa-miR827.

Transcriptome-wide identification of microRNA targets in rice

MicroRNA (miRNA)-guided target RNA expression is vital for a wide variety of biological processes in eukaryotes. Currently, miRBase (version 13) lists 142 and 353 miRNAs from Arabidopsis and rice (Oryza sativa), respectively. The integration of miRNAs in diverse biological networks relies upon the confirmation of their RNA targets. In contrast with the well-characterized miRNA targets that are cleaved in Arabidopsis, only a few such targets have been confirmed in rice. To identify small RNA targets in rice, we applied the 'degradome sequencing' approach, which globally identifies the remnants of small RNA-directed target cleavage by sequencing the 5' ends of uncapped RNAs. One hundred and sixty targets of 53 miRNA families (24 conserved and 29 rice-specific) and five targets of TAS3-small interfering RNAs (siRNAs) were identified. Surprisingly, an additional conserved target for miR398, which has not been reported so far, has been validated. Besides conserved homologous transcripts, 23 non-conserved genes for nine conserved miRNAs and 56 genes for 29 rice-specific miRNAs were also identified as targets. Besides miRNA targets, the rice degradome contained fragments derived from MIRNA precursors. A closer inspection of these fragments revealed a unique pattern distinct from siRNA-producing loci. This attribute can serve as one of the ancillary criteria for separating miRNAs from siRNAs in plants.

Sequence variation and selection of small RNAs in domesticated rice

Background

Endogenous non-coding small RNAs (21-24 nt) play an important role in post-transcriptional gene regulation in plants. Domestication selection is the most important evolutionary force in shaping crop genomes. The extent of polymorphism at small RNA loci in domesticated rice and whether small RNA loci are targets of domestication selection have not yet been determined.

Results

A polymorphism survey of 94 small RNA loci (88 MIRNAs, four TAS3 loci and two miRNA-like long hairpins) was conducted in domesticated rice, generating 2 Mb of sequence data. Many mutations (substitution or insertion/deletion) were observed at small RNA loci in domesticated rice, e.g. 12 mutation sites were observed in the mature miRNA sequences of 11 MIRNAs (12.5% of the investigated MIRNAs). Several small RNA loci showed significant signals for positive selection and/or potential domestication selection.

Conclusions

Sequence variation at miRNAs and other small RNAs is higher than expected in domesticated rice. Like protein-coding genes, non-coding small RNA loci could be targets of domestication selection and play an important role in rice domestication and improvement.

DNA methylation mediated by a microRNA pathway

In plants, the known microRNAs (miRNAs) are produced as approximately 21 nucleotide (nt) duplexes from their precursors by Dicer-like 1 (DCL1). They are incorporated into Argonaute 1 (AGO1) protein to regulate target gene expression primarily through mRNA cleavage. We report here the discovery of a class of miRNAs in the model monocot rice (Oryza sativa). These are 24 nt in length and require another member of the Dicer family, DCL3, for their biogenesis. The 24 nt long miRNAs (lmiRNAs) are loaded into AGO4 clade proteins according to hierarchical rules, depending on the upstream biogenesis machinery and the 5'-terminal nucleotide. We demonstrated that lmiRNAs direct DNA methylation at loci from which they are produced as well as in trans at their target genes and play roles in gene regulation. Considered together, our findings define a miRNA pathway that mediates DNA methylation.

MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana

MicroRNAs (miRNAs) are short regulatory RNAs processed from partially self-complementary foldbacks within longer MIRNA primary transcripts. Several MIRNA families are conserved deeply through land plants, but many are present only in closely related species or are species specific. The finding of numerous evolutionarily young MIRNA, many with low expression and few if any targets, supports a rapid birth-death model for MIRNA evolution. A systematic analysis of MIRNA genes and families in the close relatives, Arabidopsis thaliana and Arabidopsis lyrata, was conducted using both whole-genome comparisons and high-throughput sequencing of small RNAs. Orthologs of 143 A. thaliana MIRNA genes were identified in A. lyrata, with nine having significant sequence or processing changes that likely alter function. In addition, at least 13% of MIRNA genes in each species are unique, despite their relatively recent speciation (approximately 10 million years ago). Alignment of MIRNA foldbacks to the Arabidopsis genomes revealed evidence for recent origins of 32 families by inverted or direct duplication of mostly protein-coding gene sequences, but less than half of these yield miRNA that are predicted to target transcripts from the originating gene family. miRNA nucleotide divergence between A. lyrata and A. thaliana orthologs was higher for young MIRNA genes, consistent with reduced purifying selection compared with deeply conserved MIRNA genes. Additionally, target sites of younger miRNA were lost more frequently than for deeply conserved families. In summary, our systematic analyses emphasize the dynamic nature of the MIRNA complement of plant genomes.