Stars and symbiosis: microRNA- and microRNA*-mediated transcript cleavage involved in arbuscular mycorrhizal symbiosis

A dicistronic precursor encoding miR398 and the legume-specific miR2119 coregulates CSD1 and ADH1 mRNAs in response to water deficit

Plant microRNAs are commonly encoded in transcripts containing a single microRNA precursor. Processing by DICER-LIKE 1 and associated factors results in the production of a small RNA, followed by its incorporation into an AGO-containing protein complex to guide silencing of an mRNA possessing a complementary target sequence. Certain microRNA loci contain more than one precursor stem-loop structure, thus encoding more than one microRNA in the same transcript. Here, we describe a unique case where the evolutionary conserved miR398a is encoded in the same transcript as the legume-specific miR2119. The dicistronic arrangement found in common bean was also observed in other legumes. In Phaseolus vulgaris, mature miR398 and miR2119 are repressed in response to water deficit, and we demonstrate that both are functional as they target the mRNAs for CSD1 and ADH1, respectively. Our results indicate that the repression of miR398 and miR2119 leads to coordinated up-regulation of CSD1 and ADH1 mRNAs in response to water deficit in common bean and possibly in other legumes. Furthermore, we show that miRNA directed CSD1 and ADH1 mRNAs up-regulation also occurs when common bean plants are exposed to flooding, suggesting that plant redox status and fermentation metabolism must be closely coordinated under different adverse conditions.

Identification and characterisation of microRNAs and their target genes in phosphate-starved Nicotiana benthamiana by small RNA deep sequencing and 5'RACE analysis

BACKGROUND

Phosphorus is an important macronutrient that is severely lacking in soils. In plants, specific microRNAs (miRNAs) essential for nutrient management and the regulation of stress responses are responsible for the control of many phosphate starvation responses. Further understanding of conserved and species-specific microRNA species has potential implications for the development of crops tolerant to soils with low phosphate.

RESULTS

This study identified and characterised phosphate starvation-responsive miRNAs in the native Australian tobacco Nicotiana benthamiana. Small RNA libraries were constructed and sequenced from phosphate-starved plant leaves, stems and roots. Twenty-four conserved miRNA families and 36 species-specific miRNAs were identified. The majority of highly phosphate starvation-responsive miRNAs were highly conserved, comprising of members from the miR399, miR827, and miR2111 families. In addition, two miRNA-star species were identified to be phosphate starvation-responsive. A total of seven miRNA targets were confirmed using RLM-5'RACE to be cleaved by five miRNA families, including two confirmed cleavage targets for Nbe-miR399 species, one for Nbe-miR2111, and two for Nbe-miR398. A number of N. benthamiana-specific features for conserved miRNAs were identified, including species-specific miRNA targets predicted or confirmed for miR399, miR827, and miR398.

CONCLUSIONS

Our results give an insight into the phosphate starvation-responsive miRNAs of Nicotiana benthamiana, and indicate that the phosphate starvation response pathways in N. benthamiana contain both highly conserved and species-specific components.

Complex regulation of microRNAs in roots of competitively-grown isogenic Nicotiana attenuata plants with different capacities to interact with arbuscular mycorrhizal fungi

BACKGROUND

Nicotiana attenuata is an ecological model plant whose 2.57 Gb genome has recently been sequenced and assembled and for which miRNAs and their genomic locations have been identified. To understand how this plant's miRNAs are reconfigured during plant-arbuscular mycorrhizal fungal (AMF) interactions and whether hostplant calcium- and calmodulin dependent protein kinase (CCaMK) expression which regulates the AMF interaction also modulates miRNAs levels and regulation, we performed a large-scale miRNA analysis of this plant-AMF interaction.

RESULTS

Next generation sequencing of miRNAs in roots of empty vector (EV) N. attenuata plants and an isogenic line silenced in CCaMK expression (irCCaMK) impaired in AMF-interactions grown under competitive conditions with and without AMF inoculum revealed a total of 149 unique miRNAs: 67 conserved and 82 novel ones. The majority of the miRNAs had a length of 21 nucleotides. MiRNA abundances were highly variable ranging from 400 to more than 25,000 reads per million. The miRNA profile of irCCaMK plants impaired in AMF colonization was distinct from fully AMF-functional EV plants grown in the same pot. Six conserved miRNAs were present in all conditions and accumulated differentially depending on treatment and genotype; five (miR6153, miR403a-3p, miR7122a, miR167-5p and miR482d, but not miR399a-3p) showed the highest accumulation in AMF inoculated EV plants compared to inoculated irCCaMK plants. Furthermore, the accumulation patterns of sequence variants of selected conserved miRNAs showed a very distinct pattern related to AMF colonization - one variant of miR473-5p specifically accumulated in AMF-inoculated plants. Also abundances of miR403a-3p, miR171a-3p and one of the sequence variants of miR172a-3p increased in AMF-inoculated EV compared to inoculated irCCaMK plants and to non-inoculated EV plants, while miR399a-3p was most strongly enriched in AMF inoculated irCCaMK plants grown in competition with EV. The analysis of putative targets of selected miRNAs revealed an involvement in P starvation (miR399), phytohormone signaling (Nat-R-PN59, miR172, miR393) and defense (e.g. miR482, miR8667, Nat-R-PN-47).

CONCLUSIONS

Our study demonstrates (1) a large-scale reprograming of miRNAs induced by AMF colonization and (2) that the impaired AMF signaling due to CCaMK silencing and the resulting reduced competitive ability of irCCaMK plants play a role in modulating signal-dependent miRNA accumulation.

Genome-wide identification and profiling of microRNAs in Paulownia tomentosa cambial tissues in response to seasonal changes

MicroRNAs (miRNAs), a group of endogenous small non-coding RNAs, have been shown to play essential roles in the regulation of gene expression at the post-transcriptional level. Although Paulownia tomentosa is an ecologically and economically important timber species due to its rapid growth, few efforts have focused on small RNAs (sRNAs) in the cambial tissues during winter and summer transition. In the present study, we identified 33 known miRNA families and 29 novel miRNAs which include 20 putative novel miRNAs* in P. tomentosa cambial tissues during winter and summer transition. Through differential expression analysis, we showed that 15 known miRNAs and 8 novel miRNAs were preferentially abundant in certain stage of cambial tissues. Based on the P. tomentosa mRNA transcriptome database, 1667 and 78 potential targets were predicted for 29 known and 20 novel miRNAs, respectively and the predicted targets are mostly transcription factors and functional genes. The targets of these miRNAs were enriched in "metabolic process" and "transcription regulation" by using Gene Ontology enrichment analysis. In addition, KEGG pathway analyses revealed the involvement of miRNAs in starch and sucrose metabolism and plant-pathogen interaction metabolism pathways. Noticeably, qRT-PCR expression analysis demonstrated that 9 miRNAs and their targets were existed a negative correlation in P. tomentosa cambial tissues. This study is the first to examine known and novel miRNAs and their potential targets in P. tomentosa cambial tissues during winter and summer transition and identify several candidate genes potentially regulating cambial phase transition, and thus provide a framework for further understanding of miRNAs functions in the regulation of cambial phase transition and wood formation in trees.

The Polycistronic miR166k-166h Positively Regulates Rice Immunity via Post-transcriptional Control of EIN2

