MicroRNA-directed regulation: to cleave or not to cleave

Comparative analysis of the carrot miRNAome in response to salt stress

Soil salinity adversely affects the yield and quality of crops, including carrot. During salt stress, plant growth and development are impaired by restricted water uptake and ion cytotoxicity, leading to nutrient imbalance and oxidative burst. However, the molecular mechanisms of the carrot plant response to salt stress remain unclear. The occurrence and expression of miRNAs that are potentially involved in the regulation of carrot tolerance to salinity stress were investigated. The results of small RNA sequencing revealed that salt-sensitive (DH1) and salt-tolerant (DLBA) carrot varieties had different miRNA expression profiles. A total of 95 miRNAs were identified, including 71 novel miRNAs, of which 30 and 23 were unique to DH1 and DLBA, respectively. The comparison of NGS and qPCR results allowed identification of two conserved and five novel miRNA involved in carrot response to salt stress, and which differentiated the salt-tolerant and salt-sensitive varieties. Degradome analysis supported by in silico-based predictions and followed by expression analysis of exemplary target genes pointed at genes related to proline, glutathione, and glutamate metabolism pathways as potential miRNA targets involved in salt tolerance, and indicated that the regulation of osmoprotection and antioxidant protection, earlier identified as being more efficient in the tolerant variety, may be controlled by miRNAs. Furthermore, potential miRNA target genes involved in chloroplast protection, signal transduction and the synthesis and modification of cell wall components were indicated in plants growing in saline soil.

MicroRNA2871b of Dongxiang Wild Rice (Oryza rufipogon Griff.) Negatively Regulates Cold and Salt Stress Tolerance in Transgenic Rice Plants

Cold and salt stresses are major environmental factors that constrain rice production. Understanding their mechanisms is important to enhance cold and salt stress tolerance in rice. MicroRNAs (miRNAs) are a class of non-coding RNAs with only 21-24 nucleotides that are gene regulators in plants and animals. Previously, miR2871b expression was suppressed by cold stress in Dongxiang wild rice (DXWR, Oryza rufipogon Griff.). However, its biological functions in abiotic stress responses remain elusive. In the present study, miR2871b of DWXR was overexpressed to investigate its function under stress conditions. When miR2871b of DWXR was introduced into rice plants, the transgenic lines were more sensitive to cold and salt stresses, and their tolerance to cold and salt stress decreased. The increased expression of miR2871b in rice plants also increased the levels of reactive oxygen species (ROS) and malondialdehyde (MDA); however, it markedly decreased the activities of peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) and the contents of proline (Pro) and soluble sugar (SS). These data suggested that miR2871b of DXWR has negative regulatory effects on cold and salt stress tolerance. Meanwhile, 412 differentially expressed genes (DEGs) were found in rice transgenic plants using transcriptome sequencing, among which 266 genes were up-regulated and 146 genes were down-regulated. Furthermore, the upstream cis-acting elements and downstream targets of miR2871b were predicted and analyzed, and several critical acting elements (ABRE and TC-rich repeats) and potential target genes (LOC_Os03g41200, LOC_Os07g47620, and LOC_Os04g30260) were obtained. Collectively, these results generated herein further elucidate the vital roles of miR2871b in regulating cold and salt responses of DXWR.

PeGRF6-PeGIF1 complex regulates cell proliferation in the leaf of Phalaenopsis equestris

Phalaenopsis equestris is an ornamental plant with very large leaves. In this study, we identified genes related to the regulation of leaf development in Phalaenopsis and explored their mechanism of action. Sequence alignment and phylogenetic analyses revealed that PeGRF6 in the PeGRF family of P. equestris has similarities with the Arabidopsis genes AtGRF1 and AtGRF2, which are known to be involved in the regulation of leaf development. Among the PeGRFs, PeGRF6 was continuously and stably expressed at various stages of leaf development. The functions of PeGRF6 and of its complex formed with PeGIF1 in leaf development were verified by virus-induced gene silencing (VIGS) technology. The results show that the PeGRF6-PeGIF1 complex forms in the nucleus and positively regulates leaf cell proliferation via influencing cell size. Interestingly, VIGS suppression of PeGRF6 resulted in anthocyanin accumulation in Phalaenopsis leaves. Analyses of the regulatory mechanism of the miR396-PeGRF6 model based on the P. equestris small RNA library constructed here suggested that PeGRF6 transcripts are cleaved by Peq-miR396. These results show that, compared with PeGRF6 or PeGIF1 alone, the PeGRF6-PeGIF1 complex plays a more important role in the leaf development of Phalaenopsis, possibly by regulating the expression of cell cycle-related genes.

Expression Profile of Selected Genes Involved in Na+ Homeostasis and In Silico miRNA Identification in Medicago sativa and Medicago arborea under Salinity Stress

The accumulation of ions due to increased salinity in the soil is one of the major abiotic stressors of cultivated plants that negatively affect their productivity. The model plant, Medicago truncatula, is the only Medicago species that has been extensively studied, whereas research into increased salinity adaptation of two important forage legumes, M. sativa and M. arborea, has been limited. In the present study, the expression of six genes, namely SOS1, SOS3, NHX2, AKT, AVP and HKT1 was monitored to investigate the manner in which sodium ions are blocked and transferred to the various plant parts. In addition, in silico miRNA analysis was performed to identify miRNAs that possibly control the expression of the genes studied. The following treatments were applied: (1) salt stress, with initial treatment of 50 mM NaCl and gradual acclimatization every 10 days, (2) salt shock, with continuous application of 100 mM NaCl concentration and (3) no application of NaCl. Results showed that M. arborea appeared to overexpress and activate all available mechanisms of resistance in conditions of increased salinity, while M. sativa acted in a more targeted way, overexpressing the HKT1 and AKT genes that contribute to the accumulation of sodium ions, particularly in the root. Regarding miRNA in silico analysis, five miRNAs with significant complementarity to putative target genes, AKT1, AVP and SOS3 were identified and served as a first step in investigating miRNA regulatory networks. Further miRNA expression studies will validate these results. Our findings contribute to the understanding of the molecular mechanisms underlying salt-responsiveness in Medicago and could be used in the future for generating salt-tolerant genotypes in crop improvement programs.

SimiR396d targets SiGRF1 to regulate drought tolerance and root growth in foxtail millet

MicroRNAs play critical roles in growth, development and abiotic stress responses. SimR396d is a miRNA whose expression level is much higher in foxtail millet roots than other tissues. Whether SimR396d is involved in foxtail millet root growth and response to abiotic stress is still unknown. Here, we demonstrate that SimiR396d modulates both drought response and root growth in foxtail millet. The expression of SimiR396d is induced by PEG treatment. Overexpression of SimiR396d enhances drought tolerance and root length, while knockdown SimiR396d expression using target mimics of SimiR396d (MIM396) resulted in reduced drought tolerance and shortened root length. Furthermore, we identified and confirmed a plant-specific transcription factor, growth-regulating factor 1 (SiGRF1), as a direct target of SimiR396d. Overexpression of SiGRF1 in foxtail millet resulted in suppressed root growth and reduced sensitivity to drought stress. Moreover, ethylene signaling is necessary for SimiR396d and SiGRF1 to participate in the regulation of plant root growth. These results revealed a pivotal role of SimiR396d in drought tolerance and root growth in foxtail millet. SimiR396d-SiGRF1 regulatory module provides a strategy to improve drought-stress resistance of crop.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Identification of GROWTH-REGULATING FACTOR transcription factors in lettuce (Lactuca sativa) genome and functional analysis of LsaGRF5 in leaf size regulation

