Ethylene-responsive miRNAs in roots of Medicago truncatula identified by high-throughput sequencing at whole genome level

Integrative Transcriptome, miRNAs, Degradome, and Phytohormone Analysis of Brassica rapa L. in Response to Plasmodiophora brassicae

Clubroot is an infectious root disease caused by Plasmodiophora brassicae in Brassica crops, which can cause immeasurable losses. We analyzed integrative transcriptome, small RNAs, degradome, and phytohormone comprehensively to explore the infection mechanism of P. brassicae. In this study, root samples of Brassica rapa resistant line material BrT24 (R-line) and susceptible line material Y510-9 (S-line) were collected at four different time points for cytological, transcriptome, miRNA, and degradome analyses. We found the critical period of disease resistance and infection were at 0-3 DAI (days after inoculation) and 9-20 DAI, respectively. Based on our finding, we further analyzed the data of 9 DAI vs. 20 DAI of S-line and predicted the key genes ARF8, NAC1, NAC4, TCP10, SPL14, REV, and AtHB, which were related to clubroot disease development and regulating disease resistance mechanisms. These genes are mainly related to auxin, cytokinin, jasmonic acid, and ethylene cycles. We proposed a regulatory model of plant hormones under the mRNA-miRNA regulation in the critical period of P. brassicae infection by using the present data of the integrative transcriptome, small RNAs, degradome, and phytohormone with our previously published results. Our integrative analysis provided new insights into the regulation relationship of miRNAs and plant hormones during the process of disease infection with P. brassicae.

Glutathione imparts stress tolerance against Alternaria brassicicola infection via miRNA mediated gene regulation

Glutathione (GSH) is well known to play a crucial role in imparting resistance against various pathogen invasions. Nevertheless, the role of GSH in regulating miRNA-mediated defense response is yet to be explored. To decipher the GSH-mediated regulation of miRNA expression during necrotrophic infection in Arabidopsis thaliana, wild-type Col-0 and AtECS1, the transgenic line exhibiting enhanced GSH content, were infected with necrotrophic pathogen Alternaria brassicicola. AtECS1 plants exhibited enhanced resistance as compared to wild-type. MiRNA next-generation sequencing (NGS) was performed to compare the miRNA expression in Col-0 and AtECS1 leaves. Under control condition, differentially expressed 96 known miRNAs and 17 novel miRNAs viz. ath-miR8167f, ath-miR1886.3, ath-miR3932b-5p, etc. were identified. However, under infected condition, 73 known and 43 novel differentially expressed miRNAs viz. ath-miR5652, ath-miR160b, ath-miR865-5p, etc. were identified. Functional annotation and enrichment analysis revealed that several miRNAs that target defense-related genes like leucine-rich repeat protein kinase, MYB transcription factors, TCP8, etc. were down regulated in the AtECS1 line, which, in turn, relieves the repression of their target gene expression, leading to resistance against infection. Together, the present investigation suggests that GSH plays a decisive role in modulating the miRNA-mediated regulation of defense-related genes during pathogen invasion.

Identification of miRNAs Mediating Seed Storability of Maize during Germination Stage by High-Throughput Sequencing, Transcriptome and Degradome Sequencing

Seed storability is an important trait for improving grain quality and germplasm conservation, but little is known about the regulatory mechanisms and gene networks involved. MicroRNAs (miRNAs) are small non-coding RNAs regulating the translation and accumulation of their target mRNAs by means of sequence complementarity and have recently emerged as critical regulators of seed germination. Here, we used the germinating embryos of two maize inbred lines with significant differences in seed storability to identify the miRNAs and target genes involved. We identified a total of 218 previously known and 448 novel miRNAs by miRNA sequencing and degradome analysis, of which 27 known and 11 newly predicted miRNAs are differentially expressed in two maize inbred lines, as measured by Gene Ontology (GO) enrichment analysis. We then combined transcriptome sequencing and real-time quantitative polymerase chain reaction (RT-PCR) to screen and confirm six pairs of differentially expressed miRNAs associated with seed storability, along with their negative regulatory target genes. The enrichment analysis suggested that the miRNAs/target gene mediation of seed storability occurs via the ethylene activation signaling pathway, hormone synthesis and signal transduction, as well as plant organ morphogenesis. Our results should help elucidate the mechanisms through which miRNAs are involved in seed storability in maize.