MicroRNAs (miRNAs) are small RNAs acting as regulators of gene expression at the post-transcriptional level. In plants, most miRNAs are generated from independent transcriptional units, and only a few polycistronic miRNAs have been described. miR166 is a conserved miRNA in plants targeting the HD-ZIP III transcription factor genes. Here, we show that a polycistronic miRNA comprising two miR166 family members, miR166k and miR166h, functions as a positive regulator of rice immunity. Rice plants with activated MIR166k-166h expression showed enhanced resistance to infection by the fungal pathogens Magnaporthe oryzae and Fusarium fujikuroi, the causal agents of the rice blast and bakanae disease, respectively. Disease resistance in rice plants with activated MIR166k-166h expression was associated with a stronger expression of defense responses during pathogen infection. Stronger induction of MIR166k-166h expression occurred in resistant but not susceptible rice cultivars. Notably, the ethylene-insensitive 2 (EIN2) gene was identified as a novel target gene for miR166k. The regulatory role of the miR166h-166k polycistron on the newly identified target gene results from the activity of the miR166k-5p specie generated from the miR166k-166h precursor. Collectively, our findings support a role for miR166k-5p in rice immunity by controlling EIN2 expression. Because rice blast is one of the most destructive diseases of cultivated rice worldwide, unraveling miR166k-166h-mediated mechanisms underlying blast resistance could ultimately help in designing appropriate strategies for rice protection.

microRNA regulation in an ancient obligate endosymbiosis

Although many insects are associated with obligate bacterial endosymbionts, the mechanisms by which these host/endosymbiont associations are regulated remain mysterious. While microRNAs (miRNAs) have been recently identified as regulators of host/microbe interactions, including host/pathogen and host/facultative endosymbiont interactions, the role miRNAs may play in mediating host/obligate endosymbiont interactions is virtually unknown. Here, we identified conserved miRNAs that potentially mediate symbiotic interactions between aphids and their obligate endosymbiont, Buchnera aphidicola. Using small RNA sequence data from Myzus persicae and Acyrthosiphon pisum, we annotated 93 M. persicae and 89 A. pisum miRNAs, among which 69 were shared. We found 14 miRNAs that were either highly expressed in aphid bacteriome, the Buchnera-housing tissue, or differentially expressed in bacteriome vs. gut, a non-Buchnera-housing tissue. Strikingly, 10 of these 14 miRNAs have been implicated previously in other host/microbe interaction studies. Investigating the interaction networks of these miRNAs using a custom computational pipeline, we identified 103 miRNA::mRNA interactions shared between M. persicae and A. pisum. Functional annotation of the shared mRNA targets revealed only two over-represented cluster of orthologous group categories: amino acid transport and metabolism, and signal transduction mechanisms. Our work supports a role for miRNAs in mediating host/symbiont interactions between aphids and their obligate endosymbiont Buchnera. In addition, our results highlight the probable importance of signal transduction mechanisms to host/endosymbiont coevolution.

Regulation of low temperature stress in plants by microRNAs

Low temperature is one of the most common environmental stresses that seriously affect the growth and development of plants. However, plants have the plasticity in their defence mechanisms enabling them to tolerate and, sometimes, even survive adverse environmental conditions. MicroRNAs (miRNAs) are small non-coding RNAs, approximately 18-24 nucleotides in length, and are being increasingly recognized as regulators of gene expression at the post-transcriptional level and have the ability to influence a broad range of biological processes. There is growing evidence in the literature that reprogramming of gene expression mediated through miRNAs is a major defence mechanism in plants enabling them to respond to stresses. To date, numerous studies have established the importance of miRNA-based regulation of gene expression under low temperature stress. Individual miRNAs can modulate the expression of multiple mRNA targets, and, therefore, the manipulation of a single miRNA has the potential to affect multiple biological processes. Numerous functional studies have attempted to identify the miRNA-target interactions and have elaborated the role of several miRNAs in cold-stress regulation. This review summarizes the current understanding of miRNA-mediated modulation of the expression of key genes as well as genetic and regulatory pathways, involved in low temperature stress responses in plants.

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.

Evolution of the Symbiosis-Specific GRAS Regulatory Network in Bryophytes

Arbuscular mycorrhiza is one of the most common plant symbiotic interactions observed today. Due to their nearly ubiquitous occurrence and their beneficial impact on both partners it was suggested that this mutualistic interaction was crucial for plants to colonize the terrestrial habitat approximately 500 Ma ago. On the plant side the association is established via the common symbiotic pathway (CSP). This pathway allows the recognition of the fungal symbiotic partner, subsequent signaling to the nucleus, and initiation of the symbiotic program with respect to specific gene expression and cellular re-organization. The downstream part of the CSP is a regulatory network that coordinates the transcription of genes necessary to establish the symbiosis, comprising multiple GRAS transcription factors (TFs). These regulate their own expression as an intricate transcriptional network. Deduced from non-host genome data the loss of genes encoding CSP components coincides with the loss of the interaction itself. Here, we analyzed bryophyte species with special emphasis on the moss Physcomitrella patens, supposed to be a non-host, for the composition of the GRAS regulatory network components. We show lineage specific losses and expansions of several of these factors in bryophytes, potentially coinciding with the proposed host/non-host status of the lineages. We evaluate losses and expansions and infer clade-specific evolution of GRAS TFs.

Salicylic Acid Perturbs sRNA-Gibberellin Regulatory Network in Immune Response of Potato to Potato virus Y Infection

Potato virus Y is the most economically important potato viral pathogen. We aimed at unraveling the roles of small RNAs (sRNAs) in the complex immune signaling network controlling the establishment of tolerant response of potato cv. Désirée to the virus. We constructed a sRNA regulatory network connecting sRNAs and their targets to link sRNA level responses to physiological processes. We discovered an interesting novel sRNAs-gibberellin regulatory circuit being activated as early as 3 days post inoculation (dpi) before viral multiplication can be detected. Two endogenous sRNAs, miR167 and phasiRNA931 were predicted to regulate gibberellin biosynthesis genes GA20-oxidase and GA3-oxidase. The increased expression of phasiRNA931 was also reflected in decreased levels of GA3-oxidase transcripts. Moreover, decreased concentration of gibberellin confirmed this regulation. The functional relation between lower activity of gibberellin signaling and reduced disease severity was previously confirmed in Arabidopsis-virus interaction using knockout mutants. We further showed that this regulation is salicylic acid-dependent as the response of sRNA network was attenuated in salicylic acid-depleted transgenic counterpart NahG-Désirée expressing severe disease symptoms. Besides downregulation of gibberellin signaling, regulation of immune receptor transcripts by miR6022 as well as upregulation of miR164, miR167, miR169, miR171, miR319, miR390, and miR393 in tolerant Désirée, revealed striking similarities to responses observed in mutualistic symbiotic interactions. The intertwining of different regulatory networks revealed, shows how developmental signaling, disease symptom development, and stress signaling can be balanced.

ghr-miR5272a-mediated regulation of GhMKK6 gene transcription contributes to the immune response in cotton

Fusarium wilt is a major biotic stress affecting the productivity of cotton (Gossypium hirsutum). Although mitogen-activated protein kinase (MAPK) cascades play critical roles in plant disease resistance, their intricate regulation under fungal stress remains unclear, especially with regards to microRNA-mediated regulation of MAPK gene expression. In this study, we report that the MAPK kinase gene GhMKK6 and ghr-miR5272a work together in cotton resistance to Fusarium wilt. Silencing GhMKK6 in cotton decreased resistance to F. oxysporum by repressing the expression of known disease-resistance genes. Furthermore, although GhMKK6 played a positive role in disease resistance, excessive GhMKK6 activation caused an excessive hypersensitive response. ghr-miR5272a, a major regulator, prevents this excessive response by regulating GhMKK6 expression. ghr-miR5272a targets the GhMKK6 3'-untranslated region in cotton. Overexpressing miR5272a decreased the expression of GhMKK6 and disease-resistance genes, and increased sensitivity to F. oxysporum, yielding a similar phenotype to GhMKK6-silenced cotton. Overall, these results demonstrate that the ghr-miR5272a-mediated regulation of GhMKK6 expression contributes to the immune response in cotton, and reveal a new feedback loop mechanism in plant disease response.