BACKGROUND

GROWTH-REGULATING FACTORs (GRFs), a type of plant-specific transcription factors, play important roles in regulating plant growth and development. Although GRF gene family has been identified in various plant species, a genome-wide analysis of this family in lettuce (Lactuca sativa L.) has not been reported yet.

RESULTS

Here we identified 15 GRF genes in lettuce and performed comprehensive analysis of them, including chromosomal locations, gene structures, and conserved motifs. Through phylogenic analysis, we divided LsaGRFs into six groups. Transactivation assays and subcellular localization of LsaGRF5 showed that this protein is likely to act as a transcriptional factor in the cell nucleus. Furthermore, transgenic lettuce lines overexpressing LsaGRF5 exhibited larger leaves, while smaller leaves were observed in LsaMIR396a overexpression lines, in which LsaGRF5 was down-regulated.

CONCLUSIONS

These results in lettuce provide insight into the molecular mechanism of GRF gene family in regulating leaf growth and development and foundational information for genetic improvement of the lettuce variations specialized in leaf character.

Oomycete small RNAs bind to the plant RNA-induced silencing complex for virulence

The exchange of small RNAs (sRNAs) between hosts and pathogens can lead to gene silencing in the recipient organism, a mechanism termed cross-kingdom RNAi (ck-RNAi). While fungal sRNAs promoting virulence are established, the significance of ck-RNAi in distinct plant pathogens is not clear. Here, we describe that sRNAs of the pathogen Hyaloperonospora arabidopsidis, which represents the kingdom of oomycetes and is phylogenetically distant from fungi, employ the host plant’s Argonaute (AGO)/RNA-induced silencing complex for virulence. To demonstrate H. arabidopsidis sRNA (HpasRNA) functionality in ck-RNAi, we designed a novel CRISPR endoribonuclease Csy4/GUS reporter that enabled in situ visualization of HpasRNA-induced target suppression in Arabidopsis. The significant role of HpasRNAs together with AtAGO1 in virulence was revealed in plant atago1 mutants and by transgenic Arabidopsis expressing a short-tandem-target-mimic to block HpasRNAs, that both exhibited enhanced resistance. HpasRNA-targeted plant genes contributed to host immunity, as Arabidopsis gene knockout mutants displayed quantitatively enhanced susceptibility.

Conserved miR396b-GRF Regulation Is Involved in Abiotic Stress Responses in Pitaya (Hylocereus polyrhizus)

MicroRNA396 (miR396) is a conserved microRNA family that targets growth-regulating factors (GRFs), which play significant roles in plant growth and stress responses. Available evidence justifies the idea that miR396-targeted GRFs have important functions in many plant species; however, no genome-wide analysis of the pitaya (Hylocereus polyrhizus) miR396 gene has yet been reported. Further, its biological functions remain elusive. To uncover the regulatory roles of miR396 and its targets, the hairpin sequence of pitaya miR396b and the open reading frame (ORF) of its target, HpGRF6, were isolated from pitaya. Phylogenetic analysis showed that the precursor miR396b (MIR396b) gene of plants might be clustered into three major groups, and, generally, a more recent evolutionary relationship in the intra-family has been demonstrated. The sequence analysis indicated that the binding site of hpo-miR396b in HpGRF6 is located at the conserved motif which codes the conserved "RSRKPVE" amino acid in the Trp-Arg-Cys (WRC) region. In addition, degradome sequencing analysis confirmed that four GRFs (GRF1, c56908.graph_c0; GRF4, c52862.graph_c0; GRF6, c39378.graph_c0 and GRF9, c54658.graph_c0) are hpo-miR396b targets that are regulated by specific cleavage at the binding site between the 10th and 11th nucleotides from the 5' terminus of hpo-miR396b. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that hpo-miR396b is down-regulated when confronted with drought stress (15% polyethylene glycol, PEG), and its expression fluctuates under other abiotic stresses, i.e., low temperature (4 ± 1 °C), high temperature (42 ± 1 °C), NaCl (100 mM), and abscisic acid (ABA; 0.38 mM). Conversely, the expression of HpGRF6 showed the opposite trend to exposure to these abiotic stresses. Taken together, hpo-miR396b plays a regulatory role in the control of HpGRF6, which might influence the abiotic stress response of pitaya. This is the first documentation of this role in pitaya and improves the understanding of the molecular mechanisms underlying the tolerance to drought stress in this fruit.

DRB1, DRB2 and DRB4 Are Required for Appropriate Regulation of the microRNA399/PHOSPHATE2 Expression Module in Arabidopsis thaliana

Adequate phosphorous (P) is essential to plant cells to ensure normal plant growth and development. Therefore, plants employ elegant mechanisms to regulate P abundance across their developmentally distinct tissues. One such mechanism is PHOSPHATE2 (PHO2)-directed ubiquitin-mediated degradation of a cohort of phosphate (PO4) transporters. PHO2 is itself under tight regulation by the PO4 responsive microRNA (miRNA), miR399. The DOUBLE-STRANDED RNA BINDING (DRB) proteins, DRB1, DRB2 and DRB4, have each been assigned a specific functional role in the Arabidopsis thaliana (Arabidopsis) miRNA pathway. Here, we assessed the requirement of DRB1, DRB2 and DRB4 to regulate the miR399/PHO2 expression module under PO4 starvations conditions. Via the phenotypic and molecular assessment of the knockout mutant plant lines, drb1, drb2 and drb4, we show here that; (1) DRB1 and DRB2 are required to maintain P homeostasis in Arabidopsis shoot and root tissues; (2) DRB1 is the primary DRB required for miR399 production; (3) DRB2 and DRB4 play secondary roles in regulating miR399 production, and; (4) miR399 appears to direct expression regulation of the PHO2 transcript via both an mRNA cleavage and translational repression mode of RNA silencing. Together, the hierarchical contribution of DRB1, DRB2 and DRB4 demonstrated here to be required for the appropriate regulation of the miR399/PHO2 expression module identifies the extreme importance of P homeostasis maintenance in Arabidopsis to ensure that numerous vital cellular processes are maintained across Arabidopsis tissues under a changing cellular environment.

Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA-target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.

MicroRNA‑758 inhibits tumorous behavior in tongue squamous cell carcinoma by directly targeting metadherin

Numerous microRNAs (miRNAs) are dysregulated in tongue squamous cell carcinoma (TSCC), and their dysregulation has been demonstrated to have a strong correlation with TSCC progression via regulation of their targets. Therefore, miRNAs have potential use in the diagnosis and treatment of patients with TSCC. In the present study, miRNA‑758 (miR‑758) expression in TSCC tissues and cell lines was detected through reverse transcription‑quantitative polymerase chain reaction, and the effects of miR‑758 on TSCC cell proliferation and invasion were investigated by using Cell Counting kit‑8 and Transwell invasion assays. A luciferase reporter assay was performed to determine the target interaction between miR‑758 and metadherin (MTDH) in TSCC cells. The results revealed that miR‑758 was downregulated in TSCC tissues and cell lines. miR‑758 overexpression restricted the proliferation and invasion of TSCC cells. Additionally, MTDH was verified as a direct target gene of miR‑758 in TSCC cells. Furthermore, MTDH was observed to be upregulated in TSCC tissues, and the upregulation of MTDH was inversely correlated with miR‑758 expression. Moreover, restored MTDH expression significantly counteracted the suppressive effects of miR‑758 overexpression on TSCC cells. These results suggested that miR‑758 may prevent TSCC progression and development by directly targeting MTDH, thereby providing evidence that miR‑758 is a novel therapeutic target for the treatment of patients with TSCC.