Evolutionary and functional analysis of two-component system in chickpea reveals CaRR13, a TypeB RR, as positive regulator of symbiosis

The critical role of cytokinin in early nodulation in legumes is well known. In our study, exogenous cytokinin application to roots of the important crop legume, chickpea (Cicer arietinum L.), led to the formation of pseudo-nodules even in the absence of rhizobia. Hence, a genome-wide analysis of the cytokinin signalling, two-component system (TCS) genes, was conducted in chickpea, Medicago and Cajanus cajan. The integrated phylogenetic, evolutionary and expression analysis of the TCS genes was carried out, which revealed that histidine kinases (HKs) were highly conserved, whereas there was diversification leading to neofunctionalization at the level of response regulators (RRs) especially the TypeB RRs. Further, the functional role of the CaHKs in nodulation was established by complementation of the sln1Δ mutant of yeast and cre1 mutants of (Medicago) which led to restoration of the nodule-deficient phenotype. Additionally, the highest expressing TypeB RR of chickpea, CaRR13, was functionally characterized. Its localization in the nucleus and its Y1H assay-based interaction with the promoter of the early nodulation gene CaNSP2 indicated its role as a transcription factor regulating early nodulation. Overexpression, RNAi lines and complementation of cre1 mutants with CaRR13 revealed its critical involvement as an important signalling molecule regulating early events of nodule organogenesis in chickpea.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

Identification and characterization of heat-responsive microRNAs at the booting stage in two rice varieties, 9311 and Nagina 22

MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs that play important roles in abiotic stress responses in plants, but their regulatory roles in the adaptive response to heat stress at the booting stage in two rice varieties, 9311 and Nagina 22, remain largely unknown. In this study, 464 known miRNAs and 123 potential novel miRNAs were identified. Of these miRNAs, a total of 90 differentially expressed miRNAs were obtained with 9311 libraries as the control group, of which 54 were upregulated and 36 were downregulated. To gain insight into functional significance, 2773 potential target genes of these 90 differentially expressed miRNAs were predicted. GO enrichment analysis showed that the predicted target genes of differentially expressed miRNAs included NACs, LACs, CSD, and Hsp40. KEGG pathway analysis showed that the target genes of these differentially expressed miRNAs were significantly enriched in the plant hormone signal transduction pathway. The expression levels of 10 differentially expressed miRNAs and their target genes obtained by qRT-PCR were largely consistent with the sequencing results. This study lays a foundation for the elucidation of the miRNA-mediated regulatory mechanisms in rice at elevated temperatures.

A brief appraisal of ethylene signaling under abiotic stress in plants

For years, ethylene has been known to humankind as the plant hormone responsible for fruit ripening. However, the multitasking aspect of ethylene is still being investigated as ever. It is one of the most diversified signaling molecules which acclimatize plant under adverse conditions. It promotes adventitious root formation, stem and petiole elongation, opening and closing of stomatal aperture, reduces salinity and metal stress, etc. Presence of ethylene checks the production and scavenging of reactive oxygen species by strengthening the antioxidant machinery. Meanwhile, it interacts with other signaling molecules and initiates a cascade of adaptive responses. In the present mini review, the biosynthesis and sources of ethylene production, interaction with other signaling molecules, and its exogenous application under different abiotic stresses have been discussed.

Emerging Connections between Small RNAs and Phytohormones

Small RNAs (sRNAs), mainly including miRNAs and siRNAs, are ubiquitous in eukaryotes. sRNAs mostly negatively regulate gene expression via (post-)transcriptional gene silencing through DNA methylation, mRNA cleavage, or translation inhibition. The mechanisms of sRNA biogenesis and function in diverse biological processes, as well as the interactions between sRNAs and environmental factors, like (a)biotic stress, have been deeply explored. Phytohormones are central in the plant's response to stress, and multiple recent studies highlight an emerging role for sRNAs in the direct response to, or the regulation of, plant hormonal pathways. In this review, we discuss recent progress on the unraveling of crossregulation between sRNAs and nine plant hormones.

