Genome-wide identification of Medicago truncatula microRNAs and their targets reveals their differential regulation by heavy metal

Root Development in Medicago truncatula: Lessons from Genetics to Functional Genomics

This decade introduced "omics" approaches, such as genomics, transcriptomics, proteomics, and metabolomics in association with reverse and forward genetic approaches, developed earlier, to try to identify molecular pathways involved in the development or in the response to environmental conditions as well as in animals and plants. This review summarizes studies that utilized "omics" strategies to unravel the root development in the model legume Medicago truncatula and how external factors such as soil mineral status or the presence of bacteria and fungi affect root system architecture in this species. We also compare these "omics" data to the knowledges concerning the Arabidopsis thaliana root development, nowadays considered as the model of allorhiz root systems. However, unlike legumes, this species is unable to interact with soil nitrogen-fixing rhizobia and arbuscular-mycorrhizal (AM) fungi to develop novel root-derived symbiotic structures. Differences in root organization, development, and regulatory pathways between these two model species have been highlighted.

Annotation and characterization of Cd-responsive metal transporter genes in rapeseed (Brassica napus)

In higher plants, heavy metal transporters are responsible for metal uptake, translocation and homeostasis. These metals include essential metals such as zinc (Zn) or manganese (Mn) and non-essential metals like cadmium (Cd) or lead (Pb). Although a few heavy metal transporters have been well identified in model plants (e.g. Arabidopsis and rice), little is known about their functionality in rapeseed (Brassica napus). B. napus is an important oil crop ranking the third largest sources of vegetable oil over the world. Importantly, B. napus has long been considered as a desirable candidate for phytoremediation owning to its massive dry weight productivity and moderate to high Cd accumulation. In this study, 270 metal transporter genes (MTGs) from B. napus genome were identified and annotated using bioinformatics and high-throughput sequencing. Most of the MTGs (74.8%, 202/270) were validated by RNA-sequencing (RNA-seq) the seedling libraries. Based on the sequence identity, nine superfamilies including YSL, OPT, NRAMP, COPT, ZIP, CDF/MTP, HMA, MRP and PDR have been classified. RNA-sequencing profiled 202 non-redundant MTGs from B. napus seedlings, of which, 108 MTGs were differentially expressed and 62 genes were significantly induced under Cd stress. These differentially expressed genes (DEGs) are dispersed in the rapeseed genome. Some of the genes were well confirmed by qRT-PCR. Analysis of the genomic distribution of MTGs on B. napus chromosomes revealed that their evolutional expansion was probably through localized allele duplications.

Phenanthrene-responsive microRNAs and their targets in wheat roots

MicroRNAs (miRNAs) play key roles in plant growth, development and responses to abiotic stress. Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental pollutants. However, it is yet unknown how miRNAs work during PAH uptake by plant roots. Thus, in this study we ascertain phenanthrene (a model PAH)-responsive miRNAs using small RNA high-throughput deep sequencing and their target genes in wheat roots. We identified 108 conserved and non-conserved miRNA members belonging to 82 miRNA families and found 11 differentially expressed miRNAs, among which four miRNAs (miR156, miR164, miR171a and miR9678-3p) were up-regulated and the other seven miRNAs (miR398, miR531, miR1121, miR5048-5p, miR9653b, miR9773 and miR9778) were down-regulated. ABC-transporter-related Gene CA704421 and CA697226 did not respond to phenanthrene exposure. miR156 and miR164 might regulate directly the growth and development of wheat roots by targeting SPL and NAC, respectively. miR398 and miR1121 could regulate oxidative reactions to respond to phenanthrene stress. Additionally, miR9773 might involve phenanthrene metabolism through acting on CYP450. Therefore, it is concluded that phenanthrene triggers variation in miRNA expression, which is associated with uptake of and response to phenanthrene. These findings are of significance for further understanding miRNA regulation mechanisms on PAH uptake, and providing guidance for screening of resistant cultivars in crop production and phytoremediation of PAH-contaminated soils or water at genetic level.

Genome-wide characterization and expression analyses of superoxide dismutase (SOD) genes in Gossypium hirsutum

BACKGROUND

Superoxide dismutases (SODs) are a key antioxidant enzyme family, which have been implicated in protecting plants against the toxic effects of reactive oxygen species. Despite current studies have shown that the gene family are involved in plant growth and developmental processes and biotic and abiotic stress responses, little is known about its functional role in upland cotton.

RESULTS

In the present study, we comprehensively analyzed the characteristics of the SOD gene family in upland cotton (Gossypium hirsutum). Based on their conserved motifs, 18 GhSOD genes were identified and phylogenetically classified into five subgroups which corroborated their classifications based on gene-structure patterns and subcellular localizations. The GhSOD sequences were distributed at different densities across 12 of the 26 chromosomes. The conserved domains, gene family evolution cis-acting elements of promoter regions and miRNA-mediated posttranscriptional regulation were predicted and analyzed. In addition, the expression pattern of 18 GhSOD genes were tested in different tissues/organs and developmental stages, and different abiotic stresses and abscisic acid, which indicated that the SOD gene family possessed temporal and spatial specificity expression specificity and may play important roles in reactive oxygen species scavenging caused by various stresses in upland cotton.

CONCLUSIONS

This study describes the first genome-wide analysis of the upland cotton SOD gene family, and the results will help establish a foundation for the further cloning and functional verification of the GhSOD gene family during stress responses, leading to crop 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.

Transcriptomic and physiological analyses of Medicago sativa L. roots in response to lead stress

Lead (Pb) is one of the nonessential and toxic metals that threaten the environment and human health. Medicago sativa L. is a legume with high salt tolerance and high biomass production. It is not only a globally important forage crop but is also an ideal plant for phytoremediation. However, the biological and molecular mechanisms that respond to heavy metals are still not well defined in M. sativa. In this study, de novo and strand-specific RNA-sequencing was performed to identify genes involved in the Pb stress response in M. sativa roots. A total of 415,350 unigenes were obtained from the assembled cDNA libraries, among which 5,416 were identified as significantly differentially expressed genes (DEGs) (false discovery rate < 0.005) between cDNA libraries from control and Pb-treated plants. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs showed they mainly clustered with terms associated with binding, transport, membranes, and the pathways related to signal and energy metabolism. Moreover, a number of candidate genes included antioxidant enzymes, metal transporters, and transcription factors involved in heavy metal response were upregulated under Pb stress. Quantitative real-time PCR(qRT-PCR) validation of the expression patterns of 10 randomly selected candidate DEGs were consistent with the transcriptome analysis results. Thus, this study offers new information towards the investigation of biological changes and molecular mechanisms related to Pb stress response in plants.

Gene Expression Characteristics and Regulation Mechanisms of Superoxide Dismutase and Its Physiological Roles in Plants under Stress

Superoxide dismutases (SODs) are key enzymes functioning as the first line of antioxidant defense by virtue of the ability to convert highly reactive superoxide radicals to hydrogen peroxide and molecular oxygen. SOD plays a central role in protecting plants against the toxic effects of reactive oxygen species generated during normal cellular metabolic activity or as a result of various environmental stresses. Our review focuses on the characteristics of expression of SOD genes, the mechanisms regulating expression of SOD genes at transcriptional, posttranscriptional, and translation levels, and their functional role(s) during development and in response to biotic or abiotic stresses. We propose two important research directions: studying SOD at the genome-wide or proteome-wide level, and improving plant stress tolerances by selecting varieties using transgenic technology.