Characterization of Vv-miR156: Vv-SPL pairs involved in the modulation of grape berry development and ripening

SPL is a plant-specific transcription factor family. Many researchers reported that SPL members targeted by miR156s could play crucial roles in the modulation of plant growth and development. Although there are similar reports on grapes, till now little is known about grape berry development and ripening. To gain more insight into how grape miR156s (Vv-miR156s) modulated the above given processes of grape berries by mediating their target gene Vv-SPLs, here we identified the precise sequences of Vv-miR156s in 'Giant Rose' grape berries, predicted their potential targets, and revealed that the matching degree of various Vv-miR156: Vv-SPL pairs exhibited some discrepancy, implying the divergence of their interaction. Subsequently, we also discovered similar motifs such as ABRE, CGTCA and ERE, which are more specific to berry development and ripening, within the promoters of both Vv-MIR156s and Vv-SPLs. With berry development and ripening, meanwhile, Vv-miR156a, b/c/d, e and f/g/i exhibited an overall increasing expression trend, while their targets showed opposite trends at the corresponding stages. Additionally, exogenous ABA and NAA application promoted or curbed the expression of Vv-miR156s to some extent, before grape berry ripening stage. The cleavage products, sites and frequencies of Vv-miR156a, b/c/d, e, f/g/i and their respective targets (Vv-SPL2, 9, 10, 16) during grape berry development and ripening process were validated by our developed PPM-RACE and modified RLM-RACE together with qRT-PCR, which demonstrated that Vv-miR156s can be involved in the modulation of grape berry development and ripening process by mediating the expression of Vv-SPL2, 9, 10, 16. Our findings lay an important foundation for further recognizing their functions in grape berries, and enrich the knowledge of the regulatory mechanism of miRNA-mediated grape berry development and ripening.

microRNA-338-3p inhibits proliferation, migration, invasion, and EMT in osteosarcoma cells by targeting activator of 90 kDa heat shock protein ATPase homolog 1

Background

Osteosarcoma (OS) is a rare, malignant bone tumor that primarily affects adolescents and has a high degree of malignancy and high incidence of recurrence and metastasis. Our study aimed to explore the role of miR-338-3p in OS cells.

Methods

qRT-qPCR was performed to quantify miR-338-3p expression levels in OS tissue samples and in three common OS cell lines. MG-63 and Saos2 cells were separately transfected with miR-338-3p or NC mimics and miR-338-3p expression levels was determined by qRT-PCR. Cell proliferation was monitored using the Cell Counting Kit-8. Flow cytometer analysis was carried out to determine the distribution of cell cycle stages and apoptosis. Transwell assay was performed to measure the migratory and invasive capacities of MG-63 and Saos2 cells. The expression of Vimentin and E-cadherin was detected by western blot. Luciferase reporter assay, qRT-PCR and western blotting were performed to confirm the target of miR-338-3p.

Results

Analysis by qRT-PCR revealed that miR-338-3p was downregulated in the tissue samples of 20 OS patients when compared with that in their matched adjacent non-tumor tissues. Furthermore, miR-338-3p was significantly downregulated in three common OS cell lines, namely, MG-63, Saos2, and HOS, when compared with that in the human osteoblast cell line hFOB1.19. Analysis by luciferase reporter assay, qRT-PCR, and western blotting revealed that activator of 90 kDa heat shock protein ATPase homolog 1 (AHSA1) is a direct target of miR-338-3p. miR-338-3p overexpression led to significant reduction in AHSA1 protein levels in MG63 and Saos2 cells. miR-338-3p overexpression reduced cell viability and migration and invasion behavior of MG63 and Saos2 cells. In addition, miR-338-3p overexpression suppressed epithelial–mesenchymal transition (EMT), induced a significant G1-phase arrest and did not affect the apoptosis in both MG-63 and Saos2 cells. Moreover, overexpression of AHSA1 reversed the inhibitory effect of miR-338-3p overexpression on proliferation, cell cycle, apoptosis, EMT, migration, and invasion of MG63 and Saos2 cells, thereby suggesting that miR-338-3p acts as a tumor suppressor in OS cells by targeting AHSA1.

Conclusions

miR-338-3p/AHSA1 can serve as a potential therapeutic target for OS therapy.

Peroxisomal Targeting as a Sensitive Tool to Detect Protein-Small RNA Interactions through in Vivo Piggybacking

Peroxisomes are organelles that play key roles in eukaryotic metabolism. Their protein complement is entirely imported from the cytoplasm thanks to a unique pathway that is able to translocate folded proteins and protein complexes across the peroxisomal membrane. The import of molecules bound to a protein targeted to peroxisomes is an active process known as 'piggybacking' and we have recently shown that P15, a virus-encoded protein possessing a peroxisomal targeting sequence, is able to piggyback siRNAs into peroxisomes. Here, we extend this observation by analyzing the small RNA repertoire found in peroxisomes of P15-expressing plants. A direct comparison with the P15-associated small RNA retrieved during immunoprecipitation (IP) experiments, revealed that in vivo piggybacking coupled to peroxisome isolation could be a more sensitive means to determine the various small RNA species bound by a given protein. This increased sensitivity of peroxisome isolation as opposed to IP experiments was also striking when we analyzed the small RNA population bound by the Tomato bushy stunt virus-encoded P19, one of the best characterized viral suppressors of RNA silencing (VSR), artificially targeted to peroxisomes. These results support that peroxisomal targeting should be considered as a novel/alternative experimental approach to assess in vivo interactions that allows detection of labile binding events. The advantages and limitations of this approach are discussed.

Peroxisomal Targeting as a Sensitive Tool to Detect Protein-Small RNA Interactions through in Vivo Piggybacking

Peroxisomes are organelles that play key roles in eukaryotic metabolism. Their protein complement is entirely imported from the cytoplasm thanks to a unique pathway that is able to translocate folded proteins and protein complexes across the peroxisomal membrane. The import of molecules bound to a protein targeted to peroxisomes is an active process known as 'piggybacking' and we have recently shown that P15, a virus-encoded protein possessing a peroxisomal targeting sequence, is able to piggyback siRNAs into peroxisomes. Here, we extend this observation by analyzing the small RNA repertoire found in peroxisomes of P15-expressing plants. A direct comparison with the P15-associated small RNA retrieved during immunoprecipitation (IP) experiments, revealed that in vivo piggybacking coupled to peroxisome isolation could be a more sensitive means to determine the various small RNA species bound by a given protein. This increased sensitivity of peroxisome isolation as opposed to IP experiments was also striking when we analyzed the small RNA population bound by the Tomato bushy stunt virus-encoded P19, one of the best characterized viral suppressors of RNA silencing (VSR), artificially targeted to peroxisomes. These results support that peroxisomal targeting should be considered as a novel/alternative experimental approach to assess in vivo interactions that allows detection of labile binding events. The advantages and limitations of this approach are discussed.

MicroRNA-like RNAs from the same miRNA precursors play a role in cassava chilling responses