Small RNA sequencing revealed various microRNAs involved in ethylene-triggered flowering process in Aechmea fasciata

Ethylene-triggered flowering is a common phenomenon in plants of the family Bromeliaceae, but its molecular mechanism remains unclear. As a classical group of small RNAs, microRNAs play an essential role in the regulation of flowering. In this study, we found that various miRNAs participate in the ethylene-triggered flowering process in Aechmea fasciata via small RNA sequencing using juvenile and adult plants treated with ethylene for 24 hours. Finally, 63 known miRNAs, 52 novel miRNAs and 1721 target genes were identified or predicted. Expression changes of specific miRNAs were validated by qRT-PCR and northern blotting. Some predicted targets, including SPL, GAMYB and ARF, were verified in RLM-RACE experiments. Gene Ontology (GO) and KEGG analysis showed that numerous developmental and RNA-related processes were enriched. Integrated analysis of the transcriptomic data with small RNA sequencing revealed that numerous miRNAs and targets involved in ethylene-triggered flowering in A. fasciata. Our study is helpful for illuminating the molecular basis of the ethylene-triggered flowering phenomenon in Bromeliaceae.

Analyses of microRNA166 gene structure, expression, and function during the early stage of somatic embryogenesis in Dimocarpus longan Lour

MicroRNA166 (miR166) contributes to post-transcriptional regulation by binding the mRNAs of HD-ZIP III genes, which affects plant growth and development. The structural characteristics, expression, and functions of miR166 genes during the early somatic embryogenesis stage in Dimocarpus longan remain unknown. We isolated the transcripts of pri-miR166 S78 with two transcription initiation sites (TSSs) and pri-miR166 S338 with one TSS. These sequences contain potential smORFs and encode different miRNA peptides (miPEPs). Additionally, their promoters contain cis-acting elements responsive to diverse stimuli. The pre-miR166 S78 and pre-miR166 S338 expression levels were up-regulated in response to 2,4-D, abscisic acid, and ethylene. Although the expression patterns induced by hormones were similar, there were differences in the extent of the response, with pre-miR166 S338 more responsive than pre-miR166 S78. Thus, miRNA transcription and maturation are not simply linearly correlated. Moreover, pre-miR166 S78 and pre-miR166 S338 expression levels were down-regulated, whereas ATHB15 (target gene) expression was up-regulated, from the longan embryonic callus to the globular embryo stages. These results are indicative of a negative regulatory relationship between miR166 and ATHB15 during the early somatic embryogenesis stage in longan. At the same stages, miR166a.2-agomir, miR166a.2-antagomir, and miPEP166 S338 increased or decreased the expression of miR166a.2 and ATHB15, but with no consistent patterns or linear synchronization, from which we've found some reasons for it.

MiR319 mediated salt tolerance by ethylene

Salinity-induced accumulation of certain microRNAs accompanied by gaseous phytohormone ethylene production has been recognized as a mechanism of plant salt tolerance. MicroRNA319 (miR319) has been characterized as an important player in abiotic stress resistance in some C3 plants, such as Arabidopsis thaliana and rice. However, its role in the dedicated biomass plant switchgrass (Panicum virgatum L.), a C4 plant, has not been reported. Here, we show crosstalk between miR319 and ethylene (ET) for increasing salt tolerance. By overexpressing Osa-MIR319b and a target mimicry form of miR319 (MIM319), we showed that miR319 positively regulated ET synthesis and salt tolerance in switchgrass. By experimental treatments, we demonstrated that ET-mediated salt tolerance in switchgrass was dose-dependent, and miR319 regulated the switchgrass salt response by fine-tuning ET synthesis. Further experiments showed that the repression of a miR319 target, PvPCF5, in switchgrass also led to enhanced ethylene accumulation and salt tolerance in transgenic plants. Genome-wide transcriptome analysis demonstrated that overexpression of miR319 (OE-miR319) down-regulated the expression of key genes in the methionine (Met) cycle but promoted the expression of genes in ethylene synthesis. The results enrich our understanding of the synergistic effects of the miR319-PvPCF5 module and ethylene synthesis in the salt tolerance of switchgrass, a C4 bioenergy plant.