An expression atlas of miRNAs in Arabidopsis thaliana

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. MiRNA expression often exhibits spatial and temporal specificity. However, genome-wide miRNA expression patterns in different organs during development of Arabidopsis thaliana have not yet been systemically investigated. In this study, we sequenced small RNA libraries generated from 27 different organ/tissue types, which cover the entire life cycle of Arabidopsis. Analysis of the sequencing data revealed that most miRNAs are ubiquitously expressed, whereas a small set of miRNAs display highly specific expression patterns. In addition, different miRNA members within the same family have distinct spatial and temporal expression patterns. Moreover, we found that some miRNAs are produced from different arms of their hairpin precursors at different developmental stages. This work provides new insights into the regulation of miRNA biogenesis and a rich resource for future investigation of miRNA functions in Arabidopsis.

Identification of tomato miRNAs responsive to root colonization by endophytic Pochonia chlamydosporia

The molecular mechanisms active during the endophytic phase of the fungus Pochonia chlamydosporia are still poorly understood. In particular, few data are available on the links between the endophyte and the root response, as modulated by noncoding small RNAs. In this study, we describe the microRNAs (miRNAs) that are differentially expressed (DE) in the roots of tomato, colonized by P. chlamydosporia. A genome-wide NGS expression profiling of small RNAs in roots, either colonized or not by the fungus, showed 26 miRNAs upregulated in inoculated roots. Their predicted target genes are involved in the plant information processing system, which recognizes, percepts, and transmits signals, with higher representations in processes such as apoptosis and plant defense regulation. RNAseq data showed that predicted miRNA target genes were downregulated in tomato roots after 4, 7, 10, and 21 days post P. chlamydosporia inoculation. The differential expression of four miRNAs was further validated using qPCR analysis. The P. chlamydosporia endophytic lifestyle in tomato roots included an intricate network of miRNAs and targets. Data provide a first platform of DE tomato miRNAs after P. chlamydosporia colonization. They indicated that several miRNAs are involved in the host response to the fungus, playing important roles for its recognition as a symbiotic microorganism, allowing endophytism by modulating the host defense reaction. Data also indicated that endophytism affects tRNA fragmentation. This is the first study on miRNAs induced by P. chlamydosporia endophytism and related development regulation effects in Solanum lycopersicum.

Characterization of microRNAs of Beta macrocarpa and their responses to Beet necrotic yellow vein virus infection

Plant microRNAs (miRNAs) are a class of non-coding RNAs that play important roles in plant development, defense, and symptom development. Here, 547 known miRNAs representing 129 miRNA families, and 282 potential novel miRNAs were identified in Beta macrocarpa using small RNA deep sequencing. A phylogenetic analysis was performed, and 8 Beta lineage-specific miRNAs were identified. Through a differential expression analysis, miRNAs associated with Beet necrotic yellow vein virus (BNYVV) infection were identified and confirmed using a microarray analysis and stem-loop RT-qPCR. In total, 103 known miRNAs representing 38 miRNA families, and 45 potential novel miRNAs were differentially regulated, with at least a two-fold change, in BNYVV-infected plants compared with that of the mock-inoculated control. Targets of these differentially expressed miRNAs were also predicted by degradome sequencing. These differentially expressed miRNAs were involved in hormone biosynthesis and signal transduction pathways, and enhanced axillary bud development and plant defenses. This work is the first to describe miRNAs of the plant genus Beta and may offer a reference for miRNA research in other species in the genus. It provides valuable information on the pathogenicity mechanisms of BNYVV.

Long-distance movement of phosphate starvation-responsive microRNAs in Arabidopsis

Plant microRNAs are small RNAs that are important for genetic regulation of processes such as plant development or environmental responses. Specific microRNAs accumulate in the phloem during phosphate starvation, and may act as long-distance signalling molecules. We performed quantitative PCR on Arabidopsis hypocotyl micrograft tissues of wild-type and hen1-6 mutants to assess the mobility of several phosphate starvation-responsive microRNA species. In addition to the previously confirmed mobile species miR399d, the corresponding microRNA* (miR399d*) was identified for the first time as mobile between shoots and roots. Translocation by phosphate-responsive microRNAs miR827 and miR2111a between shoots and roots during phosphate starvation was evident, while their respective microRNA*s were not mobile. The results suggest that long-distance mobility of microRNA species is selective and can occur without the corresponding duplex strand. Movement of miR399d* and root-localised accumulation of miR2111a* opens the potential for persisting microRNA*s to be mobile and functional in novel pathways during phosphate starvation responses.

Genome wide identification of microRNAs involved in fatty acid and lipid metabolism of Brassica napus by small RNA and degradome sequencing

Rapeseed (Brassica napus) is an important cash crop considered as the third largest oil crop worldwide. Rapeseed oil contains various saturation or unsaturation fatty acids, these fatty acids, whose could incorporation with TAG form into lipids stored in seeds play various roles in the metabolic activity. The different fatty acids in B. napus seeds determine oil quality, define if the oil is edible or must be used as industrial material. miRNAs are kind of non-coding sRNAs that could regulate gene expressions through post-transcriptional modification to their target transcripts playing important roles in plant metabolic activities. We employed high-throughput sequencing to identify the miRNAs and their target transcripts involved in fatty acids and lipids metabolism in different development of B. napus seeds. As a result, we identified 826 miRNA sequences, including 523 conserved and 303 newly miRNAs. From the degradome sequencing, we found 589 mRNA could be targeted by 236 miRNAs, it includes 49 novel miRNAs and 187 conserved miRNAs. The miRNA-target couple suggests that bna-5p-163957_18, bna-5p-396192_7, miR9563a-p3, miR9563b-p5, miR838-p3, miR156e-p3, miR159c and miR1134 could target PDP, LACS9, MFPA, ADSL1, ACO32, C0401, GDL73, PlCD6, OLEO3 and WSD1. These target transcripts are involving in acetyl-CoA generate and carbon chain desaturase, regulating the levels of very long chain fatty acids, β-oxidation and lipids transport and metabolism process. At the same, we employed the q-PCR to valid the expression of miRNAs and their target transcripts that involve in fatty acid and lipid metabolism, the result suggested that the miRNA and their transcript expression are negative correlation, which in accord with the expression of miRNA and its target transcript. The study findings suggest that the identified miRNA may play important role in the fatty acids and lipids metabolism in seeds of B. napus.

Unraveling the microRNA of Caragana korshinskii along a precipitation gradient on the Loess Plateau, China, using high-throughput sequencing

Drought remains one of the main factors that negatively affect plant growth and development. Caragana korshinskii is widely distributed on the Loess Plateau, China, where it mediates soil and water loss and helps prevent desertification. However, little is known about the stress response mechanisms of C. korshinskii in water-starved environments. MicroRNAs (miRNAs) have been implicated in the regulation of plant responses to several types of biotic and abiotic stress. Here, we describe the miRNAs of wild C. korshinskii from Huangling, Yulin, and Dalad Banner, which occur along a precipitation gradient. Using next-generation sequencing technology, we obtained a total of 13 710 681, 15 048 945, and 15 198 442 reads for each location, respectively; after filtering and BLAST analysis, 490 conserved miRNAs and 96 novel miRNAs were characterized from the sRNAome data, with key functions determined using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses. We also designed stem-loop qRT-PCR to validate the expression patterns of 5 conserved miRNAs (miR390, miR398, miR530, miR2119, and miR5559) that obviously responded to water stress in plants grown both under natural and experimental water stress conditions and found that the expression levels of miR2119 and miR5559 were negatively correlated with their predicted target genes. This study is the first to identify miRNAs from wild C. korshinskii and provides a basis for future studies of miRNA-mediated gene regulation of stress responses in C. korshinskii.

Characterization and Function of MicroRNA∗s in Plants

MicroRNAs, a group of non-coding RNA molecules, play essential roles in a wide range of cellular processes in different molecules, cells, and organisms. In plants, microRNAs are a class of 20- to 24-nucleotides endogenous small RNAs that repress gene expression. The microRNA guide strand (miRNA) and its complementary strand (miRNA^∗) both originate from the miRNA/miRNA^∗ duplex. Generally, the guide strands act as post-transcriptional regulators that suppress gene expression by cleaving their target mRNA transcripts, whereas the complementary strands were thought to be degraded as 'passenger strands.' However, the complementary strand has been confirmed to possess significant biological functionality in recent reports. In this review, we summarized the binding characteristics of the miRNA^∗ strands with ARGONAUTE proteins, their tissue-specific accumulations and their biological functions, illustrating the essential roles of miRNA^∗s in biological processes and therefore providing directions for further exploration.