MicroRNAs (miRNAs) are known to play important roles in various cellular processes and stress responses. MiRNAs can be identified by analyzing reads from high-throughput deep sequencing. The reads realigned to miRNA precursors besides canonical miRNAs were initially considered as sequencing noise and ignored from further analysis. Here we reported a small-RNA species of phased and half-phased miRNA-like RNAs different from canonical miRNAs from cassava miRNA precursors detected under four distinct chilling conditions. They can form abundant multiple small RNAs arranged along precursors in a tandem and phased or half-phased fashion. Some of these miRNA-like RNAs were experimentally confirmed by re-amplification and re-sequencing, and have a similar qRT-PCR detection ratio as their cognate canonical miRNAs. The target genes of those phased and half-phased miRNA-like RNAs function in process of cell growth metabolism and play roles in protein kinase. Half-phased miR171d.3 was confirmed to have cleavage activities on its target gene P-glycoprotein 11, a broad substrate efflux pump across cellular membranes, which is thought to provide protection for tropical cassava during sharp temperature decease. Our results showed that the RNAs from miRNA precursors are miRNA-like small RNAs that are viable negative gene regulators and may have potential functions in cassava chilling responses.

Novel Insights into miRNA Regulation of Storage Protein Biosynthesis during Wheat Caryopsis Development under Drought Stress

Drought stress is a significant abiotic stress factor that affects wheat yield and quality. MicroRNA (miRNA) plays an important role in regulating caryopsis development in response to drought stress. However, little is known about the expression characteristics of miRNAs and how they regulate protein accumulation in wheat caryopsis under drought stress. To address this, two small RNA libraries of wheat caryopsis under control and drought stress conditions were constructed and sequenced. A total of 125 miRNAs were identified in the two samples, of which 110 were known and 15 were novel. A total of 1,981 miRNA target genes were predicted and functional annotations were obtained from various databases for 1,641 of them. Four miRNAs were identified as differential expression under drought stress, and the expression patterns of three of them were consistent with results obtained by reverse transcription polymerase chain reaction (RT-PCR) and reverse transcription quantitative polymerase chain reaction (RT-qPCR). Moreover, three miRNA-target pairs showed negative regulation tendency, as revealed by RT-qPCR. Functional enrichment and pathway analysis revealed that four pathways might be involved in storage protein biosynthesis. Furthermore, drought stress significantly increased the accumulation of protein bodies and protein content in wheat endosperm. In summary, our findings suggest that drought stress may enhance storage protein by regulating the expression of miRNAs and their target genes.

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.

MicroRNAs and Their Regulatory Role in Sugarcane

Sugarcane, one of the most photosynthetically efficient crops, is an important source of sugar and feedstock for green energy and co-generation. The high level of polyploidy and genomic peculiarities in this crop point towards a complex mechanism of regulation for the economically important traits like sugar content, cane yield related traits, resistance to biotic and abiotic stresses etc. The regulatory pathways for these traits comprise of a number of genes, transcription factors and different categories of RNAs like small interference RNAs (siRNAs), and Micro RNAs (miRNAs). MicroRNAs (miRNAs) are found to play an important regulatory role in many crops. As in other crops, several miRNAs have been identified in sugarcane too and these are speculated to have a role in regulating the various metabolic processes. Role of miRNAs in relation to drought tolerance has been studied to a great extent in this crop. miRNAs have been predicted to be linked to expression of other traits like disease resistance, salinity tolerance, waterlogging and axillary bud growth in sugarcane. miRNAs can have a significant role in biomass production in sugarcane, as reported in several biofuel crops. Till now, miRNAs linked to sugar accumulation have not been identified in sugarcane, but studies suggest an important role for miRNAs in sugar metabolic pathway in crops like Sorghum and switch grass. It is presumed that in sugarcane too, sugar accumulation as well as the other important metabolic pathways might be regulated to some extent by the miRNAs. The review examines the progress made in understanding the miRNA regulation in sugarcane and the extent to which miRNA mediated regulation can be utilized in sugarcane improvement.

Identification of Drought-Responsive MicroRNAs from Roots and Leaves of Alfalfa by High-Throughput Sequencing

Alfalfa, an important forage legume, is an ideal crop for sustainable agriculture and a potential crop for bioenergy resources. Drought, one of the most common environmental stresses, substantially affects plant growth, development, and productivity. MicroRNAs (miRNAs) are newly discovered gene expression regulators that have been linked to several plant stress responses. To elucidate the role of miRNAs in drought stress regulation of alfalfa, a high-throughput sequencing approach was used to analyze 12 small RNA libraries comprising of four samples, each with three biological replicates. From the 12 libraries, we identified 348 known miRNAs belonging to 80 miRNA families, and 281 novel miRNAs, using Mireap software. Eighteen known miRNAs in roots and 12 known miRNAs in leaves were screened as drought-responsive miRNAs. With the exception of miR319d and miR157a which were upregulated under drought stress, the expression pattern of drought-responsive miRNAs was different between roots and leaves in alfalfa. This is the first study that has identified miR3512, miR3630, miR5213, miR5294, miR5368 and miR6173 as drought-responsive miRNAs. Target transcripts of drought-responsive miRNAs were computationally predicted. All 447 target genes for the known miRNAs were predicted using an online tool. This study provides a significant insight on understanding drought-responsive mechanisms of alfalfa.

Characterization of miR061 and its target genes in grapevine responding to exogenous gibberellic acid

MicroRNAs (miRNAs), as an important growth regulator, are also involved in gibberellic acid (GA) signaling, revealing much relationship between miRNAs and GA in various plant responses. Grape is highly sensitive to GA₃, which plays a significant regulatory role in regulation of flower development, berry expansion, berry set, berry ripening, and seedlessness induction; further, it was found that grapevine miR061 (VvmiR061) is a GA₃ responsive miRNA. In this study, grapevine REV (VvREV) and HOX32 (VvHOX32), two target genes of VvmiR061, were predicted, verified, and cloned; homologous conservation was analyzed in various plants. The expression profiles of both VvmiR061 and its target genes (VvREV and VvHOX32) under GA₃ treatment were detected by qRT-PCR during grapevine flower and berry development. Results revealed that GA₃ treatment has upregulated the transcription of VvREV and VvHOX32, while it downregulated the expression of VvmiR061. The function of VvmiR061 in cleaving target genes VvREV and VvHOX32 was diminished by GA₃ treatment during flower developmental process. The results of this study exhibited the importance of VvmiR061 in regulating flower development and GA₃ signaling pathway and also contributed some to the knowledge of small RNA-mediated regulation in grape.

Identification of microRNAs in Response to Drought in Common Wild Rice (Oryza rufipogon Griff.) Shoots and Roots

Background

Drought is the most important factor that limits rice production in drought-prone environments. Plant microRNAs (miRNAs) are involved in biotic and abiotic stress responses. Common wild rice (Oryza rufipogon Griff.) contains abundant drought-resistant genes, which provide an opportunity to explore these excellent resources as contributors to improve rice resistance, productivity, and quality.

Results

In this study, we constructed four small RNA libraries, called CL and CR from PEG6000-free samples and DL and DR from PEG6000-treated samples, where ‘R’ indicates the root tissue and ‘L’ indicates the shoot tissue. A total of 200 miRNAs were identified to be differentially expressed under the drought-treated conditions (16% PEG6000 for 24 h), and the changes in the miRNA expression profile of the shoot were distinct from those of the root. At the miRNA level, 77 known miRNAs, which belong to 23 families, including 40 up-regulated and 37 down-regulated in the shoot, and 85 known miRNAs in 46 families, including 65 up-regulated and 20 down-regulated in the root, were identified as differentially expressed. In addition, we predicted 26 new miRNA candidates from the shoot and 43 from the root that were differentially expressed during the drought stress. The quantitative real-time PCR analysis results were consistent with high-throughput sequencing data. Moreover, 88 miRNAs that were differentially-expressed were predicted to match with 197 targets for drought-stress.