The miR396-GRF Regulatory Module Controls the Embryogenic Response in Arabidopsis via an Auxin-Related Pathway

In plants, microRNAs have been indicated to control various developmental processes, including somatic embryogenesis (SE), which is triggered in the in vitro cultured somatic cells of plants. Although a transcriptomic analysis has indicated that numerous MIRNAs are differentially expressed in the SE of different plants, the role of specific miRNAs in the embryogenic reprogramming of the somatic cell transcriptome is still poorly understood. In this study, we focused on performing a functional analysis of miR396 in SE given that the transcripts of MIR396 genes and the mature molecules of miR396 were found to be increased during an SE culture of Arabidopsis. In terms of miR396 in embryogenic induction, we observed the SE-associated expression pattern of MIR396b in explants of the β-glucuronidase (GUS) reporter line. In order to gain insight into the miR396-controlled mechanism that is involved in SE induction, the embryogenic response of mir396 mutants and the 35S:MIR396b overexpressor line to media with different 2,4-Dichlorophenoxyacetic acid (2,4-D) concentrations was evaluated. The results suggested that miR396 might contribute to SE induction by controlling the sensitivity of tissues to auxin treatment. Within the targets of miR396 that are associated with SE induction, we identified genes encoding the GROWTH-REGULATING FACTOR (GRF) transcription factors, including GRF1, GRF4, GRF7, GRF8, and GRF9. Moreover, the study suggested a regulatory relationship between miR396, GRF, and the PLETHORA (PLT1 and PLT2) genes during SE induction. A complex regulatory relationship within the miR396–GRF1/4/8/9–PLT1/2 module that involves the negative and positive control of GRFs and PLT (respectively) by miR396 might be assumed.

Identification of Ethylene Responsive miRNAs and Their Targets from Newly Harvested Banana Fruits Using High-Throughput Sequencing

The roles of microRNAs (miRNAs) related to ethylene response in banana fruits remain unknown because many miRNAs are differentially expressed as the fruit ripens, making the identification of ethylene-responsive miRNAs difficult. Using newly harvested banana fruits (within 5 h after harvest) as material, we found that these fruit did not ripen when treated with 5 μL/L of ethylene for 12 h at 22 °C. Two miRNA libraries were generated from newly harvested banana fruits with and without ethylene treatment and sequenced. In total, 128 known miRNAs belonging to 42 miRNA families were obtained, and 12 novel miRNAs were identified. Among them, 22 were differentially expressed in response to ethylene treatment, among which 6 known miRNAs and their putative targets were validated using qRT-PCR. These putative targets encoded proteins including GATA, ARF, DLC, and AGO, etc. KEGG and GO analyses showed that miRNAs differentially expressed in response to ethylene mainly function in the molecular and biological processes.

The MicroRNA319d/TCP10 Node Regulates the Common Bean - Rhizobia Nitrogen-Fixing Symbiosis