Variation of DNA methylation patterns associated with gene expression in rice (Oryza sativa) exposed to cadmium

We report genome-wide single-base resolution maps of methylated cytosines and transcriptome change in Cd-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between Cd-exposed and Cd-free rice genomes. There are 2320 non-redundant differentially methylated regions detected in the genome. RNA sequencing revealed 2092 DNA methylation-modified genes differentially expressed under Cd exposure. More genes were found hypermethylated than those hypomethylated in CG, CHH and CHG (where H is A, C or T) contexts in upstream, gene body and downstream regions. Many of the genes were involved in stress response, metal transport and transcription factors. Most of the DNA methylation-modified genes were transcriptionally altered under Cd stress. A subset of loss of function mutants defective in DNA methylation and histone modification activities was used to identify transcript abundance of selected genes. Compared with wide type, mutation of MET1 and DRM2 resulted in general lower transcript levels of the genes under Cd stress. Transcripts of OsIRO2, OsPR1b and Os09g02214 in drm2 were significantly reduced. A commonly used DNA methylation inhibitor 5-azacytidine was employed to investigate whether DNA demethylation affected physiological consequences. 5-azacytidine provision decreased general DNA methylation levels of selected genes, but promoted growth of rice seedlings and Cd accumulation in rice plant.

Identification of a group of XTHs genes responding to heavy metal mercury, salinity and drought stresses in Medicago truncatula

Xyloglucan endotransglucosylase/hydrolases (XTH) are one of the key enzymes regulating cell wall construction, extension and metabolism. In the study, 44 XTH protein genes from Medicago truncatula genome were identified using bioinformatics, microarray and RT-PCR. Each XTH was showed to possess a highly conserved domain ((D/N)-E-(I/L/F/V)-D-(F/I/L)-E-(F/L)-L-G), and most of XTHs possess four Cys in the C terminal region, which suggests the potential for generating disulfide bonds. Based on the XTH protein sequences, these XTHscan be classified into three major families and each family can be subdivided into more groups. Examination of the genomic location of XTH genes on M. truncatula chromosomes showed that the evolutional expansion of the genes was possibly attributed to localized gene duplications. To investigate the possible involvement of the XTHs responding to heavy metals and other abiotic stresses, the XTH genes were exposed to heavy metal (Hg or Cu), salt and drought stresses. There were 28, 21 and 21 MtXTH genes found to respond to HgCl2, salt and drought stresses, respectively, but their expression were different under the stresses. Some of the XTH genes were well confirmed by quantitative RT-PCR (qRT-PCR). We further specified expression of a XTH gene Medtr4g128580 (MtXTH3) under different environmental stresses, and showed that MtXTH3 was induced by Hg exposure. These results indicated that a group of MtXTHs could be differentially expressed under the environmental stresses.

Rice (Oryza sativa) Laccases Involved in Modification and Detoxification of Herbicides Atrazine and Isoproturon Residues in Plants

Atrazine (ATR) and isoproturon (IPU) as herbicides have become serious environmental contaminants due to their overuse in crop production. Although ATR and IPU in soils are easily absorbed by many crops, the mechanisms for their degradation or detoxification in plants are poorly understood. This study identified a group of novel genes encoding laccases (EC 1.10.3.2) that are possibly involved in catabolism or detoxification of ATR and IPU residues in rice. Transcriptome profiling shows at least 22 differentially expressed laccase genes in ATR/IPU-exposed rice. Some of the laccase genes were validated by RT-PCR analysis. The biochemical properties of the laccases were analyzed, and their activities in rice were induced under ATR/IPU exposure. To investigate the roles of laccases in degrading or detoxifying ATR/IPU in rice, transgenic yeast cells (Pichia pastoris X-33) expressing two rice laccase genes (LOC_Os01g63180 and LOC_Os12g15680) were generated. Both transformants were found to accumulate less ATR/IPU compared to the control. The ATR/IPU-degraded products in the transformed yeast cells using UPLC-TOF-MS/MS were further characterized. Two metabolites, hydroxy-dehydrogenated atrazine (HDHA) and 2-OH-isopropyl-IPU, catalyzed by laccases were detected in the eukaryotic cells. These results indicate that the laccase-coding genes identified here could confer degradation or detoxification of the herbicides and suggest that the laccases could be one of the important enzymatic pathways responsible for ATR/IPU degradation/detoxification in rice.

Identification of transcriptome involved in atrazine detoxification and degradation in alfalfa (Medicago sativa) exposed to realistic environmental contamination

Plants are constantly exposed to a variety of toxic compounds (or xenobiotics) such as pesticides (or herbicides). Atrazine (ATZ) as herbicide has become one of the environmental contaminants due to its intensive use during crop production. Plants have evolved strategies to cope with the adverse impact of ATZ. However, the mechanism for ATZ degradation and detoxification in plants is largely unknown. Here we employed a global RNA-sequencing (RNA-Seq) strategy to dissect transcriptome variation in alfalfa (Medicago sativa) exposed to ATZ. Four libraries were constructed including Root-ATZ (root control, ATZ-free), Shoot-ATZ, Root+ATZ (root treated with ATZ) and Shoot+ATZ. Hierarchical clustering was performed to display the expression patterns for all differentially expressed genes (DEGs) under ATZ exposure. Transcripts involved in ATZ detoxification, stress responses (e.g. oxidation and reduction, conjugation and hydrolytic reactions), and regulations of cysteine biosynthesis were identified. Several genes encoding glycosyltransferases, glutathione S-transferases or ABC transporters were up-regulated notably. Also, many other genes involved in oxidation-reduction, conjugation, and hydrolysis for herbicide degradation were differentially expressed. These results suggest that ATZ in alfalfa can be detoxified or degraded through different pathways. The expression patterns of some DEGs by high-throughput sequencing were well confirmed by qRT-PCR. Our results not only highlight the transcriptional complexity in alfalfa exposed to ATZ but represent a major improvement for analyzing transcriptional changes on a large scale as well.

Toxicity responses of Cu and Cd: the involvement of miRNAs and the transcription factor SPL7

BACKGROUND

MicroRNAs are important posttranscriptional regulators of gene expression playing a role in developmental processes as well as in stress responses, including metal stress responses. Despite the identification of several metal-responsive miRNAs, the regulation and the role of these miRNAs and their targets remain to be explored. In this study, miRNAs involved in the response to Cd and Cu excess in Arabidopsis thaliana are identified. In addition, the involvement of the transcription factor SPL7, namely the key regulator of Cu homeostasis, in these metal stress responses is demonstrated by the use of an spl7 knockout mutant. Furthermore, more insight is given in the Cd-induced Cu deficiency response through determining the effects of adding supplemental Cu to Cd-exposed plants.

RESULTS

Thirteen miRNAs were identified in response to Cu and Cd excess in A. thaliana. Several of these miRNAs (miR397a, miR398b/c and miR857) were oppositely affected under Cu and Cd exposure. The induced expression of these miRNAs after Cd exposure was totally abolished in the spl7 mutant (SQUAMOSA promoter binding protein like7), indicating a major role for SPL7 in the Cd response. Plants exposed to Cd showed a higher Cu content in the roots, whereas the Cu content in the leaves of the spl7 mutant was reduced. Furthermore, the Cd-induced Cu deficiency response disappeared when supplemental Cu was added.

CONCLUSIONS