Symbiont modulates expression of specific gene categories in Angomonas deanei

Trypanosomatids are parasites that cause disease in humans, animals, and plants. Most are non-pathogenic and some harbor a symbiotic bacterium. Endosymbiosis is part of the evolutionary process of vital cell functions such as respiration and photosynthesis. Angomonas deanei is an example of a symbiont-containing trypanosomatid. In this paper, we sought to investigate how symbionts influence host cells by characterising and comparing the transcriptomes of the symbiont-containing A. deanei (wild type) and the symbiont-free aposymbiotic strains. The comparison revealed that the presence of the symbiont modulates several differentially expressed genes. Empirical analysis of differential gene expression showed that 216 of the 7625 modulated genes were significantly changed. Finally, gene set enrichment analysis revealed that the largest categories of genes that downregulated in the absence of the symbiont were those involved in oxidation-reduction process, ATP hydrolysis coupled proton transport and glycolysis. In contrast, among the upregulated gene categories were those involved in proteolysis, microtubule-based movement, and cellular metabolic process. Our results provide valuable information for dissecting the mechanism of endosymbiosis in A. deanei.

Plant microRNAs: key regulators of root architecture and biotic interactions

Contents 22 I. 22 II. 24 III. 25 IV. 27 V. 29 VI. 10 31 References 32 SUMMARY: Plants have evolved a remarkable faculty of adaptation to deal with various and changing environmental conditions. In this context, the roots have taken over nutritional aspects and the root system architecture can be modulated in response to nutrient availability or biotic interactions with soil microorganisms. This adaptability requires a fine tuning of gene expression. Indeed, root specification and development are highly complex processes requiring gene regulatory networks involved in hormonal regulations and cell identity. Among the different molecular partners governing root development, microRNAs (miRNAs) are key players for the fast regulation of gene expression. miRNAs are small RNAs involved in most developmental processes and are required for the normal growth of organisms, by the negative regulation of key genes, such as transcription factors and hormone receptors. Here, we review the known roles of miRNAs in root specification and development, from the embryonic roots to the establishment of root symbioses, highlighting the major roles of miRNAs in these processes.

Identification of novel miRNAs from drought tolerant rice variety Nagina 22

MicroRNAs regulate a spectrum of developmental and biochemical processes in plants and animals. Thus, knowledge of the entire miRNome is essential to understand the complete regulatory schema of any organism. The current study attempts to unravel yet undiscovered miRNA genes in rice. Analysis of small RNA libraries from various tissues of drought-tolerant ‘aus’ rice variety Nagina 22 (N22) identified 71 novel miRNAs. These were validated based on precursor hairpin structure, small RNA mapping pattern, ‘star’ sequence, conservation and identification of targets based on degradome data. While some novel miRNAs were conserved in other monocots and dicots, most appear to be lineage-specific. They were segregated into two different classes based on the closeness to the classical miRNA definition. Interestingly, evidence of a miRNA-like cleavage was found even for miRNAs that lie beyond the classical definition. Several novel miRNAs displayed tissue-enriched and/or drought responsive expression. Generation and analysis of the degradome data from N22 along with publicly available degradome identified several high confidence targets implicated in regulation of fundamental processes such as flowering and stress response. Thus, discovery of these novel miRNAs considerably expands the dimension of the miRNA-mediated regulation in rice.

Validation of Suitable Reference Genes for Assessing Gene Expression of MicroRNAs in Lonicera japonica

MicroRNAs (miRNAs), which play crucial regulatory roles in plant secondary metabolism and responses to the environment, could be developed as promising biomarkers for different varieties and production areas of herbal medicines. However, limited information is available for miRNAs from Lonicera japonica, which is widely used in East Asian countries owing to various pharmaceutically active secondary metabolites. Selection of suitable reference genes for quantification of target miRNA expression through quantitative real-time (qRT)-PCR is important for elucidating the molecular mechanisms of secondary metabolic regulation in different tissues and varieties of L. japonica. For precise normalization of gene expression data in L. japonica, 16 candidate miRNAs were examined in three tissues, as well as 21 cultivated varieties collected from 16 production areas, using GeNorm, NormFinder, and RefFinder algorithms. Our results revealed combination of u534122 and u3868172 as the best reference genes across all samples. Their specificity was confirmed by detecting the cycling threshold (C t) value ranges in different varieties of L. japonica collected from diverse production areas, suggesting the use of these two reference miRNAs is sufficient for accurate transcript normalization with different tissues, varieties, and production areas. To our knowledge, this is the first report on validation of reference miRNAs in honeysuckle (Lonicera spp.). Restuls from this study can further facilitate discovery of functional regulatory miRNAs in different varieties of L. japonica.

Identification of Arbuscular Mycorrhiza (AM)-Responsive microRNAs in Tomato

A majority of land plants can form symbiosis with arbuscular mycorrhizal (AM) fungi. MicroRNAs (miRNAs) have been implicated to regulate this process in legumes, but their involvement in non-legume species is largely unknown. In this study, by performing deep sequencing of sRNA libraries in tomato roots and comparing with tomato genome, a total of 700 potential miRNAs were predicted, among them, 187 are known plant miRNAs that have been previously deposited in miRBase. Unlike the profiles in other plants such as rice and Arabidopsis, a large proportion of predicted tomato miRNAs was 24 nt in length. A similar pattern was observed in the potato genome but not in tobacco, indicating a Solanum genus-specific expansion of 24-nt miRNAs. About 40% identified tomato miRNAs showed significantly altered expressions upon Rhizophagus irregularis inoculation, suggesting the potential roles of these novel miRNAs in AM symbiosis. The differential expression of five known and six novel miRNAs were further validated using qPCR analysis. Interestingly, three up-regulated known tomato miRNAs belong to a known miR171 family, a member of which has been reported in Medicago truncatula to regulate AM symbiosis. Thus, the miR171 family likely regulates AM symbiosis conservatively across different plant lineages. More than 1000 genes targeted by potential AM-responsive miRNAs were provided and their roles in AM symbiosis are worth further exploring.

The use of high-throughput sequencing methods for plant microRNA research

MicroRNA (miRNA) acts as a critical regulator of gene expression at post-transcriptional and occasionally transcriptional levels in plants. Identification of reliable miRNA genes, monitoring the procedures of transcription, processing and maturation of the miRNAs, quantification of the accumulation levels of the miRNAs in specific biological samples, and validation of miRNA-target interactions become the basis for thoroughly understanding of the miRNA-mediated regulatory networks and the underlying mechanisms. Great progresses have been achieved for sequencing technology. Based on the high degree of sequencing depth and coverage, the high-throughput sequencing (HTS, also called next-generation sequencing) technology provides unprecedentedly efficient way for genome-wide or transcriptome-wide studies. In this review, we will introduce several HTS platform-based methods useful for plant miRNA research, including RNA-seq (RNA sequencing), RNA-PET-seq (paired end tag sequencing of RNAs), sRNA-seq (small RNA sequencing), dsRNA-seq (double-stranded RNA sequencing), ssRNA-seq (single-stranded RNA sequencing) and degradome-seq (degradome sequencing). In particular, we will provide some special cases to illustrate the novel use of HTS methods for investigation of the processing modes of the miRNA precursors, identification of the RNA editing sites on miRNA precursors, mature miRNAs and target transcripts, re-examination of the current miRNA registries, and discovery of novel miRNA species and novel miRNA-target interactions. Summarily, we opinioned that integrative use of the above mentioned HTS methods could make the studies on miRNAs more efficient.

Noncoding RNAs, Emerging Regulators in Root Endosymbioses