Micro-RNAs from legume plants are emerging as relevant regulators of the rhizobia nitrogen-fixing symbiosis. In this work we functionally characterized the role of the node conformed by micro-RNA319 (miR319) - TEOSINTE BRANCHED/CYCLOIDEA/PCF (TCP) transcription factor in the common bean (Phaseolus vulgaris) - Rhizobium tropici symbiosis. The miR319d, one of nine miR319 isoforms from common bean, was highly expressed in root and nodules from inoculated plants as compared to roots from fertilized plants. The miR319d targets TCP10 (Phvul.005G067950), identified by degradome analysis, whose expression showed a negative correlation with miR319d expression. The phenotypic analysis of R. tropici-inoculated composite plants with transgenic roots/nodules overexpressing or silencing the function of miR319d demonstrated the relevant role of the miR319d/TCP10 node in the common bean rhizobia symbiosis. Increased miR319d resulted in reduced root length/width ratio, increased rhizobial infection evidenced by more deformed root hairs and infection threads, and decreased nodule formation and nitrogenase activity per plant. In addition, these plants with lower TCP10 levels showed decreased expression level of the jasmonic acid (JA) biosynthetic gene: LOX2. The transcription of LOX2 by TCPs has been demonstrated for Arabidopsis and in several plants LOX2 level and JA content have been associate with TCP levels. On this basis, we propose that in roots/nodules of inoculated common bean plants TCP10 could be the transcriptional regulator of LOX2 and the miR319d/TCP10 node could affect nodulation through JA signaling. However, given the complexity of nodulation, the participation of other signaling pathways in the phenotypes observed cannot be ruled out.

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.

MicroRNAs Are Intensively Regulated during Induction of Somatic Embryogenesis in Arabidopsis

Several genes encoding transcription factors (TFs) were indicated to have a key role in the induction of somatic embryogenesis (SE), which is triggered in the somatic cells of plants. In order to further explore the genetic regulatory network that is involved in the embryogenic transition induced in plant somatic cells, micro-RNA (miRNAs) molecules, the products of MIRNA (MIR) genes and the common regulators of TF transcripts, were analyzed in an embryogenic culture of Arabidopsis thaliana. In total, the expression of 190 genes of the 114 MIRNA families was monitored during SE induction and the levels of the primary (pri-miRNAs) transcripts vs. the mature miRNAs were investigated. The results revealed that the majority (98%) of the MIR genes were active and that most of them (64%) were differentially expressed during SE. A distinct attribute of the MIR expression in SE was the strong repression of MIR transcripts at the early stage of SE followed by their significant up-regulation in the advanced stage of SE. Comparison of the mature miRNAs vs. pri-miRNAs suggested that the extensive post-transcriptional regulation of miRNA is associated with SE induction. Candidate miRNA molecules of the assumed function in the embryogenic response were identified among the mature miRNAs that had a differential expression in SE, including miR156, miR157, miR159, miR160, miR164, miR166, miR169, miR319, miR390, miR393, miR396, and miR398. Consistent with the central role of phytohormones and stress factors in SE induction, the functions of the candidate miRNAs were annotated to phytohormone and stress responses. To confirm the functions of the candidate miRNAs in SE, the expression patterns of the mature miRNAs and their presumed targets were compared and regulatory relation during SE was indicated for most of the analyzed miRNA-target pairs. The results of the study contribute to the refinement of the miRNA-controlled regulatory pathways that operate during embryogenic induction in plants and provide a valuable platform for the identification of the genes that are targeted by the candidate miRNAs in SE induction.

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.

Small RNA and degradome sequencing reveals microRNAs and their targets involved in tomato pedicel abscission

We constructed small RNA and degradome sequencing libraries to identify miRNAs and targets involved in tomato pedicel abscission, and confirmed their roles via quantitative real-time PCR. MicroRNAs (miRNAs) are endogenous small RNAs which play crucial negatively regulatory roles at both the transcriptional and post-transcriptional levels in plants; however, limited knowledge is available on the expression profiles of miRNAs and their target genes during tomato pedicel abscission. Taking advantage of small RNA (sRNA) and degradome sequencing technology, a total of 56 known and 11 novel candidate miRNAs targeting 223 mRNA genes were confirmed during pedicel abscission. Gene ontology annotation and KEGG pathway analysis showed that these target genes were significantly enriched in intracellular, membrane-bounded organelle-related biological processes as well as in metabolic, plant-pathogen interaction and hormone signaling pathways. We screened 17 miRNA/target pairs for further analysis and performed quantitative real-time PCR to identify the roles. Cluster analysis of selected miRNAs revealed that the expression profiles of miRNAs varied in different stages of abscission and could be impacted by ethylene treatment. In the present study, the correlations between miRNAs and targets suggested a complex regulatory network of miRNA-mediated target interaction during pedicel abscission. Additionally, the expression profiles of miRNAs and their targets changed by ethylene might be a considerable reason why ethylene promotes pedicel abscission. Our study provides new insights into the expression and regulatory profiles of miRNAs during tomato pedicel abscission.