Copper- and Cd-responsive miRNAs were identified and several of them are SPL7-dependently regulated. SPL7 seems to be a shared component between both the Cu toxicity and the Cd toxicity response, yet oppositely regulated, that is inactivated after Cu exposure and activated after Cd exposure. Since SPL7 is the key regulator of Cu homeostasis, and Cd affects the Cu homeostasis, we hypothesize that SPL7 is activated in response to Cd possibly due to a Cd-induced Cu deficiency. Since adding additional Cu to Cd-exposed plants resulted in the disappearance of the Cu deficiency response, Cd possibly provokes Cu deficiency, thereby activating SPL7 and inducing subsequently the Cu deficiency response.

Identification of microRNAs and their targets in Paulownia fortunei plants free from phytoplasma pathogen after methyl methane sulfonate treatment

MicroRNAs (miRNAs) play major roles in plant responses to various biotic and abiotic stresses by regulating gene expression at the transcriptional and post-transcriptional levels. Paulownia witches' broom (PaWB) disease caused by phytoplasmas reduces Paulownia production worldwide. In this study, we investigated the miRNA-mediated plant response to PaWB phytoplasma by Illumina sequencing and degradome analysis of Paulownia fortunei small RNAs (sRNAs). The sRNA and degradome libraries were constructed from healthy and diseased P. fortunei plants and the plants free from phytoplasma pathogen after 60 mg L(-1) methyl methane sulfonate treatment. A total of 96 P. fortunei-conserved miRNAs and 83 putative novel miRNAs were identified. Among them, 37 miRNAs (17 conserved, 20 novel) were found to be differentially expressed in response to PaWB phytoplasma infection. In addition, 114 target genes for 18 of the conserved miRNA families and 33 target genes for 15 of the novel miRNAs in P. fortunei were detected. The expression patterns of 14 of the PaWB phytoplasma-responsive miRNAs and 12 target genes were determined by quantitative real-time polymerase chain reaction (qPCR) experiments. A functional analysis of the miRNA targets indicated that these targeted genes may regulate transcription, stress response, nitrogen metabolism, and various other activities. Our results will help identify the potential roles of miRNAs involved in protecting P. fortunei from diseases.

The Small-RNA Profiles of Almond (Prunus dulcis Mill.) Reproductive Tissues in Response to Cold Stress

Spring frost is an important environmental stress that threatens the production of Prunus trees. However, little information is available regarding molecular response of these plants to the frost stress. Using high throughput sequencing, this study was conducted to identify differentially expressed miRNAs, both the conserved and the non-conserved ones, in the reproductive tissues of almond tolerant H genotype under cold stress. Analysis of 50 to 58 million raw reads led to identification of 174 unique conserved and 59 novel microRNAs (miRNAs). Differential expression pattern analysis showed that 50 miRNA families were expressed differentially in one or both of almond reproductive tissues (anther and ovary). Out of these 50 miRNA families, 12 and 15 displayed up-regulation and down-regulation, respectively. The distribution of conserved miRNA families indicated that miR482f harbor the highest number of members. Confirmation of miRNAs expression patterns by quantitative real- time PCR (qPCR) was performed in cold tolerant (H genotype) alongside a sensitive variety (Sh12 genotype). Our analysis revealed differential expression for 9 miRNAs in anther and 3 miRNAs in ovary between these two varieties. Target prediction of miRNAs followed by differential expression analysis resulted in identification of 83 target genes, mostly transcription factors. This study comprehensively catalogued expressed miRNAs under different temperatures in two reproductive tissues (anther and ovary). Results of current study and the previous RNA-seq study, which was conducted in the same tissues by our group, provide a unique opportunity to understand the molecular basis of responses of almond to cold stress. The results can also enhance the possibility for gene manipulation to develop cold tolerant plants.

Integration of small RNAs, degradome and transcriptome sequencing in hyperaccumulator Sedum alfredii uncovers a complex regulatory network and provides insights into cadmium phytoremediation

The hyperaccumulating ecotype of Sedum alfredii Hance is a cadmium (Cd)/zinc/lead co-hyperaccumulating species of Crassulaceae. It is a promising phytoremediation candidate accumulating substantial heavy metal ions without obvious signs of poisoning. However, few studies have focused on the regulatory roles of miRNAs and their targets in the hyperaccumulating ecotype of S. alfredii. Here, we combined analyses of the transcriptomics, sRNAs and the degradome to generate a comprehensive resource focused on identifying key regulatory miRNA-target circuits under Cd stress. A total of 87 721 unigenes and 356 miRNAs were identified by deep sequencing, and 79 miRNAs were differentially expressed under Cd stress. Furthermore, 754 target genes of 194 miRNAs were validated by degradome sequencing. A gene ontology (GO) enrichment analysis of differential miRNA targets revealed that auxin, redox-related secondary metabolism and metal transport pathways responded to Cd stress. An integrated analysis uncovered 39 pairs of miRNA targets that displayed negatively correlated expression profiles. Ten miRNA-target pairs also exhibited negative correlations according to a real-time quantitative PCR analysis. Moreover, a coexpression regulatory network was constructed based on profiles of differentially expressed genes. Two hub genes, ARF4 (auxin response factor 4) and AAP3 (amino acid permease 3), which might play central roles in the regulation of Cd-responsive genes, were uncovered. These results suggest that comprehensive analyses of the transcriptomics, sRNAs and the degradome provided a useful platform for investigating Cd hyperaccumulation in S. alfredii, and may provide new insights into the genetic engineering of phytoremediation.

Dehydration-responsive miRNAs in foxtail millet: genome-wide identification, characterization and expression profiling

A set of novel and known dehydration-responsive miRNAs have been identified in foxtail millet. These findings provide new insights into understanding the functional role of miRNAs and their respective targets in regulating plant response to dehydration stress. MicroRNAs perform significant regulatory roles in growth, development and stress response of plants. Though the miRNA-mediated gene regulatory networks under dehydration stress remain largely unexplored in plant including foxtail millet (Setaria italica), which is a natural abiotic stress tolerant crop. To find out the dehydration-responsive miRNAs at the global level, four small RNA libraries were constructed from control and dehydration stress treated seedlings of two foxtail millet cultivars showing contrasting tolerance behavior towards dehydration stress. Using Illumina sequencing technology, 55 known and 136 novel miRNAs were identified, representing 22 and 48 miRNA families, respectively. Eighteen known and 33 novel miRNAs were differentially expressed during dehydration stress. After the stress treatment, 32 dehydration-responsive miRNAs were up-regulated in tolerant cultivar and 22 miRNAs were down-regulated in sensitive cultivar, suggesting that miRNA-mediated molecular regulation might play important roles in providing contrasting characteristics to these cultivars. Predicted targets of identified miRNAs were found to encode various transcription factors and functional enzymes, indicating their involvement in broad spectrum regulatory functions and biological processes. Further, differential expression patterns of seven known miRNAs were validated by northern blot and expression of ten novel dehydration-responsive miRNAs were confirmed by SL-qRT PCR. Differential expression behavior of five miRNA-target genes was verified under dehydration stress treatment and two of them also validated by RLM RACE. Overall, the present study highlights the importance of dehydration stress-associated post-transcriptional regulation governed by miRNAs and their targets in a naturally stress-tolerant model crop.

Small RNA and degradome sequencing reveals important microRNA function in Astragalus chrysochlorus response to selenium stimuli