Conclusion

Our results suggest that the miRNAs of O. rufipogon are responsive to drought stress. The differentially expressed miRNAs that are tissue-specific under drought conditions could play different roles in the regulation of the auxin pathway, the flowering pathway, the drought pathway, and lateral root formation. Thus, the present study provides an account of tissue-specific miRNAs that are involved in the drought adaption of O. rufipogon.

Identification of Salt Tolerance-related microRNAs and Their Targets in Maize (Zea mays L.) Using High-throughput Sequencing and Degradome Analysis

To identify the known and novel microRNAs (miRNAs) and their targets that are involved in the response and adaptation of maize (Zea mays) to salt stress, miRNAs and their targets were identified by a combined analysis of the deep sequencing of small RNAs (sRNA) and degradome libraries. The identities were confirmed by a quantitative expression analysis with over 100 million raw reads of sRNA and degradome sequences. A total of 1040 previously known miRNAs were identified from four maize libraries, with 762 and 726 miRNAs derived from leaves and roots, respectively, and 448 miRNAs that were common between the leaves and roots. A total of 37 potential new miRNAs were selected based on the same criteria in response to salt stress. In addition to known miR167 and miR164 species, novel putative miR167 and miR164 species were also identified. Deep sequencing of miRNAs and the degradome [with quantitative reverse transcription polymerase chain reaction (qRT-PCR) analyses of their targets] showed that more than one species of novel miRNA may play key roles in the response to salinity in maize. Furthermore, the interaction between miRNAs and their targets may play various roles in different parts of maize in response to salinity.

Analysis of high iron rice lines reveals new miRNAs that target iron transporters in roots

The present study highlights the molecular regulation of iron transport in soyFER1-overexpressing transgenic rice. Accumulation of iron in three different seed developmental stages, milk, dough, and mature, has been examined. The transgenic seeds of the milk stage showed significant augmentation of iron and zinc levels compared with wild-type seeds, and similar results were observed throughout the dough and mature stages. To investigate the regulation of iron transport, the role of miRNAs was studied in roots of transgenic rice. Sequencing of small RNA libraries revealed 153 known and 41 novel miRNAs in roots. Among them, 59 known and 14 novel miRNAs were found to be significantly expressed. miR166, miR399, and miR408 were identified as playing a vital role in iron uptake in roots of transgenic plants . Most importantly, four putative novel miRNAs, namely miR11, miR26, miR30, and miR31, were found to be down-regulated in roots of transgenic plants. For all these four novel miRNAs, natural resistance-associated macrophage protein 4 (NRAMP4), encoding a metal transporter, was predicted as a target gene. It is hypothesized that the NRAMP4 transporter is activated in roots of transgenic plants due to the lower abundance of its corresponding putative novel miRNAs. The relative transcript level of the NRAMP4 transcript was increased from 0.107 in the wild type to 65.24 and 55.39 in transgenic plants, which demonstrates the elevated amount of iron transport in transgenic plants. In addition, up-regulation of OsYSL15, OsFRO2, and OsIRT1 in roots also facilitates iron loading in transgenic seeds.

Study on Expression Modes and Cleavage Role of miR156b/c/d and its Target Gene Vv-SPL9 During the Whole Growth Stage of Grapevine

miR156 regulates the expression of its target SPL (PROMOTER BINDING-LIKE) genes during flower and fruit development, diverse developmental stage transitions, especially from vegetative to reproductive growth phases, by cleaving the target mRNA SPL of one plant-specific transcription factor. However, systematic reports on grapevine have yet to be presented. Here, the precise sequence of miR156 (vvi-miR156b/c/d) in grapevine "Takatsuma" was cloned with a previously cloned grapevine SPL (Vv-SPL9). Expression profiles in 18 grapevine tissues were identified through stem-loop RT-PCR. The interaction mode between vvi-miR156b/c/d and Vv-SPL9 was further validated by detecting the cleavage site and cleavage products of 3'- and 5'-ends via an integrated approach of 5'-RLM-RACE (RNA ligase-mediated 5'-rapid amplification of cDNA ends), 3'-PPM-RACE (poly(A) polymerase-mediated 3'-rapid amplification of cDNA ends), and qRT-PCR (real time reverse transcriptase-polymerase chain reaction). The variation in their cleavage roles in the whole growth stage of grapevine was also systematically investigated. Results showed that vvi-miR156b/c/d exhibited typical temporal-spatial-specific expression levels. The expression levels were higher in vegetative organs, such as leaf, than in reproductive organs, such as tendrils, flowers, and berries. A significant variation was observed during vegetative-to-reproductive transition. The expression patterns of Vv-SPL9 showed the opposite trends with those of vvi-miR156b. We confirmed that the cleavage site was at the 10th site of vvi-miR156b/c/d complementary to Vv-SPL9 in "Takatsuma" grapevine. We also identified the temporal-spatial variation of the cleavage products. This variation can indicate the regulatory function of miR156 on SPL in grapevines. Our findings provide further insights into the functions of vvi-miR156b/c/d and its target Vv-SPL9, and also help enrich our knowledge of small RNA-mediated regulation in grapevine.

MicroRNAs in Control of Plant Development

In the long evolutionary history, plant has evolved elaborate regulatory network to control functional gene expression for surviving and thriving, such as transcription factor-regulated transcriptional programming. However, plenty of evidences from the past decade studies demonstrate that the 21-24 nucleotides small RNA molecules, majorly microRNAs (miRNAs) play dominant roles in post-transcriptional gene regulation through base pairing with their complementary mRNA targets, especially prefer to target transcription factors in plants. Here, we review current progresses on miRNA-controlled plant development, from miRNA biogenesis dysregulation-caused pleiotropic developmental defects to specific developmental processes, such as SAM regulation, leaf and root system regulation, and plant floral transition. We also summarize some miRNAs that are experimentally proved to greatly affect crop plant productivity and quality. In addition, recent reports show that a single miRNA usually displays multiple regulatory roles, such as organ development, phase transition, and stresses responses. Thus, we infer that miRNA may act as a node molecule to coordinate the balance between plant development and environmental clues, which may shed the light on finding key regulator or regulatory pathway for uncovering the mysterious molecular network.

Genome-Wide Investigation Using sRNA-Seq, Degradome-Seq and Transcriptome-Seq Reveals Regulatory Networks of microRNAs and Their Target Genes in Soybean during Soybean mosaic virus Infection