Identification of nitrogen starvation-responsive microRNAs in Chrysanthemum nankingense

MicroRNA (miRNA) is involved in many developmental processes and various abiotic stress responses in plants. As nitrogen is a limited element for plant growth, comparative analyses of miRNAs responding to low nitrogen stress is important for improving the nitrogen use efficiency (NUE). We used high-throughput sequencing to detect the response of miRNAs to low nitrogen stress in the roots and leaves of Chrysanthemum nankingense. Compared with the control, the differential expression was more than 2-fold in 81 miRNAs in roots and 101 miRNAs in leaves. The identified miRNAs showed overlapping or unique response to nitrate limitation in roots and leaves, including several members of known miRNA families with low nitrogen stress response, such as miR156, miR169, and miR393. The potential target genes of these miRNAs were also identified. The total amount of predicted target genes was 219, and the corresponding amount of matched miRNAs was 37 in roots and 44 in leaves. Moreover, we used 5' RLM-RACE to map the cleavage sites in four predicted target genes. The differential expression level of miRNAs and target genes was verified by quantitative real-time polymerase chain reaction (qRT-PCR). According to the functional characteristics of the predicted target genes, they were divided into three main categories: transcription factors, kinases, and metabolism.

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.

Involvement of ethylene-responsive microRNAs and their targets in increased latex yield in the rubber tree in response to ethylene treatment

The rubber tree is an economically important plant that produces natural rubber for various industrial uses. The application of ethylene contributes to increased latex production in rubber trees; however, the molecular biology behind the effects of ethylene on latex yield remains to be elucidated. Recently, the intersection between microRNA (miRNA) regulation and phytohormone responses has been revealed. Insight into the regulation of miRNAs and their target genes should help to determine the functional importance of miRNAs as well as the role of miRNAs in signaling under ethylene stimulation in the rubber tree. In this study, hbr-miR159 and hbr-miR166 were down-regulated in bark under ethylene treatment. The ethylene also down-regulated ATHB15-like (Class III Homeodomain Leucine Zipper, HD-ZIP III) which have been extensively implicated in the regulation of primary and secondary vascular tissue pattern formation. The strong negative-regulation of ARF6/ARF8 caused by hbr-miR167 involved in an attenuation of vascular development and may gradually lead to bark dryness syndrome in the long term ethylene treatment. The negative correlation of hbr-miR172 and its target REF3 in the inner soft bark under ethylene treatment results in dramatic increases in latex yield in the ethylene-sensitive clone of the rubber tree. The overall results suggested that the differential expression of HD-ZIP III, miR167/ARF6, ARF8, and miR172/REF3 and related genes may play possible roles in the response to ethylene treatment, resulting in longer latex flow and increased latex yield.

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.

Characterization of rubber tree microRNA in phytohormone response using large genomic DNA libraries, promoter sequence and gene expression analysis

The para rubber tree is the most widely cultivated tree species for producing natural rubber (NR) latex. Unfortunately, rubber tree characteristics such as a long life cycle, heterozygous genetic backgrounds, and poorly understood genetic profiles are the obstacles to breeding new rubber tree varieties, such as those with improved NR yields. Recent evidence has revealed the potential importance of controlling microRNA (miRNA) decay in some aspects of NR regulation. To gain a better understanding of miRNAs and their relationship with rubber tree gene regulation networks, large genomic DNA insert-containing libraries were generated to complement the incomplete draft genome sequence and applied as a new powerful tool to predict a function of interested genes. Bacterial artificial chromosome and fosmid libraries, containing a total of 120,576 clones with an average insert size of 43.35 kb, provided approximately 2.42 haploid genome equivalents of coverage based on the estimated 2.15 gb rubber tree genome. Based on these library sequences, the precursors of 1 member of rubber tree-specific miRNAs and 12 members of conserved miRNAs were successfully identified. A panel of miRNAs was characterized for phytohormone response by precisely identifying phytohormone-responsive motifs in their promoter sequences. Furthermore, the quantitative real-time PCR on ethylene stimulation of rubber trees was performed to demonstrate that the miR2118, miR159, miR164 and miR166 are responsive to ethylene, thus confirmed the prediction by genomic DNA analysis. The cis-regulatory elements identified in the promoter regions of these miRNA genes help augment our understanding of miRNA gene regulation and provide a foundation for further investigation of the regulation of rubber tree miRNAs.