Selenium (Se), an essential element, plays important roles in human health as well as environmental sustainability. Se hyperaccumulating plants are thought as an alternative selenium resource, recently. Astragalus species are known as hyperaccumulator of Se by converting it to nonaminoacid compounds. However, Se-metabolism-related hyperaccumulation is not elucidated in plants yet. MicroRNAs (miRNAs) are key molecules in many biological and metabolic processes via targeting mRNAs, which may also play an important role in Se accumulation in plants. In this study, we identified 418 known miRNAs, belonging to 380 families, and 151 novel miRNAs induced by Se exposure in Astragalus chyrsochlorus callus. Among known miRNAs, the expression of 287 families was common in both libraries, besides 71 families were expressed only in Se-treated sample, whereas 60 conserved families were expressed in control tissue. miR1507a, miR1869 and miR2867-3p were mostly up-regulated, whereas miR1507-5p and miR8781b were significantly down-regulated by Se exposure. Computational analysis shows that the targets of miRNAs are involved in different types of biological mechanisms including 47 types of cellular component, 103 types of molecular function and 144 types of biological process. Degradome analysis shows that 1256 mRNAs were targeted by 499 miRNAs. We conclude that some known and novel miRNAs such as miR167a, miR319, miR1507a, miR4346, miR7767-3p, miR7800, miR9748 and miR-n93 target transcription factors, disease resistance proteins and some specific genes like cysteine synthase and might be related to plant hormone signal transduction, plant-pathogen interaction and sulphur metabolism pathways.

Genome-wide identification of DNA methylation provides insights into the association of gene expression in rice exposed to pesticide atrazine

Atrazine (ATR) is a pesticide widely used for controlling weeds for crop production. Crop contamination with ATR negatively affects crop growth and development. This study presents the first genome-wide single-base-resolution maps of DNA methylation in ATR-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between the ATR-exposed and ATR-free (control) rice. Most of DNA methyltransferases, histone methyltransferases and DNA demethylase were differentially regulated by ATR. We found more genes hypermethylated than those hypomethylated in the regions of upstream, genebody and downstream under ATR exposure. A stringent group of 674 genes (p < 0.05, two-fold change) with a strong preference of differential expression in ATR-exposed rice was identified. Some of the genes were identified in a subset of loss of function mutants defective in DNA methylation/demethylation. Provision of 5-azacytidine (AZA, inhibitor of DNA methylation) promoted the rice growth and reduced ATR content. By UPLC/Q-TOF-MS/MS, 8 degraded products and 9 conjugates of ATR in AZA-treated rice were characterized. Two of them has been newly identified in this study. Our data show that ATR-induced changes in DNA methylation marks are possibly involved in an epigenetic mechanism associated with activation of specific genes responsible for ATR degradation and detoxification.

MicroRNA390 Is Involved in Cadmium Tolerance and Accumulation in Rice

Cadmium (Cd) is a non-essential heavy metal that is toxic to plants. microRNAs (miRNAs) are 21-nucleotide RNAs that are ubiquitous regulators of gene expression at the post-transcriptional level. Several plant miRNAs, such as miR390, have vital roles in plant growth, development and responses to environmental stresses including heavy metal stress. In this study, the expression of mature miR390 was significantly down-regulated under Cd stress in rice. Consequently, the target gene of miR390, OsSRK was dramatically induced by Cd treatment. Transgenic rice plants overexpressing miR390 displayed reduced Cd tolerance and higher Cd accumulation compared with wild-type plants. Simultaneously, expression of OsSRK was less pronounced in 35S:MIR390 plants than in wild-type. These results indicate that miR390 was a negative regulator involved in Cd stress tolerance in rice.

MicroRNAs As Potential Targets for Abiotic Stress Tolerance in Plants

The microRNAs (miRNAs) are small (20-24 nt) sized, non-coding, single stranded riboregulator RNAs abundant in higher organisms. Recent findings have established that plants assign miRNAs as critical post-transcriptional regulators of gene expression in sequence-specific manner to respond to numerous abiotic stresses they face during their growth cycle. These small RNAs regulate gene expression via translational inhibition. Usually, stress induced miRNAs downregulate their target mRNAs, whereas, their downregulation leads to accumulation and function of positive regulators. In the past decade, investigations were mainly aimed to identify plant miRNAs, responsive to individual or multiple environmental factors, profiling their expression patterns and recognizing their roles in stress responses and tolerance. Altered expressions of miRNAs implicated in plant growth and development have been reported in several plant species subjected to abiotic stress conditions such as drought, salinity, extreme temperatures, nutrient deprivation, and heavy metals. These findings indicate that miRNAs may hold the key as potential targets for genetic manipulations to engineer abiotic stress tolerance in crop plants. This review is aimed to provide recent updates on plant miRNAs, their biogenesis and functions, target prediction and identification, computational tools and databases available for plant miRNAs, and their roles in abiotic stress-responses and adaptive mechanisms in major crop plants. Besides, the recent case studies for overexpressing the selected miRNAs for miRNA-mediated enhanced abiotic stress tolerance of transgenic plants have been discussed.

RNA-sequencing Oryza sativa transcriptome in response to herbicide isoprotruon and characterization of genes involved in IPU detoxification

The comprehensive analysis of transcriptome and UPLC-MS/MS in rice was performed to explore the regulatory mechanism of mRNA level and chemical metabolism in response to herbicide isoproturon.

Genome-wide identification of abiotic stress-regulated and novel microRNAs in mulberry leaf

As the most important food plant for sericultural industry, mulberry trees have to suffer from a wide range of abiotic and biotic stresses, such as drought and high salinity. MicroRNAs (miRNAs) have been proved to play important roles in abiotic stresses regulation in many plants. However, there are seldom reports on the miRNAs expression profiles upon abiotic challenges in mulberry. In this study, three small RNA libraries from mulberry leaf tissue with or without drought or salt treatment were constructed and deep sequenced. Total of 48 conserved miRNAs (including miRNA*) and 162 novel miRNAs were identified (processing precision value>0.1). A total of 270 and 1963 target genes were predicted for conserved miRNAs and novel miRNAs, respectively. 13 differentially expressed miRNAs were detected under drought or salt stresses by deep sequencing and qRT-PCR. 5' RLM-RACE validated Morus 013341 to be the target gene of miR-395a. Our results provided initial clue to further study molecular mechanism on abiotic stresses regulation in mulberry.

Transcriptome-wide analysis of chromium-stress responsive microRNAs to explore miRNA-mediated regulatory networks in radish (Raphanus sativus L.)

MicroRNAs (miRNAs) are small noncoding RNAs that play pivotal roles in plant growth, development and stress response. Chromium (Cr) is one of common environmental contaminants possessing potential health hazards to living organisms. To date, little is known about the regulatory roles of miRNAs in response to Cr stress in radish. To systematically identify Cr-responsive miRNAs and their targets in radish, two sRNA libraries derived from Cr-free (CK) and Cr-treated (Cr200) roots were constructed. With Solexa sequencing, 81 known and 72 novel miRNAs were identified, from which 54 known and 16 novel miRNAs were significantly differentially expressed under Cr stress. Several target genes for Cr-responsive miRNAs encode different transcription factor (TF) families, including SPLs, MYBs, ERFs and bZIPs, might regulate corresponding HM-related transcriptional processes in plants. Notably, a few key responsive enzymes or proteins, including HMA, YSL1 and ABC transporter protein were involved in Cr uptake and homeostasis process. Furthermore, the expression patterns of some Cr-responsive miRNAs and their targets were validated by RT-qPCR. This study represents the first characterization of Cr-responsive miRNAs and their targets in radish. The outcomes of this study could provide novel insights into miRNA-mediated regulatory mechanisms underlying plant response to Cr stress in root vegetable crops.