MicroRNAs (miRNAs) play key roles in a variety of cellular processes through regulation of their target gene expression. Accumulated experimental evidence has demonstrated that infections by viruses are associated with the altered expression profile of miRNAs and their mRNA targets in the host. However, the regulatory network of miRNA-mRNA interactions during viral infection remains largely unknown. In this study, we performed small RNA (sRNA)-seq, degradome-seq and as well as a genome-wide transcriptome analysis to profile the global gene and miRNA expression in soybean following infections by three different Soybean mosaic virus (SMV) isolates, L (G2 strain), LRB (G2 strain) and G7 (G7 strain). sRNA-seq analyses revealed a total of 253 soybean miRNAs with a two-fold or greater change in abundance compared with the mock-inoculated control. 125 transcripts were identified as the potential cleavage targets of 105 miRNAs and validated by degradome-seq analyses. Genome-wide transcriptome analysis showed that total 2679 genes are differentially expressed in response to SMV infection including 71 genes predicted as involved in defense response. Finally, complex miRNA-mRNA regulatory networks were derived using the RNAseq, small RNAseq and degradome data. This work represents a comprehensive, global approach to examining virus-host interactions. Genes responsive to SMV infection are identified as are their potential miRNA regulators. Additionally, regulatory changes of the miRNAs themselves are described and the regulatory relationships were supported with degradome data. Taken together these data provide new insights into molecular SMV-soybean interactions and offer candidate miRNAs and their targets for further elucidation of the SMV infection process.

A characterization of grapevine of GRAS domain transcription factor gene family

GRAS domain genes are a group of important plant-specific transcription factors that have been reported to be involved in plant development. In order to know the roles of GRAS genes in grapevine, a widely cultivated fruit crop, the study on grapevine GRAS (VvGRAS) was carried out, and from which, 43 were identified from 12× assemble grapevine genomic sequences. Further, the genomic structures, synteny, phylogeny, expression profiles in different tissues of these genes, and their roles in response to stress were investigated. Among the genes, two potential target genes (VvSCL15 and VvSCL22) for VvmiR171 were experimentally verified by PPM-RACE and RLM-RACE, in that not only the cleavage sites of miR171 on the target mRNA were mapped but also the cleaved fragments and their expressing patterns were detected. Transgenic Arabidopsis plants over expression VvSCL15 showed lower tolerance to drought and salt treatments.

miR-RACE: an effective approach to accurately determine the sequence of computationally identified miRNAs

Computational prediction of microRNAs (miRNAs) is one of the most important approaches in microRNA studies. While validation of the predicted microRNAs' precise sequences is essential for further studies on their biogenesis, evolution, and functions, computational miRNA prediction methods, however, often fail to predict the accurate sequence of the mature miRNA within the precursor at the nucleotide precision level. Here, we depict a highly efficient method for determining the precise sequences of computationally predicted miRNAs. The method combines the generation of miRNA-enriched libraries, with 5'- and 3'-end adaptors being linked to the miRNA molecules, the reverse transcription of small RNAs with an oligo-d(T) anchor primer, two specific 5'- and 3'-miRNA-RACE (miR-RACE) PCR reactions and sequence-directed cloning. The efficiency of this method was demonstrated by the precise sequence validation of computationally predicted miRNAs in citrus, apple, and other fruit crops. Our ongoing research indicates that miR-RACE is also very useful to verify the sequences of putative miRNAs obtained by deep sequencing of small RNA libraries. The protocol of miR-RACE is rapid and can be completed within 2-3 days. miR-RACE should make the bioinformatic prediction of miRNAs more powerful and accurate.

Determination of the precise sequences of computationally predicted miRNAs in Citrus reticulata by miR-RACE and characterization of the related target genes using RLM-RACE

MicroRNAs play vital roles in various biological and metabolic processes by regulating the expression of their target genes in model plants. Since there are limited reports on miRNAs in Citrus reticulata (Crt-miRNAs), the determination of precise sequences of miRNAs is essential to further analyze the functions of miRNAs in Citrus reticulata. Here, miR-RACE, a recently developed technique for determination of the potential miRNAs computationally, was employed to identify the precise sequences of Crt-miRNAs. Tissue- and development-specific expression of nine miRNAs were identified by quantitative RT-PCR in the leaves, stems, flowers and fruits Subsequently, 10 potential target genes were predicated for the eight Crt-miRNAs, most of which were transcription factors and disease resistance proteins. Four target genes were experimentally validated by Poly (A) polymerase-mediated 3′ rapid amplification of cDNA ends and RNA ligase-mediated 5′ rapid amplification of cDNA ends (PPM-RACE and RLM-RACE). Our findings showed that regulatory miRNAs in C. reticulata may play a key role in regulating growth, development, and response to disease. Future work is required to study the functions of miRNAs and their targets of C. reticulata.

Reprogramming of plant cells induced by 6b oncoproteins from the plant pathogen Agrobacterium

Reprogramming of plant cells is an event characterized by dedifferentiation, reacquisition of totipotency, and enhanced cell proliferation, and is typically observed during formation of the callus, which is dependent on plant hormones. The callus-like cell mass, called a crown gall tumor, is induced at the sites of infection by Agrobacterium species through the expression of hormone-synthesizing genes encoded in the T-DNA region, which probably involves a similar reprogramming process. One of the T-DNA genes, 6b, can also by itself induce reprogramming of differentiated cells to generate tumors and is therefore recognized as an oncogene acting in plant cells. The 6b genes belong to a group of Agrobacterium T-DNA genes, which include rolB, rolC, and orf13. These genes encode proteins with weakly conserved sequences and may be derived from a common evolutionary origin. Most of these members can modify plant growth and morphogenesis in various ways, in most cases without affecting the levels of plant hormones. Recent studies have suggested that the molecular function of 6b might be to modify the patterns of transcription in the host nuclei, particularly by directly targeting the host transcription factors or by changing the epigenetic status of the host chromatin through intrinsic histone chaperone activity. In light of the recent findings on zygotic resetting of nucleosomal histone variants in Arabidopsis thaliana, one attractive idea is that acquisition of totipotency might be facilitated by global changes of epigenetic status, which might be induced by replacement of histone variants in the zygote after fertilization and in differentiated cells upon stimulation by plant hormones as well as by expression of the 6b gene.

A Novel Arabidopsis microRNA promotes IAA biosynthesis via the indole-3-acetaldoxime pathway by suppressing superroot1

IAA is a plant hormone that plays important roles in regulating growth and responses to environmental changes. Indole-3-acetaldoxime (IAOx) has been proposed as an important intermediate in the biosynthesis of IAA and two other indole compounds, indole glucosinolates and camalexin. Although the IAOx-dependent IAA biosynthesis pathway has been well studied, the mechanisms of its regulation remain elusive. Here, we report the identification of a novel microRNA, miR10515, which targets superroot1 (SUR1), the gene encoding an indole glucosinolate biosynthetic enzyme. miR10515 was induced by high temperature. Overexpression of MIR10515 resulted in a high-IAA phenotype, while the loss of function of miR10515 resulted in a low-IAA phenotype; these phenotypes were more severe at high temperature. Our results further demonstrated that miR10515 promoted IAA biosynthesis via the IAOx pathway by blocking the indole glucosinolate and camalexin biosynthetic pathways. Phytochrome interacting factor4 (PIF4), a dominant regulator of plant development in response to high temperature, was not required for miR10515 expression. These results provide information on the IAOx metabolic branching point and its biological importance.

Pathogenesis of Soybean mosaic virus in soybean carrying Rsv1 gene is associated with miRNA and siRNA pathways, and breakdown of AGO1 homeostasis