Epigenetic control of mobile DNA as an interface between experience and genome change

Mobile DNA in the genome is subject to RNA-targeted epigenetic control. This control regulates the activity of transposons, retrotransposons and genomic proviruses. Many different life history experiences alter the activities of mobile DNA and the expression of genetic loci regulated by nearby insertions. The same experiences induce alterations in epigenetic formatting and lead to trans-generational modifications of genome expression and stability. These observations lead to the hypothesis that epigenetic formatting directed by non-coding RNA provides a molecular interface between life history events and genome alteration.

miRNAs in the crosstalk between phytohormone signalling pathways

Phytohormones are signal molecules produced within the plant that control its growth and development through the regulation of gene expression. Interaction between different phytohormone pathways is essential in coordinating tissue outgrowth in response to environmental changes, such as the adaptation of root development to water deficit or the initiation of seed germination during imbibition. Recently, microRNAs (miRNAs) have emerged as key regulators of phytohormone response pathways in planta by affecting their metabolism, distribution, and perception. Here we review current knowledge on the miRNA-mediated regulations involved in phytohormone crosstalk. We focus on the miRNAs exhibiting regulatory links with more than one phytohormone pathway and discuss their possible implication in coordinating multiple phytohormone responses during specific developmental processes.

Genome-wide identification of microRNAs in Medicago truncatula by high-throughput sequencing

MicroRNAs (miRNAs) are small, endogenous RNAs that play important regulatory roles in development and stress response in plants by negatively regulating gene expression post-transcriptionally. Medicago truncatula has been used as a model plant to study functional genomics of legume plants. It has also been widely used to functionally study miRNAs. Identification of miRNAs at the whole-genome level is essential for functional characterization of miRNAs in plants. High-throughput sequencing is a powerful technology to identify miRNAs. In this chapter, the methods used for construction of a small RNA library and high-throughput sequencing involving total RNA isolation, small RNA purification, adapter ligation, reverse transcription, PCR amplification, and Solexa sequencing are described. Bioinformatics and analysis of differential expression of miRNAs including primary disposal, miRNA identification, target prediction, and expression analysis are also discussed. These methodologies associated with identification and functional characterization of miRNAs may provide useful tools for readers to study miRNAs in plants in general and Medicago truncatula in particular.

Transcriptomics of the interaction between the monopartite phloem-limited geminivirus tomato yellow leaf curl Sardinia virus and Solanum lycopersicum highlights a role for plant hormones, autophagy and plant immune system fine tuning during infection

Tomato yellow leaf curl Sardinia virus (TYLCSV), a DNA virus belonging to the genus Begomovirus, causes severe losses in tomato crops. It infects only a limited number of cells in the vascular tissues, making difficult to detect changes in host gene expression linked to its presence. Here we present the first microarray study of transcriptional changes induced by the phloem-limited geminivirus TYLCSV infecting tomato, its natural host. The analysis was performed on the midrib of mature leaves, a material naturally enriched in vascular tissues. A total of 2206 genes were up-regulated and 1398 were down-regulated in infected plants, with an overrepresentation of genes involved in hormone metabolism and responses, nucleic acid metabolism, regulation of transcription, ubiquitin-proteasome pathway and autophagy among those up-regulated, and in primary and secondary metabolism, phosphorylation, transcription and methylation-dependent chromatin silencing among those down-regulated. Our analysis showed a series of responses, such as the induction of GA- and ABA-responsive genes, the activation of the autophagic process and the fine tuning of the plant immune system, observed only in TYLCSV-tomato compatible interaction so far. On the other hand, comparisons with transcriptional changes observed in other geminivirus-plant interactions highlighted common host responses consisting in the deregulation of biotic stress responsive genes, key enzymes in the ethylene biosynthesis and methylation cycle, components of the ubiquitin proteasome system and DNA polymerases II. The involvement of conserved miRNAs and of solanaceous- and tomato-specific miRNAs in geminivirus infection, investigated by integrating differential gene expression data with miRNA targeting data, is discussed.