A collection of cytochrome P450 monooxygenase genes involved in modification and detoxification of herbicide atrazine in rice (Oryza sativa) plants

Plant cytochrome P450 monooxygenases constitute one of the largest families of protein genes involved in plant growth, development and acclimation to biotic and abiotic stresses. However, whether these genes respond to organic toxic compounds and their biological functions for detoxifying toxic compounds such as herbicides in rice are poorly understood. The present study identified 201 genes encoding cytochrome P450s from an atrazine-exposed rice transcriptome through high-throughput sequencing. Of these, 69 cytochrome P450 genes were validated by microarray and some of them were confirmed by real time PCR. Activities of NADPH-cytochrome P450 reductase (CPR) and p-nitroanisole O-demethylase (PNOD) related to toxicity were determined and significantly induced by atrazine exposure. To dissect the mechanism underlying atrazine modification and detoxification by P450, metabolites (or derivatives) of atrazine in plants were analyzed by ultra performance liquid chromatography mass spectrometry (UPLC/MS). Major metabolites comprised desmethylatrazine (DMA), desethylatrazine (DEA), desisopropylatrazine (DIA), hydroxyatrazine (HA), hydroxyethylatrazine (HEA) and hydroxyisopropylatrazine (HIA). All of them were chemically modified by P450s. Furthermore, two specific inhibitors of piperonyl butoxide (PBO) and malathion (MAL) were used to assess the correlation between the P450s activity and rice responses including accumulation of atrazine in tissues, shoot and root growth and detoxification.

Discovery of microRNAs and transcript targets related to witches' broom disease in Paulownia fortunei by high-throughput sequencing and degradome approach

Paulownia witches' broom (PaWB) caused by the phytoplasma is a devastating disease of Paulownia trees. It has caused heavy yield losses to Paulownia production worldwide. However, knowledge of the transcriptional and post-transcriptional regulation of gene expression by microRNAs (miRNAs), especially miRNAs responsive to PaWB disease stress, is still rudimentary. In this study, to identify miRNAs and their transcript targets that are responsive to PaWB disease stress, six sequencing libraries were constructed from healthy (PF), PaWB-infected (PFI), and PaWB-infected, 20 mg L(-1) methyl methane sulfonate-treated (PFI20) P. fortunei seedlings. As a result, 95 conserved miRNAs belonging to 18 miRNA families, as well as 122 potential novel miRNAs, were identified. Most of them were found to be a response to PaWB disease-induced stress, and the expression levels of these miRNAs were validated by quantitative real-time PCR analysis. The study simultaneously identified 109 target genes from the P. fortunei for 14 conserved miRNA families and 24 novel miRNAs by degradome sequencing. Furthermore, the functions of the miRNA targets were annotated based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. The results presented here provide the groundwork for further analysis of miRNAs and target genes responsive to the PaWB disease stress, and could be also useful for addressing new questions to better understand the mechanisms of plant infection by phytoplasma in the future.

Plant miRNAs: biogenesis, organization and origins

MicroRNAs, or miRNAs, are posttranscriptional regulators of gene expression. A wealth of observations and findings suggest highly complex, multicomponent, and intermingled pathways governing miRNA biogenesis and miRNA-mediated gene silencing. Plant miRNA genes are usually found as individual entities scattered around the intergenic and-to a much lesser extent-intragenic space, while miRNA gene clusters, formed by tandem or segmental duplications, also exist in plant genomes. Genome duplications are proposed to contribute to miRNA family expansions, as well. Evolutionarily young miRNAs retaining extensive homology to their loci of origin deliver important clues into miRNA origins and evolution. Additionally, imprecisely processed miRNAs evidence noncanonical routes of biogenesis, which may affect miRNA expression levels or targeting capabilities. Majority of the knowledge regarding miRNAs comes from model plant species. As ongoing research progressively expands into nonmodel systems, our understanding of miRNAs and miRNA-related pathways changes which opens up new perspectives and frontiers in miRNA research.

Identification and characterization of long non-coding RNAs involved in osmotic and salt stress in Medicago truncatula using genome-wide high-throughput sequencing

BACKGROUND

Long non-coding RNAs (lncRNAs) have been shown to play crucially regulatory roles in diverse biological processes involving complex mechanisms. However, information regarding the number, sequences, characteristics and potential functions of lncRNAs in plants is so far overly limited.

RESULTS

Using high-throughput sequencing and bioinformatics analysis, we identified a total of 23,324 putative lncRNAs from control, osmotic stress- and salt stress-treated leaf and root samples of Medicago truncatula, a model legume species. Out of these lncRNAs, 7,863 and 5,561 lncRNAs were identified from osmotic stress-treated leaf and root samples, respectively. While, 7,361 and 7,874 lncRNAs were identified from salt stress-treated leaf and root samples, respectively. To reveal their potential functions, we analyzed Gene Ontology (GO) terms of genes that overlap with or are neighbors of the stress-responsive lncRNAs. Enrichments in GO terms in biological processes such as signal transduction, energy synthesis, molecule metabolism, detoxification, transcription and translation were found.

CONCLUSIONS

LncRNAs are likely involved in regulating plant's responses and adaptation to osmotic and salt stresses in complex regulatory networks with protein-coding genes. These findings are of importance for our understanding of the potential roles of lncRNAs in responses of plants in general and M. truncatula in particular to abiotic stresses.

Genome-wide analysis reveals dynamic changes in expression of microRNAs during vascular cambium development in Chinese fir, Cunninghamia lanceolata

MicroRNAs (miRNAs) are small noncoding regulatory RNAs that play key roles in the process of plant development. To date, extensive studies of miRNAs have been performed in a few model plants, but few efforts have focused on small RNAs (sRNAs) in conifers because of the lack of reference sequences for their enormous genomes. In this study, Solexa sequencing of three sRNA libraries obtained from dormant, reactivating, and active vascular cambium in Chinese fir (Cunninghamia lanceolata) using tangential cryosectioning identified 20 known miRNA families and 18 novel potential miRNAs, of which nine novel miRNA precursors were validated by RT-PCR and sequencing. More than half of these novel miRNAs displayed stage-specific expression patterns in the vascular cambium. Furthermore, analysing the 103 miRNAs and their predicted targets indicated that about 70% appeared to negatively regulate their targets, of which two target genes involved in the regulation of cambial cell division were validated via RNA ligase-mediated rapid amplification of 5'-cDNA ends (RLM 5'-RACE) and transient co-expression in Nicotiana benthamiana leaves. Interestingly, miRNA156 and miRNA172 may regulate the phase transition in vascular cambium from dormancy to active growth. These results provide new insights into the important regulatory functions of miRNAs in vascular cambium development and wood formation in conifers.

Cesium Toxicity Alters MicroRNA Processing and AGO1 Expressions in Arabidopsis thaliana

MicroRNAs (miRNAs) are short RNA fragments that play important roles in controlled gene silencing, thus regulating many biological processes in plants. Recent studies have indicated that plants modulate miRNAs to sustain their survival in response to a variety of environmental stimuli, such as biotic stresses, cold, drought, nutritional starvation, and toxic heavy metals. Cesium and radio-cesium contaminations have arisen as serious problems that both impede plant growth and enter the food chain through contaminated plants. Many studies have been performed to define plant responses against cesium intoxication. However, the complete profile of miRNAs in plants during cesium intoxication has not been established. Here we show the differential expression of the miRNAs that are mostly down-regulated during cesium intoxication. Furthermore, we found that cesium toxicity disrupts both the processing of pri-miRNAs and AGONOUTE 1 (AGO1)-mediated gene silencing. AGO 1 seems to be especially destabilized by cesium toxicity, possibly through a proteolytic regulatory pathway. Our study presents a comprehensive profile of cesium-responsive miRNAs, which is distinct from that of potassium, and suggests two possible mechanisms underlying the cesium toxicity on miRNA metabolism.