Profiling small RNAs in soybean Williams 82 (rsv), susceptible to Soybean mosaic virus (SMV, the genus Potyvirus, family Potyviridae) strains G2 and G7, and soybean PI96983 (Rsv1), resistant to G2 but susceptible to G7, identified the microRNA miR168 that was highly overexpressed only in G7-infected PI96983 showing a lethal systemic hypersensitive response (LSHR). Overexpression of miR168 was in parallel with the high-level expression of AGO1 mRNA, high-level accumulation of miR168-mediated AGO1 mRNA cleavage products but with severely repressed AGO1 protein. In contrast, AGO1 mRNA, degradation products and protein remained without significant changes in G2- and G7-infected Williams 82. Moreover, knock-down of SGS3, an essential component in RNA silencing, suppressed AGO1 siRNA, partially recovered repressed AGO1 protein, and alleviated LSHR severity in G7-infected Rsv1 soybean. These results suggest that both miRNA and siRNA pathways are involved in G7 infection of Rsv1 soybean, and LSHR is associated with breakdown of AGO1 homeostasis.

Small RNAs as important regulators for the hybrid vigour of super-hybrid rice

Heterosis is an important biological phenomenon; however, the role of small RNA (sRNA) in heterosis of hybrid rice remains poorly described. Here, we performed sRNA profiling of F1 super-hybrid rice LYP9 and its parents using high-throughput sequencing technology, and identified 355 distinct mature microRNAs and trans-acting small interfering RNAs, 69 of which were differentially expressed sRNAs (DES) between the hybrid and the mid-parental value. Among these, 34 DES were predicted to target 176 transcripts, of which 112 encoded 94 transcription factors. Further analysis showed that 67.6% of DES expression levels were negatively correlated with their target mRNAs either in flag leaves or panicles. The target genes of DES were significantly enriched in some important biological processes, including the auxin signalling pathway, in which existed a regulatory network mediated by DES and their targets, closely associated with plant growth and development. Overall, 20.8% of DES and their target genes were significantly enriched in quantitative trait loci of small intervals related to important rice agronomic traits including growth vigour, grain yield, and plant architecture, suggesting that the interaction between sRNAs and their targets contributes to the heterotic phenotypes of hybrid rice. Our findings revealed that sRNAs might play important roles in hybrid vigour of super-hybrid rice by regulating their target genes, especially in controlling the auxin signalling pathway. The above finding provides a novel insight into the molecular mechanism of heterosis.

System analysis of microRNAs in the development and aluminium stress responses of the maize root system

MicroRNAs (miRNAs) are a class of regulatory small RNAs (sRNAs) that down-regulate target genes through mRNA cleavage or translational inhibition. miRNA is known to play an important role in the root development and environmental responses in both the Arabidopsis and rice. However, little information is available to form a complete view of miRNAs in the development of the maize root system and Al stress responses in maize. Four sRNA libraries were generated and sequenced from the early developmental stage of primary roots (PRY), the later developmental stage of maize primary roots (PRO), seminal roots (SR) and crown roots (CR). Through integrative analysis, we identified 278 miRNAs (246 conserved and 32 novel ones) and found that the expression patterns of miRNAs differed dramatically in different maize roots. The potential targets of the identified conserved and novel miRNAs were also predicted. In addition, our data showed that CR is more resistant to Al stress compared with PR and SR, and the differentially expressed miRNAs are likely to play significant roles in different roots in response to environmental stress such as Al stress. Here, we demonstrate that the expression patterns of miRNAs are highly diversified in different maize roots. The differentially expressed miRNAs are correlated with both the development and environmental responses in the maize root. This study not only improves our knowledge about the roles of miRNAs in maize root development but also reveals the potential role of miRNAs in the environmental responses of different maize roots.

A versatile system for functional analysis of genes and microRNAs in cotton

Cotton is an important economic crop worldwide. Due to its long growth cycle, large genome size and recalcitrance to stable transformation, traditional methods for the analysis of gene function in this crop are difficult and labour intensive. Here, we report a cotton leaf crumple virus (CLCrV)-based vector and its application in gene function analysis through virus-induced gene silencing (VIGS) and overexpression of microRNAs (miRNAs), small tandem target mimic (STTM) and artificial miRNA (amiRNA) in cotton via an Agrobacterium-mediated infiltration approach. Using this system, we were able to efficiently silence two endogenous genes, magnesium chelatase subunit I (CHLI) and elongation factor-1α (EF-1α), in Gossypium species and the Bacillus thuringiensis cry1A gene in transgenic cotton. Furthermore, our results show that this vector can be used to ectopically express endogenous miR156 in G. hirsutum, causing a reduction in miR156-targeted RNA transcripts resulting in the development of abnormal leaf phenotypes. Ectopic expression of miR165/166 STTM with this vector led to downward curling and crumpled leaves, and a significant increase in the miR165/166 target mRNAs. This versatile system is easy to use and can provide more uniform and persistent gene silencing in cotton, thereby providing a powerful approach for gene discovery in cotton.

Mla- and Rom1-mediated control of microRNA398 and chloroplast copper/zinc superoxide dismutase regulates cell death in response to the barley powdery mildew fungus

• Barley (Hordeum vulgare L.) Mildew resistance locus a (Mla) confers allele-specific interactions with natural variants of the ascomycete fungus Blumeria graminis f. sp. hordei (Bgh), the causal agent of powdery mildew disease. Significant reprogramming of Mla-mediated gene expression occurs upon infection by this obligate biotrophic pathogen. • We utilized a proteomics-based approach, combined with barley mla, required for Mla12 resistance1 (rar1), and restoration of Mla resistance1 (rom1) mutants, to identify components of Mla-directed signaling. • Loss-of-function mutations in Mla and Rar1 both resulted in the reduced accumulation of chloroplast copper/zinc superoxide dismutase 1 (HvSOD1), whereas loss of function in Rom1 re-established HvSOD1 levels. In addition, both Mla and Rom1 negatively regulated hvu-microRNA398 (hvu-miR398), and up-regulation of miR398 was coupled to reduced HvSOD1 expression. Barley stripe mosaic virus (BSMV)-mediated over-expression of both barley and Arabidopsis miR398 repressed accumulation of HvSOD1, and BSMV-induced gene silencing of HvSod1 impeded Mla-triggered H₂O₂ and hypersensitive reaction (HR) at barley-Bgh interaction sites. • These data indicate that Mla- and Rom1-regulated hvu-miR398 represses HvSOD1 accumulation, influencing effector-induced HR in response to the powdery mildew fungus.

Regulation of MIR genes in response to abiotic stress in Hevea brasiliensis

Increasing demand for natural rubber (NR) calls for an increase in latex yield and also an extension of rubber plantations in marginal zones. Both harvesting and abiotic stresses lead to tapping panel dryness through the production of reactive oxygen species. Many microRNAs regulated during abiotic stress modulate growth and development. The objective of this paper was to study the regulation of microRNAs in response to different types of abiotic stress and hormone treatments in Hevea. Regulation of MIR genes differs depending on the tissue and abiotic stress applied. A negative co-regulation between HbMIR398b with its chloroplastic HbCuZnSOD target messenger is observed in response to salinity. The involvement of MIR gene regulation during latex harvesting and tapping panel dryness (TPD) occurrence is further discussed.

Genome-wide identification and analysis of Japonica and Indica cultivar-preferred transcripts in rice using 983 Affymetrix array data

BACKGROUND

Accumulation of genome-wide transcriptome data provides new insight on a genomic scale which cannot be gained by analyses of individual data. The majority of rice (O. sativa) species are japonica and indica cultivars. Genome-wide identification of genes differentially expressed between japonica and indica cultivars will be very useful in understanding the domestication and evolution of rice species.

RESULTS