Endosymbiosis interactions allow plants to grow in nutrient-deficient soil environments. The arbuscular mycorrhizal (AM) symbiosis is an ancestral interaction between land plants and fungi, whereas nitrogen-fixing symbioses are highly specific for certain plants, notably major crop legumes. The signaling pathways triggered by specific lipochitooligosaccharide molecules involved in these interactions have common components that also overlap with plant root development. These pathways include receptor-like kinases, transcription factors (TFs), and various intermediate signaling effectors, including noncoding (nc)RNAs. These latter molecules have emerged as major regulators of gene expression and small ncRNAs, composed of micro (mi)RNAs and small interfering (si)RNAs, are known to control gene expression at transcriptional (chromatin) or posttranscriptional levels. In this review, we describe exciting recent data connecting variants of conserved si/miRNAs with the regulation of TFs, such as NSP2, NFY-A1, auxin-response factors, and AP2-like proteins, known to be involved in symbiosis. The link between hormonal regulations and these si- and miRNA-TF nodes is proposed in a model in which different feedback loops or regulations controlling endosymbiosis signaling are integrated. The diversity and emerging regulatory networks of young legume miRNAs are also highlighted.

Identification and Expression Analyses of miRNAs from Two Contrasting Flower Color Cultivars of Canna by Deep Sequencing

miRNAs are endogenous small RNA (sRNA) that play critical roles in plant development processes. Canna is an ornamental plant belonging to family Cannaceae. Here, we report for the first time the identification and differential expression of miRNAs in two contrasting flower color cultivars of Canna, Tropical sunrise and Red president. A total of 313 known miRNAs belonging to 78 miRNA families were identified from both the cultivars. Thirty one miRNAs (17 miRNA families) were specific to Tropical sunrise and 43 miRNAs (10 miRNA families) were specific to Red president. Thirty two and 18 putative new miRNAs were identified from Tropical sunrise and Red president, respectively. One hundred and nine miRNAs were differentially expressed in the two cultivars targeting 1343 genes. Among these, 16 miRNAs families targeting60 genes were involved in flower development related traits and five miRNA families targeting five genes were involved in phenyl propanoid and pigment metabolic processes. We further validated the expression analysis of a few miRNA and their target genes by qRT-PCR. Transcription factors were the major miRNA targets identified. Target validation of a few randomly selected miRNAs by RLM-RACE was performed but was successful with only miR162. These findings will help in understanding flower development processes, particularly the color development in Canna.

Novel and conserved microRNAs in soybean floral whorls

MicroRNAs (miRNAs) correspond to a class of endogenous small non-coding RNAs (19-24 nt) that regulates the gene expression, through mRNA target cleavage or translation inhibition. In plants, miRNAs have been shown to play pivotal roles in a wide variety of metabolic and biological processes like plant growth, development, and response to biotic and abiotic stress. Soybean is one of the most important crops worldwide, due to the production of oil and its high protein content. The reproductive phase is considered the most important for soybean yield, which is mainly intended to produce the grains. The identification of miRNAs is not yet saturated in soybean, and there are no studies linking them to the different floral organs. In this study, three different mature soybean floral whorls were used in the construction of sRNA libraries. The sequencing of petal, carpel and stamen libraries generated a total of 10,165,661 sequences. Subsequent analyses identified 200 miRNAs sequences, among which, 41 were novel miRNAs, 80 were conserved soybean miRNAs, 31 were new antisense conserved soybean miRNAs and 46 were soybean miRNAs isoforms. We also found a new miRNA conserved in other plant species, and finally one miRNA-sibling of a soybean conserved miRNA. Conserved and novel miRNAs were evaluated by RT-qPCR. We observed a differential expression across the three whorls for six miRNAs. Computational predicted targets for miRNAs analyzed by RT-qPCR were identified and present functions related to reproductive process in plants. In summary, the increased accumulation of specific and novel miRNAs in different whorls indicates that miRNAs are an important part of the regulatory network in soybean flower.

A transcriptome-wide, organ-specific regulatory map of Dendrobium officinale, an important traditional Chinese orchid herb

Dendrobium officinale is an important traditional Chinese herb. Here, we did a transcriptome-wide, organ-specific study on this valuable plant by combining RNA, small RNA (sRNA) and degradome sequencing. RNA sequencing of four organs (flower, root, leaf and stem) of Dendrobium officinale enabled us to obtain 536,558 assembled transcripts, from which 2,645, 256, 42 and 54 were identified to be highly expressed in the four organs respectively. Based on sRNA sequencing, 2,038, 2, 21 and 24 sRNAs were identified to be specifically accumulated in the four organs respectively. A total of 1,047 mature microRNA (miRNA) candidates were detected. Based on secondary structure predictions and sequencing, tens of potential miRNA precursors were identified from the assembled transcripts. Interestingly, phase-distributed sRNAs with degradome-based processing evidences were discovered on the long-stem structures of two precursors. Target identification was performed for the 1,047 miRNA candidates, resulting in the discovery of 1,257 miRNA--target pairs. Finally, some biological meaningful subnetworks involving hormone signaling, development, secondary metabolism and Argonaute 1-related regulation were established. All of the sequencing data sets are available at NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra/). Summarily, our study provides a valuable resource for the in-depth molecular and functional studies on this important Chinese orchid herb.

Dissecting miRNAs in Wheat D Genome Progenitor, Aegilops tauschii

As the post-transcriptional regulators of gene expression, microRNAs or miRNAs comprise an integral part of understanding how genomes function. Although miRNAs have been a major focus of recent efforts, miRNA research is still in its infancy in most plant species. Aegilops tauschii, the D genome progenitor of bread wheat, is a wild diploid grass exhibiting remarkable population diversity. Due to the direct ancestry and the diverse gene pool, A. tauschii is a promising source for bread wheat improvement. In this study, a total of 87 Aegilops miRNA families, including 51 previously unknown, were computationally identified both at the subgenomic level, using flow-sorted A. tauschii 5D chromosome, and at the whole genome level. Predictions at the genomic and subgenomic levels suggested A. tauschii 5D chromosome as rich in pre-miRNAs that are highly associated with Class II DNA transposons. In order to gain insights into miRNA evolution, putative 5D chromosome miRNAs were compared to its modern ortholog, Triticum aestivum 5D chromosome, revealing that 48 of the 58 A. tauschii 5D miRNAs were conserved in orthologous T. aestivum 5D chromosome. The expression profiles of selected miRNAs (miR167, miR5205, miR5175, miR5523) provided the first experimental evidence for miR5175, miR5205 and miR5523, and revealed differential expressional changes in response to drought in different genetic backgrounds for miR167 and miR5175. Interestingly, while miR5523 coding regions were present and expressed as pre-miR5523 in both T. aestivum and A. tauschii, the expression of mature miR5523 was observed only in A. tauschii under normal conditions, pointing out to an interference at the downstream processing of pre-miR5523 in T. aestivum. Overall, this study expands our knowledge on the miRNA catalog of A. tauschii, locating a subset specifically to the 5D chromosome, with ample functional and comparative insight which should contribute to and complement efforts to develop drought tolerant wheat varieties.

Characterization of drought- and heat-responsive microRNAs in switchgrass

Recent investigations revealed that microRNAs (miRNAs) play crucial roles in plant acclimation to stress conditions. Switchgrass, one of the important biofuel crop species can withstand hot and dry climates but the molecular basis of stress tolerance is relatively unknown. To identify miRNAs that are important for tolerating drought or heat, small RNAs were profiled in leaves of adult plants exposed to drought or heat. Sequence analysis enabled the identification of 29 conserved and 62 novel miRNA families. Notably, the abundances of several conserved and novel miRNAs were dramatically altered following drought or heat. Using at least one fold (log2) change as cut off, we observed that 13 conserved miRNA families were differentially regulated by both stresses, and, five and four families were specifically regulated by drought and heat, respectively. Similarly, using a more stringent cut off of two fold (log2) regulation, we found 5 and 16 novel miRNA families were upregulated but 6 and 7 families were downregulated under drought and heat, respectively. The stress-altered expression of a subset of miRNAs and their targets was confirmed using quantitative PCR. Overall, the switchgrass plants exposed to drought or heat revealed similarities as well as differences with respect to miRNA regulation, which could be important for enduring different stress conditions.