Transcriptome-Wide Identification of miRNAs and Their Targets from Typha angustifolia by RNA-Seq and Their Response to Cadmium Stress

MicroRNAs (miRNAs) play important roles in plant responses to environmental stress. In this work, we used high-throughput sequencing to analyze transcriptome and small RNAs (sRNAs) in Typha angustifolia under cadmium (Cd) stress. 57,608,230 raw reads were obtained from deep sequencing of a pooled cDNA library. Sequence assembly and analysis yielded 102,473 unigenes. We subsequently sequenced two sRNA libraries from T. angustifolia with or without Cd exposure respectively. Based on transcriptome data of T. angustifolia, we catalogued and analyzed the sRNAs, resulting in the identification of 114 conserved miRNAs and 41 novel candidate miRNAs in both small RNA libraries. In silico analysis revealed 764 targets for 89 conserved miRNAs and 21 novel miRNAs. Statistical analysis on sequencing reads abundance and experimental validation revealed that 4 conserved and 6 novel miRNAs showed specific expression. Combined with function of target genes, these results suggested that miRNAs might play a role in plant Cd stress response. This study provided the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in T. angustifolia, which provide a framework for further analysis of miRNAs and their role in regulating plant responses to Cd stress.

The F-box family genes as key elements in response to salt, heavy mental, and drought stresses in Medicago truncatula

F-box protein is a subunit of Skp1-Rbx1-Cul1-F-box protein (SCF) complex with typically conserved F-box motifs of approximately 40 amino acids and is one of the largest protein families in eukaryotes. F-box proteins play critical roles in selective and specific protein degradation through the 26S proteasome. In this study, we bioinformatically identified 972 putative F-box proteins from Medicago truncatula genome. Our analysis showed that in addition to the conserved motif, the F-box proteins have several other functional domains in their C-terminal regions (e.g., LRRs, Kelch, FBA, and PP2), some of which were found to be M. truncatula species-specific. By phylogenetic analysis of the F-box motifs, these proteins can be classified into three major families, and each family can be further grouped into more subgroups. Analysis of the genomic distribution of F-box genes on M. truncatula chromosomes revealed that the evolutional expansion of F-box genes in M. truncatula was probably due to localized gene duplications. To investigate the possible response of the F-box genes to abiotic stresses, both publicly available and customer-prepared microarrays were analyzed. Most of the F-box protein genes can be responding to salt and heavy metal stresses. Real-time PCR analysis confirmed that some of the F-box protein genes containing heat, drought, salicylic acid, and abscisic acid responsive cis-elements were able to respond to the abiotic stresses.

MicroRNA390-Directed TAS3 Cleavage Leads to the Production of tasiRNA-ARF3/4 During Somatic Embryogenesis in Dimocarpus longan Lour

Trans-acting short-interfering RNAs (tasiRNAs) originate from TAS3 families through microRNA (miRNA) 390-guided cleavage of primary transcripts and target auxin response factors (ARF3/-4), which are involved in the normal development of lateral roots and flowers in plants. However, their roles in embryo development are still unclear. Here, the pathway miR390-TAS3-ARF3/-4 was identified systematically for the first time during somatic embryo development in Dimocarpus longan. We identified the miR390 primary transcript and promoter. The promoter contained cis-acting elements responsive to stimuli such as light, salicylic acid, anaerobic induction, fungal elicitor, circadian control, and heat stress. The longan TAS3 transcript, containing two miR390-binding sites, was isolated; the miR390- guided cleavage site located near the 3' end of the TAS3 transcript was verified. Eight TAS3-tasiRNAs with the 21-nucleotides phase were found among longan small RNA data, further confirming that miR390-directed TAS3 cleavage leads to the production of tasiRNA in longan. Among them, TAS3_5'D5+ and 5'D6+ tasiRNAs were highly abundant, and verified to target ARF3 and -4, implying that miR390-guided TAS3 cleavage with 21-nucleotides phase leading to the production of tasiRNA-ARF is conserved in plants. Pri-miR390 was highly expressed in friable-embryogenic callus (EC), and less expressed in incomplete compact pro-embryogenic cultures, while miR390 showed its lowest expression in EC and highest expression in torpedo-shaped embryos (TEs). DlTAS3 and DlARF4 both exhibited their lowest expressions in EC, and reached their peaks in the globular embryos stage, which were mainly inversely proportional to the expression of miR390, especially at the globular embryos to cotyledonary embryos (CEs) stages. While DlARF3 showed little variation from the EC to TEs stages, and exhibited its lowest expression in the CEs stage. There was a general lack of correlation between the expressions of DlARF3 and miR390. In addition, pri-miR390, DlTAS3, DlARF3 and -4 were up-regulated by 2,4-D in a concentration-dependent manner. They were also preferentially expressed in roots, pulp, and seeds of 'Sijimi' longan, implying their extended roles in the development of longan roots and fruit. This study provided insights into a possible role of miR390-tasiRNAs-ARF in plant somatic embryo development.

miRNA plays a role in the antagonistic effect of selenium on arsenic stress in rice seedlings

We report that miRNA plays an important role in the selenium mediated antagonistic effect on arsenic stress in rice seedlings.

Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play versatile roles in plants, functioning in processes such as growth, development and stress responses. Chilling is a common abiotic stress that seriously affects plants growth and development. Recently, chilling-responsive miRNAs have been detected in several plant species. However, little is known about the miRNAs in the model plant tomato. 'LA1777' (Solanum habrochaites) has been shown to survive chilling stress due to its various characteristics.

RESULTS

Here, two small RNA libraries and two degradome libraries were produced from chilling-treated (CT) and non-chilling-treated (NT) leaves of S. habrochaites seedlings. Following high-throughput sequencing and filtering, 161 conserved and 236 novel miRNAs were identified in the two libraries. Of these miRNAs, 192 increased in the response to chilling stress while 205 decreased. Furthermore, the target genes of the miRNAs were predicted using a degradome sequencing approach. It was found that 62 target genes were cleaved by 42 conserved miRNAs, while nine target genes were cleaved by nine novel miRNAs. Additionally, nine miRNAs and six target genes were validated by quantitative real-time PCR (qRT-PCR). Target gene functional analysis showed that most target genes played positive roles in the chilling response, primarily by regulating the expression of anti-stress proteins, antioxidant enzyme and genes involved in cell wall formation.

CONCLUSIONS

Tomato is an important model plant for basic biological research. In this study, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, and the target genes were analyzed by degradome sequencing. The work helps identify chilling-responsive miRNAs in tomato and increases the number of identified miRNAs involved in chilling stress. Furthermore, the work provides a foundation for further study of the regulation of miRNAs in the plant response to chilling stress.

Genome-wide identification of different dormant Medicago sativa L. MicroRNAs in response to fall dormancy

BACKGROUND

MicroRNAs (miRNAs) are a class of regulatory small RNAs (sRNAs) that regulate gene post-transcriptional expression in plants and animals. High-throughput sequencing technology is capable of identifying small RNAs in plant species. Alfalfa (Medicago sativa L.) is one of the most widely cultivated perennial forage legumes worldwide, and fall dormancy is an adaptive characteristic related to the biomass production and winter survival in alfalfa. Here, we applied high-throughput sRNA sequencing to identify some miRNAs that were responsive to fall dormancy in standard variety (Maverick and CUF101) of alfalfa.