In this study, we analyzed 983 of the 1866 entries in the Affymetrix array data in the public database: 595 generated from indica and 388 from japonica rice cultivars. To discover differentially expressed genes in each cultivar, we performed significance analysis of microarrays for normalized data, and identified 490 genes preferentially expressed in japonica and 104 genes in indica. Gene Ontology analyses revealed that defense response-related genes are significantly enriched in both cultivars, indicating that japonica and indica might be under strong selection pressure for these traits during domestication. In addition, 36 (34.6%) of 104 genes preferentially expressed in indica and 256 (52.2%) of 490 genes preferentially expressed in japonica were annotated as genes of unknown function. Biotic stress overview in the MapMan toolkit revealed key elements of the signaling pathway for defense response in japonica or indica eQTLs.

CONCLUSIONS

The percentage of screened genes preferentially expressed in indica was 4-fold higher (34.6%) and that in japonica was 5-fold (52.2%) higher than expected (11.1%), suggesting that genes of unknown function are responsible for the novel traits that distinguish japonica and indica cultivars. The identification of 10 functionally characterized genes expressed preferentially in either japonica or indica highlights the significance of our candidate genes during the domestication of rice species. Functional analysis of the roles of individual components of stress-mediated signaling pathways will shed light on potential molecular mechanisms to improve disease resistance in rice.

Unrelated viral suppressors of RNA silencing mediate the control of ARGONAUTE1 level

Various plant viruses ubiquitously mediate the induction of miR168, resulting in the control of ARGONAUTE 1 (AGO1), which is the pivotal component of the microRNA (miRNA) regulation pathway and can also exhibit antiviral function. Here, we demonstrate that miR168-driven control of AGO1 can persist for a long time in virus-infected plants and can be an important component of symptom development. We also show that infection of RNA viruses belonging to various genera is associated with the transcriptional induction of the MIR168 precursor gene. Moreover, in a transient expression study, we reveal that different unrelated viral suppressors of RNA silencing (VSRs) are responsible for the enhanced accumulation of miR168. The induction of miR168 accumulation is an early function of VSRs and this activity is associated with the control of the endogenous AGO1 protein level. The common ability of unrelated VSRs to induce the miR168 level implies that this activity might be a component of the host defence suppression in plant-virus interactions.

Characterization of target mRNAs for grapevine microRNAs with an integrated strategy of modified RLM-RACE, newly developed PPM-RACE and qPCRs

MicroRNAs (miRNAs) regulate target gene expression by mediating target gene cleavage or inhibition of translation at transcriptional and post-transcriptional levels in higher plants. Until now, many grapevine microRNAs (Vv-miRNAs) have been identified and quite a number of miRNA target genes were also verified by various analysis. However, global interaction of miRNAs with their target genes still remained to perform more research. We reported experimental validation of a number of miRNA target genes in table grapevine that had been previously identified by bioinformatics in our earlier studies. To verify more predicted target genes of Vv-miRNAs and elucidate the modes by which these Vv-miRNAs work on their target genes, 31 unverified potential target genes for 18 Vv-miRNAs were experimentally verified by a new integrated strategy employing a modified 5'-RLM-RACE (RNA ligase-mediated 5' rapid amplification of cDNA ends), 3'-PPM-RACE (poly(A) polymerase-mediated 3' rapid amplification of cDNA ends) and qRT-PCRs of cleavage products. The results showed that these Vv-miRNAs negatively regulated expression of their target messenger RNAs (mRNAs) through guiding corresponding target mRNA cleavage, of which about 94.4% Vv-miRNAs cleaved their target mRNAs mainly at the tenth nucleotide of 5'-end of miRNAs. Expression levels of both miRNAs and their target mRNAs in eight tissues exhibited inverse relationships, and expressions both of cleaved targets and miRNAs indicated a cleavage mode of Vv-miRNAs on their target genes. Our results confirm the importance of Vv-miRNAs in grapevine growth and development, and suggest more study on Vv-miRNAs and targets can enrich the knowledge of miRNA mediated-regulation in grapevine.

The hot pepper (Capsicum annuum) microRNA transcriptome reveals novel and conserved targets: a foundation for understanding MicroRNA functional roles in hot pepper

MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 21 nt in length which play important roles in regulating gene expression in plants. Although many miRNA studies have focused on a few model plants, miRNAs and their target genes remain largely unknown in hot pepper (Capsicum annuum), one of the most important crops cultivated worldwide. Here, we employed high-throughput sequencing technology to identify miRNAs in pepper extensively from 10 different libraries, including leaf, stem, root, flower, and six developmental stage fruits. Based on a bioinformatics pipeline, we successfully identified 29 and 35 families of conserved and novel miRNAs, respectively. Northern blot analysis was used to validate further the expression of representative miRNAs and to analyze their tissue-specific or developmental stage-specific expression patterns. Moreover, we computationally predicted miRNA targets, many of which were experimentally confirmed using 5′ rapid amplification of cDNA ends analysis. One of the validated novel targets of miR-396 was a domain rearranged methyltransferase, the major de novo methylation enzyme, involved in RNA-directed DNA methylation in plants. This work provides the first reliable draft of the pepper miRNA transcriptome. It offers an expanded picture of pepper miRNAs in relation to other plants, providing a basis for understanding the functional roles of miRNAs in pepper.

Plant nucleotide binding site-leucine-rich repeat (NBS-LRR) genes: active guardians in host defense responses

The most represented group of resistance genes are those of the nucleotide binding site-leucine-rich repeat (NBS-LRR) class. These genes are very numerous in the plant genome, and they often occur in clusters at specific loci following gene duplication and amplification events. To date, hundreds of resistance genes and relatively few quantitative trait loci for plant resistance to pathogens have been mapped in different species, with some also cloned. When these NBS-LRR genes have been physically or genetically mapped, many cases have shown co-localization between resistance loci and NBS-LRR genes. This has allowed the identification of candidate genes for resistance, and the development of molecular markers linked to R genes. This review is focused on recent genomics studies that have described the abundance, distribution and evolution of NBS-LRR genes in plant genomes. Furthermore, in terms of their expression, NBS-LRR genes are under fine regulation by cis- and trans-acting elements. Recent findings have provided insights into the roles of alternative splicing, the ubiquitin/ proteasome system, and miRNAs and secondary siRNAs in the regulation of NBS-LRR gene expression at the post-transcriptional, post-translational and epigenetic levels. The possibility to use this knowledge for genetic improvement of plant resistance to pathogens is discussed.

Widespread long noncoding RNAs as endogenous target mimics for microRNAs in plants

Target mimicry is a recently identified regulatory mechanism for microRNA (miRNA) functions in plants in which the decoy RNAs bind to miRNAs via complementary sequences and therefore block the interaction between miRNAs and their authentic targets. Both endogenous decoy RNAs (miRNA target mimics) and engineered artificial RNAs can induce target mimicry effects. Yet until now, only the Induced by Phosphate Starvation1 RNA has been proven to be a functional endogenous microRNA target mimic (eTM). In this work, we developed a computational method and systematically identified intergenic or noncoding gene-originated eTMs for 20 conserved miRNAs in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). The predicted miRNA binding sites were well conserved among eTMs of the same miRNA, whereas sequences outside of the binding sites varied a lot. We proved that the eTMs of miR160 and miR166 are functional target mimics and identified their roles in the regulation of plant development. The effectiveness of eTMs for three other miRNAs was also confirmed by transient agroinfiltration assay.