Identification and functional characterization of soybean root hair microRNAs expressed in response to Bradyrhizobium japonicum infection

Three soybean [Glycine max (L) Merr.] small RNA libraries were generated and sequenced using the Illumina platform to examine the role of miRNAs during soybean nodulation. The small RNA libraries were derived from root hairs inoculated with Bradyrhizobium japonicum (In_RH) or mock-inoculated with water (Un_RH), as well as from the comparable inoculated stripped root samples (i.e. inoculated roots with the root hairs removed). Sequencing of these libraries identified a total of 114 miRNAs, including 22 novel miRNAs. A comparison of miRNA abundance among the 114 miRNAs identified 66 miRNAs that were differentially expressed between root hairs and stripped roots, and 48 miRNAs that were differentially regulated in infected root hairs in response to B. japonicum when compared to uninfected root hairs (P ≤ 0.05). A parallel analysis of RNA ends (PARE) library was constructed and sequenced to reveal a total of 405 soybean miRNA targets, with most predicted to encode transcription factors or proteins involved in protein modification, protein degradation and hormone pathways. The roles of gma-miR4416 and gma-miR2606b during nodulation were further analysed. Ectopic expression of these two miRNAs in soybean roots resulted in significant changes in nodule numbers. miRNA target information suggested that gma-miR2606b regulates a Mannosyl-oligosaccharide 1, 2-alpha-mannosidase gene, while gma-miR4416 regulates the expression of a rhizobium-induced peroxidase 1 (RIP1)-like peroxidase gene, GmRIP1, during nodulation.

Regulation of Non-coding RNAs in Heat Stress Responses of Plants

Heat stress is an important factor limiting plant growth, development, and productivity; thus, plants have evolved special adaptive mechanisms to cope with high-temperature stress. Non-coding RNAs (ncRNAs) are a class of regulatory RNAs that play an important role in many biological processes. Recently developed advanced technologies, such as genome-wide transcriptomic analysis, have revealed that abundant ncRNAs are expressed under heat stress. Although this area of research is still in its infancy, an increasing number of several classes of regulatory ncRNA (i.e., miRNA, siRNA, and lncRNA) related to heat stress responses have been reported. In this mini-review, we discuss our current understanding of the role of ncRNAs in heat stress responses in plants, especially miRNAs, siRNAs, and their targets. For example, the miR398-CSD/CCS-HSF, miR396-WRKY6, miR159-GAMYB, and TAS1-HTT-HSF pathways regulate plant heat tolerance. We highlight the hormone/development-related miRNAs involved in heat stress, and discuss the regulatory networks of miRNA-targets. We also note that DNA methylation and alternative splicing could affect miRNA expression under heat stress, and some lncRNAs could respond to heat stress. Finally, we briefly discuss future prospects concerning the ncRNA-related mechanisms of heat stress responses in plants.

Hyphal Branching during Arbuscule Development Requires Reduced Arbuscular Mycorrhiza1

During arbuscular mycorrhizal symbiosis, arbuscule development in the root cortical cell and simultaneous deposition of the plant periarbuscular membrane generate the interface for symbiotic nutrient exchange. The transcriptional changes that accompany arbuscule development are extensive and well documented. By contrast, the transcriptional regulators that control these programs are largely unknown. Here, we provide a detailed characterization of an insertion allele of Medicago truncatula Reduced Arbuscular Mycorrhiza1 (RAM1), ram1-3, which reveals that RAM1 is not necessary to enable hyphopodium formation or hyphal entry into the root but is essential to support arbuscule branching. In ram1-3, arbuscules consist only of the arbuscule trunk and in some cases, a few initial thick hyphal branches. ram1-3 is also insensitive to phosphate-mediated regulation of the symbiosis. Transcript analysis of ram1-3 and ectopic expression of RAM1 indicate that RAM1 regulates expression of EXO70I and Stunted Arbuscule, two genes whose loss of function impacts arbuscule branching. Furthermore, RAM1 regulates expression of a transcription factor Required for Arbuscule Development (RAD1). RAD1 is also required for arbuscular mycorrhizal symbiosis, and rad1 mutants show reduced colonization. RAM1 itself is induced in colonized root cortical cells, and expression of RAM1 and RAD1 is modulated by DELLAs. Thus, the data suggest that DELLAs regulate arbuscule development through modulation of RAM1 and RAD1 and that the precise transcriptional control essential to place proteins in the periarbuscular membrane is controlled, at least in part, by RAM1.

Comprehensive analysis of small RNA-seq data reveals that combination of miRNA with its isomiRs increase the accuracy of target prediction in Arabidopsis thaliana

Along with the canonical miRNA, distinct miRNA-like sequences called sibling miRNAs (sib-miRs) are generated from the same pre-miRNA. Among them, isomeric sequences featuring slight variations at the terminals, relative to the canonical miRNA, constitute a pool of isomeric sibling miRNAs (isomiRs). Despite the high prevalence of isomiRs in eukaryotes, their features and relevance remain elusive. In this study, we performed a comprehensive analysis of mature precursor miRNA (pre-miRNA) sequences from Arabidopsis to understand their features and regulatory targets. The influence of isomiR terminal heterogeneity in target binding was examined comprehensively. Our comprehensive analyses suggested a novel computational strategy that utilizes miRNA and its isomiRs to enhance the accuracy of their regulatory target prediction in Arabidopsis. A few targets are shared by several members of isomiRs; however, this phenomenon was not typical. Gene Ontology (GO) enrichment analysis showed that commonly targeted mRNAs were enriched for certain GO terms. Moreover, comparison of these commonly targeted genes with validated targets from published data demonstrated that the validated targets are bound by most isomiRs and not only the canonical miRNA. Furthermore, the biological role of isomiRs in target cleavage was supported by degradome data. Incorporating this finding, we predicted potential target genes of several miRNAs and confirmed them by experimental assays. This study proposes a novel strategy to improve the accuracy of predicting miRNA targets through combined use of miRNA with its isomiRs.

micro RNA 172 (miR172) signals epidermal infection and is expressed in cells primed for bacterial invasion in Lotus japonicus roots and nodules

Legumes interact with rhizobial bacteria to form nitrogen-fixing root nodules. Host signalling following mutual recognition ensures a specific response, but is only partially understood. Focusing on the stage of epidermal infection with Mesorhizobium loti, we analysed endogenous small RNAs (sRNAs) of the model legume Lotus japonicus to investigate their involvement in host response regulation. We used Illumina sequencing to annotate the L. japonicus sRNA-ome and isolate infection-responsive sRNAs, followed by candidate-based functional characterization. Sequences from four libraries revealed 219 novel L. japonicus micro RNAs (miRNAs) from 114 newly assigned families, and 76 infection-responsive sRNAs. Unlike infection-associated coding genes such as NODULE INCEPTION (NIN), a micro RNA 172 (miR172) isoform showed strong accumulation in dependency of both Nodulation (Nod) factor and compatible rhizobia. The genetics of miR172 induction support the existence of distinct epidermal and cortical signalling events. MIR172a promoter activity followed a previously unseen pattern preceding infection thread progression in epidermal and cortical cells. Nodule-associated miR172a expression was infection-independent, representing the second of two genetically separable activity waves. The combined data provide a valuable resource for further study, and identify miR172 as an sRNA marking successful epidermal infection. We show that miR172 acts upstream of several APETALA2-type (AP2) transcription factors, and suggest that it has a role in fine-tuning AP2 levels during bacterial symbiosis.

Combined genetic and transcriptomic analysis reveals three major signalling pathways activated by Myc-LCOs in Medicago truncatula