Identification of MiRNA from eggplant (Solanum melongena L.) by small RNA deep sequencing and their response to Verticillium dahliae infection

MiRNAs are a class of non-coding small RNAs that play important roles in the regulation of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data, including eggplant (Solanum melongena L.). To identify miRNAs in eggplant and their response to Verticillium dahliae infection, a fungal pathogen for which clear understanding of infection mechanisms and effective cure methods are currently lacking, we deep-sequenced two small RNA (sRNA) libraries prepared from mock-infected and infected seedlings of eggplants. Specifically, 30,830,792 reads produced 7,716,328 unique miRNAs representing 99 known miRNA families that have been identified in other plant species. Two novel putative miRNAs were predicted with eggplant ESTs. The potential targets of the identified known and novel miRNAs were also predicted based on sequence homology search. It was observed that the length distribution of obtained sRNAs and the expression of 6 miRNA families were obviously different between the two libraries. These results provide a framework for further analysis of miRNAs and their role in regulating plant response to fungal infection and Verticillium wilt in particular.

Integrative analysis of miRNA and mRNA profiles in response to ethylene in rose petals during flower opening

MicroRNAs play an important role in plant development and plant responses to various biotic and abiotic stimuli. As one of the most important ornamental crops, rose (Rosa hybrida) possesses several specific morphological and physiological features, including recurrent flowering, highly divergent flower shapes, colors and volatiles. Ethylene plays an important role in regulating petal cell expansion during rose flower opening. Here, we report the population and expression profiles of miRNAs in rose petals during flower opening and in response to ethylene based on high throughput sequencing. We identified a total of 33 conserved miRNAs, as well as 47 putative novel miRNAs were identified from rose petals. The conserved and novel targets to those miRNAs were predicted using the rose floral transcriptome database. Expression profiling revealed that expression of 28 known (84.8% of known miRNAs) and 39 novel (83.0% of novel miRNAs) miRNAs was substantially changed in rose petals during the earlier opening period. We also found that 28 known and 22 novel miRNAs showed expression changes in response to ethylene treatment. Furthermore, we performed integrative analysis of expression profiles of miRNAs and their targets. We found that ethylene-caused expression changes of five miRNAs (miR156, miR164, miR166, miR5139 and rhy-miRC1) were inversely correlated to those of their seven target genes. These results indicate that these miRNA/target modules might be regulated by ethylene and were involved in ethylene-regulated petal growth.

Stable transformation of Medicago truncatula cv. Jemalong for gene analysis using Agrobacterium tumefaciens

Medicago truncatula is a model legume that has all the genomic resources to be suitable as a model for functional genomics. Transformation to produce transgenic plants is part of this toolkit, enabling a spectrum of approaches to study gene function: including knockdown, overexpression, reporter genes for gene expression, and proteins tagged with fluorescent proteins such as GFP. A special genetic line is necessary for transformation and Jemalong 2HA derived from cv. Jemalong is used in the methods described. Leaf explants can be used for the transformation of the embryonic stem cells to produce the transgenic somatic embryos for regeneration. An auxin and a cytokinin are the key hormone requirements for regeneration by somatic embryogenesis but other hormones such as abscisic acid can be used to augment the system. As the explants used in this system are from leaves, rather than immature embryos or meristematic tissues often used in other species, it is a quite straightforward system. Agrobacterium tumefaciens containing a binary vector suitable for the particular objectives is used to deliver the transgene of interest.