RESULTS

Four sRNA libraries were generated and sequenced from alfalfa leaves in two typical varieties at distinct seasons. Through integrative analysis, we identified 51 novel miRNA candidates of 206 families. Additionally, we identified 28 miRNAs associated with fall dormancy in standard variety (Maverick and CUF101), including 20 known miRNAs and eight novel miRNAs. Both high-throughput sequencing and RT-qPCR confirmed that eight known miRNA members were up-regulated and six known miRNA members were down-regulated in response to fall dormancy in standard variety (Maverick and CUF101). Among the 51 novel miRNA candidates, five miRNAs were up-regulated and three miRNAs were down-regulated in response to fall dormancy in standard variety (Maverick and CUF101), and five of them were confirmed by Northern blot analysis.

CONCLUSION

We identified 20 known miRNAs and eight new miRNA candidates that were responsive to fall dormancy in standard variety (Maverick and CUF101) by high-throughput sequencing of small RNAs from Medicago sativa. Our data provide a useful resource for investigating miRNA-mediated regulatory mechanisms of fall dormancy in alfalfa, and these findings are important for our understanding of the roles played by miRNAs in the response of plants to abiotic stress in general and fall dormancy in alfalfa.

Genome-wide fungal stress responsive miRNA expression in wheat

MicroRNAs (miRNAs) are small non-coding class of RNAs. They were identified in many plants with their diverse regulatory roles in several cellular and metabolic processes. A number of miRNAs were involved in biotic and abiotic stress responses. Here, fungal stress responsive wheat miRNAs were analyzed by using miRNA-microarray strategy. Two different fungi (Fusarium culmorum and Bipolaris sorokiniana) were inoculated on resistant and sensitive wheat cultivars. A total of 87 differentially regulated miRNAs were detected in the 8 × 15 K array including all of the available plant miRNAs. Using bioinformatics tools, the target transcripts of responsive miRNAs were predicted, and related biological processes and mechanisms were assessed. A number of the miRNAs such as miR2592s, miR869.1, miR169b were highly differentially regulated showing more than 200-fold change upon fungal-inoculation. Some of the miRNAs were identified as fungal-inoculation responsive for the first time. The analyses showed that some of the differentially regulated miRNAs targeted resistance-related genes such as LRR, glucuronosyl transferase, peroxidase and Pto kinase. The comparison of the two miRNA-microarray analyses indicated that fungal-responsive wheat miRNAs were differentially regulated in pathogen- and cultivar-specific manners.

Orthologous plant microRNAs: microregulators with great potential for improving stress tolerance in plants

Small RNAs that are highly conserved across many plant species are involved in stress responses. Plants are exposed to many types of unfavorable conditions during their life cycle that result in some degree of stress. Recent studies on microRNAs (miRNAs) have highlighted their great potential as regulators of stress tolerance in plants. One of the possible ways in which plants counter environmental stresses is by altering their gene expression by the action of miRNAs. miRNAs regulate the expression of target genes by hybridizing to their nascent reverse complementary sequences marking them for cleavage in the nucleus or translational repression in the cytoplasm. Some miRNAs have been reported to be key regulators in biotic as well as abiotic stress responses across many species. The present review highlights some of the regulatory roles of orthologous plant miRNAs in response to various types of stress conditions.

Differential profiling analysis of miRNAs reveals a regulatory role in low N stress response of Populus

Nitrogen (N) is an essential mineral element for plant growth processes, and its availability severely affects the productivity of plants, especially trees. MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 21 nucleotides in length that play important roles in plant growth, development and stress responses. To identify Populus miRNAs and their functions in response to nutrition stress, high-throughput sequencing was performed using Populus tomentosa plantlets treated with or without low concentrations of N. We identified 160 conserved miRNAs, 15 known but non-conserved miRNAs, 2 candidate novel miRNAs and 71 corresponding miRNA*s. Differential expression analysis showed that expression of the 21 conserved miRNA families was significantly altered. Real-time quantitative PCR (qPCR) was used to further validate and analyze the dynamic expression of the identified miRNAs. A total of 218 target genes from the low-N-responsive miRNAs were predicted, and their functions were further annotated in combination with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. These results suggest that miRNAs play important roles in the response of Populus to low N stress. Furthermore, this study provides the first identification and profiles of N stress-responsive miRNAs from trees.

The role of small RNAs in wide hybridisation and allopolyploidisation between Brassica rapa and Brassica nigra

BACKGROUND

An allopolyploid formation consists of the two processes of hybridisation and chromosome doubling. Hybridisation makes a different genome combined in the same cell, and genome "shock" and instability occur during this process, whereas chromosome doubling results in doubling and reconstructing the genome dosage. Recent studies have demonstrated that small RNAs, play an important role in maintaining the genome reconstruction and stability. However, to date, little is known regarding the role of small RNAs during the process of wide hybridisation and chromosome doubling, which is essential to elucidate the mechanism of polyploidisation. Therefore, the genetic and DNA methylation alterations and changes in the siRNA and miRNA were assessed during the formation of an allodiploid and its allotetraploid between Brassica rapa and Brassica nigra in the present study.

RESULTS

The phenotypic analysis exhibited that the allotetraploid had high heterosis compared with their parents and the allodiploid. The methylation-sensitive amplification polymorphism (MSAP) analysis indicated that the proportion of changes in the methylation pattern of the allodiploid was significantly higher than that found in the allotetraploid, while the DNA methylation ratio was higher in the parents than the allodiploid and allotetraploid. The small RNAs results showed that the expression levels of miRNAs increased in the allodiploid and allotetraploid compared with the parents, and the expression levels of siRNAs increased and decreased compared with the parents B. rapa and B. nigra, respectively. Moreover, the percentages of miRNAs increased with an increase in the polyploidy levels, but the percentages of siRNAs and DNA methylation alterations decreased with an increase in the polyploidy levels. Furthermore, qRT-PCR analysis showed that the expression levels of the target genes were negatively corrected with the expressed miRNAs.

CONCLUSIONS

The study showed that siRNAs and DNA methylation play an important role in maintaining the genome stability in the formation of an allotetraploid. The miRNAs regulate gene expression and induce the phenotype variation, which may play an important role in the occurrence of heterosis in the allotetraploid. The findings of this study may provide new information for elucidating that the allotetraploids have a growth advantage over the parents and the allodiploids.

Integrated analysis of miRNA and mRNA expression profiles in response to Cd exposure in rice seedlings

BACKGROUND

Independent transcriptome profile analyses of miRNAs or mRNAs under conditions of cadmium (Cd) stress have been widely reported in plants. However, a combined analysis of sRNA sequencing expression data with miRNA target expression data to infer the relative activities of miRNAs that regulate gene expression changes resulting from Cd stress has not been reported in rice. To elucidate the roles played by miRNAs in the regulation of changes in gene expression in response to Cd stress in rice (Oryza sativa L.), we simultaneously characterized changes in the miRNA and mRNA profiles following treatment with Cd.

RESULTS