Myc-LCOs are newly identified symbiotic signals produced by arbuscular mycorrhizal (AM) fungi. Like rhizobial Nod factors, they are lipo-chitooligosaccharides that activate the common symbiotic signalling pathway (CSSP) in plants. To increase our limited understanding of the roles of Myc-LCOs we aimed to analyse Myc-LCO-induced transcriptional changes and their genetic control. Whole genome RNA sequencing (RNA-seq) was performed on roots of Medicago truncatula wild-type plants, and dmi3 and nsp1 symbiotic mutants affected in nodulation and mycorrhizal signalling. Plants were treated separately with the two major types of Myc-LCOs, sulphated and nonsulphated. Generalized linear model analysis identified 2201 differentially expressed genes and classified them according to genotype and/or treatment effects. Three genetic pathways for Myc-LCO-regulation of transcriptomic reprogramming were highlighted: DMI3- and NSP1-dependent; DMI3-dependent and NSP1-independent; and DMI3- and NSP1-independent. Comprehensive analysis revealed overlaps with previous AM studies, and highlighted certain functions, especially signalling components and transcription factors. These data provide new insights into mycorrhizal signalling mechanisms, supporting a role for NSP1, and specialisation for NSP1-dependent and -independent pathways downstream of DMI3. Our data also indicate significant Myc-LCO-activated signalling upstream of DMI3 and/or parallel to the CSSP and some constitutive activity of the CSSP.

MicroRNA156 as a promising tool for alfalfa improvement

A precursor of miR156 (MsmiR156d) was cloned and overexpressed in alfalfa (Medicago sativa L.) as a means to enhance alfalfa biomass yield. Of the five predicted SPL genes encoded by the alfalfa genome, three (SPL6, SPL12 and SPL13) contain miR156 cleavage sites and their expression was down-regulated in transgenic alfalfa plants overexpressing miR156. These transgenic plants had reduced internode length and stem thickness, enhanced shoot branching, increased trichome density, a delay in flowering time and elevated biomass production. Minor effects on sugar, starch, lignin and cellulose contents were also observed. Moreover, transgenic alfalfa plants had increased root length, while nodulation was maintained. The multitude of traits affected by miR156 may be due to the network of genes regulated by the three target SPLs. Our results show that the miR156/SPL system has strong potential as a tool to substantially improve quality and yield traits in alfalfa.

Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j-3p in soybean nodulation

Plant microRNAs (miRNAs) play important regulatory roles in a number of developmental processes. The present work investigated the roles of miRNAs during nodule development in the crop legume soybean (Glycine max). Fifteen soybean small RNA libraries were sequenced from different stages of nodule development, including young nodules, mature nodules and senescent nodules. In order to identify the regulatory targets of the miRNAs, five parallel analysis of RNA ends (PARE) libraries were also sequenced from the same stages of nodule development. Sequencing identified 284 miRNAs, including 178 novel soybean miRNAs. Analysis of miRNA abundance identified 139 miRNAs whose expression was significantly regulated during nodule development, including 12 miRNAs whose expression changed > 10-fold. Analysis of the PARE libraries identified 533 miRNA targets, including three nodulation-related genes and eight nodule-specific genes. miR393j-3p was selected for detailed analysis as its expression was significantly regulated during nodule formation, and it targeted a nodulin gene, Early Nodulin 93 (ENOD93). Strong, ectopic expression of miR393j-3p, as well as RNAi silencing of ENOD93 expression, significantly reduced nodule formation. The data indicate that miR393j-3p regulation of ENOD93 mRNA abundance is a key control point for soybean nodule formation.

Discovery of Conservation and Diversification of miR171 Genes by Phylogenetic Analysis based on Global Genomes

The microRNA171 (miR171) family is widely distributed and highly conserved in a range of species and plays critical roles in regulating plant growth and development through repressing expression of SCARECROW-LIKE (SCL) transcription factors. However, information on the evolutionary conservation and functional diversification of the miRNA171 family members remains scanty. We reconstructed the phylogenetic relationships among miR171 precursor and mature sequences so as to investigate the extent and degree of evolutionary conservation of miR171 in Arabidopsis thaliana (L.) Heynh. (ath), grape (Vitis vinifera L.) (vvi), poplar (Populus trichocarpa Torr. & A.Gray ex Hook.) (ptc), and rice (Oryza sativa L.) (osa). Despite strong conservation of over 80%, some mature miR171 sequences, such as ptc-miR171j -l, and -m and osa-miR171g, -h, and -i, have undergone critical sequence variation, leading to functional diversification, since they target non-SCL gene transcript(s). Phylogenetic analyses revealed a combination of old ancestral relationships and recent lineage-specific diversification in the miR171 family within the four model plants. The cis-regulatory motifs on the upstream promoter sequences of miR171 genes were highly divergent and shared some similar elements, indicating their possible contribution to the functional variation observed within the miR171 family. This study will buttress our understanding of the functional differentiation of miRNAs and the relationships of miRNA-target pairs based on the evolutionary history of miRNA genes.

High throughput sequencing of small RNA component of leaves and inflorescence revealed conserved and novel miRNAs as well as phasiRNA loci in chickpea

Among legumes, chickpea (Cicer arietinum L.) is the second most important crop after soybean. MicroRNAs (miRNAs) play important roles by regulating target gene expression important for plant development and tolerance to stress conditions. Additionally, recently discovered phased siRNAs (phasiRNAs), a new class of small RNAs, are abundantly produced in legumes. Nevertheless, little is known about these regulatory molecules in chickpea. The small RNA population was sequenced from leaves and flowers of chickpea to identify conserved and novel miRNAs as well as phasiRNAs/phasiRNA loci. Bioinformatics analysis revealed 157 miRNA loci for the 96 highly conserved and known miRNA homologs belonging to 38 miRNA families in chickpea. Furthermore, 20 novel miRNAs belonging to 17 miRNA families were identified. Sequence analysis revealed approximately 60 phasiRNA loci. Potential target genes likely to be regulated by these miRNAs were predicted and some were confirmed by modified 5' RACE assay. Predicted targets are mostly transcription factors that might be important for developmental processes, and others include superoxide dismutases, plantacyanin, laccases and F-box proteins that could participate in stress responses and protein degradation. Overall, this study provides an inventory of miRNA-target gene interactions for chickpea, useful for the comparative analysis of small RNAs among legumes.

miRVine: a microRNA expression atlas of grapevine based on small RNA sequencing

Background

miRNAs are the most abundant class of small non-coding RNAs, and they are involved in post-transcriptional regulations, playing a crucial role in the refinement of genetic programming during plant development. Here we present a comprehensive picture of miRNA regulation in Vitis vinifera L. plant during its complete life cycle. Furthering our knowledge about the post-transcriptional regulation of plant development is fundamental to understand the biology of such an important crop.

Results

We analyzed 70 small RNA libraries, prepared from berries, inflorescences, tendrils, buds, carpels, stamens and other samples at different developmental stages. One-hundred and ten known and 175 novel miRNAs have been identified and a wide grapevine expression atlas has been described. The distribution of miRNA abundance reveals that 22 novel miRNAs are specific to stamen, and two of them are, interestingly, involved in ethylene biosynthesis, while only few miRNAs are highly specific to other organs. Thirty-eight miRNAs are present in all our samples, suggesting a role in key regulatory circuit. On the basis of miRNAs abundance and distribution across samples and on the estimated correlation, we suggest that miRNA expression define organ identity. We performed target prediction analysis and focused on miRNA expression analysis in berries and inflorescence during their development, providing an initial functional description of the identified miRNAs.

Conclusions

Our findings represent a very extensive miRNA expression atlas in grapevine, allowing the definition of how the spatio-temporal distribution of miRNAs defines organ identity. We describe miRNAs abundance in specific tissues not previously described in grapevine and contribute to future targeted functional analyses. Finally, we present a deep characterization of miRNA involvement in berry and inflorescence development, suggesting a role for miRNA-driven hormonal regulation.