A total of 163 miRNAs and 2,574 mRNAs were identified to be differentially expressed under Cd stress, and the changes in the gene expression profile in the shoot were distinct from those in the root. At the miRNA level, 141 known miRNAs belonging to 48 families, and 39 known miRNAs in 23 families were identified to be differentially expressed in the root and shoot, respectively. In addition, we identified eight new miRNA candidates from the root and five from the shoot that were differentially expressed in response to Cd treatment. For the mRNAs, we identified 1,044 genes in the root and 448 genes in the shoot that were up-regulated, while 572 and 645 genes were down-regulated in the root and shoot, respectively. GO and KEGG enrichment analyses showed that genes encoding secondary, metabolite synthases, signaling molecules, and ABC transporters were significantly enriched in the root, while only ribosomal protein and carotenoid biosynthesis genes were significantly enriched in the shoot. Then 10 known miRNA-mRNA interaction pairs and six new candidate ones, that showed the opposite expression patterns, were identified by aligning our two datasets against online databases and by using the UEA sRNA toolkit respectively.

CONCLUSIONS

This study is the first to use high throughput DNA sequencing to simultaneously detect changes in miRNA and mRNA expression patterns in the root and shoot in response to Cd treatment. These integrated high-throughput expression data provide a valuable resource to examine global genome expression changes in response to Cd treatment and how these are regulated by miRNAs.

The small RNA diversity from Medicago truncatula roots under biotic interactions evidences the environmental plasticity of the miRNAome

BACKGROUND

Legume roots show a remarkable plasticity to adapt their architecture to biotic and abiotic constraints, including symbiotic interactions. However, global analysis of miRNA regulation in roots is limited, and a global view of the evolution of miRNA-mediated diversification in different ecotypes is lacking.

RESULTS

In the model legume Medicago truncatula, we analyze the small RNA transcriptome of roots submitted to symbiotic and pathogenic interactions. Genome mapping and a computational pipeline identify 416 miRNA candidates, including known and novel variants of 78 miRNA families present in miRBase. Stringent criteria of pre-miRNA prediction yield 52 new mtr-miRNAs, including 27 miRtrons. Analyzing miRNA precursor polymorphisms in 26 M. truncatula ecotypes identifies higher sequence polymorphism in conserved rather than Medicago-specific miRNA precursors. An average of 19 targets, mainly involved in environmental responses and signalling, is predicted per novel miRNA. We identify miRNAs responsive to bacterial and fungal pathogens or symbionts as well as their related Nod and Myc-LCO symbiotic signals. Network analyses reveal modules of new and conserved co-expressed miRNAs that regulate distinct sets of targets, highlighting potential miRNA-regulated biological pathways relevant to pathogenic and symbiotic interactions.

CONCLUSIONS

We identify 52 novel genuine miRNAs and large plasticity of the root miRNAome in response to the environment, and also in response to purified Myc/Nod signaling molecules. The new miRNAs identified and their sequence variation across M. truncatula ecotypes may be crucial to understand the adaptation of root growth to the soil environment, notably in the agriculturally important legume crops.

Transcriptome/Degradome-wide discovery of microRNAs and transcript targets in two Paulownia australis genotypes

MicroRNAs (miRNAs) are involved in plant growth, development, and response to biotic and abiotic stresses. Most of the miRNAs that have been identified in model plants are well characterized, but till now, no reports have previously been published concerning miRNAs in Paulownia australis. In order to investigate miRNA-guided transcript target regulation in P. australis, small RNA libraries from two P. australis (diploids, PA2; and autotetraploids, PA4) genotypes were subjected to Solexa sequencing. As a result, 10,691,271 (PA2) and 10,712,733 (PA4) clean reads were obtained, and 45 conserved miRNAs belonging to 15 families, and 31 potential novel miRNAs candidates were identified. Compared with their expression levels in the PA2 plants, 26 miRNAs were up-regulated and 15 miRNAs were down-regulated in the PA4 plants. The relative expressions of 12 miRNAs were validated by quantitative real-time polymerase chain reaction. Using the degradome approach, 53 transcript targets were identified and annotated based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. These targets were associated with development, stimulus response, metabolism, signaling transduction and biological regulation. Among them, 11 targets, including TCP transcription factors, auxin response factors, squamosa promoter-binding-like proteins, scarecrow-like proteins, L-type lectin-domain containing receptor kinases and zinc finger CCCH domain-containing protein, cleaved by four known miRNA family and two potentially novel miRNAs were found to be involved in regulating plant development, biotic and abiotic stresses. The findings will be helpful to facilitate studies on the functions of miRNAs and their transcript targets in Paulownia.

Role of microRNAs in aluminum stress in plants

Aluminum (Al) stress is a major factor limiting crop production. The primary symptom of Al toxicity is to inhibit root growth. Plant responses to Al require precise regulation of gene expression at transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) are 20-23 nucleotides length non-coding RNAs, which promote the cleavage of target mRNAs. We have summarized some Al-responsive miRNAs identified, especially proposed the regulatory roles of miR319, miR390, miR393, miR319a.2, and miR398 in Al stress signaling network. The cross-talk between miRNAs and signaling pathways also has been discussed.

Bioinformatic identification and experimental validation of miRNAs from foxtail millet (Setaria italica)

MiRNAs are a novel group of non-coding small RNAs that negatively regulate gene expression. Many miRNAs have been identified and investigated extensively in plant species with sequenced genomes. However, few miRNAs have been identified in foxtail millet (Setaria italica), which is an ancient cereal crop of great importance for dry land agriculture. In this study, 271 foxtail millet miRNAs belonging to 44 families were identified using a bioinformatics approach. Twenty-three pairs of sense/antisense miRNAs belonging to 13 families, and 18 miRNA clusters containing members of 8 families were discovered in foxtail millet. We identified 432 potential targets for 38 miRNA families, most of which were predicted to be involved in plant development, signal transduction, metabolic pathways, disease resistance, and environmental stress responses. Gene ontology (GO) analysis revealed that 101, 56, and 23 target genes were involved in molecular functions, biological processes, and cellular components, respectively. We investigated the expression patterns of 43 selected miRNAs using qRT-PCR analysis. All of the miRNAs were expressed ubiquitously with many exhibiting different expression levels in different tissues. We validated five predicted targets of four miRNAs using the RNA ligase mediated rapid amplification of cDNA end (5'-RLM-RACE) method.

Identification and comparative analysis of the Tegillarca granosa haemocytes microRNA transcriptome in response to Cd using a deep sequencing approach

Background

MicroRNAs (miRNAs) are endogenous non-coding small RNAs (sRNAs) that can base pair with their target mRNAs, which represses their translation or induces their degradation in various biological processes. To identify miRNAs regulated by heavy metal stress, we constructed two sRNA libraries for the blood clam Tegillarca granosa: one for organisms exposed to toxic levels of cadmium (Cd) and one for a control group.

Results

Sequencing of the two libraries and subsequent analysis revealed 215 conserved and 39 new miRNAs. Most of the new miRNAs in T. granosa were up- or down-regulated in response to Cd exposure. There were significant differences in expression between the Cd and control groups for 16 miRNAs. Of these, five miRNAs were significantly up-regulated and 11 were significantly down-regulated in the Cd stress library. Potential targets were predicted for the 16 differential miRNAs in pre-miRNAs identified according to sequence homology. Some of the predicted miRNA targets are associated with regulation of the response to stress induced by heavy metals. Five differentially expressed miRNAs (Tgr-nmiR-8, Tgr-nmiR-21, Tgr-miR-2a, Tgr-miR-10a-5p, and Tgr-miR-184b) were validated by qRT-PCR.

Conclusion

Our study is the first large-scale identification of miRNAs in T. granosa haemocytes. Our findings suggest that some miRNAs and their target genes and pathways may play critical roles in the responses of this species to environmental heavy metal stresses.