Electronic supplementary material

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

Network of GRAS transcription factors involved in the control of arbuscule development in Lotus japonicus

Arbuscular mycorrhizal (AM) fungi, in symbiosis with plants, facilitate acquisition of nutrients from the soil to their host. After penetration, intracellular hyphae form fine-branched structures in cortical cells termed arbuscules, representing the major site where bidirectional nutrient exchange takes place between the host plant and fungus. Transcriptional mechanisms underlying this cellular reprogramming are still poorly understood. GRAS proteins are an important family of transcriptional regulators in plants, named after the first three members: GIBBERELLIC ACID-INSENSITIVE, REPRESSOR of GAI, and SCARECROW. Here, we show that among 45 transcription factors up-regulated in mycorrhizal roots of the legume Lotus japonicus, expression of a unique GRAS protein particularly increases in arbuscule-containing cells under low phosphate conditions and displays a phylogenetic pattern characteristic of symbiotic genes. Allelic rad1 mutants display a strongly reduced number of arbuscules, which undergo accelerated degeneration. In further studies, two RAD1-interacting proteins were identified. One of them is the closest homolog of Medicago truncatula, REDUCED ARBUSCULAR MYCORRHIZATION1 (RAM1), which was reported to regulate a glycerol-3-phosphate acyl transferase that promotes cutin biosynthesis to enhance hyphopodia formation. As in M. truncatula, the L. japonicus ram1 mutant lines show compromised AM colonization and stunted arbuscules. Our findings provide, to our knowledge, new insight into the transcriptional program underlying the host's response to AM colonization and propose a function of GRAS transcription factors including RAD1 and RAM1 during arbuscule development.

Rapid divergence and high diversity of miRNAs and miRNA targets in the Camelineae

MicroRNAs (miRNAs) are short RNAs involved in gene regulation through translational inhibition and transcript cleavage. After processing from imperfect fold-back structures, miRNAs are incorporated into RNA-induced silencing complexes (RISCs) before targeting transcripts with varying degrees of complementarity. Some miRNAs are evolutionarily deep-rooted, and sequence complementarity with their targets is maintained through purifying selection. Both Arabidopsis and Capsella belong to the tribe Camelineae in the Brassicaceae, with Capsella rubella serving as an outgroup to the genus Arabidopsis. The genome sequence of C. rubella has recently been released, which allows characterization of its miRNA complement in comparison with Arabidopsis thaliana and Arabidopsis lyrata. Through next-generation sequencing, we identify high-confidence miRNA candidates specific to the C. rubella lineage. Only a few lineage-specific miRNAs have been studied for evolutionary constraints, and there have been no systematic studies of miRNA target diversity within or divergence between closely related plant species. Therefore we contrast sequence variation in miRNAs and their targets within A. thaliana, and between A. thaliana, A. lyrata and C. rubella. We document a surprising amount of small-scale variation in miRNA-target pairs, where many miRNAs are predicted to have species-specific targets in addition to ones that are shared between species. Our results emphasize that the transitive nature of many miRNA-target pairs can be observed even on a relatively short evolutionary time-scale, with non-random occurrences of differences in miRNAs and their complements in the miRNA precursors, the miRNA* sequences.

Profiling microRNA expression during multi-staged date palm (Phoenix dactylifera L.) fruit development

MicroRNAs (miRNAs) play crucial roles in multiple stages of plant development and regulate gene expression at posttranscriptional and translational levels. In this study, we first identified 238 conserved miRNAs in date palm (Phoenix dactylifera) based on a high-quality genome assembly and defined 78 fruit-development-associated (FDA) miRNAs, whose expression profiles are variable at different fruit development stages. Using experimental data, we subsequently detected 276 novel P. dactylifera-specific FDA miRNAs and predicted their targets. We also revealed that FDA miRNAs function mainly in regulating genes involved in starch/sucrose metabolisms and other carbon metabolic pathways; among them, 221 FDA miRNAs exhibit negative correlation with their corresponding targets, which suggests their direct regulatory roles on mRNA targets. Our data define a comprehensive set of conserved and novel FDA miRNAs along with their expression profiles, which provide a basis for further experimentation in assigning discrete functions of these miRNAs in P. dactylifera fruit development.

Identification of microRNAs differentially expressed involved in male flower development

Hickory (Carya cathayensis Sarg.) is one of the most economically important woody trees in eastern China, but its long flowering phase delays yield. Our understanding of the regulatory roles of microRNAs (miRNAs) in male flower development in hickory remains poor. Using high-throughput sequencing technology, we have pyrosequenced two small RNA libraries from two male flower differentiation stages in hickory. Analysis of the sequencing data identified 114 conserved miRNAs that belonged to 23 miRNA families, five novel miRNAs including their corresponding miRNA*s, and 22 plausible miRNA candidates. Differential expression analysis revealed 12 miRNA sequences that were upregulated in the later (reproductive) stage of male flower development. Quantitative real-time PCR showed similar expression trends as that of the deep sequencing. Novel miRNAs and plausible miRNA candidates were predicted using bioinformatic analysis methods. The miRNAs newly identified in this study have increased the number of known miRNAs in hickory, and the identification of differentially expressed miRNAs will provide new avenues for studies into miRNAs involved in the process of male flower development in hickory and other related trees.

NPK macronutrients and microRNA homeostasis

Macronutrients are essential elements for plant growth and development. In natural, non-cultivated systems, the availability of macronutrients is not a limiting factor of growth, due to fast recycling mechanisms. However, their availability might be an issue in modern agricultural practices, since soil has been frequently over exploited. From a crop management perspective, the nitrogen (N), phosphorus (P), and potassium (K) are three important limiting factors and therefore frequently added as fertilizers. NPK are among the nutrients that have been reported to alter post-embryonic root developmental processes and consequently, impairs crop yield. To cope with nutrients scarcity, plants have evolved several mechanisms involved in metabolic, physiological, and developmental adaptations. In this scenario, microRNAs (miRNAs) have emerged as additional key regulators of nutrients uptake and assimilation. Some studies have demonstrated the intrinsic relation between miRNAs and their targets, and how they can modulate plants to deal with the NPK availability. In this review, we focus on miRNAs and their regulation of targets involved in NPK metabolism. In general, NPK starvation is related with miRNAs that are involved in root-architectural changes and uptake activity modulation. We further show that several miRNAs were discovered to be involved in plant-microbe symbiosis during N and P uptake, and in this way we present a global view of some studies that were conducted in the last years. The integration of current knowledge about miRNA-NPK signaling may help future studies to focus in good candidates genes for the development of important tools for plant nutritional breeding.

Soybean miR172c targets the repressive AP2 transcription factor NNC1 to activate ENOD40 expression and regulate nodule initiation

MicroRNAs are noncoding RNAs that act as master regulators to modulate various biological processes by posttranscriptionally repressing their target genes. Repression of their target mRNA(s) can modulate signaling cascades and subsequent cellular events. Recently, a role for miR172 in soybean (Glycine max) nodulation has been described; however, the molecular mechanism through which miR172 acts to regulate nodulation has yet to be explored. Here, we demonstrate that soybean miR172c modulates both rhizobium infection and nodule organogenesis. miR172c was induced in soybean roots inoculated with either compatible Bradyrhizobium japonicum or lipooligosaccharide Nod factor and was highly upregulated during nodule development. Reduced activity and overexpression of miR172c caused dramatic changes in nodule initiation and nodule number. We show that soybean miR172c regulates nodule formation by repressing its target gene, Nodule Number Control1, which encodes a protein that directly targets the promoter of the early nodulin gene, ENOD40. Interestingly, transcriptional levels of miR172c were regulated by both Nod Factor Receptor1α/5α-mediated activation and by autoregulation of nodulation-mediated inhibition. Thus, we established a direct link between miR172c and the Nod factor signaling pathway in addition to adding a new layer to the precise nodulation regulation mechanism of soybean.